生态学报  2020, Vol. 40 Issue (8): 2802-2811

文章信息

张倩, 杨晶, 姚宝辉, 蔡志远, 孙小妹, 王缠, 郭怀亮, 谭宇尘, 苏军虎
ZHANG Qian, YANG Jing, YAO Baohui, CAI Zhiyuan, SUN Xiaomei, WANG Chan, GUO Huailiang, TAN Yuchen, SU Junhu
放牧管理模式对高寒草甸鼢鼠鼠丘群落演替的影响
Effects of different grazing strategies on zokor mounds plant community succession in alpine meadow
生态学报. 2020, 40(8): 2802-2811
Acta Ecologica Sinica. 2020, 40(8): 2802-2811
http://dx.doi.org/10.5846/stxb201902120249

文章历史

收稿日期: 2019-02-12
网络出版日期: 2019-12-31
放牧管理模式对高寒草甸鼢鼠鼠丘群落演替的影响
张倩1,2 , 杨晶1,2 , 姚宝辉1,2 , 蔡志远1,2 , 孙小妹2,3 , 王缠1,2 , 郭怀亮1,2 , 谭宇尘1,2 , 苏军虎1,2     
1. 甘肃农业大学草业学院/草业生态系统教育部重点实验室/甘肃省草业工程实验室/中美草地畜牧业可持续发展研究中心, 兰州 730070;
2. 甘肃农业大学-新西兰梅西大学草地生物多样性研究中心, 兰州 730070;
3. 甘肃农业大学资源与环境学院, 兰州 730070
摘要: 高寒草甸放牧利用下高原鼢鼠(Eospalax baileyi)等危害的发生是草地管理的关键难题,分析放牧管理模式对鼢鼠鼠丘植被群落演替的影响能为草地管理提供重要依据。研究选择划区轮牧(RG)、生长季休牧(GSG)、连续放牧(CG)和禁牧(PG)4种放牧管理模式,以及各模式下不同年限鼠丘(一年(ZM1)、两年(ZM2)、三年(ZM3)和多年鼠丘(ZMM))与对照(CM)草地。分析不同放牧管理模式对鼠丘植被群落特征和生物量等的影响,结果发现:PG和GSG下所有年限鼠丘的植被高度、地上生物量均高于RG和CG;RG和GSG下ZM1和ZM2物种Shannon-wiener指数均高于对照样地物种Shannon-wiener指数。主成分分析表明:RG下地上生物量和物种丰富度指数是影响鼠丘植被群落演替的重要因子,PG下地上生物量、盖度、Shannon-wiener指数和均匀度指数是鼠丘植被群落演替的重要因子,CG下物种丰富度和重要值是影响鼠丘植被群落演替的重要指标,GSG下Shannon-wiener指数、盖度和高度是影响鼠丘植被群落演替的重要指标。可见,不同放牧制度对鼠丘植被群落演替的影响不同,禁牧和生长季休牧管理模式能够较好地恢复鼠丘植被群落演替。
关键词: 放牧管理    高原鼢鼠鼠丘    群落特征    演替    
Effects of different grazing strategies on zokor mounds plant community succession in alpine meadow
ZHANG Qian1,2 , YANG Jing1,2 , YAO Baohui1,2 , CAI Zhiyuan1,2 , SUN Xiaomei2,3 , WANG Chan1,2 , GUO Huailiang1,2 , TAN Yuchen1,2 , SU Junhu1,2     
1. College of Grassland Science, Key Laboratory of Grassland Ecosystem(Ministry of Education), Pratacultural Engineering Laboratory of Gansu Province, Sino-U. S. Centers for Grazing land Ecosystem Sustainability, Gansu Agricultural University, Lanzhou 730070, China;
2. Gansu Agricultural University-Massey University Research Centre for Grassland Biodiversity, Gansu Agricultural University, Lanzhou 730070, China;
3. College of Resource and Environmental Science, Gansu Agricultural University, Lanzhou 730070, China
Abstract: The rodent damage by plateau zokor in the process of alpine meadow grazing is a key problem in its utilization and management. Grazing is one of the most important ways of using grassland. The grazing system is used as a utilization system in grazing management. It regulates the use of livestock in the time and space of grazing land, and adjusts the grazing intensity and grazing frequency to balance vegetation grow and livestock. A reasonable grazing system can improve the diversity of grassland plant communities and the stability of grassland ecosystems, and is of great significance for promoting the healthy development and sustainable use of grassland. Plateau zokor is one of the dominant mounds of alpine grassland on the Qinghai-Tibet Plateau. It mainly feeds 3-20 cm of plant roots and stems below the surface. Excavation and hill-building activities of plateau zokor is one of the important driving factors for grassland degradation. The zokor mound is one of the important ways to affect grassland. It is mainly reflected in the coverage of grassland original vegetation, the reduction of species diversity in the community and affecting the productivity of grassland vegetation and the composition of aboveground vegetation. Vegetation community succession of zokor mounds under different grazing management models can provide important basis for grassland management. To analyze the effects of different grazing management models on the vegetation community characteristics and biomass of the rodent mound, four different grazing management models were selected in this study: rotational grazing (RG), growing season grazing (GSG), continuous grazing (CG) and prohibition grazing (PG), with different ages of zokor mounds in each model including one-year (ZM1), two-year (ZM2), three-year (ZM3), multi-year (ZMM) zokor mounds, and control meadow (CM). The results showed that under PG and GSG, the vegetation height and aboveground biomass of zokor mounds in all ages were higher than those of RG and CG. The Shannon-Werner index of ZM1 and ZM2 species under RG and GSG was higher than the species of control meadow. The results of principal component analysis (PCA) showed that aboveground biomass and species richness under RG, while aboveground biomass, coverage, Shannon-Wiener diversity index, and uniformity under PG were important factors of vegetation community succession on zokor mounds. Species richness and important value under CG and Shannon-Wiener diversity index, coverage and height under GSG were critical indicators affecting the vegetation community succession of zokor mounds. In conclusion, different grazing modes have different effects on vegetation community succession of zokor mounds. It can be effectively restored under PG and GSG.
Key Words: grazing management    plateau zokor mounds    community characteristics    succession    

