_黑土滩_退化草地_高寒湿地及其交错区植物群落结构多样性

三江源区“黑土型”退化草地形成机理与恢复模式研究三江源区“黑土型”退化草地形成机理与恢复模式研究摘要：三江源区作为中国重要的生态保护区，其草地生态系统的健康与恢复至关重要。

然而，近年来，三江源区的“黑土型”退化草地面积不断扩大，严重威胁到该区生态环境的可持续发展。

为了深入了解“黑土型”退化草地的形成机理和制定恢复模式，本文通过野外调查、实验室分析及文献研究等方法，对三江源区的“黑土型”退化草地进行综合研究。

一、引言三江源区是我国重要的自然保护区，以其独特的高寒草甸和湿地生态系统而著称。

然而，近年来，该区的草地面积不断减少，退化面积日益扩大，对生态系统的可持续发展产生了极大的威胁。

尤其是“黑土型”退化草地是三江源区面临的一大难题。

二、“黑土型”退化草地的形成机理1. 水土流失土地利用变化、气候变化等原因导致了三江源区土壤的侵蚀和水土流失，加速了“黑土型”退化草地的形成。

2. 过度放牧过度放牧对“黑土型”退化草地的形成起到了重要的推动作用。

过度放牧导致植被被过度踩踏和食用，加剧了草地的退化。

3. 土壤养分失衡长期过度放牧会导致草地土壤中养分的流失，使得土壤变得贫瘠，草地植被无法得到充分的养分供应，进而引发“黑土型”退化草地。

三、“黑土型”退化草地的恢复模式1. 放牧管理措施科学合理的放牧管理有助于改善“黑土型”退化草地。

比如，合理规划放牧牲畜的饲养数量和放牧时间，通过轮牧、休牧等方式实现草地的恢复。

2. 种草补播草地进行种草补播是恢复“黑土型”退化草地的有效手段之一。

通过选择适应当地环境的草种进行补播，可以快速恢复草地的植被覆盖，减缓土壤侵蚀的速度。

3. 抚育措施抚育是指通过人工干预，对退化的草地进行保护和修复。

比如，通过人工修筑草坝、草沟等措施，增加草地的保水性和保肥性，帮助草地恢复。

四、研究结果与讨论通过对三江源区“黑土型”退化草地进行实地调查，发现该区域的退化草地主要集中在放牧区和水源区，且空间分布呈现明显的梯度。

