大熊猫栖息地森林群落乔木层与林下泥巴山筇竹关系研究.pdf

研究报告大熊猫栖息地森林群落乔木层与林下泥巴山大熊猫栖息地森林群落乔木层与林下泥巴山筇竹关系研究筇竹关系研究秦炜锐1ab,2,刘景怡2,宋心强2,潘翰3,程勇2,付明霞2,张远彬4,周材权1b,杨彪1a*1.西华师范大学 a.栗子坪大熊猫生态与保护四川省野外科学观测研究站,b.生态研究院，四川南充637009；2.荥经县大相岭自然保护区管护中心，四川雅安625000；3.北京市企业家环保基金会，北京100020；4.中国科学院、水利部成都山地灾害与环境研究所，成都6100412023-02-06收稿,2023-04-18网络版发表开放科学(资源服务)标识码(OSID):基金项目四川省科学技术厅重大科技专项项目（2021YFS0280，22ZDYF2289）；西华师范大学科研启动项目（18Q038）；国家自然科学基金项目资助（32270551）；秦炜锐（1997），女，硕士研究生，*通信作者摘要摘要为探索森林群落乔木层与林下大熊猫取食竹的相关关系，利用样线法和样方法，对大相岭大熊猫栖息地森林群落乔木层结构特征与林下泥巴山筇竹（Qiongzhuea multigemmia）生长指标进行调查，构建广义线性混合模型（GLMM）和线性混合模型（LMM）进行分析。结果表明：（1）乔木郁闭度对一年生、多年生泥巴山筇竹笋的基径，以及泥巴山筇竹活竹株数、死竹株数和笋高度的影响差异均显著；乔木株数对泥巴山筇竹发笋量、死竹株数和活竹株数的影响差异均显著；乔木平均高度对泥巴山筇竹发笋量、一年生和多年生泥巴山筇竹的高度及基径的影响差异均显著。（2）泥巴山筇竹发笋量随乔木株数和平均高度增加而减少；活竹株数和死竹株数均随乔木层郁闭度和乔木株数增加而减少；一年生泥巴山筇竹平均高度随乔木平均高度增加而减少；一年生泥巴山筇竹平均基径均随乔木层郁闭度和乔木平均高度增加而减少；多年生泥巴山筇竹平均高度和基径均随乔木平均高度增加而增加；泥巴山筇竹笋高度随乔木层郁闭度增加而增加；多年生泥巴山筇竹平均基径和笋基径均随乔木层郁闭度增加而减少。本研究能为该区域大熊猫栖息地的恢复提供量化依据。关键词关键词大熊猫栖息地;栖息地恢复;林冠乔木层;泥巴山筇竹;广义线性混合模型大熊猫（Ailuropoda melanoleuca）是全球生物多样性保护的旗舰物种与伞护物种。过去几十年大熊猫保护取得了举世瞩目的成就，但该物种仍然面临繁殖力弱、存活率低、栖息地破碎化等问题1。根据全国第四次大熊猫调查的结果，仍有 30%的大熊猫栖息地质量有问题2。由于自然隔离和人为干扰，大熊猫栖息地被隔离成 33 个斑块，面积小于 10000hm2的斑块有 9 个，栖息地破碎化将给野生大熊猫种群的生存、繁衍带来威胁3。大相岭山系有大熊猫38 只2，是六大山系中种群数量最少的4。该区域栖息地面积小、破碎化程度高、人类活动干扰严重，亟待加强保护4。伴随大熊猫国家公园的建立，栖息地恢复工作是未来大熊猫国家公园保护管理工作的重点。大熊猫是森林类型的物种，竹类资源、森林起源以及遮蔽条件是影响大熊猫生境利用的主要因素5。森林中的乔木为大熊猫提供躲避、休息以及巢穴。林下竹子作为大熊猫主要的食物来源，是大熊猫生存的基础。已有研究表明，森林群落乔木层郁闭度对林下植物的生长起重要作用6-9，乔木是森林群落引用格式引用格式：秦炜锐,刘景怡,宋心强,等.大熊猫栖息地森林群落乔木层与林下泥巴山筇竹关系研究J.四川林业科技,2023,44（5）:714.QIN W R,SONG X Q,PAN H,et al.Study on the relationship between tree layer and understory Qiongzhuea multigemmia of forestcommunityingiantpandahabitatJ.JournalofSichuanForestryScienceandTechnology,2023,44(5):714.doi:10.12172/202302060001.第44卷第5期四川林业科技Vol.44,No.52023年10月JournalofSichuanForestryScienceandTechnologyOct.,20232023四川林业科技四川林业科技编辑部编辑部恢复过程中重要的生物因子10，乔木层通过对光的透射、反射和吸收，直接影响林下光照的强度和分布，而光照对于绝大多数森林树种而言，是决定更新个体能否生存和生长的关键11。乔木层郁闭度对竹子的高度、基径和密度等生长发育具有重要的影响7,12。因此，研究乔木层结构特征与林下竹子生长指标之间的关系是栖息地恢复的基础，对提升大熊猫栖息地质量、保护大熊猫及全球生物多样性具有至关重要的作用。泥巴山筇竹是非常典型的大熊猫取食竹，是大相岭山系泥巴山区域大熊猫最主要的取食来源。而有关泥巴山筇竹的研究少之又少。基于此，本研究选择大相岭山系的大熊猫栖息地作为研究对象，针对大熊猫栖息地森林群落乔木层结构特征与林下泥巴山筇竹生长指标开展调查，定量分析两者之间的关系，并建立大熊猫栖息地森林群落乔木层特征与林下泥巴山筇竹生长关系的模型，以期为该区域大熊猫栖息地恢复提供量化依据。1研究区概况研究区概况研究区选择在大熊猫国家公园荥经片区泥巴山廊道区域原生森林生态系统林下分布有泥巴山筇竹（Qiongzhuea multigemmia）的区域（图 1）。该区域海拔 17242445m，属亚热带高山山地气候。植被类型以冷杉（Abies fabri）、云杉（Picea asperata）、槭树（Acerspp.）、中华木荷（Schima sinensis）、刺榛（Corylus ferox）、天师栗（Aesculus chinensisvar.wilsonii）、皂柳（Salix wallichiana）、青冈（Quercusglauca）等为主要树种的针阔混交林，林下以泥巴山筇竹为主，是理想的大熊猫栖息地。2研究方法研究方法2.1野外调查野外调查调查于 2022 年 56 月进行。设置样线，在大熊猫国家公园荥经片区泥巴山廊道区域大熊猫种群动态调查监测网格内布设 20m20m 的乔木样方共68 个，记录样方的海拔（m）、坡度（）、坡向（）、坡位（背坡、山肩、山顶、麓坡、趾坡）、坡型（均匀坡、复合坡、凹坡、凸坡），以及乔木层郁闭度、乔木株数和乔木平均高度等。每个乔木样方内分别按“品”字形获取 3 个 2m2m 的泥巴山筇竹样方，共计 204 个泥巴山筇竹样方。计数每个样方中泥巴山筇竹发笋量、活竹株数和死竹株数，测不同年份泥巴山筇竹的基径和高度、笋基径和笋高度。2.2数据处理方法数据处理方法本研究以野外调查所得的 204 个泥巴山筇竹样方的样本数据为建模样本。依据乔木层数（3 层），将乔木平均高度划分为 3 个等级：矮：12m；依据不同样方内的乔木株数之间的差异，将乔木株数划分为 4 个等级：A：020 株，B：2140 株，C：4160 株，D：60 株；由于竹子在不同乔木层郁闭度下有差异，而通常乔木层郁闭度又被分为 4 个等级，所以将乔木层郁闭度划分为 4 个等级：A：00.25，B：0.260.5，C：0.510.75，D：0.76-1.00。Yij=a0+a1Xij+bi+eijYija0a1Xijbieijg(ij)=ij=Xij+Zijig(.)ijijXijZiji利用 R 软件进行单因素方差分析，以检验乔木层结构特征对泥巴山筇竹生长指标的影响，然后进行逐步回归分析，选取显著性最高的因子，建立估测模型，根据模型评价指标，选择最优回归模型。对于不同年份泥巴山筇竹高度和基径、笋高度和基径数据，考虑环境因子（海拔、坡度、坡向、坡位和坡型）的随机效应，构建线性混合效应模型13：，其中，表示第 i 类的第 j 个观测的反应变量，表示固定截距，表示固定斜率，表示固定效应，表示随机效应（只影响截距），表示残差；对于泥巴山筇竹无性系种群数量（包括发笋量、活竹株数和死竹株数）数据，选择泊松分布形式，考虑环境因子（海拔、坡度、坡向、坡位和坡型）的随机效应，构建广义线性混合模型14-16：，其中，表示联结函数，表示条件线性预测值，表示条件均数（考虑了海拔、坡度、坡向、坡位和坡型的随机效应），表示解释变量，表示固定效应参数矢量，表示与随机效应相关的解释变量，表示随机效应。3结果与分析结果与分析3.1泥巴山筇竹及乔木层调查结果泥巴山筇竹及乔木层调查结果结果显示（表 1）：每个泥巴山筇竹样方内泥巴山筇竹平均发笋量为（3.064.58）株，死竹平均株数为（10.838.05）株，活竹平均株数为（41.7126.28）株，一年生泥巴山筇竹平均基径为（7.53四川林业科技2023年第 44卷第 5期81.86）mm，一年生泥巴山筇竹平均高度为（175.9049.72）cm，多年生泥巴山筇竹平均基径为（8.411.77）mm，多年生泥巴山筇竹平均高度为（199.7855.67），笋平均基径为（8.352.05），笋平均高度为（134.0957.12）cm。每个植物样方内乔木平均高度为（7.422.61）m，乔木平均株数为（36.6620.86）株，乔木平均郁闭度为（0.420.20）。3.2乔木层结构特征与泥巴山筇竹生长指标的乔木层结构特征与泥巴山筇竹生长指标的关系关系单因素方差分析的结果表明，乔木郁闭度对一年生、多年生和泥巴山筇竹笋的基径的影响差异显著（p0.05），对泥巴山筇竹活竹株数、死竹株数和笋高度的影响差异极显著（p0.01）。乔木株数对泥巴山筇竹死竹株数的影响差异显著（p0.05），对泥巴山筇竹发笋量和活竹株数的影响差异极显著（p0.01）。乔木平均高度对泥巴山筇竹发笋量和多年生竹基径的影响差异显著（p0.05），对一年生泥巴山筇竹高度和基径、多年生泥巴山筇竹高度的影响差异极显著（p0.01）（表 2）。3.3模型构建模型构建根据方差分析结果，挑选对泥巴山筇竹生长指标有显著性影响的因子，将其作为模型的固定效表1泥巴山筇竹及乔木层调查结果（平均值标准差）Tab.1SurveyresultsofQiongzhuea multigemmiaandtreelayers(MeanSD)乔木平均高度/mAverageheightoftrees(m)乔木株数/株Numberoftrees乔木层郁闭度Treecanopydensity发笋量/株Shootingquantity死竹株数/株Numberofdeadbambooplants活竹株数/株Numberoflivebambooplants7.422.6136.6620.860.420.203.064.5810.838.0541.7126.28一年生竹基径/mmBasaldiameterofannualbamboo(mm)一年生竹高度/cmHeightofannualbamboo(cm)多年生竹基径/mmBasaldiameterofperennialbamboo(mm)多年生竹高度/cmHeightofperennialbamboo(cm)笋基径/mmBasaldiameterofbambooshoots(mm)笋高度/cmBambooshootheight(cm)7.531.86175.9049.728.411.77199.7855.678.352.05134.0957.12雅安片区在大熊猫国家公园的位置雅安片区荥经片区在大熊猫国家公园雅安片区的位置荥经片区图例研究地点大熊猫国家公园荥经片区-核心保护区大熊猫国家公园荥经片区-一般控制区Miles015306090120N图1研究区域Fig.1Studyarea秦炜锐，等：大熊猫栖息地森林群落乔木层与林下泥巴山筇竹关系研究研究报告9应。因此，可将乔木株数和乔木平均高度 2 个因子，保留在泥巴山筇竹发笋量模型中；将乔木层郁闭度和乔木株数 2 个因子，均保留在活竹株数和死竹株数模型中；将乔木平均高度因子保留在一年生泥巴山筇竹平均高度模型中；将乔木层郁闭度和乔木平均高度 2 个因子，均保留在一年生和多年生泥巴山筇竹平均基径模型中；将乔木平均高度因子保留在多年生泥巴山筇竹平均高度模型中；将乔木层郁闭度因子保留在泥巴山筇竹笋高度和笋基径模型中。采用赤池信息准则（AIC）指标来比较不同模型间的模拟效果。AIC 值越小，表明模拟效果越好。结果表明，所挑选的主要因子模拟效果均好于没有显著性影响的因子的模型拟合效果。使用不同的有显著性的因子作为固定效应，其混合模型的精度不同。在泥巴山筇竹发笋量、活竹株数和死竹株数模型中，以乔木株数作为固定效应的模型 AIC 值均最小，拟合精度均最高；在一年生平均基径模型中，以乔木平均高度作为固定效应的模型 AIC 值最小，拟合精度最高；在多年生平均基径模型中，以乔木层郁闭度作为固定效应的模型 AIC 值最小，拟合精度最高（表 3）。结果显示，泥巴山筇竹活竹株数、死竹株数、一年生平均高度、一年生平均基径、多年生平均高表2乔木层结构特征不同组别的泥巴山筇竹生长指标方差分析Tab.2VarianceanalysisofgrowthindexofQiongzhuea multigemmiaindifferentgroupsoftreelayerstructurecharacteristics乔木层指标Treelayerindex泥巴山筇竹指标Qiongzhuea multigemmiaindexdfFP乔木层郁闭度Treelayerdensity发笋量Shootingquantitydf1=3,df2=1980.8310.478活竹株数Numberoflivebambooplantsdf1=3,df2=20013.590.001*死竹株数Numberofdeadbambooplantsdf1=3,df2=1784.540.004*一年生竹高度Heightofannualbamboodf1=3,df2=122.4730.064一年生竹基径Basaldiameterofannualbamboodf1=3,df2=1522.6970.048*多年生竹高度Heightofperennialbamboodf1=3,df2=1971.5430.205多年生竹基径Basaldiameterofperennialbamboodf1=3,df2=1973.6840.013*笋高度Bambooshootheightdf1=3,df2=1004.2890.007*笋基径Basaldiameterofbambooshootsdf1=3,df2=1003.3260.023*乔木株数Numberoftrees发笋量Shootingquantitydf1=3,df2=1988.4320.001*活竹株数Numberoflivebambooplantsdf1=3,df2=2008.4470.001*死竹株数Numberofdeadbambooplantsdf1=3,df2=1783.1590.026*一年生竹高度Heightofannualbamboodf1=3,df2=1521.1240.341一年生竹基径Basaldiameterofannualbamboodf1=3,df2=1520.4930.688多年生竹高度Heightofperennialbamboodf1=3,df2=1971.3970.245多年生竹基径Basaldiameterofperennialbamboodf1=3,df2=1971.7890.15笋高度Bambooshootheightdf1=3,df2=1001.8790.138笋基径Basaldiameterofbambooshootsdf1=3,df2=1000.5840.627乔木平均高度Averageheightoftrees发笋量Shootingquantitydf1=2,df2=1094.6250.011*活竹株数Numberoflivebambooplantsdf1=2,df2=2010.7660.466死竹株数Numberofdeadbambooplantsdf1=2,df2=2012.0560.131一年生竹高度Heightofannualbamboodf1=2,df2=1537.510.001*一年生竹基径Basaldiameterofannualbamboodf1=2,df2=1537.2610.001*多年生竹高度Heightofperennialbamboodf1=2,df2=1987.5560.001*多年生竹基径Basaldiameterofperennialbamboodf1=2,df2=1983.7470.025*笋高度Bambooshootheightdf1=2,df2=1011.0960.338笋基径Basaldiameterofbambooshootsdf1=2,df2=1010.1360.873注：*为P0.05，*为P0.01，*为P0.001。Note:*meansP0.05,*meansP0.01and*meansP0.001.四川林业科技2023年第 44卷第 5期10度、多年生平均基径、笋高度和笋基径模型的固定效应截距均达到极显著水平（P0.001），说明乔木层结构特征各参数对模型变化影响极显著。乔木株数的固定效应系数对泥巴山筇竹发笋量模型变化影响显著（P0.05），乔木株数的固定效应系数对泥巴山筇竹活竹和死竹株数模型变化影响均极显著（P0.01），乔木平均高度的固定效应系数对多年生泥巴山筇竹平均高度模型变化影响显著（P0.05）（表 3）。随着乔木株数和乔木平均高度的增加，泥巴山筇竹发笋量均呈减少的趋势；随着乔木层郁闭度和乔木株数的增加，泥巴山筇竹活竹株数和死竹株数均呈减少的趋势；随着乔木平均高度的增加，一年生泥巴山筇竹平均高度呈减少的趋势；随着乔木层郁闭度和乔木平均高度的增加，一年生泥巴山筇竹平均基径均呈减少的趋势；随着乔木平均高度的增加，多年生泥巴山筇竹平均高度和基径均呈增加的趋势；随着乔木层郁闭度的增加，泥巴山筇竹笋高度呈增加的趋势，多年生平均基径、笋基径呈减少的趋势（表 3）。4结论与讨论结论与讨论竹子和木本植物是大熊猫栖息地的重要组成部分17，竹子生长状况与木本植物之间具有关系18。大熊猫、森林和竹子三者协同进化19。本研究在大熊猫国家公园荥经片区泥巴山廊道区域原生森林生态系统林下分布有泥巴山筇竹的区域，利用样线法和样方法，开展森林群落乔木层与林下泥巴山筇竹调查，并利用 GLMM 和 LMM 模型，构建了大熊猫栖息地森林群落乔木层结构特征与林下泥巴山筇竹生长指标的关系，这为大熊猫栖息地恢复提供了定量的数据支撑。修复效果、修复要求和修复考虑因素是大熊猫栖息地修复的主要研究主题20。竹子的空间分布对表3泥巴山筇竹生长指标模型拟合结果Tab.3ModelfittingresultsofQiongzhuea multigemmia八月竹指标Bambooindex乔木层指标Treelayerindex固定效应截距Fixedeffectintercept固定效应系数Fixedeffectcoefficient随机效应截距Randomeffectintercept残差ResidualAIC发笋量Shootingquantity乔木株数Numberoftrees0.95270.0191*4.5231842.5乔木平均高度Averageheightoftrees0.51940.20516.9233844.3活竹株数Numberoflivebambooplants乔木层郁闭度Treecanopydensity4.0594*0.53170.78872392.1乔木株数Numberoftrees3.9996*0.0061*0.30522388.5死竹株数Numberofdeadbambooplants乔木层郁闭度Treecanopydensity2.3415*0.42070.29061362.1乔木株数Numberoftrees2.5597*0.0115*0.25651354.7一年生平均高度Averageannualheight乔木平均高度Averageheightoftrees170.3135*0.23131923.8400867.54001625.2一年生平均基径Annualaveragebasaldiameter乔木层郁闭度Treecanopydensity7.8847*0.87542.87441.1910598.7乔木平均高度Averageheightoftrees8.4845*0.12482.74241.1820597.3多年生平均高度Perennialaverageheight乔木平均高度Averageheightoftrees154.7610*5.8200*1863.38301065.00002105.2多年生平均基径Perennialaveragebasaldiameter乔木层郁闭度Treecanopydensity8.4666*0.87162.22881.0310717.5乔木平均高度Averageheightoftrees7.5236*0.06512.23961.0293717.7笋高度Bambooshootheight乔木层郁闭度Treecanopydensity192.3600*11.88002109.62001070.56002109.7笋基径Basaldiameterofbambooshoots乔木层郁闭度Treecanopydensity8.7150*1.80431.73552.4800439.2注：*为P0.05，*为P0.01，*为P0.001。Note:*meansP0.05,*meansP0.01and*meansP0.001.秦炜锐，等：大熊猫栖息地森林群落乔木层与林下泥巴山筇竹关系研究研究报告11于大熊猫栖息地的分布具有重要影响21，竹子密度、基径、高度和幼竹所占比例是衡量森林中竹子资源的重要指标22，也是表征竹子生长状况的指标18。大熊猫的栖息地选择对获得高可用性的美味竹子和良好的庇护所具有重要的影响23。幼竹的可利用性是竹林生境选择的重要驱动因素24，竹林密度、高度、基径、竹笋的生长、出笋量、无性系根茎等，均受温度、湿度、光照强度等环境因子的显著影响25。乔木作为森林群落的最高层，其郁闭度的大小会对林下环境（如光照强度、温度和湿度等）产生影响，从而影响竹子的生长和天然更新。冠层郁闭度和层片结构影响同一群落内缺苞箭竹（Fargesia denudata）天然更新26；林冠条件影响矮竹秆动态变化27。有研究表明，八月竹（Chimonobambusa szech-uanensis）发笋量随乔木层郁闭度增加而减少28，峨热竹（Arundinaria spanostachya）发笋率、老笋比例和成竹死亡率均随上层乔木郁闭度增加而降低29，峨热竹发笋量随乔木平均高度增加而减少30-32。本研究中，乔木层郁闭度对泥巴山筇竹发笋量的影响差异不显著，但乔木株数和乔木平均高度对泥巴山筇竹发笋量的影响差异显著，这可以解释为乔木层郁闭度对泥巴山筇竹发笋量的影响可能与乔木株数和乔木平均高度有关。随着乔木株数和乔木平均高度的增加，乔木树冠增加，冠幅增大，乔木层郁闭度也增加，泥巴山筇竹发笋量减少。次生林竹子生长密度大，营养质量低31；不同郁闭度对不同竹种的生长更新的影响不同，如全光照下细枝箭竹（F.stenoclada）有较高的密度、基径和高度33；茂密的矮竹会通过争夺光、水、营养，对树苗的再生产生负面影响，影响冠层树木种子的保留34。泥巴山筇竹活竹株数和死竹株数均随着乔木层郁闭度和乔木株数的增加，这可能是因为乔木株数与乔木层郁闭度增加，乔木根系增多，环境也越来越隐蔽，对林下泥巴山筇竹根系的发育产生影响，使竹子幼苗的根系生长脆弱，根系竞争能力差，影响幼苗对水分和养分的吸收，造成幼苗死亡35。同时也说明了过密的泥巴山筇竹纯林会影响冠层树木种子的保留。解蕊等36的研究结果表明缺苞箭竹的高度随乔木郁闭度增加而增加；宋利霞37和张聪等38的研究结果表明华西箭竹（F.nitida）的高度与乔木郁闭度呈负相关关系。而本研究结果表明：一年生泥巴山筇竹平均高度和基径均随乔木平均高度的增加而减少，而多年生泥巴山筇竹平均高度随乔木平均高度的增加而增加，多年生泥巴山筇竹平均基径随乔木层郁闭度的增加而减少。这可能是因为后期光照强度降低导致了泥巴山筇竹无性系分株的光合作用下降，积累的光合产物和能量降低39，从而引起多年生泥巴山筇竹平均基径的变小；也可能是因为随着竹龄的增长，泥巴山筇竹根系吸收能力下降，以至于基径变小，高度增加。同时，泥巴山筇竹笋高度随着乔木层郁闭度的增加而增加，说明郁闭度过高会使泥巴山筇竹笋高度增加，也体现了不同克隆植物对应环境异质性策略的差异。总的来说，乔木株数和乔木平均高度、乔木层郁闭度对泥巴山筇竹的生长发育均有着不同的影响。大熊猫喜欢活动于上层有乔木、下层有竹林的环境中，不同的山系大熊猫喜欢的乔木层郁闭度不同，但基本不利用竹子密度过高的纯竹林40。大熊猫国家公园荥经片区是大相岭山系大熊猫种群的核心分布和集中分布区，也是连接大相岭山系和邛崃山山系大熊猫种群基因交流的关键区。区内有大片上层无乔木层覆盖的泥巴山筇竹纯林，以及林下无竹或极少有竹子生长的人工次生林，不属于大熊猫理想的栖息地。对大熊猫国家公园荥经片区大熊猫栖息地进行恢复是恢复大相岭山系大熊猫栖息地的重要基础，是保护、维持大熊猫种群数量的有效途径。因此，要为该区域大熊猫创造良好的栖息地环境，可以将森林和竹子之间的关系，作为大熊猫栖息地人工恢复的参考标准。通过研究了解了大相岭山系大熊猫栖息地森林群落乔木层结构特征与林下泥巴山筇竹的生长指标之间的关系，可以科学合理的采取经营措施，如通过栽种本土阔叶树种，提高乔木层郁闭度，从而控制林下泥巴山筇竹无性系种群数量，降低泥巴山筇竹林密度，提高泥巴山筇竹笋质量，以确保合理的乔木层结构特征促进大熊猫取食竹生长更新。参考文献1WEIFW,SWAISGOODR,HUYB,etal.ProgressintheecologyandconservationofgiantpandasJ.ConservationBiology,2015,29(6):14971507.四川林业科技2023年第 44卷第 5期122国家林业和草原局.全国第四次大熊猫调查报告M.北京:科学出版社,2021.3唐小平，贾建生，王志臣，等.全国第四次大熊猫调查方案设计及主要结果分析J.林业资源管理，2015（1）：1116.4冉江洪，曾宗永，刘世昌，等.四川大相岭大熊猫种群及栖息地调查J.四川大学学报(自然科学版)，2006（4）：889893.5康东伟.大熊猫的生境选择研究D.北京:北京林业大学,2015.6TAYLOYAH,QINZS.CulmDynamicsandDryMatterProductionofBamboosintheWolongandTangjiaheGiantPandaReserves,Sichuan,ChinaJ.JournalofAppliedEcology,1987,24(2):419433.7REID,DG,TAYLOYAH,HUJC,etal.EnvironmentalinfluencesonbambooBashania fangianagrowthandimplicationsforgiantpanadaconservationJ.JournalofAppliedEcology,1991,28(3):855868.8王金锡，马自贵，刘长祥，等.缺苞箭竹生长发育规律初步研究J.竹子研究汇刊，1991（3）：3848.9秦自生，艾伦泰勒，蔡绪慎，等.冷箭竹生物学特性研究J.西华师范大学学报(自然科学版)，1994，15（2）：107113.10李媛良，汪思龙，宿秀江，等.天然林重建过程中单优箬叶竹灌丛对树木更新的影响J.生态学报，2009，29（12）：66156621.11游莉.伏牛山自然保护区森林群落冠层结构及光环境特征研究D.开封:河南大学,2009.12VEBLEN T T.Growth patterns of Chusquea bamboos in the understory of Chilean Nothofagus forests and their influences in forestdynamicsJ.BulletinoftheTorreyBotanicalClub,1982,109(4):474487.13许崇华，崔珺，黄兴召，等.基于线性混合效应模型的杉木树高-胸径模型J.西北农林科技大学学报:自然科学版，2017，45（6）：5360.14李丽霞，郜艳晖，张丕德，等.广义线性混合效应模型及其应用J.现代预防医学，2007，34（11）：21032104.15李春明，付卓.基于广义混合效应模型的云冷杉林天然更新计数方法研究J.西南林业大学学报:自然科学，2020，40（5）：108114.16李春明.基于广义线性混合效应模型的蒙古栎林单木枯损建模及影响因子分析J.林业科学研究，2020，33（6）：105113.17KANGDW,WANGXR,LIS,etal.ComparingtheplantdiversitybetweenartificialforestandnaturegrowthforestinagiantpandahabitatJ.ScientificReports,2017,7(1):17.18KANGD,LVJ,LIS,etal.RelationshipbetweenbamboogrowthstatusandwoodyplantsinagiantpandahabitatJ.Ecologicalindicators,2019,98:840843.19申国珍.大熊猫栖息地恢复研究D.北京:北京林业大学,2002.20KangDW.Areviewofthehabitatrestorationofgiantpandasfrom2012to2021:ResearchtopicsandadvancesJ.ScienceofTheTotalEnvironment,2022,852:158207.21LINDERMANM,BEARERS,ANL,etal.Theeffectsofunderstorybambooonbroad-scaleestimatesofgiantpandahabitatJ.BiologicalConservation,2005,121(3):383390.22申国珍，李俊清，任艳林，等.大熊猫适宜栖息地恢复指标研究J.北京林业大学学报，2002，24（4）：15.23FuM,PanH,SongX,etal.Back-and-forthshiftsinhabitatselectionbygiantpandasoverthepasttwodecadesintheDaxianglingMountains,southwesternChinaJ.JournalforNatureConservation,2022,66:126129.24KangDW,YangHW,LiJQ,etal.HabitatusebygiantpandasAiluropodamelanoleucaintheWanglangNatureReserve,Sichuan,ChinaJ.ZoologicalStudies,2013,52(1):154159.25潘红丽，李迈和，田雨，等.卧龙自然保护区油竹子形态学特征及地上部生物量对海拔梯度的响应J.四川林业科技，2010（3）：3440.26贾炅，武吉华.四川王朗自然保护区大熊猫主食竹天然更新J.北京师范大学学报:自然科学版，1991，027（002）：250256.27TomohiroItabashi,ShinjiAkada,KiyoshiIshida,etal.Culmdynamicsofdwarfbamboo(SasakurilensisMakino&Shibata)inrelationtoforestcanopyconditionsinbeechforestsJ.AdvancesinBambooScience.2023,228黄荣澄，刘香东，冉江洪，等.大熊猫主食竹八月竹笋期生长发育规律初步研究J.四川大学学报:自然科学版，2011，48（2）：469473.29魏辅文，冯祚建.相岭山系大、小熊猫主食竹类峨热竹的生长发育与环境因子间的相互关系J.生态学报，1999，19（5）：710714.30HongMS,YuanSB,YangZS,etal.ComparisonofmicrohabitatselectionandtraceabundanceofgiantpandasbetweenprimaryandsecondaryforestsinLizipingNaturReserve,China:Effectsofselectivelogging.MammalianBiology,2015,80,373-379.31Hong,MS,WeiW,YangZS,etal.EffectsoftimberharvestingonArundinariaspanostachyabambooandfeeding-siteselectionbygiantpandasinlizipingnaturereserve,ChinaJ.ForestEcologyandManagement,2016,373:7480.32XieJM,HongMS.EffectsofecologicalfactorsongrowthofArundinariaspanostachyashootsinLizipingNationalNatureReserve,ChinaJ.GlobalEcologyandConservation,2020,23:e01121.33彭科，陈旭，张宣，等.四川白水河自然保护区大熊猫主食竹种群特征及其影响因素J.四川林业科技，2020，41（5）：94104.34QianF,ZhangTD,GuoQX,etal.DenseunderstorydwarfbambooalterstheretentionofcanopytreeseedsJ.ActaOecologica,2016,73:3844.35马志贵，王金锡.缺苞箭竹更新与生态因子关系的研究J.四川林业科技，1989，10（4）：3439+75.36解蕊，李俊清，赵雪，等.林冠环境对亚高山针叶林下缺苞箭竹生物量分配和克隆形态的影响J.植物生态学报，2010，34（6）：753760.37宋利霞.华西箭竹克隆种群对卧龙自然保护区暗针叶林不同林冠环境的响应研究D.重庆:西南大学,2006.38张聪，曾涛，唐明坤，等.九寨沟自然保护区华西箭竹生长研究J.四川大学学报(自然科学版)，2010，47（5）：11371143.秦炜锐，等：大熊猫栖息地森林群落乔木层与林下泥巴山筇竹关系研究研究报告1339杨在娟，岳春雷，汪奎宏.光照强度对雷竹无性系生长的影响J.浙江林业科技，2002（3）：7476+83.40胡锦矗.大熊猫生物学研究与进展M.四川科学技术出版社,1990.Summary for“大熊猫栖息地森林群落乔木层与林下泥巴山筇竹关系研究”StudyontheRelationshipbetweenTreeLayerandUnderstoryQiongzhueamultigemmiaofForestCommunityinGiantPandaHabitatQINWeirui1ab,2,LIUJingyi2,SONGXinqiang2,PANHan3,CHENGYong2,FUMingxia2,ZHANGYuanbin4,ZHOUCaiquan1b,YANGBiao1a*1ChinaWestNormalUniversitya.SichuanLizipingGiantP