放牧是草地最基本且尤为重要的利用方式之一, 放牧制度作为放牧管理中的利用体系, 从时间和空间角度出发, 规定牛、羊等对草地的使用强度和频率, 保持植物生长和家畜数量之间达到平衡[1]。合理有效的放牧制度对增加草地物种多样性具有重要意义, 也能促进草地健康可持续发展[2]。青藏高原草地退化日益加剧, 过度放牧是导致草地退化的重要原因之一[3]。而家畜主要通过采食、践踏和排泄粪便影响草地植物群落, 以及对物种竞争等产生不同程度的干扰[4]。放牧优化假说认为, 当放牧强度处于中等程度时, 草地初级生产力最大, 即通过优化的放牧制度能够增加草地初级生产力[5], 现有的放牧制度如划区轮牧是放牧强度在空间和时间上调整, 而生长季休牧的放牧制度是放牧强度在时间上的调整, 而连续放牧的管理制度并没有考虑放牧强度的时间和空间效应[6]。不同的放牧制度对群落的优势种、植被特征和草地生产力等产生不同的影响[7], 禁牧有利于调节食草动物与植物的关系, 也是草地管理最常用的手段, 可以在短期内恢复草地建群种的盖度及促进草地生态系统良性循环[8], 但长期禁牧会导致草地生产力降低等[9], 不同的放牧制度, 如划区轮牧、生长季休牧和连续放牧等对草地的影响不同, 客观准确评价放牧制度对生态、生产之间的权衡, 具有重要的现实指导意义。

高寒草甸放牧利用下高原鼢鼠(Eospalax baileyi)等危害的发生是草地管理的关键难题之一, 高原鼢鼠是青藏高原高寒草地优势害鼠之一, 主要采食地表以下3—20 cm的植物根、茎[10]。高原鼢鼠的挖掘、造丘活动是草地退化的重要驱动因子之一, 其所形成的鼠丘是影响草地的重要方式之一, 主要表现在对草地原有植被的覆盖、降低群落中物种多样性和改变植被空间结构, 同时也影响草地植物生产力和地上植被组成等[11]。Koide等[12]研究表明, 鼢鼠扰动会降低优势禾草的地上和地下生物量。不同放牧制度下, 家畜的践踏、排泄物及摄食方式均会对啮齿动物的种群数量产生不同影响[13]。中度干扰假说认为, 中等程度的干扰, 一方面会加速竞争力较弱的物种获得资源, 提高草地物种多样性[14];另一方面会导致时间和空间的异质性, 有利于生态位分化[15]。研究发现, 休牧和轮牧下显著降低禾本科和杂类草的生物量和丰富度, 而高原鼢鼠干扰对草地植物的物种丰富度、Shannon-wiener指数和Simpson指数等均无显著影响[16], 但轻度或中度高原鼢鼠干扰能够提高群落物种多样性, 以及物种多样性指数(Simpson指数和Shannon-wiener指数)随鼠丘群落演替的增加而增加[17-18]。不同放牧制度对高原鼢鼠鼠丘的影响不同, 不同的放牧制度也会导致不同的生态机遇, 给鼠丘植被演替增加机会, 可能会有不同的演替情况, 以及使群落中物种间的竞争发生变化。而鼠丘植被的演替变化, 直接影响草地的利用以及高原鼢鼠危害的发生情况等, 从放牧制度对鼠丘植被演替影响的认知中能获得更有效的草地管理模式[19]。然而迄今未见该方面的报道, 在草地持续退化的态势下, 从鼠丘植被特征以及优化放牧管理策略的角度入手, 对退化草地的综合治理和有效管理具有重要意义。基于此, 本研究通过放牧管理模式对高原鼢鼠鼠丘群落结构影响的研究, 探讨放牧和高原鼢鼠扰动下鼠丘植被生物量(地上和地下)、物种多样性和群落演替, 旨在为高原鼢鼠放牧管理, 以及高寒草甸资源的合理利用和退化草地的恢复与重建提供理论依据。

1 研究区域与方法 1.1 研究区域概况

研究区位于甘肃省天祝县抓喜秀龙镇甘肃农业大学天祝高山草原试验站(37°11′24″—37°12′ N, 102°45′—102°46′48″ E), 平均海拔在2900—3100 m, 年均气温1℃, 年均日照时数2500—2700 h, 降水主要集中在6—8月, 年均降雨量400 mm, 年均蒸发量1600 mm, 无绝对无霜期, 亚高山土壤类型[20]。草地类型为高寒草甸, 主要优势种为冰草(Agropyron cristatum)、苔草(Carex supina)、黄花棘豆(Medicago falcata)、二裂萎菱菜(Potentilla bifurca)、鹅绒萎菱菜(Potentilla anserina)、密花香薷(Elsholtzia densa)等杂类草为主。

1.2 样地设置