中国环境科学 2021,41(3)：1405~1414 China Environmental Science 甘南高寒草甸植物群落物种多度分布特征刘旻霞*,张娅娅,李全弟,李博文,孙瑞弟,宋佳颖(西北师范大学地理与环境科学学院,甘肃兰州 730070)摘要：为分析甘南高寒草甸植物群落物种多度分布格局对海拔梯度的响应,探讨物种多样性的形成规律及维持机制,采用典型样地调查法采集数据信息,并通过RAD软件程序包对实验数据进行拟合分析.结果表明:植物群落物种多样性随海拔的升高呈现先增加后降低的“中间膨胀”模式,海拔3500m处群落物种多样性达到最大.利用5个模型对海拔梯度上全部物种多度分布进行拟合发现,甘南高寒草甸植物群落物种多度分布符合生态位模型,尤其是geo 模型,利用该模型可以很好地拟合群落物种多度分布格局,同时也说明生态位分化对高寒草甸植物群落物种多度分布起主要作用.常见种多度分布的最优拟合模型为geo模型,其与全部物种的最优拟合模型一致;而稀有种多度分布的最优拟合模型为rane模型,说明常见种是影响不同海拔梯度群落物种多度分布格局的主要原因,其对群落生产力和稳定性的维持有不可替代的作用,但稀有种因其独特的功能性状也在一定程度上影响着群落结构,二者以各自不同的方式共同维系着高寒草甸的物种多样性.关键词：高寒草甸；物种多度分布；生态位模型；常见种；稀有种；生态位分化中图分类号：X176文献标识码：A 文章编号：1000-6923(2021)03-1405-10Species abundance distribution characteristics of alpine meadow plant community in Gannan. LIU Min-xia*, ZHANG Ya-ya, LI Quan-di, LI Bo-wen, SUN Rui-di, SONG Jia-ying (College of Geography and Environmental Sciences, Northwest Normal University, Lanzhou 730070, China). China Environmental Science, 2021,41(3)：1405~1414Abstract：The response of species abundance distribution (SAD) pattern of alpine meadow plant community to altitude gradient was analyzed, and the formation rule and maintenance mechanism of species diversity were discussed. In this paper, the typical sample plot survey method was used to collect data, and the experimental data were fitted and analyzed by RAD (rank abundance distribution) software package. T he results showed that the species diversity of plant community showed a “Hump relationship” pattern with the increase of altitude, and reached the maximum at 3500m. Five models were used to fit the SAD on the altitude, it was found that the overall SAD of alpine meadow plant community was in line with the niche model, especially the geo. The model could be used to fit the distribution pattern of community species abundance, and it also showed that niche differentiation plays a major role in the SAD of alpine meadow plant community. The best fitting model of common SAD was geo, which was consistent with the best fitting model of overall species, while the best fitting model of rare SAD was rane. T he above results showed that common species were the main factors affecting the SAD pattern of communities at different altitudes, and they played an irreplaceable role in maintaining the productivity and stability of communities. However, rare species also affected the community structure to a certain extent because of their unique functional characteristics, and they maintain the stability of species diversity of alpine meadow communities in different ways.Key words：alpine meadow；species abundance；niche model；common species；rare species；niche differentiation物种多度分布格局分析对于理解植物群落物种多样性的形成和维持机制具有重要意义.多度可以反映出一个物种在植物群落中占用资源的能力,不同的植物群落具有不同的多度分布格局.物种多度分布格局是多个物种不断地相互作用、相互影响的结果,通过研究多度分布格局可揭示植物群落的组成,进而了解群落中不同物种间的关系和作用机制[1].物种多度分布和物种多样性指数都是研究群落物种多样性的重要方法,物种多度通过群落内不同物种个体数量的分布情况来反映群落的物种多样性,对于认识一个群落来说,多度分布比多样性指数更能说明问题,因而成为描述群落结构及格局的有力工具[2].物种多度分布研究始于20世纪30年代,最早由Motomura[3]提出几何级数模型用于描述水生动物种多度关系.Fisher等[4]和Preston[5]分别提出对数级数模型与对数正态模型,在对鸟类、昆虫等群落的调查研究中取得了不错的拟合效果,但此两种模型在解释群落生态学过程的形成方面出现了困难,因此,基于生态位分化理论,MacArthur[6]和Sugihara[7]先后提出了断棍模型和连续断裂模型,为对数级数模型收稿日期：2020-07-30基金项目：国家自然科学基金资助项目(31760135)* 责任作者, 教授,***************1406 中 国 环 境 科 学 41卷和对数正态模型提供了生态学解释.之后,Tokeshi [8]对前人的研究做了回顾与总结,进一步提出了幂分割模型并得到了广泛的应用.以Hubbell [9-10]为代表的生态学家为解释生物多样性的分布模式在不同空间尺度上的形成机制,建立了群落中性理论模型,且在热带雨林群落物种多度的研究方面得到了很好的印证.总体而言,多度格局的拟合模型自发展以来大致出现了3种类型,即统计模型、生态位理论模型和中性理论模型;由于统计模型缺乏明确的生态学意义,难以对实际数据做出明确的解释,因而生态位理论和中性理论逐渐成为解释物种多度分布格局的两大理论模型[11-12].近年来,生态位理论模型因其明确的生态学意义而在物种多度分布模式的研究中得到了广泛的应用[13-15].青藏高原是全球海拔最高的一个独特地域单元,受长期的自然演替过程和地质活动的影响,形成了世界上最大的高山生态系统(高寒草甸和草原生态系统),其中青藏高原东部的甘南高寒草甸生物资源和物种多样性相当丰富,其对青藏高原高寒草甸生态系统有着巨大的影响.但由于全球气候变化和频繁的人类活动,甘南高寒草甸草场不断退化,生物多样性锐减,生态系统结构和功能受损[16],因此,保护和恢复甘南高寒草甸生态系统物种多样性刻不容缓.纵观目前有关青藏高原植物群落物种多度分布的研究较多,但大多都集中于不同坡向[15,17]或不同取样面积[18]上,而关于垂直带谱的研究依然比较鲜见,尤其是沿海拔梯度上植物群落物种多度分布格局的相关研究.海拔作为一个重要的环境因子会影响植物的生长和分布,虽然植物的生长并不受海拔的直接影响,但海拔的变化会使温度、光照、水分等因子发生变化从而影响植物的生长及生存,植物为了适应这种改变,会通过其生理结构、形态变化等进行响应[19].基于上述背景,选取青藏高原东北部的甘南高寒草甸区植物群落作为研究对象,采用RAD 软件程序对群落物种多度分布进行模型拟合,以期回答以下3个科学问题:(1)物种多样性随海拔梯度是如何变化的?(2)群落全部物种多度分布格局随海拔的上升是如何变化的?是否可以利用生态位模型进行拟合?(3)常见种和稀有种的多度分布格局有何不同?对以上科学问题的回答可以填补关于垂直带谱上高寒草甸植物群落物种多度分布格局研究的空缺,同时也有助于我们以这样的方法探索高寒草甸植物群落物种多样性的分布规律及形成机制,揭示该地区植物群落构建过程中的资源分配模式,从而为青藏高原高寒草甸的植被恢复和生物多样性的保护提供技术支撑和理论依据. 1 材料与方法 1.1 研究区概况图1 研究区位置示意 Fig.1 Study area location map研究区地处甘肃省西南部的甘南高寒草甸区(33°06′~36º10´N,100º46′~104º44′E)(图1),是典型的高寒草甸生态系统.该地区处于青藏高原与黄土高原的过渡带,年均温2.3℃,气温南高北低,年日照时数2000~2500h 左右,年平均降雨量约600~800mm,多集中在6~9月,降雨变率大,年蒸发量约1238mm.甘南州草地总面积2.72106hm 2,占甘南总面积的70.28%,草地类型多样,其中高寒草甸类面积最大,植被以耐寒的多年生草本植物为主,主要有金露梅(Potentilla fruticosa )、矮嵩草(Kobresia humilis )、嵩3期 刘旻霞等：甘南高寒草甸植物群落物种多度分布特征 1407草(Kobresia Willd )、火绒草(Leontopodium nanum )、垂穗披碱草(Elymusnu t ans Griseb )、条叶银莲花(Anemone t rullifolia )、乳白香青(Anaphalislactea maxim )、珠芽蓼(Polygonum vivipurum )、圆穗蓼(Polygonum macrophyllum )、麻花艽(Gentiana straminea )、蛇含委陵菜(Potentilla kleiniana )、羊茅(Festuca ovina )、米口袋(Gueldenstaedtia multiflora )、莓叶委陵菜(Po t en t illa fragarioides )、甘肃棘豆(Oxy t ropis kansuensis )、马先蒿(Pedicularis resupinata )、华丽龙胆(Gentiana sino -ornata )等.土壤表层有致密紧实的草皮层,腐殖质和有机质含量较高,属于典型的高寒草甸土. 1.2 样地设置2019年7月中旬在甘南高寒草甸进行野外群落学调查并采集数据.在研究区试验点内沿海拔3000m 、3250m 、3500m 、3750m 、4000m 各选3座山体,于坡度接近、坡向相同的山坡分别设置研究样地(每个海拔梯度共设置3个样地作为重复),在每个样地随机布设8个50cm×50cm 的样方,调查并记录每个样方中物种的密度、盖度、高度等基本特征并采集生物量.采用梅花五点法在每个样方内取0~20cm 深处的土样,取得土样,带回实验室进行土壤养分的测定,同时测量该海拔处的土壤温度(ST)、大气温度、坡度(土温、大气温度从8:00~ 18:00每隔1h 进行测量,每组5个重复).样地基本信息见表1.表1 样地基本信息Table 1 Basic information of sample plots海拔(m)经度(E)纬度(N)坡向 坡度(°) 放牧强度 是否围栏3000 101°32′28″ 33°39′44″ 北 18. 3 ± 1.6a 相同 否 3250 101°57′34″ 33°47′39″ 北 18. 6 ± 2.3a 相同 否 3500 101°46′26″ 33°58′50″ 北 19.2 ± 2.8b 相同 否 3750 102°06′38″ 34°27′45″ 北 18.8 ± 1.7a 相同 否 4000 102°52′55″ 34°5′31″ 北19.9 ± 1.7b相同否注:不同字母表示差异显著(平均值±标准误).1.3 研究方法1.3.1 全部物种多度分布模型(1)生态位重叠模型(overlapping niche model,简称over),该模型假设群落总生态位为一条棒,每个物种的多度为棒上随机两点之间的距离,且各物种多度间不存在联系,物种都是按照需求获取生存资源,这样物种间就会出现生态位重叠[20].计算公式如下:121[1]12i i i i P P +⎡⎤=−×−⎢⎥+⎣⎦ (1) (2)几何级数模型(geometric series,简称geo),该模型定义群落中第一个优势种优先占领一部分资源,第二个优势种占领剩余资源的部分,这样依次类推,直到最后一个物种没有剩余资源可占领,体现出了生态位优势占领假说.计算公式如下:()1[]1i i E P P P −=− (2) 式中:P i 为第i 个物种的多度,P 为各物种所占资源占总资源的比例,E 为群落资源总量[20].(3)断棍模型(broken stick model,简称bro),该模型认为所有物种会在群落中同时出现,且其竞争能力相近,其多度反映了竞争物种间资源的随机分配是沿着一个一维梯度进行的.其第i 个物种的多度比例期望值为:11s i x i P s x==∑ (i =1,2,3,4,…,s) (3) 式中:i 为观察群落的物种数,S 为物种总数[20].(4)随机分配模型(random assignment model,简称rane),该模型假定群落中各物种的多度之间没有联系,不存在种间竞争,一般情况下物种不能全部占用其生态位.其第i 个种的期望多度为[20]:0.510.5i s iP =− (4) (5)对数正态分布模型(log -normal distribution model,简称norm),该模型认为群落中总个体数的对数符合正态分布,则第i 个物种的多度为:()()log log i N E μσθ⎡⎤+⎣⎦= (5) 式中:N 表示物种总数,μ和σ分别表示正态分布的均值和方差,θ表示正态偏差[20].1.3.2 常见种、稀有种多度分布模型(1)几何级数模型(geometric series,简称geo),公式同上.1408 中 国 环 境 科 学 41卷(2)随机分配模型(random assignment model,简称rane),公式同上.1.3.3 数据处理及模型拟合优势度检验 使用RAD 软件进行物种多度拟合,采用Excel 2010进行数据处理,用Origin 2017进行绘图.以r (用模型的n 个副本的平均密度进行计算,采用最小二乘法中加入密度最大差异和物种数量的校正因子,从而增强了模型拟合优度的区分能力)、Oc (经过Ulrich 校正过的Preston 倍程分组后的χ2检验法)、CL (置信区间)值作为优势度检验的3个指标(r 、Oc 、CL 值都由RAD 拟合输出),r 和Oc 值大于0小于100,当其越接近10或等于10时拟合效果更好;置信区间CL 值一般默认为是大于或等于95%,等于100%时拟合效果最好.1.3.4 物种重要值 计算物种重要值(Importance value),将其结果作为评价各物种的群落优势度、划分常见种(Z ≥0.01)与稀有种(Z <0.01)的指标[21].其计算如下:()4Z =+++相对多度相对盖度相对频度相对高度式中:Z 为每个样地每种植物的重要值. 1.3.5 物种多样性指数计算(1)Margalef 丰富度指数[22]:()()1ln D S N =− (6)(2)Simpson 指数[23]:21si i N D N =⎛⎞=⎜⎟⎝⎠∑ (7)(3)Shannon -Weiner 指数[24]:()1ln si i i H P P ==−∑ (8)式中:S 为物种数;N 为个体总数;N i 为第i 种物种个体数,i =1,2,3… S;P i 是第i 种比例多度,给定为:P i = N i /N .2 结果与分析2.1 群落物种多样性随海拔梯度的变化越复杂的群落结构,其多样性指数、丰富度指数等就越高,这不仅在水平方向上有所体现,且在垂直方向上也有较大差异.从图2可以看出,随海拔的升高,物种丰富度、Margalef 丰富度指数与Simpson 指数的变化幅度较大且变化趋势基本一致,均随海拔的升高先增大后减小;Shannon -Weiner 指数海拔3500m 最大,其次是海拔3000m,其余3个海拔的变化幅度基本相同.由此表明,群落物种多样性沿海拔梯度出现了较大的差异,中间海拔梯度(3500m)物种多样性最大.3000325035003750 40004.55.05.56.06.57.07.58.08.59.09.5海拔(m)海拔(m)M a r g a l e f 丰富度指数3000 3250 3500 375040000.70.8 0.9 1.0 1.1S i m p s o n 指数3000325035003750 40001.891.921.951.982.01S h a n n o n -W e i n e r 指数3000 3250 3500 37504000物种丰富度图2 物种多样性随海拔梯度的变化Fig.2 Variation of species diversity with altitude gradients2.2 不同海拔全部物种多度分布拟合曲线的变化全部物种多度分布曲线斜率表明(图3),随海拔的升高,曲线斜率呈先减小后增大的趋势,海拔3000m 处曲线斜率为0.059,3500m 处为0.041,4000m3期 刘旻霞等：甘南高寒草甸植物群落物种多度分布特征 1409处为0.112.由此发现在海拔3500m 曲线斜率最小,曲线最为平缓,说明此群落物种分布较为均匀,各物种相对多度差异较小,且此处全部物种相对多度之和在5个海拔中最大、物种丰富度最高、多样性最大.而海拔4000m 恰恰相反.0 20 4060 8010-310-210-1100相对多度物种序列图3 高寒草甸各海拔物种多度分布曲线Fig.3 SAD curve of alpine meadow at different altitudes2.3 全部物种多度分布的模型拟合优势度检验采用over 、geo 、bro 、rane 和norm 模型对不同海拔群落的全部物种多度分布进行拟合优势度检验,结果见表2.海拔3000m 、3250m 、3750m 、4000m 全部物种多度分布均符合geo 、bro 、rane 及over 4种模型,其中geo 模型的拟合效果优于其他3种模型;norm 模型在此4个海拔处都无法拟合,其r 值、Oc 值、CL 值与检验标准差异较大.海拔3500m 处,物种多度分布符合全部5种模型,其中拟合效果最好的是geo 模型(r =10.28,Oc =9.10,CL =1);其次是bro 模型(r =14.49,Oc =10.38,CL =1),其他3个模型的拟合效果一般.表2 高寒草甸不同海拔物种多度分布模型拟合优势度检验 Table 2 Dominance test of SAD model fitting at differentaltitudes in alpine meadow海拔(m)模型rOc Clgeo 21.18 13.71 1.00 bro 41.33 15.90 1.00 3000 over 29.71 12.18 1.00 norm 35.89 43.35 0 rane 32.07 10.92 0.95 geo 12.36 16.01 1.00 bro 28.56 19.90 1.00 3250 over 24.71 12.31 1.00 norm 105.63 37.29 0 rane 54.38 32.56 0.95geo 10.28 9.10 1.00bro 14.49 10.38 1.00 3500 over 48.09 14.33 0.97 norm 28.69 5.44 1.00 rane 18.39 34.84 0.96 geo 30.18 25.37 1.00 bro 35.32 28.72 0.98 3750 over 54.68 28.42 0.96 norm 85.31 112.73 0.69 rane 47.53 29.83 1.00 geo 22.21 32.82 1.00 bro 70.42 18.20 0.96 4000 over 84.33 32.89 0.96 norm 78.18 107.33 0.73 rane 44.68 20.06 0.97 注:r 为采用最小二乘法中加入了密度最大差异和物种数量的校正因子;Oc 为Preston 倍程分组后的χ2检验法;CL 为置信区间.加黑数据为每个海拔处的最优模型拟合结果.2.4 全部物种多度拟合值和实测值的比较0 10 2030 4001020304010-310-210-110001020 30 40 50 3000m相对多度3250m相对多度3500m相对多度 物种序列物种序列物种序列10-3 10-210-1 100 10-310-210-11000 1020 30 40010203010-310-210-11003750m相对多度实测值 拟合值geo bro over norm rane物种序列物种序列4000m相对多度10-310-210-1 100图4 甘南高寒草甸不同海拔实测值和拟合值的对比Fig.4 Comparison of measured and fitted values of alpine meadow at different altitudes1410 中国环境科学 41卷图4表明,全部物种多度拟合值和实测值的变化趋势相似.海拔3000m、3250m、3750m和4000m 处geo、bro、over和rane 4个模型的拟合值和实测值的变化趋势基本一致,均随物种序列的增加逐渐减小;norm模型对此4个海拔群落物种多度的拟合结果与实测值的变化趋势有很大差异,实测值随物种序列的增加逐渐减小,而拟合值随物种序列的增加基本无变化,说明其拟合结果对物种多度分布只有统计学意义,而在资源分配上并无实际意义.海拔3500m处,5种模型的拟合值均随物种序列的增加逐渐减小,和实测值的变化趋势一致,说明5种模型对于该海拔群落物种的资源分配模式都具有一定的指导意义.而在5个模型中,geo模型对各海拔物种多度的拟合效果最好,其均表现出随着物种序列的增加,物种多度的拟合值逐渐小于实测值的特征,说明在资源有限的情况下,不能使每个物种都维持较高的多度,这种逐渐减少的特征也符合geo模型的分布模式.2.5 常见种和稀有种多度的模型拟合优势度检验常见种的模型拟合优势度检验结果表明,geo 模型对各个海拔常见种的拟合效果整体较好,其r 值和Oc值均小于100接近于10,置信区间CL值均接近1;但相比之下,geo模型对海拔3500m常见种多度的拟合效果明显优于其他海拔(r=10.43, Oc=10.59, CL=1).rane模型只对海拔3500m处常见种的多度分布可以拟合,但拟合效果一般(r= 20.44,Oc=59.83, CL=0.96),而对其他4个海拔无法拟合(表3).稀有种的模型拟合优势度检验结果表明,各个海拔稀有种的多度分布均符合rane模型,其中该模型对3500m的拟合效果最好(r=12.32,Oc=8.70, CL=1).geo模型只能拟合海拔3250m和3500m处稀有种的多度分布,但拟合效果都一般,而在其余海拔处无法拟合.表3 常见种和稀有种的模型拟合优势度检验Table 3 Model fitting dominance test of common and rare species物种常见种稀有种海拔(m) 模型r Oc CL r Oc CL 3000 geo 24.58 8.39 1 55.41 4.70 0 rane87.3256.17 0 36.783.13 13250 geo 31.26 18.71 1 47.25 15.72 0.98 rane65.8294.22 0 43.286.71 0.963500 geo 10.43 10.59 1 43.59 4.20 1 rane 20.44 59.83 0.96 12.32 8.70 1 3750 geo 42.58 27.13 0.98 89.05 101.23 0.87 rane116.3254.750.87 8.8337.26 14000 geo 23.46 36.59 1 33.09 115.06 0.79 rane65.28104.820.8532.494.83 1 注:r为采用最小二乘法中加入了密度最大差异和物种数量的校正因子;Oc为Preston倍程分组后的x2检验法;CL为置信区间.2.6 常见种和稀有种多度拟合值和实测值的比较就常见种来说,geo模型对其多度的拟合结果与实测值的变化趋势基本吻合,均随物种序列的增加逐渐减小(图5(a));rane模型只对海拔3500m 群落常见种多度分布有较好的拟合效果,其他4个海拔常见种多度实测值随物种序列的增加逐渐减小,而拟合值随物种序列的增加无任何变化(图5(c)),其拟合结果对常见种多度分布没有任何实际意义.由此说明常见种在环境中资源的获取方式偏向固定分配.就稀有种来说,geo模型只在海拔3250m和3500m处物种多度的拟合值和实测值的变化趋势一致,在其余海拔处拟合值与实测值差异较大,其实测值均随物种序列的增加逐渐减小,而拟合值基本无变化或变化甚微(图5(b));rane模型对各个海拔稀有种的拟合结果与实测值的变化趋势基本一致(图5(d)),其拟合效果明显优于geo模型,说明稀有种在各个海拔的资源分配模式偏向随机分配.3期刘旻霞等：甘南高寒草甸植物群落物种多度分布特征 14110 1020 30 40 50607010-310-210-11000510152025 30 35 10-210-1100(a)常见种geo 3500m4000m 3750m 3250m 3000m 相对多度(b)稀有种geo4000m 3750m 3500m 3250m 3000m相对多度实测值geo510152025 30 3510-210-11000 1020 30 40 50 60708010-310-210-1100(c)常见种rane4000m3750m (d)稀有种rane4000m 3750m 3500m 3250m 3000m 相对多度物种序列3500m 3250m 3000m 相对多度物种序列实测值rane图5 不同海拔常见种、稀有种实测值和拟合值比较Fig.5 Comparison of measured and fitted values of common and rare species at different altitudes3 讨论海拔作为重要的地形因子在很大程度上会影响植物的生长发育,不同海拔常常因其存在异质性变化的水热环境因子而引起植物性状发生适应性变化,形成具有不同生活策略的植物群落[25],因此从海拔梯度方面探讨其对群落物种多度分布的影响对于研究高寒草甸植物群落的物种多样性变化具有重要意义.3.1 物种多样性随海拔梯度的变化物种多度分布曲线的斜率变化可以判断群落的物种丰富度,进而判断群落物种多样性的变化,多度分布曲线斜率越小,则群落物种多样性越大;反之,物种多样性越小[26].本文通过研究海拔梯度上植物群落物种多样性(图2)及物种多度分布(图3)发现,甘南高寒草甸植物群落物种多样性随海拔的变化出现“中间膨胀”模式,这与其他研究结果一致[27],即物种多样性分布存在随海拔的升高先增大后减小的规律,最大值往往出现在中间海拔.其可能的解释是低海拔梯度植物群落在发育过程中,充分的环境资源使其维持着较高的物种多样性,但由于低海拔梯度受外界干扰较大,尤其是人类活动的影响会导致某些物种丧失[19],对植被的空间分布产生较大的影响,这在一定程度上会降低群落物种多样性.在中间海拔3500m 处,外界干扰变小,水热条件及土壤养分达到最优,丰富的环境资源足以支撑更多不同物种的生长发育,其物种多样性也会随之增大.当海拔超过3500m 向更高海拔过渡时,温度逐渐降低,生境越来越严酷,环境的筛选作用降低了群落物种多样性.3.2 全部物种多度分布对海拔梯度的响应由于自然群落结构的复杂性,不同物种多度分布模型表示的生态学意义更能够真实有效地反映群落的内在特征[20].甘南高寒草甸植物群落全部物种多度分布随海拔升高发生了明显的变化,导致与其符合的多度分布模型也随之改变,同时也改变了其群落构建的生态学过程(表2,图4).海拔3000m 和3250m 处,植物群落物种多度分布符合geo 、bro 、over 和rane 4种模型,其中,geo 模型的拟合效果较好,其次是over 模型,这可能是因为低海拔地区环境条件相对优越,植物赖以生存的资源较为充足,物种间的竞争作用较小,确定性的生态位划分在群落构建过程中占据主要地位;但同一群落各物种相似的生态适应性会导致产生较高的生态位重叠,尤其对于物种数量多且差异较小、物种分布比较均匀的群落来说更为显著[28-29].海拔3500m 处群落物种多度分布符合全部5种模型,其中geo 模型的拟合效果甚优,其次是bro 模型,而geo 模型和bro 模型都代表了资1412 中国环境科学 41卷源分配的确定性特征.这可能是因为此处属中海拔地段,水热条件适合更多的物种生存,因此物种丰富度较高,根据生态位补偿效应可知不同物种能更有效地获得较为充足的生存资源,相比竞争,资源共享在此处可能更为突出,所以其资源分配很大程度上属于固定分配模式,这与其它研究结果一致[30].海拔3750m处,除norm外其他4种模型均可拟合群落物种多度分布,但其拟合优势度都较低,这可能是因为该海拔处于3500m和4000m的过渡地带,逐渐减少的环境资源在一定程度上无法满足群落所有物种的发展需求,植物的生长处于不稳定的状态,导致群落物种分布格局受制于多个生态学过程.在海拔4000m处,同样是geo模型的拟合效果相对较好,可能的解释是高海拔地区环境条件恶劣,其生境因子对物种的生长要求严苛,且在此处物种赖以生存的资源有限,物种之间存在较为激烈的竞争,受这2大因素的影响,一般的物种难以在此处存活,尤其是敏感度高的稀有种[31],在长期的资源竞争中,一些弱势群体被淘汰,未被淘汰的物种各自占据着所需生态位稳定共存.由此说明,生态位模型可以很好地拟合青藏高原高寒草甸植物群落物种多度分布,尤其是geo模型,但这并不说明其群落多度分布格局一定符合geo模型所代表的生态学过程,由于geo模型划分资源的方式以确定性模式为主,因此至少可以推测出高寒草甸群落物种多度分布主要受生态位分化作用的影响,这与其他研究结果一致[15,18].3.3 常见种和稀有种多度分布随海拔梯度的变化通过常见种和稀有种的多度分布模式来探讨植物群落物种多度分布的变化是研究植物群落物种多度分布格局的一种重要且有效的方法[32].常见种多度分布模式通常趋向于群落整体多度分布模式或与之持有较高的一致性,而稀有种因其独特的功能性状对于群落物种多样性的增加及生态系统功能的维持也有突出的贡献.王世雄等[33]通过研究黄土高原辽东栎群落稀有种和常见种对物种多样性的贡献发现,稀有种和常见种对物种多样性的贡献存在空间尺度差异性,并且这种差异不依赖于稀有种和常见种的相对比例.刘旻霞等[34]研究了甘南亚高寒草甸稀有种对物种多样性和物种多度分布格局的贡献,结果表明稀有种在甘南亚高寒草甸物种多样性中的相对贡献高于非稀有种.本研究发现,常见种的最优拟合模型为geo模型,其与全部物种多度分布的最优模型一致(表3,图5a,图5c),说明其多度分布模式是引起群落物种多度分布形成的主要原因,同时也说明常见种在生境中获取资源的方式可能以固定分配为主.这是因为作为常见种的矮嵩草、嵩草等植物,它们也是高寒草甸的优势种或亚优势种,这些物种具有较强的适应严酷环境的功能性状,生存能力较强,它们对群落类型的形成和物种丰富度的主导作用受自身比例及稀有种比例等因素的影响甚微[35].就稀有种来说,rane模型可以拟合所有海拔梯度上稀有种的多度分布;geo模型除了在海拔3250m和3500m处可以拟合外,其他海拔皆无法拟合(表3,图5b,图5d).其中,rane模型的拟合效果优于geo模型,说明稀有种的资源分配模式更倾向于随机分配.这可能是因为稀有种本身对各类因素的干扰更为敏感,竞争力无法与其他常见种或者优势物种相抗衡,在激烈的竞争过程中会逐渐减少甚至丧失,其生长状态缺乏稳定性,所以对于稀有种来说,其资源获取模式可能更偏向随机分配,这与相关研究[15][36][37]结果一致.由此表明,群落内物种分布不只遵循一种规律,利用不同模型拟合群落内常见种和稀有种多度分布效果可能会更好[8].4结论4.1甘南高寒草甸植物群落物种多样性随海拔的不同出现“中间膨胀”模式.4.2物种多度分布随海拔的不同出现了较为显著的变化,海拔3500m处群落多度拟合符合geo、bro、over、rane及norm 5种模型;其余4个海拔除norm 模型外,其他4种模型均可拟合;在5个海拔梯度上,geo模型的拟合效果最佳;生态位分化在高寒草甸植物群落的构建过程中起主导作用,利用生态位模型(尤其是geo模型)可以很好地拟合高寒草甸植物群落物种多度分布.4.3常见种多度分布的最优拟合模型为geo模型,与全部物种多度分布的最优拟合模型一致,稀有种多度分布的最优拟合模型为rane模型;常见种对不同海拔梯度群落物种多度分布的影响高于稀有种,其对群落生产力和稳定性的维持有不可替代的作用,但稀有种因其独特的功能性状也在一定程度上影响着群落结构,二者以各自不同的方式共同维系。

