用微信扫码二维码

分享至好友和朋友圈

WE ARE COMMITTED TO REPORTING THE LATEST FORESTRY ACADEMIC ACHIEVEMENTS

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

四川栗子坪国家级自然保护区大熊猫主食竹形态特征及其立地环境初探

袁圆 杨毅 代勤龙 龙珏洁 王佳 洪明生

袁圆, 杨毅, 代勤龙, 等. 四川栗子坪国家级自然保护区大熊猫主食竹形态特征及其立地环境初探[J]. 四川林业科技, 2021, 42(2): 8−14 doi: 10.12172/202010070002
引用本文: 袁圆, 杨毅, 代勤龙, 等. 四川栗子坪国家级自然保护区大熊猫主食竹形态特征及其立地环境初探[J]. 四川林业科技, 2021, 42(2): 8−14 doi: 10.12172/202010070002
Yuan Y, Yang Y, Dai Q L, et al. Study on morphological characteristics of giant panda’s staple food bamboo and its environmental features in Liziping National Nature Reserve, Sichuan[J]. Journal of Sichuan Forestry Science and Technology, 2021, 42(2): 8−14 doi: 10.12172/202010070002
Citation: Yuan Y, Yang Y, Dai Q L, et al. Study on morphological characteristics of giant panda’s staple food bamboo and its environmental features in Liziping National Nature Reserve, Sichuan[J]. Journal of Sichuan Forestry Science and Technology, 2021, 42(2): 8−14 doi: 10.12172/202010070002

四川栗子坪国家级自然保护区大熊猫主食竹形态特征及其立地环境初探


doi: 10.12172/202010070002
详细信息
    作者简介:

    袁圆(1995—),女,硕士,1453378516@qq.com

    杨毅(1981—),男,工程师,硕士,651843296@qq.com #共同第一作者

    通讯作者: mingshenghong@cqu.edu.cn
  • 基金项目:  2019年国家重点生态功能区转移支付禁止开发区补助资金项目;国家自然科学基金(31900337);西华师范大学基本科研业务费(18B025,18Q040);四川省科技计划项目资助(2019YFG0460)

Study on Morphological Characteristics of Giant Panda’s Staple Food Bamboo and Its Environmental Features in Liziping National Nature Reserve, Sichuan