2018年8月初在试验站附近选择划区轮牧(Rotational grazing, RG, 面积0.95 hm2, 放牧时间从2012年到2018年, 轮牧时间为180 d, 轮牧次数为3次, 轮牧强度为5.38羊hm2/a, 轮牧起始时间为5月15日, 结束时间为11月15日)、连续放牧(Continuous grazing, CG, 面积1.44 hm2, 全年连续放牧, 放牧率为8.37羊hm2/a)、生长季休牧(Growing season grazing, GSG, 面积1.47hm2, 时间为2015年到2018年, 生长季放牧率为0, 非生长季放牧率为4.53羊hm2/a, 每年在4月15日到9月15日间休牧, 其余时间自由放牧)和禁牧(Prohibition grazing, PG, 面积5.2 hm2, 时间为2011年到2018年, 全年禁牧), 每个草地利用模式样地仅设置1个样地。4种放牧制度样地, 在每个样地中选择高原鼢鼠鼠丘样地以及在鼢鼠鼠丘周围未干扰的草地中选定对照(CM)植物群落, 以空间序列代替演替时间序列的方法[17], 研究不同放牧制度对不同演替阶段高原鼢鼠土丘植物群落结构的影响。结合土壤性状、植被盖度及物种组成等判断鼠丘年限, 具体为一年鼠丘(ZM1), 即1次越冬, 土壤较疏松, 植被总盖度大于20%, 一年植物占90%以上;两年鼠丘(ZM2), 即2次越冬, 土壤较紧实, 植被总盖度大于60%, 多年禾本科植物小于15%;三年鼠丘(ZM3), 即3次越冬, 土壤紧实, 植被总盖度大于75%, 一年植物小于10%;多年鼠丘(ZMM), 即4次以上越冬, 土壤紧实, 植被总盖度大于85%, CM样地, 土壤紧实, 植被总盖度大于90%[21]

1.3 取样及计算方法 1.3.1 群落特征和生物量取样

分别在每个样地中选择12个鼠丘(每个演替阶段选择3个, 鼠丘可近似看成椭圆形, 各年限鼠丘的长轴×短轴均在45 cm×45 cm — 65 cm×68 cm范围之内)样方(50 cm×50 cm), 以及在鼢鼠鼠丘周围未干扰的草地(距离鼠丘至少50 cm以外)中随机选取3个50 cm×50 cm的自然植被样方作为对照(CM), 用针刺法测定物种总盖度;用卷尺测量物种的高度;用样圆测定物种频度, 样圆面积0.25 m2。用最大值标准化法计算物种高度、盖度和频度的相对值。在每个50 cm×50 cm样方内, 齐地面刈割样方内所有植物放入信封。在烘箱中恒温70℃烘干48 h至恒量后, 称重记录。地下生物量测定采用20 cm×20 cm的铁丝框, 在地上部分刈割完的样方内, 分0—10 cm、10—20 cm和20—30 cm层取出地下土壤样品装入60目的网袋, 抖落松散土壤并用清水冲洗干净, 在70℃的烘箱内烘干48 h至恒量, 并称重记录, 生物量都以g/m2表示。

1.3.2 物种重要值、丰富度指数、Pielou均匀度指数和Shannon-wiener指数计算方法

式中, IV为物种的重要值, A为相对高度, B为相对盖度, C为相对频度;O为物种丰富度指数, S为样方中物种总数, N代表样方各物种重要值之和;Pi为物种i的相对重要值, Jsw为物种Pielou均匀度指数;H为物种Shannon-wiener指数[22]

1.4 数据分析

应用Excel 2007对所得数据进行简单记录和处理并制作相关图表;采用SPSS 19.0软件对同一放牧方式下不同年限鼠丘的群落特征进行单因素方差分析。

对不同放牧制度下群落演替的驱动因子进行主成分分析, 主要步骤为:(1)将选取的指标变量建立原始数据矩阵;(2)原始数据矩阵的标准化;(3)依据特征值≥ 1的原则提取公因子, 确定公因子数, 并运用Varimax旋转法对初始因子进行旋转, 以便给予明确的因子解释[23]

2 结果与分析 2.1 两种因素影响下群落特征变化

在交互作用下, 群落的盖度、高度、重要值、丰富度指数O、Shannon-wiener指数H、地上生物量均有显著影响(P < 0.05), Pielou均匀度指数Jsw和地下生物量交互作用不显著(P>0.05)(表 1)。

表 1 两种因素影响下群落特征变化 Table 1 Changes in community characteristics under the influence of two factors
指标
Index
因素A(年限鼠丘)
Factor A (different zokor mounds)
因素B(放牧方式)
Factor B (grazing method)
交互作用(A/B)
Interaction (A/B)
df F P df F P df F P
盖度Coverage 4 84.678 0.000 3 14.665 0.000 12 2.863 0.006
高度Height 4 4.310 0.005 3 40.665 0.000 12 3.073 0.004
重要值Important value 4 5.342 0.002 3 2.995 0.042 12 2.163 0.034
丰富度指数Richness index 4 0.652 0.629 3 2.963 0.043 12 3.269 0.002
均匀度指数Uniformity index 4 1.152 0.346 3 2.663 0.061 12 1.493 0.167
Shannon-wiener指数 4 1.327 0.277 3 4.416 0.009 12 4.519 0.000
地上生物量AGB 4 86.400 0.000 3 21.169 0.000 12 3.181 0.003
地下生物量UGB 4 2.420 0.064 3 1.842 0.155 12 0.128 1.000
AGB:地上生物量, Above ground biomass; UGB:地下生物量, Under ground biomass
2.2 不同放牧管理模式对鼠丘群落特征的影响