黄河源区高寒草原的植被退化与土壤退化特征一、本文概述《黄河源区高寒草原的植被退化与土壤退化特征》一文旨在深入探讨黄河源区高寒草原的植被和土壤退化现象，分析其原因，并提出相应的应对策略。

文章首先概述了黄河源区高寒草原的地理特征和生态环境，然后详细分析了植被退化和土壤退化的现状及其特征，接着探讨了导致退化的主要因素，包括气候变化、人类活动、过度放牧等。

在此基础上，文章进一步分析了植被和土壤退化对当地生态系统和人类生产生活的影响，最后提出了针对性的建议和措施，以期为保护和恢复黄河源区高寒草原的生态环境提供参考。

本文的研究不仅对于理解高寒草原生态系统的脆弱性和稳定性具有重要意义，也为全球气候变化背景下的生态保护提供了有益借鉴。

二、研究区域与方法黄河源区位于中国青海省境内，地理位置介于北纬32°30′至36°20′，东经95°30′至102°30′之间。

该区域是黄河的发源地，拥有广阔的高寒草原生态系统，海拔在3000米以上。

黄河源区的气候属于高原寒带半湿润气候，冬季漫长且寒冷，夏季短暂而温和。

年降水量较低，且主要集中在夏季。

黄河源区的高寒草原生态系统对维护区域生态平衡和黄河水源涵养具有极其重要的意义。