More Information
    Corresponding author: mingshenghong@cqu.edu.cn
  • 摘要: 大熊猫在进化过程中形成了一系列策略适应竹子食物,如在不同季节可采食所处不同海拔的竹类,而不同海拔的竹子具有特定的形态特征和环境特点。本研究以四川栗子坪国家级自然保护区为研究区域,通过野外调查揭示小相岭山系大熊猫种群不同主食竹的形态特征和立地环境差异。结果表明,峨热竹、石棉玉山竹和丰实箭竹的基经和高度具有显著差异,低海拔的丰实箭竹具有最大的基经和高度,而高海拔峨热竹的基经和高度最小。不同竹种的立地环境不同,主要表现在海拔和乔木郁闭度上,不同竹种所需求的最适宜生长的乔木郁闭度不同。本研究结果可为该区域大熊猫栖息地恢复和管理,特别是大熊猫主食竹的精细化管理提供参考依据。
    #共同第一作者
  • 图  1  不同竹种生境变量判别的散点图。

    DF1,判别函数1,(F(2, 47)=456.20,P<0.001);DF2,判别函数2,(F(2, 47)=54.74,P<0.001)

    Fig.  1  Discriminant scatter diagram of different bamboo species' habitat variables

    DF1, Discriminant function 1 (F(2, 47)=456.20,P<0.001); DF2, Discriminant function 2 (F(2, 47)=54.74,P<0.001)

    表  1  各变量的定义和测量

    Tab.  1  Definition and measurement of variables

    变量(单位)Variables(Units)描述与测量Definition and Measurement
    海拔(m)生境样方(20 m×20 m)中心点的海拔
    坡度(°)生境样方(20 m×20 m)中心点的坡度
    水源距离(m)生境样方(20 m×20 m)中心点到可看见或听见水流的距离
    乔木高(m)生境样方(20 m×20 m)内高于5 m的乔木的平均高度
    乔木胸径(cm)生境样方(20 m×20 m)内高于5 m的乔木平均胸径
    乔木数量生境样方(20 m×20 m)内高于5 m的乔木数
    乔木郁闭度(%)生境样方(20 m×20 m)内高于5 m的乔木郁闭度(%)
    灌木高(m)生境样方(20 m×20 m)内高≤5 m的灌木的平均高度
    灌木盖度(%)生境样方(20 m×20 m)内高≤5 m的灌木平均胸径
    灌木数量生境样方(20 m×20 m)内高≤5 m的灌木数
    竹子盖度(%)3个竹子样方(1 m×1 m)平均竹子盖度(%)
    活竹总数3个竹子样方(1 m×1 m)平均活竹子数量
    竹子死亡数3个竹子样方(1 m×1 m)平均死亡竹子数量
    一年生竹数量3个竹子样方(1 m×1 m)平均一年生竹子数量
    一年生竹高(cm)3个竹子样方(1 m×1 m)平均一年生竹子高
    一年生竹基径(mm)3个竹子样方(1 m×1 m)平均一年生竹子基经
    二年生竹数量3个竹子样方(1 m×1 m)平均二年生竹子数量
    二年生竹高(cm)3个竹子样方(1 m×1 m)平均二年生竹子高
    二年生竹基径(mm)3个竹子样方(1 m×1 m)平均二年生竹子基经
    多年生竹数量3个竹子样方(1 m×1 m)平均多年生竹子数量
    多年生竹高(cm)3个竹子样方(1 m×1 m)平均多年生竹子高
    多年生竹基径(mm)3个竹子样方(1 m×1 m)平均多年生竹子基经
    下载: 导出CSV

    表  2  不同竹种特征变量之间的差异性

    Tab.  2  Differences of bamboo characteristics among different bamboo species

    变量Variables丰实箭竹F. ferax石棉玉山竹Y. lineolata峨热竹A. spanostachyaANOVA (F, P)
    竹子盖度52.77±20.9252.75±15.7749.14±8.380.30, 0.75
    活竹总数34.83±7.3233.59±5.4932.53±4.580.63, 0.54
    竹子死亡数5.65±1.62a6.80±0.79b5.55±0.90a6.22,<0.01
    一年生竹数量7.90±1.927.71±1.177.14±1.421.12, 0.34
    一年生竹高347.09±110.57a306.55±46.97a175.42±31.20b26.95,<0.001
    一年生竹基径20.24±3.06a17.04±1.42b9.24±1.28c125.98,<0.001
    二年生竹数量7.81±1.53a10.00±2.06b8.06±1.02a9.42,<0.001
    二年生竹高321.00±94.35a276.18±51.44b187.70±31.42c18.74,<0.001
    二年生竹基径17.84±3.48a18.16±1.77a10.24±1.21b62.39,<0.001
    多年生竹数量19.13±5.82a15.88±3.16b17.33±3.21ab2.45, 0.097
    多年生竹高371.10±97.05a310.89±56.24b193.98±24.97c31.36,<0.001
    多年生竹基径21.14±3.64a19.14±1.26b10.80±1.27c94.73,<0.001
    下载: 导出CSV

    表  3  不同竹种生境变量之间的差异性

    Tab.  3  Differences of habitat variables among different bamboo species

    变量Variables丰实箭竹F. ferax石棉玉山竹Y. lineolata峨热竹A. spanostachyaANOVA (F, P)
    海拔2064.06±165.10a2561.76±126.61b2896.35±128.31c145.81,<0.001
    坡度35.94±24.57a18.47±16.17b18.82±16.06b4.40, 0.02
    水源距离63.69±81.79a197.12±142.78b133.82±110.76ab5.54, 0.01
    乔木高20.56±8.0422.29±7.7018.24±4.931.44, 0.25
    乔木胸径20.31±10.16a24.65±10.99a38.35±11.69b12.27,<0.001
    乔木数量36.56±14.74a30.82±19.60a13.94±5.86b10.89,<0.001
    乔木郁闭度46.19±18.33a48.06±15.20a26.12±10.62b11.13,<0.001
    灌木高3.53±0.783.03±0.923.51±0.931.74, 0.19
    灌木盖度22.31±10.5217.47±14.6519.82±12.920.59, 0.56
    灌木数量20.63±10.26a14.12±12.94ab9.29±3.98b5.54, 0.01
    下载: 导出CSV

    表  4  竹子生境变量标准化的判别函数

    Tab.  4  Standardized discriminant functions of bamboo habitat variables

    变量Variables标准系数Standardized coefficients
    判别函数1
    Discriminant function 1
    判别函数2
    Discriminant function 2
    海拔0.9780.211
    乔木郁闭度−0.214 0.977
    下载: 导出CSV

    表  5  竹子生境变量的判别预测

    Tab.  5  Discriminant prediction of bamboo habitat variables

    竹种
    Bamboo species
    预测结果Predicted results
    丰实箭竹
    F. ferax
    石棉玉山竹
    Y. lineolata
    峨热竹
    A. spanostachya
    丰实箭竹F. ferax87.50 12.500.00
    石棉玉山竹
    Y. lineolata
    5.9088.205.90
    峨热竹
    A. spanostachya
    0.00 0.00100.00
    下载: 导出CSV
  • [1] 胡锦矗, Schaller GB, 潘文石, 等. 卧龙的大熊猫[M]. 四川科学技术出版社, 1985, 成都.
    [2] 邱占祥,祁国琴. 云南禄丰晚中新世的大熊猫祖先化石[J]. 古脊椎动物学报,1989(3):153−169.
    [3] Zhao S, Zheng P, Dong S, et al. Whole-genome sequencing of giant pandas provides insights into demographic history and local adaptation[J]. Nature Genetics, 2013, 45(1): 67−U99. doi: 10.1038/ng.2494
    [4] Jin C, Ciochon R L, Dong W, et al. The first skull of the earliest giant panda[J]. Proceedings of the National Academy of Sciences, 2007, 104(26): 10932−10937. doi: 10.1073/pnas.0704198104
    [5] Wei F, Feng Z, Wang Z, et al. Feeding strategy and resource partitioning between giant and red pandas[J]. Mammalia, 2009, 63(4): 417−429.
    [6] Wei F W, Hu Y B, Yan L et al. Giant Pandas Are Not an Evolutionary cul-de-sac: Evidence from Multidisciplinary Research[J]. Molecular Biology & Evolution, 2015(1): 4−12.
    [7] Wei F W, Swaisgood R R, Hu Y B, et al. Progress in the ecology and conservation of giant pandas[J]. Conservation. Biology, 2015(29): 1497−1507.
    [8] Hu Y B, Wu Q, Ma S, et al. Comparative genomics reveals convergent evolution between the bamboo-eating giant and red pandas[J]. Proceedings of the National Academy of Sciences, 2017, 114(5): 1081−1086. doi: 10.1073/pnas.1613870114
    [9] Hong M S, Yuan S B, Yang Z S et al. Comparison of microhabitat selection and trace abundance of giant pandas between primary and secondary forests in Liziping Natur Reserve, China: Effects of selective loggingMammalian Biology, 2015, 80, 373–379.
    [10] Hong, M S, Wei W, Yang Z S, et al. Effects of timber harvesting on Arundinaria spanostachya bamboo and feeding-site selection by giant pandas in liziping nature reserve, China[J]. Forest Ecology and Management, 2016, 373: 74−80. doi: 10.1016/j.foreco.2016.04.039
    [11] 潘文石, 吕植, 朱小健, 等. 2001. 继续生存的机会. 北京大学出版社, 2001, 北京.
    [12] Wu Q, Wang X, Ding Y, et al. Seasonal variation in nutrient utilization shapes gut microbiome structure and function in wild giant pandas[J]. Proceedings of the Royal Society B Biological Sciences, 2017, 284(1862): 20170955. doi: 10.1098/rspb.2017.0955
    [13] 魏辅文,冯祚建,王祖望. 相岭山系大熊猫和小熊猫对生境的选择[J]. 动物学报,1999,45(1):57−63. doi: 10.3969/j.issn.1674-5507.1999.01.009
    [14] 四川省林业厅. 2015. 四川的大熊猫: 四川省第四次大熊猫调查报告. 四川科技出版社, 2015, 成都.
    [15] 洪明生,王继成,杨旭煜,et al. 原始林与次生林中大熊猫微生境结构的比较[J]. 西华师范大学学报(自然科学版),2012(04):356−361.
    [16] Xie J M, Hong M S. Effects of ecological factors on growth of Arundinaria spanostachya shoots in Liziping National Nature Reserve, China[J]. Global Ecology and Conservation, 2020, 23: e01121. doi: 10.1016/j.gecco.2020.e01121
    [17] Hong M S, Wei W, Zhou H, et al. Creative conservation in China: releasing captive giant pandas into the wild[J]. Environmental Science and Pollution Research, 2019, 26: 31548−31549. doi: 10.1007/s11356-019-06384-4
    [18] Zhang L Y, Gan X H, Hou Z Y, et al. Grazing by wild giant pandas does not affect the regeneration of Arundinaria spanostachya[J]. Journal of Forestry Research, 2018: 1−8.
    [19] 周宏,袁施彬,杨志松,等. 四川栗子坪自然保护区夏季大熊猫食性与主食竹生物量的关系[J]. 兽类学报,2014,34(001):93−99.
    [20] Lei M W, Yuan S B, Yang Z S, et al. Comparison of microhabitats and foraging strategies between the captive-born Zhangxiang and wild giant pandas: implications for future reintroduction[J]. Environmental Science and Pollution Research, 2015, 22: 15089−15096. doi: 10.1007/s11356-015-4720-3
    [21] He K, Dai Q, Foss-Grant A, et al. Movement and activity of reintroduced giant pandas[J]. Ursus, 2019, 29(2): 163. doi: 10.2192/URSUS-D-17-00030.1
    [22] Nie Y G, Zhang Z J, Raubenheimer D, et al. Obligate herbivory in an ancestrally carnivorous lineage: the giant panda and bamboo from the perspective of nutritional geometry[J]. Functional Ecology, 2014, 29: 26−34.
    [23] Wei W F, Nie Y G, Zhang Z J, et al. Hunting bamboo: Foraging patch selection and utilization by giant pandas and implications for conservation[J]. Biological Conservation, 2015, 186: 260−267. doi: 10.1016/j.biocon.2015.03.023
    [24] Zhang M C, Zhang Z Z, Li Z, et al. Giant panda foraging and movement patterns in response to bamboo shoot growth[J]. Environmental Science and Pollution Research, 2018, 25: 8636e8643.
    [25] Huang C, Wang L, Gong X, et al. Silicon fertilizer and biochar effects on plant and soil PhytOC concentration and soil PhytOC stability and fractionation in subtropical bamboo plantations[J]. Science of The Total Environment, 2020: 136846.
  • [1] 秦炜锐, 刘景怡, 宋心强, 付明霞, 潘翰, 程勇, 薛飞, 周材权, 杨彪.  大熊猫国家公园荥经片区大熊猫可食竹生物量模型 . 四川林业科技, 2024, 45(): 1-7. doi: 10.12172/202309070001
    [2] 杨钰华, 杨桂梅, 欧阳学军, 贺握权, 黄柳菁.  鼎湖山国家级自然保护区蕨类植物资源调查及分析 . 四川林业科技, 2023, 44(6): 55-61. doi: 10.12172/202302160001
    [3] 席健, 朱玉东, 何可, 杨志松, 黎大勇.  四川栗子坪国家级自然保护区林麝野化培训个体的生境选择 . 四川林业科技, 2023, 44(2): 51-57. doi: 10.12172/202211090003
    [4] 张峻, 班英, 卓玛姐, 言迫, 索郎夺尔基, 杨志松.  四川若尔盖湿地国家级自然保护区鸟类多样性及区系分析 . 四川林业科技, 2022, 43(5): 11-17. doi: 10.12172/202204240001
    [5] 古玉, 陈明华, 曾进科, 马晓平.  四川美姑大风顶国家级自然保护区竹类资源与大熊猫活动区域间关系研究 . 四川林业科技, 2022, 43(4): 24-30. doi: 10.12172/202109060001
    [6] 郑阳, 田春洋, 李建威, 李威, 王飞, 杨毅, 马文宝.  四川栗子坪国家级自然保护区杜鹃花属植物多样性与分布特点 . 四川林业科技, 2022, 43(6): 80-86. doi: 10.12172/202105250001
    [7] 田春洋, 洪明生, 龙珏洁, 谢建妹.  大熊猫主食竹叶围细菌多样性的季节性变化 . 四川林业科技, 2021, 42(5): 1-7. doi: 10.12172/202107050001
    [8] 刘涌, 李俊清, 刘正霄, 睢雅妮, 兰飞, 熊书圣, 何敏, 罗亮, 刘月, 贺品.  大熊猫主食竹产地土壤重金属评价 . 四川林业科技, 2021, 42(1): 25-28. doi: 10.12172/202008110002
    [9] 杨毅, 杨伟, 骆伟, 鲁碧耕, 刘秉昕, 田春洋, 李威, 张翼.  四川栗子坪国家级自然保护区兰科植物物种多样性及种群现状 . 四川林业科技, 2020, 41(4): 121-128. doi: 10.12172/202003030001
    [10] 李健威, 李玉霞, 张勘, 骆伟, 兰琦, 何流洋, 唐博, 杨志松.  四川栗子坪国家级自然保护区野生鸟兽的红外相机初步监测 . 四川林业科技, 2020, 41(3): 7-13. doi: 10.12172/202004090001
    [11] 彭科, 陈旭, 张宣, 房超, 杨彪.  四川白水河自然保护区大熊猫主食竹种群特征及其影响因素 . 四川林业科技, 2020, 41(5): 94-104. doi: 10.12172/202005090001
    [12] 白文科, 董鑫, 王光磊, 杨旭煜, 古晓东, 周材权.  人为采笋对大叶筇竹生长及大熊猫采食的影响 . 四川林业科技, 2019, 40(4): 80-83,96. doi: 10.16779/j.cnki.1003-5508.2019.04.016
    [13] 的么罗英, 王晓, 张晋东.  卧龙自然保护区人类活动与大熊猫(Ailuropoda melanoleuca)空间利用特征分析 . 四川林业科技, 2019, 40(5): 60-65. doi: 10.16779/j.cnki.1003-5508.2019.05.012
    [14] 曾诗卉, 赵志芳, 陈琪.  四川米仓山国家级自然保护区生态环境变化及其驱动力分析 . 四川林业科技, 2019, 40(2): 62-68. doi: 10.16779/j.cnki.1003-5508.2019.02.014
    [15] 雷开明, 程胡椒, 旷培刚, 冉江洪, 孙鸿鸥, 张跃, 麦浪, 曾涛.  九寨沟自然保护区大熊猫取食利用的主食竹特征 . 四川林业科技, 2019, 40(5): 66-69. doi: 10.16779/j.cnki.1003-5508.2019.05.013
    [16] 梁玉喜, 王戈, 王洪荣.  基于GIS大熊猫活动区域环境特征分析——以大熊猫国家公园(四川)为例 . 四川林业科技, 2018, 39(4): 98-101,114. doi: 10.16779/j.cnki.1003-5508.2018.04.024
    [17] 侯正扬, 田忠琼, 李文英, 毛维莉, 甘小洪, 欧拉提子.  野生大熊猫采食对峨热竹粗灰分和干物质含量的影响 . 四川林业科技, 2018, 39(1): 50-53. doi: 10.16779/j.cnki.1003-5508.2018.01.010
    [18] 方懿.  四川米亚罗省级自然保护区人类活动遥感监测分析 . 四川林业科技, 2017, 38(3): 36-39. doi: 10.16779/j.cnki.1003-5508.2017.03.007
    [19] 李谦, 金贵祥, 邓浩, 杨岚, 潘开文, 蒋晔.  汶川“5·12”特大地震对四川白河自然保护区大熊猫主食竹及栖息地的影响研究 . 四川林业科技, 2014, 35(5): 80-83,44. doi: 10.16779/j.cnki.1003-5508.2014.05.017
    [20] 刘明冲, 杨晓军, 张清宇, 管晓, 唐浩, 杨森.  卧龙自然保护区2013年大熊猫主食竹监测分析报告 . 四川林业科技, 2014, 35(4): 45-47. doi: 10.16779/j.cnki.1003-5508.2014.04.010
  • 加载中
  • 图(1) / 表(5)
    计量
    • 文章访问数:  425
    • HTML全文浏览量:  197
    • PDF下载量:  38
    • 被引次数: 0
    出版历程
    • 收稿日期:  2020-10-07
    • 网络出版日期:  2021-01-18
    • 刊出日期:  2021-04-20

    四川栗子坪国家级自然保护区大熊猫主食竹形态特征及其立地环境初探

    doi: 10.12172/202010070002
      作者简介:

      袁圆(1995—),女,硕士,1453378516@qq.com

      杨毅(1981—),男,工程师,硕士,651843296@qq.com #共同第一作者

      通讯作者: mingshenghong@cqu.edu.cn
    基金项目:  2019年国家重点生态功能区转移支付禁止开发区补助资金项目;国家自然科学基金(31900337);西华师范大学基本科研业务费(18B025,18Q040);四川省科技计划项目资助(2019YFG0460)

    摘要: 大熊猫在进化过程中形成了一系列策略适应竹子食物,如在不同季节可采食所处不同海拔的竹类,而不同海拔的竹子具有特定的形态特征和环境特点。本研究以四川栗子坪国家级自然保护区为研究区域,通过野外调查揭示小相岭山系大熊猫种群不同主食竹的形态特征和立地环境差异。结果表明,峨热竹、石棉玉山竹和丰实箭竹的基经和高度具有显著差异,低海拔的丰实箭竹具有最大的基经和高度,而高海拔峨热竹的基经和高度最小。不同竹种的立地环境不同,主要表现在海拔和乔木郁闭度上,不同竹种所需求的最适宜生长的乔木郁闭度不同。本研究结果可为该区域大熊猫栖息地恢复和管理,特别是大熊猫主食竹的精细化管理提供参考依据。

    #共同第一作者

    English Abstract

    • 在长达数百万年的进化历史中,大熊猫(Ailuropoda melanoleuca)食性由食肉逐渐特化为专食竹子,竹子含量超过了食谱组成的90%;然而大熊猫仍然保留了典型食肉目动物的消化道特征[1-3]化石材料研究表明,出现于上新世末更新世初的小种大熊猫(Ailuropoda microta)在食性上已经特化为以竹子为食,这表明大熊猫-竹子之间的协同进化可能已经达200万年以上[4]。在进化过程中,大熊猫形成了一系列适应高纤维、低营养含量竹子食物的适应性策略,如形成了一系列特殊的生理形态特征、生态行为、基因和基因组甚至是肠道微生物等来适应特化的竹子食物,然而并没有形成高效消化高纤维含量食物的类似于草食动物的复杂消化道[5-8]

      大熊猫对竹子资源的适应首先表现在觅食策略上,如季节性垂直迁移、对生境和采食点的选择等[1,9-10];同时,在不同山系的不同季节,大熊猫食性具有较大差异,如在秦岭山系,大熊猫在非竹笋季节均主要以竹叶为食[11-12];而在四川的小相岭山系竹叶占大熊猫粪便组成超过一半的时间主要出现在夏秋季节[13]。据全国第四次大熊猫调查,四川省大熊猫栖息地内共有大熊猫采食的竹类7属32种,其中冷箭竹(Bashania faberi)分布面积最大,达全省大熊猫取食竹总面积的20.92%[14]。小相岭山系大熊猫栖息地内有大熊猫采食竹类3属7种,其中面积最大的为峨热竹(Arundinaria spanostachya),占该山系大熊猫栖息地内大熊猫取食竹总面积的38.08%;其次为石棉玉山竹(Yushania lineolata)和丰实箭竹(Fargesia ferax),分别占该山系大熊猫栖息地内大熊猫采食竹的28.02%和12.48%[14]

      为了解小相岭山系大熊猫主食竹形态特征和立地环境状况,本研究以四川栗子坪国家级自然保护区为研究区域,通过野外调查大熊猫主食竹的主要特征及其生境,了解该区域大熊猫不同主食竹的形态差异,进而揭示该区域大熊猫不同主食竹的主要立地环境特征,以期为该区域大熊猫主食竹的精细化管理提供参考依据。

      • 四川栗子坪国家级自然保护区位于四川盆地西南缘的大渡河中上游,贡嘎山东南面石棉县境内,地理位置为E102°10'33″—E102°29'07″,N28°51'02″—N29°08'42″[15]。该保护区海拔1 330~4 550 m,面积约有47 940 hm2,为高低不平的山脊和狭窄的山谷;年平均气温和降雨量分别是11.7~14.4 ℃和800~1 250 mm。该保护区内的植被呈现出典型的垂直分布,随着海拔的升高,植被由常绿阔叶林向常绿落叶阔叶林过渡,然后是针阔混交林、针叶林至高山流石滩附近的高山灌丛[16]

        根据第四次全国大熊猫调查,该保护区内有大熊猫22只[14]。为了加强大熊猫局域小种群的复壮,中国政府发起了圈养大熊猫放归计划(CGPRP)。目前,已经有1只野外救护的大熊猫和9只圈养大熊猫被放归在该自然保护区[17]

        该自然保护区内主要分布有3种可供大熊猫采食的竹种,分别是峨热竹、石棉玉山竹和丰实箭竹;峨热竹主要生长在海拔2 500 m以上区域,其分布面积最广,丰实箭竹和石棉玉山竹分布于海拔2 800 m以下[9-10,14]。该自然保护区内大熊猫全年喜食峨热竹,冬季采食部分石棉玉山竹,偶尔采食丰实箭竹[18-19]

      • 2020年5月9日至2020年5月29日在四川栗子坪国家级自然保护区内的公益海保护站区域开展野外调查。首先,在丰实箭竹林、石棉玉山竹林和峨热竹林内分别设置4条样线,样线由低海拔到高海拔设置,样线之间距离不低于200 m,样线之间大致平行分布。其次,在每条样线上设置3~5个调查样方(20 m×20 m),同一条样线上的相邻调查样方(20 m×20 m)的海拔距离不低于50 m。在每个调查样方(20 m×20 m)中记录和测量竹种、经纬度、海拔,以及乔木层、灌木层等的相关变量数据;在每个调查样方(20 m×20 m)内的中心点以及距中心点正东和正南5 m处各设置1竹子样方(1 m×1 m),并测量和记录竹子层相关变量数据,各变量的定义及测量见表1

        表 1  各变量的定义和测量

        Table 1.  Definition and measurement of variables

        变量(单位)Variables(Units)描述与测量Definition and Measurement
        海拔(m)生境样方(20 m×20 m)中心点的海拔
        坡度(°)生境样方(20 m×20 m)中心点的坡度
        水源距离(m)生境样方(20 m×20 m)中心点到可看见或听见水流的距离
        乔木高(m)生境样方(20 m×20 m)内高于5 m的乔木的平均高度
        乔木胸径(cm)生境样方(20 m×20 m)内高于5 m的乔木平均胸径
        乔木数量生境样方(20 m×20 m)内高于5 m的乔木数
        乔木郁闭度(%)生境样方(20 m×20 m)内高于5 m的乔木郁闭度(%)
        灌木高(m)生境样方(20 m×20 m)内高≤5 m的灌木的平均高度
        灌木盖度(%)生境样方(20 m×20 m)内高≤5 m的灌木平均胸径
        灌木数量生境样方(20 m×20 m)内高≤5 m的灌木数
        竹子盖度(%)3个竹子样方(1 m×1 m)平均竹子盖度(%)
        活竹总数3个竹子样方(1 m×1 m)平均活竹子数量
        竹子死亡数3个竹子样方(1 m×1 m)平均死亡竹子数量
        一年生竹数量3个竹子样方(1 m×1 m)平均一年生竹子数量
        一年生竹高(cm)3个竹子样方(1 m×1 m)平均一年生竹子高
        一年生竹基径(mm)3个竹子样方(1 m×1 m)平均一年生竹子基经
        二年生竹数量3个竹子样方(1 m×1 m)平均二年生竹子数量
        二年生竹高(cm)3个竹子样方(1 m×1 m)平均二年生竹子高
        二年生竹基径(mm)3个竹子样方(1 m×1 m)平均二年生竹子基经
        多年生竹数量3个竹子样方(1 m×1 m)平均多年生竹子数量
        多年生竹高(cm)3个竹子样方(1 m×1 m)平均多年生竹子高
        多年生竹基径(mm)3个竹子样方(1 m×1 m)平均多年生竹子基经
      • 首先,运用单因素方差分析(One-way ANOVA)比较3种竹种的生态特征及生境变量之间是否具有显著差异。其次,通过多重比较(Duncan)比较比较丰实箭竹、石棉玉山竹和峨热竹3种竹子的形态特征和生境变量两两之间是否具有显著性差异(见表2),进而了解该自然保护区3种大熊猫主食竹的生态特征差异及其生境差异。最后,分别对3种竹子的生境各变量运用逐步判别分析,从而找出影响3种竹子生境差异的主要变量,进而揭示3种大熊猫主食竹的主要不同立地环境特征。所有数据分析均在 SPSS17.0软件中进行,变量数值以平均值±标准差(Mean±SD)表示,显著性水平设为0.05。

        表 2  不同竹种特征变量之间的差异性

        Table 2.  Differences of bamboo characteristics among different bamboo species

        变量Variables丰实箭竹F. ferax石棉玉山竹Y. lineolata峨热竹A. spanostachyaANOVA (F, P)
        竹子盖度52.77±20.9252.75±15.7749.14±8.380.30, 0.75
        活竹总数34.83±7.3233.59±5.4932.53±4.580.63, 0.54
        竹子死亡数5.65±1.62a6.80±0.79b5.55±0.90a6.22,<0.01
        一年生竹数量7.90±1.927.71±1.177.14±1.421.12, 0.34
        一年生竹高347.09±110.57a306.55±46.97a175.42±31.20b26.95,<0.001
        一年生竹基径20.24±3.06a17.04±1.42b9.24±1.28c125.98,<0.001
        二年生竹数量7.81±1.53a10.00±2.06b8.06±1.02a9.42,<0.001
        二年生竹高321.00±94.35a276.18±51.44b187.70±31.42c18.74,<0.001
        二年生竹基径17.84±3.48a18.16±1.77a10.24±1.21b62.39,<0.001
        多年生竹数量19.13±5.82a15.88±3.16b17.33±3.21ab2.45, 0.097
        多年生竹高371.10±97.05a310.89±56.24b193.98±24.97c31.36,<0.001
        多年生竹基径21.14±3.64a19.14±1.26b10.80±1.27c94.73,<0.001
      • 在丰实箭竹林、石棉玉山竹林和峨热竹林内分别调查到16、17和17个20×20 m2的调查样方。通过单因素方差分析中(One-way ANOVA)多重比较(Duncan)分析发现,3种竹子的盖度、活竹子总数量、一年生竹子数量、多年生竹子数量差异不显著。而石棉玉山竹的死亡数量显著大于其他两种竹子(F=6.22,P<0.01)。所有年龄段的丰实箭竹的竹高和基经显著大于其他竹种,而峨热竹的这些指标最小(见表2)。

      • 首先,运用One-way ANOVA和Duncan多重比较比较了丰实箭竹、石棉玉山竹、峨热竹3种竹子的生境差异。分析发现3种竹子生长的海拔均具有显著差异,峨热竹最高,而丰实箭竹最低(见表3)。丰实箭竹生长的坡度一般位于边坡等,坡度明显大于其他两种竹子(F=4.40,P=0.02)。石棉玉山竹生长的地方离水源距离最远,最小的是丰实箭竹,主要生长在溪流和河流旁。3种竹子中,峨热竹生长的生境里乔木胸径最大(F=12.27,P<0.001),而丰实箭竹和石棉玉山竹生境里的乔木数量和郁闭度显著大于峨热竹的生境。分析发现3种竹子生境中的灌木高、灌木盖度和灌木数量没有显著差异。

        表 3  不同竹种生境变量之间的差异性

        Table 3.  Differences of habitat variables among different bamboo species

        变量Variables丰实箭竹F. ferax石棉玉山竹Y. lineolata峨热竹A. spanostachyaANOVA (F, P)
        海拔2064.06±165.10a2561.76±126.61b2896.35±128.31c145.81,<0.001
        坡度35.94±24.57a18.47±16.17b18.82±16.06b4.40, 0.02
        水源距离63.69±81.79a197.12±142.78b133.82±110.76ab5.54, 0.01
        乔木高20.56±8.0422.29±7.7018.24±4.931.44, 0.25
        乔木胸径20.31±10.16a24.65±10.99a38.35±11.69b12.27,<0.001
        乔木数量36.56±14.74a30.82±19.60a13.94±5.86b10.89,<0.001
        乔木郁闭度46.19±18.33a48.06±15.20a26.12±10.62b11.13,<0.001
        灌木高3.53±0.783.03±0.923.51±0.931.74, 0.19
        灌木盖度22.31±10.5217.47±14.6519.82±12.920.59, 0.56
        灌木数量20.63±10.26a14.12±12.94ab9.29±3.98b5.54, 0.01

        其次,分别对3种竹子生境特征运用逐步判别分析,从而找出3种竹子生境差异的主要变量。分析结果发现3种竹子生境中起主要作用的变量为海拔和乔木郁闭度(见表4),这两个变量对3种竹子生境的判别成功率分别达到了87.5%、88.2%和100%(见表5)。进一步发现判别函数1对3种竹子在空间分布上的贡献最大(见图1F(2, 47)=456.20,P<0.001;VS F(2, 47)=54.74,P<0.001),其中海拔在空间上贡献大于乔木郁闭度(0.978 VS −0.214)。

        表 4  竹子生境变量标准化的判别函数

        Table 4.  Standardized discriminant functions of bamboo habitat variables

        变量Variables标准系数Standardized coefficients
        判别函数1
        Discriminant function 1
        判别函数2
        Discriminant function 2
        海拔0.9780.211
        乔木郁闭度−0.214 0.977

        表 5  竹子生境变量的判别预测

        Table 5.  Discriminant prediction of bamboo habitat variables

        竹种
        Bamboo species
        预测结果Predicted results
        丰实箭竹
        F. ferax
        石棉玉山竹
        Y. lineolata
        峨热竹
        A. spanostachya
        丰实箭竹F. ferax87.50 12.500.00
        石棉玉山竹
        Y. lineolata
        5.9088.205.90
        峨热竹
        A. spanostachya
        0.00 0.00100.00

        图  1  不同竹种生境变量判别的散点图。

        Figure 1.  Discriminant scatter diagram of different bamboo species' habitat variables

      • 四川栗子坪国家级自然保护区位于小相岭山系大熊猫主要分布区的石灰窑和公益海两个大熊猫局域小种群区域[14],本研究的研究地点位于大熊猫公益海局域小种群的分布区,该区域植被较好,林下峨热竹分布面积广长势良好,为小相岭山系大熊猫主要分布区[14]。本研究发现该区域的峨热竹的基经和高度与前人的研究类似,且位于海拔2 800 m以上的区域,这和该保护区大熊猫主要分布海拔一致[5,9-10],表明峨热竹是该区域大熊猫的主要采食竹类。而石棉玉山竹主要分布于海拔2 500 m区域,这与前人的报到一致[20],推测可能是大熊猫冬季采食的竹类[16],另外,该竹类还是该区域圈养大熊猫放归地点麻麻地区域的主要分布竹类[21]。而丰实箭竹为高大竹类,其基经和高度远远大于石棉玉山竹和峨热竹(见表2),尽管全国第四次大熊猫调查报到了该山系大熊猫采食丰实箭[14],但有关大熊猫采食丰实箭竹的策略及其选择机制还需进一步研究。

        据调查发现丰实箭竹主要位于海拔2 000 m的区域,该海拔区域段与石棉玉山竹和峨热竹分布的海拔段相比显著较低(见表3),低海拔意味着河流或小溪等可见水流较高海拔分布更密集,因此,该区域的丰实箭竹离水源距离较近,且坡度较大。另一方面,低海拔表明与保护区周边居民居住地更近,更易受到人为干扰的破坏;据报道,该区域的森林砍伐主要位于低海拔的区域[9-10]。历史的砍伐暗示了在低海拔的砍伐区域高大乔木数量较少,因此,本研究中低海拔分布的丰实箭竹林环境中的乔木胸径显著小于高海拔石棉玉山竹和峨热竹林的乔木胸径(见表3)。然而,该保护区的垂直海拔上具有明显的植被变化[20],据调查发现该区域的丰实箭竹分布区主要为落叶阔叶林,石棉玉山竹分布区域主要是落叶阔叶林和针阔混交林,而峨热竹分布区主要为针阔混交林和针叶林分布[14]。本研究主要发生于五月时乔木开新枝长新叶的时间段,这可能是低海拔的乔木郁闭度大于高海拔区域的乔木郁闭度的原因。

        通过逐步判别分析发现,在区分峨热竹、石棉玉山竹、丰实箭竹三种竹子生境的变量为海拔和乔木郁闭度,表明海拔和乔木郁闭度是这3种竹子生境主要差异变量,在判别此3种竹子生境时起主要作用。不同的竹类在海拔分布上相异的现象较多,如最为典型的是冷箭竹分布于高海拔,而其他竹类分布于低海拔[14]。大量的研究表明在佛坪自然保护区,高海拔的主要分布竹类为秦岭箭竹(Fargesia qinlingensis),而低海拔为巴山木竹(Bashania fargesii),前者为该区域大熊猫夏秋季节采食的主要竹种,而后者为冬季大熊猫的主要食物[12,22-23]。在卧龙自然保护区,高海拔分布的冷箭竹可能是该区域大熊猫夏季的主要采食竹类,而冬季可能主要采食拐棍竹[1,24]。这些研究进一步表明了这些山系的不同海拔段分布着不同的竹种,与本研究区域类似。

        乔木郁闭度是区分该区域3种竹子生境的主要变量,暗示了不同竹类其最佳的乔木郁闭度生境可能不一样。研究表明在乔木郁闭度较低的次生林内,峨热竹的盖度、高度、密度等变量均大于乔木郁闭度较大的原始林[9-10],暗示了低郁闭度或者较强的光照下有利于峨热竹的发笋及其生长,然而过密的竹林可能导致死亡竹子数的增加[9-10,16],也不利于大熊猫在其中穿。充足的水源和低海拔适宜的温度可能有利于包括竹类和其他木本植物的生长[21],然而植物的生长还与所处的土壤肥力等因素相关[25],因此,有关乔木郁闭度对石棉玉山竹和丰实箭竹的生态影响机制还需进一步研究。

        本研究表明公益海大熊猫局域种群栖息地内分布的不同竹类包括峨热竹、石棉玉山竹和丰实箭竹,它们的形态特征具有较大差异,特别是表现在具有不同的基经和高度,低海拔的丰实箭竹具有最大的基经和高度,而高海拔峨热竹的基经和高度最小。不同竹种的立地环境不同,主要表现在海拔和乔木郁闭度上,不同竹种所需求的最适宜生长的乔木郁闭度不同。以上研究结果表明了在大熊猫栖息地恢复时,不仅要注意竹种本身的生长密度、基经和高度等特征,还需要兼顾其适宜的生长环境特征,特别是所适宜的生长海拔、乔木郁闭度、土壤肥力等等。因为不同的海拔意味着不同的温度条件,而不同的郁闭度会影响光照强度的大小。

    参考文献 (25)

    目录

      /

      返回文章
      返回