同一放牧管理模式下, 分析不同年限鼠丘植被群落总盖度、重要值、丰富度、均匀度和优势种等群落特征的变化, 结果发现:GSG和PG下鼠丘物种平均高度均高于CG和RG;RG和PG下ZM3和ZMM物种重要值均高于CM, 但差异不显著(P>0.05), CG和GSG下CM物种重要值均最高, RG和PG下ZM3物种重要值最高, RG和PG下CM与所有年限鼠丘物种重要值之间无显著性差异(P>0.05);RG和PG下ZM3物种Jsw高于CM, CG和GSG下CM物种Jsw最高, PG、CG和GSG下所有年限鼠丘物种Jsw与CM物种Jsw之间均无显著性差异(P>0.05), RG下CM物种Jsw显著低于ZM1、ZM2和ZMM物种Jsw(P < 0.05);RG下ZM1物种O显著高于CM物种O(P < 0.05), 且ZM3物种O显著高于ZM1和ZM2物种O(P < 0.05), PG和GSG下ZM2和ZM3物种O均高于CM物种O, 但差异不显著(P>0.05), CG下CM物种O显著高于ZM1物种O(P < 0.05)(表 2)。

表 2 不同放牧管理模式下鼠丘群落特征比较 Table 2 Comparison of characteristics of zokor mounds community under different grazing management modes
样地
Plots
鼠丘年限
Zokor mounds
总盖度/%
Total coverage
平均高度
Average height
重要值
Important value
均匀度指数
Uniformity index
丰富度指数
Richness index
优势种
Dominant species
RG ZM1 40.00±7.00c 3.46±0.71ab 0.96±0.05a 0.92±0.01a 3.56±0.05a 冰草Agropyron cristatum
苔草Carex supina
二裂萎菱菜Potentilla bifurca
ZM2 55.67±6.39bc 2.10±0.22b 1.07±0.03a 0.92±0.00a 3.30±0.05ab 冰草A. cristatum
兰石草Lancea tibetica
黄花棘豆Oxytropis ochrocephala
ZM3 63.00±5.86b 2.60±0.29b 1.21±0.04a 0.90±0.02ab 2.42±0.09c 黄花棘豆O. ochrocephala
苔草C. supina
二裂萎菱菜P. bifurca
ZMM 89.67±1.20a 3.49±0.24ab 1.15±0.16a 0.91±0.00a 3.11±0.35abc 冰草A. cristatum
苔草C. supina
黄花棘豆O. ochrocephala
CM 99.00±0.58a 4.61±0.41a 1.06±0.11a 0.85±0.02b 2.72±0.36bc 苔草C. supina
冰草A. cristatum
黄花棘豆O. ochrocephala
PG ZM1 46.33±3.84c 9.07±3.45b 1.38±0.23a 0.90±0.03a 2.12±0.19a 冰草A. cristatum
黄花棘豆O. ochrocephala
紫菀Aster tataricus
ZM2 76.33±11.68b 11.78±0.79b 0.98±0.04a 0.88±0.03a 3.10±0.44a 冰草A. cristatum
密花香薷Elsholtzia densa
节裂角茴香Hypecoum leptocarpum
ZM3 89.00±3.00ab 10.84±1.04b 1.47±0.19a 0.95±0.02a 2.67±0.65a 紫菀A. tataricus
苔草C. supina
二裂萎菱菜P. bifurca
ZMM 100.00±0.00a 22.37±3.94a 1.40±0.17a 0.91±0.06a 2.57±0.23a 冰草A. cristatum
苔草C. supina
臭蒿Artemisia hedinii
CM 100.00±0.00a 14.33±1.63b 1.31±0.13a 0.94±0.02a 2.64±0.20a 冰草A. cristatum
苔草C. supina
黄花棘豆O. ochrocephala
CG ZM1 33.33±6.17d 9.00±1.70a 1.04±0.13b 0.83±0.09a 2.27±0.14b 冰草A. cristatum
二裂萎菱菜P. bifurca
紫菀A. tataricus
ZM2 53.67±1.20c 7.04±0.76a 1.12±0.11b 0.86±0.04a 2.55±0.40ab 冰草A. cristatum
兰石草L. tibetica
多裂萎菱菜P. anserina
ZM3 85.00±4.04b 6.69±0.35a 1.10±0.12b 0.92±0.02a 2.58±0.24ab 冰草A. cristatum
二裂萎菱菜P. bifurca
黄花棘豆O. ochrocephala
ZMM 100.00±0.00a 8.73±0.12a 1.56±0.07a 0.96±0.00a 3.52±0.29a 冰草A. cristatum
鹅绒萎菱菜P. anserina
苔草C. supina
CM 100.00±0.00a 8.73±0.12a 1.56±0.07a 0.96±0.00a 3.52±0.29a 冰草A. cristatum
苔草C. supina
多裂萎菱菜P. mulatifida
GSG ZM1 43.67±22.51b 11.8±1.38a 0.65±0.33b 0.64±0.32a 2.46±1.26a 黄花棘豆O. ochrocephala
冰草A. cristatum
密花香薷E. densa
ZM2 82.67±9.13a 8.67±0.96a 1.03±0.08ab 0.93±0.01a 3.50±0.59a 冰草A. cristatum
黄花棘豆O. ochrocephala
二裂萎菱菜P. bifurca
ZM3 93.33±3.67a 9.70±2.00a 1.09±0.04ab 0.95±0.00a 3.77±0.30a 冰草A. cristatum
黄花棘豆O. ochrocephala
兰石草L.tibetica
ZMM 99.33±0.67a 9.14±0.80a 1.48±0.04a 0.90±0.01a 2.28±0.19a 冰草A. cristatum
黄花棘豆O. ochrocephala
苔草C. supina
CM 100.00±0.00a 11.43±0.39a 1.52±0.03a 0.95±0.01a 2.62±0.17a 苔草C. supina
冰草A. cristatum
黄花棘豆O. ochrocephala
RG:划区轮牧, Rotational grazing; PG:禁牧, Grazing prohibition; CG:连续放牧, Continuous grazing; GSG:生长季休牧, Growing season grazing; ZM1:一年鼠丘, One year zokor mounds; ZM2:二年鼠丘, Two years zokor mounds; ZM3:三年鼠丘, Three years zokor mounds; ZMM:多年鼠丘, Multi-year zokor mounds; CM:对照, Control meadow; 不同字母表示同一放牧方式下, 不同年限鼠丘之间的差异性(P < 0.05)
2.3 不同放牧管理模式对鼠丘群落物种多样性的影响