本研究采用野外调查与室内分析相结合的方法，对黄河源区高寒草原的植被退化和土壤退化特征进行深入研究。

在黄河源区内选择了具有代表性的高寒草原样地，进行植被和土壤的详细调查。

通过设立固定样方，对样方内的植被种类、数量、高度、盖度等参数进行定期观测和记录。

同时，采集土壤样品，分析土壤的物理和化学性质，如土壤质地、pH值、有机质含量、全氮、全磷、全钾等。

将采集的植被和土壤样品带回实验室进行进一步分析。

植被样品通过烘干、称重等方法，计算植被的生物量。

土壤样品则通过一系列化学和物理分析方法，如粒度分析、酸碱度测定、有机质测定、元素分析等，评估土壤的退化程度。

对野外调查和室内分析得到的数据进行统计和分析，采用描述性统计、方差分析、相关分析等方法，揭示黄河源区高寒草原植被退化和土壤退化的特征及其相互关系。

三江源自然保护区有何重要意义三江源自然保护区纪念碑三江源自然保护区主要保护对象为江河源头的生态系统、高原湿地、高原特有野生动物、高原原始地貌和世界第三极自然景观。

青海三江源地区地处青藏高原腹地，是长江、黄河、澜沧江三大河流的发源地。

历史上，三江源区曾是水草丰美、湖泊星罗棋布、野生动物种群繁多的高原草原草甸区，被称为生态“处女地”。