以α多样性(Shannon-wiener指数H)为物种多样性指标, 分析不同放牧管理模式对不同年限鼠丘物种多样性的影响, 发现CG下物种H最低, GSG下H较高;RG下ZM1和ZM2物种H显著高于CM物种H(P < 0.05);GSG下ZM1和ZM3物种H显著高于CM物种H(P < 0.05);PG下ZM3物种H高于CM物种H, 但差异不显著(P>0.05);CG下CM物种H显著高于ZM1、ZM2、ZM3和CM物种H(P < 0.05)(图 1)。

图 1 不同放牧管理模式在不同鼠丘年限下Shannon-wiener指数的比较 Fig. 1 Comparison of Shannon-Weiner index in different grazing management modes under different zokor mounds year RG:划区轮牧, Rotational grazing; PG:禁牧, Grazing prohibition; CG:连续放牧, Continuous grazing; GSG:生长季休牧, Growing season grazing; ZM1:一年鼠丘, One year zokor mounds; ZM2:二年鼠丘, Two years zokor mounds; ZM3:三年鼠丘, Three years zokor mounds; ZMM:多年鼠丘, Multi-year zokor mounds; CM:对照, Control meadow; 不同字母表示同一放牧方式下, 不同年限鼠丘之间的差异性(P < 0.05)
2.4 不同放牧管理模式下鼠丘群落生物量分配格局

PG、CG和GSG下, 随着鼠丘年限的增加地上生物量也增多, 但CG和GSG下CM地上生物量最高, PG下ZMM地上生物量最高, RG下地上生物量大小为:ZM1 < ZM3 < ZM2 < ZMM < CM, CG和RG下各年限鼠丘的地上生物量均高于PG和GSG下各鼠丘的生物量;所有样地的地下生物量在10 cm土层分布最多, 随着土层深度增加地下生物量减少, 对照样地中地下总生物量最高, 鼢鼠的扰动减少了地下生物量, 但各放牧管理模式下大小为:GSG>PG>CG>RG。GSG下CM地下总生物量接近3500 g/m2, 而PG、RG和GSG下ZM2地下总生物量均为最小, CG下ZM1地下总生物量最小(图 2)。

图 2 不同放牧管理模式在不同鼠丘年限下生物量分配 Fig. 2 Biodiversity allocation in different grazing management modes under different zokor mounds year RG:划区轮牧, Rotational grazing; PG:禁牧, Grazing prohibition; CG:连续放牧, Continuous grazing; GSG:生长季休牧, Growing season grazing; ZM1:一年鼠丘, One year zokor mounds; ZM2:二年鼠丘, Two years zokor mounds; ZM3:三年鼠丘, Three years zokor mounds; ZMM:多年鼠丘, Multi-year zokor mounds; CM:对照, Control meadow; AGB:地上生物量, Above ground biomass
2.5 不同放牧制度下影响鼠丘群落演替的因子分析

主成分分析结果表明, RG下2个公因子对样本方差的贡献率为74.91%, 即2个公因子包含了原始数据信息总量的74.91%, 可以用这2个公因子代表原来的8个指标, 第一主成分上, 具有高载荷因子的指标有地上生物量(0.939)、盖度(0.889), 第二主成分上, 高载荷因子的指标有O(0.874)、重要值(-0.829), 各指标在主成分中的载荷量, 代表与主成分关系的密切程度, 值越大则与主成分关系越密切, 越能代表鼠丘群落演替的变化;PG下3个公因子对样本方差的贡献率为77.23%, 第一主成分上, 具有高载荷因子的指标有地上生物量(0.934)、盖度(0.915);第二主成分上, 高载荷因子有O(0.914)、H(0.906), 第三主成分上, 高载荷因子有Jsw(0.889);CG下2个公因子对样本方差的贡献率为75.71%, 第一主成分上, 具有高载荷因子的指标有O(0.925)、H(0.871)、盖度(0.812);第二主成分上, 高载荷因子有高度(0.858)、重要值(0.860);GSG下3个公因子对样本方差的贡献率为85.99%, 第一主成分上, 具有高载荷因子的指标有H(0.974)、O(0.849), 第二主成分上, 高载荷因子有盖度(0.914)、地上生物量(0.901), 第三主成分上, 高载荷因子有高度(0.890)(表 3)。

表 3 旋转成分矩阵 Table 3 Rotation component matrix
指标
Index
RG PG CG GSG
PC1 PC2 PC1 PC2 PC3 PC1 PC2 PC1 PC2 PC3
地上生物量AGB 0.939 0.036 0.934 0.071 0.081 0.830 0.391 -0.303 0.901 0.083
丰富度指数Richness index 0.344 0.874 -0.032 0.914 -0.376 0.925 -0.152 0.849 -0.420 -0.207
盖度Coverage 0.889 -0.090 0.915 0.135 0.002 0.775 0.354 -0.013 0.914 -0.217
Shannon-wiener指数 -0.520 0.753 0.262 0.906 0.240 0.871 0.440 0.974 -0.166 -0.047
均匀度指数Uniformity index -0.792 0.355 0.040 0.138 0.889 0.390 0.761 0.800 -0.171 0.418
高度Height 0.613 -0.055 0.774 0.062 0.263 -0.148 0.858 0.208 -0.052 0.890
重要值Important value 0.001 -0.829 0.294 -0.384 0.716 0.305 0.860 -0.455 0.837 0.238
地下生物量UGB 0.513 0.714 0.379 -0.208 0.313 0.579 0.029 -0.415 0.107 0.580
特征值Eigenvalues 3.328 2.665 2.611 1.893 1.675 3.507 2.550 2.825 2.594 1.459
方差贡献率VCR/% 41.596 33.306 32.633 23.660 20.938 43.840 31.872 35.319 32.429 18.239
累计贡献率CCR/% 41.596 74.906 32.633 56.293 77.232 43.840 75.712 35.319 67.748 85.987
RG:划区轮牧, Rotational grazing; PG:禁牧, Grazing prohibition; CG:连续放牧, Continuous grazing; GSG:生长季休牧, Growing season grazing; PC1:第一主成分, First principal component; PC2:第二主成分, Second principal component; PC3:第三主成分, Third principal component; AGB:地上生物量, Above ground biomass; UGB:地下生物量, Under ground biomass; VCR:方差贡献率, Variance contribution rate; CCR:累计贡献率, Cumulative contribution rate
3 讨论 3.1 放牧管理模式对草地群落特征和地上地下生物量的影响

由于生物和非生物因素的综合作用, 造成草地生态系统破坏愈加严重, 而不合理的放牧制度是造成这一现象的主要原因, 其不仅降低了草地生态系统的生产力和恢复力, 也降低畜产品产量, 从而影响当地牧民的基本生存[24]。合理有效的放牧方式, 能促进草地生态系统的结构和功能向有利方向发展, 尤其对群落的基本特征(盖度、高度和丰富度等)和群落多样性等会造成直接影响[25]。中度干扰假说认为, 中等强度的干扰会阻碍生长旺盛植物对资源的利用, 有利于生长能力较弱的物种对资源的吸收利用, 使群落物种多样性增高[26]。有研究发现, 中等程度的干扰下物种丰富度指数和Shannon-wiener指数高于禁牧区[27], 长期禁牧会增加家畜喜食植物的竞争能力, 阻碍群落中某些植物幼苗的生长, 导致群落不能充分进行物种更新, 进而降低群落物种多样性[28]。有研究表明, 禁牧和全生长季休牧能增加植物密度、生物量以及改善草地植被组成[1]。本研究中, GSG、PG和RG下地下总生物量都高于CG(表 2), 符合中度干扰假说。原因可能由于家畜啃食和践踏的强度较大, 使植物自身不能进行自我调节导致浅根系植物直接死亡, 群落结构单一, 从而导致地下生物量降低。

3.2 啮齿动物对草地群落特征和生物量的影响

啮齿动物在草地生态系统中扮演着重要的角色[29]。啮食优良牧草, 挖掘活动损失牧草, 挖洞推土影响土壤肥力, 降低植被盖度, 促使土壤水分蒸发, 改变植被成分, 引起群落演替等, 同时也影响到高寒草地植物的生产力和草地地上植被组成、群落结构的复杂性等[30]。研究发现, 物种丰富度随高原鼢鼠新土丘数的增加呈先增加后降低的趋势, 也会随鼢鼠种群密度增加呈降低趋势, 出现这种差异的原因可能与其采用空间序列代替时间序列的方法有关[31]。有研究表明, 高原鼢鼠种群增长的过度放牧样地, 可以通过适度禁牧等休牧措施来提高草地群落多样性和生产力[32], 以及群落盖度随鼠丘密度的增加显著降低[18], 且植物高度随鼠丘年限的增加越来越大[33], 这些与本研究结果相似(图 1图 2)。原因可能由于过度放牧会使草地退化越加严重, 进而使草地高原鼢鼠种群数量增加, 而适度放牧会改善草地退化, 高原鼢鼠种群数量也会降低, 以及随鼠丘年限增加, 鼠丘上土壤水分、养分等趋于稳定并逐渐增大, 对植物的生长提供良好的条件。

3.3 放牧管理模式对鼠丘群落演替的影响

青藏高原高寒草甸对放牧、鼢鼠干扰等因素十分敏感, 微小波动都会使草地生态系统产生强烈的响应, 进而导致群落特征、植物多样性和生产力都发生很大的改变[34]。生态位理论认为物种能够共存是由于物种在生态位中的位置不同、生态位重叠较小和生态位分化等[35]。中性理论认为扩散限制是决定群落结构的重要因素, 且也受到空间尺度、环境因子、物种互利和密度等因子的影响[36]。Turnbull等[37]资源生态位模型研究表明, 生态位差异在物种共存时会在物种间产生互补效应, 但其不利于物种的稳定共存。有研究表明, 禁牧和自由放牧下物种重要值、多样性指数和优势度指数均小于季节性放牧地[19], 与本研究结果相同(图 1表 2), 说明放牧强度对鼠丘群落的恢复有一定的影响, 由于高强度放牧下, 家畜长时间采食和践踏导致植物的根系受到一定的影响。有研究表明, 禁牧和划区轮牧对快速生长的一年生植物有利, 对多年生植物可能会产生一定的限制作用[38], 这很好的解释了本研究中RG和PG下ZM3物种重要值最高的原因。也有研究发现, 合理的放牧制度能够减弱优势种的竞争力, 加速弱势物种的存活, 提高草地物种多样性[5]