近些年来，随着全球气候变暖，冰川、雪山逐年萎缩，直接影响高原湖泊和湿地的水源补给，众多的湖泊、湿地面积缩小甚至干涸，沼泽地消失，泥炭地干燥并裸露，沼泽低湿草甸植被向中旱生高原植被演变，生态环境已十分脆弱。

随着人口的无节制增加和人类无限度的生产经营活动，又大大加速了该地区生态环境恶化的进度。

特别是草地大规模的退化与沙化，不仅使该地区草地生产力和对土地的保护功能下降，优质牧草逐渐被毒、杂草所取代，一些草地危害动物如鼠类乘虚而入，导致草地载畜量减少，野生动物栖息环境质量减退，栖息地破碎化，生物多样性降低。

更为重要的是，随着源区植被与湿地生态系统的破坏，水源涵养能力急剧减退，导致三江中下游广大地区旱涝灾害频繁、工农业生产受到严重制约，并已直接威胁到了长江、黄河流域、乃至东南亚诸国的生态安全。

1、是亚洲和我国大部分地区的“生命之源”三江源风光青藏高原被称为“世界屋脊”、地球的“第三极”，青藏高原的隆起打乱了行星风系的临界尺度，迫使大气环流改变行径，成为一个独立的气候区域，孕育了黄河、长江、澜沧江、恒河、印度河等国内外许多著名的河流，是欧亚大陆上大江大河发育最多的区域。

三江源地处青藏高原腹地，起着各江河水文循环的初始作用。

据资料，长江总水量的25%、黄河总水量的49%、澜沧江总水量的15%来自该青海地区。

2、是世界上海拔最高面积最大的高原湿地区三江源地区湿地风光湿地与森林、海洋并称为全球三大生态系统。

湿地生态系统具有水陆过渡性、系统脆弱性、功能多样性和结构复杂性的基本特征，在水源涵养、减缓径流、蓄洪防旱、防灾抗灾、降解污染、维持生物多样性，为人类生产、生活提供多种资源、调节气候等方面有着其它生态系统不可替代的作用，具有巨大的生态功能，在生态安全体系中独具特色。