不同放牧制度下, 植物的生活型和繁殖方式也不同, 资源比率/异质性假说认为, 异质性环境中两种或多种限制性资源比率与单一限制性资源相比, 其生态位分化更丰富[39]。放牧制度的差异性决定不同的物种组成和更替, 是构成群落演替的基础[40]。物种丰富度指数指群落内植物种类的多少, Shannon-wiener指数是反映群落物种均匀度和丰富度的综合指标, 均匀度是指群落中每一物种数量的分配比例[41]。因此, Shannon-wiener指数、均匀度指数和丰富度指数能较全面地反映物种的多样性和植物个体在种群中分布的均匀程度。本研究主成分分析表明, RG和PG下, 地上生物量与第一主成分的关系最密切, CG下, 第一主成分上物种丰富度的载荷最高, GSG下, 第一主成分上物种Shannon-wiener指数的载荷最高(表 3), 说明不同放牧制度会导致鼠丘的群落演替因子的不同, 可能由于不同放牧管理模式下土壤水分、热量等的利用和维持不同, 供给植被生长、发育所吸收和利用的养分也不同, 因此, 鼠丘上植被的生长恢复也会表现出一定的差异性[27]

综上所述, 不同放牧管理模式下影响鼠丘群落演替的因素不同, RG下鼢鼠扰动增加了物种Shannon-wiener指数, 而CG下鼢鼠扰动降低了物种Shannon-wiener指数和地上生物量;ZM2物种总盖度、丰富度指数和Shannon-wiener指数, 均表现为GSG高于RG、PG和CG, 且CG下均最小;PG和GSG下物种平均高度、地上生物量均高于CG和RG。因此, 合理的放牧管理模式可以有效改善高寒草甸草地的群落特征, 从而优化高原鼢鼠的管理策略。

参考文献
[1]
彭祺, 王宁. 不同放牧制度对草地植被的影响. 农业科学研究, 2005, 26(1): 27-30. DOI:10.3969/j.issn.1673-0747.2005.01.007
[2]
Kotzé E, Sandhage-Hofmann A, Meinel J A, du Preez C C, Amelung W. Rangeland management impacts on the properties of clayey soils along grazing gradients in the semi-arid grassland biome of South Africa. Journal of Arid Environments, 2013, 97: 220-229. DOI:10.1016/j.jaridenv.2013.07.004
[3]
Miao F H, Guo Z G, Xue R, Wang X Z, Shen Y Y. Effects of grazing and precipitation on herbage biomass, herbage nutritive value, and yak performance in an alpine meadow on the Qinghai-Tibetan Plateau. PLoS One, 2015, 10(6): e0127275. DOI:10.1371/journal.pone.0127275
[4]
侯扶江, 杨中艺. 放牧对草地的作用. 生态学报, 2006, 26(1): 244-264. DOI:10.3321/j.issn:1000-0933.2006.01.031
[5]
McNaughton S J. Grazing as an optimization process:grass-ungulate relationships in the Serengeti. The American Naturalist, 1979, 113(5): 691-703. DOI:10.1086/283426
[6]
Cui X F, Graf H F. Recent land cover changes on the Tibetan Plateau:a review. Climatic Change, 2009, 94(1/2): 47-61.
[7]
蒙荣, 赛音吉雅, 张小芬. 不同放牧制度大针茅草原群落特征的比较. 内蒙古草业, 2003, 15(2): 10-12.
[8]
Meissner R A, Facelli J M. Effects of sheep exclusion on the soil seed bank and annual vegetation in chenopod shrublands of South Australia. Journal of Arid Environments, 1999, 42(2): 117-128. DOI:10.1006/jare.1999.0515
[9]
Li W, Liu Y Z, Wang J L, Shi S L, Cao W X. Six years of grazing exclusion is the optimum duration in the alpine meadow-steppe of the north-eastern Qinghai-Tibetan Plateau. Scientific Reports, 2018, 8(1): 17269. DOI:10.1038/s41598-018-35273-y
[10]
苏军虎, 南志标, 纪维红. 家畜放牧对草地啮齿动物影响的研究进展. 草业学报, 2016, 25(11): 136-148. DOI:10.11686/cyxb2015587
[11]
施大钊, 钟文勤. 2000年我国草原鼠害发生状况及防治对策. 草地学报, 2001, 9(4): 248-252.
[12]
Koide R T, Huenneke L F, Mooney H A. Gopher mound soil reduces growth and affects ion uptake of two annual grassland species. Oecologia, 1987, 72(2): 284-290. DOI:10.1007/BF00379280
[13]
Rebollo S, Milchunas D G, Stapp P, Augustine D J, Derner J D. Disproportionate effects of non-colonial small herbivores on structure and diversity of grassland dominated by large herbivores. Oikos, 2013, 122(12): 1757-1767. DOI:10.1111/j.1600-0706.2013.00403.x
[14]
Pickett S T A. Non-equilibrium coexistence of plants. Bulletin of the Torrey Botanical Club, 1980, 107(2): 238-248. DOI:10.2307/2484227
[15]
Schaffer W M, Fox G A. Plant strategies and the dynamics and structure of plant communities. Bulletin of Mathematical Biology, 1989, 51(3): 409-411.
[16]
刘双双.家畜与高原鼢鼠对高寒草甸植物多样性和生产力的影响[D].兰州: 兰州大学, 2018. http://cdmd.cnki.com.cn/Article/CDMD-10730-1018979798.htm
[17]
江小雷, 张卫国, 杨振宇, 杜国祯. 不同演替阶段鼢鼠土丘群落植物多样性变化研究. 应用生态学报, 2004, 15(5): 814-818. DOI:10.3321/j.issn:1001-9332.2004.05.017
[18]
于海玲, 樊江文, 李愈哲, 史文娇. 高原鼢鼠干扰对三江源区高寒草甸群落特征的影响. 应用生态学报, 2018, 29(6): 1902-1910.
[19]
王玉琴, 鲍根生, 王宏生, 曾辉, 李继鹏. 不同放牧制度对高原鼢鼠活动区高寒草地群落结构影响. 草地学报, 2018, 26(1): 134-141.
[20]
姬万忠, 王庆华. 补播对天祝高寒退化草地植被和土壤理化性质的影响. 草业科学, 2016, 33(5): 886-890.
[21]
杨莹博, 辛小娟, 艾得协措, 杨雪, 王绍美, 储诚进, 任正炜, 王刚. 鼢鼠土丘植被恢复演替过程中的物种多样性变化. 草业学报, 2010, 19(1): 14-20.
[22]
鲍根生, 王宏生, 曾辉. 有机硅助剂对"狼毒净"防除狼毒效果及草地植物群落多样性的影响. 草业科学, 2015, 32(2): 263-268.
[23]
陈琳, 赵陟峰, 赵廷宁, 王娟. 基于因子与聚类分析的不同林分土壤理化性质评价. 水土保持研究, 2011, 18(5): 191-196.
[24]
朱桂林, 山仑, 卫智军, 韩国栋. 放牧制度对短花针茅群落植物生长的影响研究. 中国生态农业学报, 2004, 12(4): 181-183.
[25]
杨利民, 韩梅, 周广胜, 王玉辉. 草地群落物种多样性维持机制的研究Ⅲ.物种分布格局. 吉林农业大学学报, 2002, 24(1): 58-61. DOI:10.3969/j.issn.1000-5684.2002.01.013
[26]
Connell J H. Diversity in tropical rain forests and coral reefs. Science, 1978, 199(4335): 1302-1310. DOI:10.1126/science.199.4335.1302
[27]
Proulx M, Mazumder A. Reversal of grazing impact on plant species richness in nutrient-poor vs. nutrient- rich ecosystems. Ecology, 1998, 79(8): 2581-2592. DOI:10.1890/0012-9658(1998)079[2581:ROGIOP]2.0.CO;2
[28]
闫玉春, 唐海萍. 围栏禁牧对内蒙古典型草原群落特征的影响. 西北植物学报, 2007, 27(6): 1225-1232. DOI:10.3321/j.issn:1000-4025.2007.06.026
[29]
Qin Y, Chen J J, Yi S H. Plateau pikas burrowing activity accelerates ecosystem carbon emission from alpine grassland on the Qinghai-Tibetan Plateau. Ecological Engineering, 2015, 84: 287-291. DOI:10.1016/j.ecoleng.2015.09.012
[30]
Klein J A, Harte J, Zhao X Q. Experimental warming, not grazing, decreases rangeland quality on the Tibetan Plateau. Ecological Applications, 2007, 17(2): 541-557. DOI:10.1890/05-0685
[31]
黄彬.高原鼢鼠种群密度与植被、土壤因素的关系[D].兰州: 兰州大学, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10730-1014305427.htm
[32]
张卫国, 江小蕾, 王树茂, 杨振宇. 鼢鼠的造丘活动及不同休牧方式对草地植被生产力的影响. 西北植物学报, 2004, 24(10): 1882-1887. DOI:10.3321/j.issn:1000-4025.2004.10.020
[33]
张燕堃.高寒草甸鼢鼠土丘年限对植物功能特征及生产力的影响[D].兰州: 兰州大学, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10730-1014305425.htm
[34]
赵新全, 马玉寿, 周华坤. 三江源区退化生态系统恢复及可持续管理化. 北京: 科学出版社, 2011.
[35]
Mahdi A, Law R, Willis A J. Large niche overlaps among coexisting plant species in a limestone grassland community. Journal of Ecology, 1989, 77(2): 386-400. DOI:10.2307/2260757
[36]
Hubbell S P. The Unified Neutral Theory of Biodiversity and Biogeography. Princeton: Princeton University Press, 2001.
[37]
Turnbull L A, Levine J M, Loreau M, Hector A. Coexistence, niches and biodiversity effects on ecosystem functioning. Ecology Letters, 2013, 16(S1): 116-127.
[38]
双全, 邢旗, 金玉, 松梅, 巴格那. 草甸草原主要植物种群动态特征对放牧制度的响应. 畜牧与饲料科学, 2010, 31(6/7): 298-300.
[39]
Tilman D. Resource Competition and Community Structure. Princeton: Princeton University Press, 1982.
[40]
汪诗平, 李永宏. 内蒙古典型草原退化机理的研究. 应用生态学报, 1999, 10(4): 437-441. DOI:10.3321/j.issn:1001-9332.1999.04.015
[41]
Haynes R J. The decomposition process: mineralization, immobilization, humus formation, and degradation//Haynes R J, ed. Mineral Nitrogen in the Plant-Soil System. London: Academic Press, 1986: 52-126.