随着人口增加，经济活动不断加剧，生物多样性正急剧下降，生态系统更是遭到了严重破坏，大量物种已经灭绝或处于灭绝的边缘^{[1, 2]}。目前，物种保护特别是濒危物种保护成了国际社会关注的热点。

秦岭冷杉(Abies chensiensis)是松科冷杉属常绿乔木，为我国特有珍稀濒危植物，主要分布于秦巴山地^[3]。秦岭冷杉树干通直，木质松软，富含冷杉胶，是非常优良的用材和林化产品资源树种^[4]。秦岭冷杉喜气候温凉湿润，通常生于山沟溪旁或阴坡，天然更新较差，加上过度采伐，分布面积日益缩小，现多为零散分布。秦岭冷杉的濒危状态已经引起不少学者的关注^{[5, 6]}。

秦岭冷杉死亡的现象在分布地陕西秦岭，神农架等都有报道^{[7, 8, 9]}，特别是河南省鲁山县石人山，秦岭冷杉出现大量死亡，死亡数量目前已达200多株^[10]。小秦岭国家级自然保护区主峰高2413.8 m，为河南最高峰，生物资源丰富，区系成分复杂。小秦岭区域的秦岭冷杉群落是河南省保存较好的野生秦岭冷杉群落。可经过多年的跟踪调查，近3—4a小秦岭区域的秦岭冷杉也出现了不明原因的衰败死亡现象，并有不断扩展的趋势。为避免秦岭冷杉死亡现象的进一步扩大，分析研究就地保护这一珍稀植物资源迫在眉睫。

许多学者针对秦岭冷杉已经开展了一些非常有意义的研究。朱秀红等在石人山研究了秦岭冷杉群落特征并进行了群落分类^[10]；李为民等描述了秦岭山地秦岭冷杉的群落结构特征^[3]；张文辉等探讨了秦岭冷杉的地理分布和生物生态学特征^[4]。这些研究结果有助于人们了解秦岭冷杉的群落结构特征。还有许多学者针对秦岭冷杉的结实特性^{[8, 9]} 、幼苗扦插^[11]和种子萌发^{[7, 12]}开展大量研究，在秦岭冷杉幼苗扦插、引种保护和培育等方面起到了指导作用。然而在秦岭冷杉生长状况不良并不断死亡原因研究方面并没有进展。是什么因素造成了秦岭冷杉的死亡？这方面一直没有相关报道。

自然生态系统是多个相互作用的过程组成的，物种的生长和死亡过程也是复杂的^{[13, 14]}。随着研究的深入，单变量统计方法已经不能完全满足生态学发展的需要，因为它无法提供所研究内容的完整信息^[15]。结构等式模型作为一种多变量的统计方法应运而生，为此类复杂问题的解决提供了一个合适的途径，近几年逐步为生态学家所应用，成为解决生态学中复杂问题的重要手段，对相关科学问题的解决提供了独特视角^{[16, 17, 18]}。

土壤含水率、土壤温度和土壤电导率是土壤理化性质的重要指标，可以影响到植物根系对矿质元素的吸收以及呼吸速率等生理过程^[19]；不同坡度的土壤养分含量以及不同海拔高度温度和降水等环境因素都存在很大差异^{[20, 21, 22]}；而秦岭冷杉的垂直分布跨越很大的海拔范围，受海拔影响明显^{[8, 20]}。定量分析这些环境因子对秦岭冷杉生长的影响对了解其种群更新和死亡具有很大意义。

本研究在对小秦岭全部秦岭冷杉定位监测的基础上，通过对活树、死树分布和结构的对比分析，并采用相关分析和结构等式模型(structural equation model,SEM)结合环境因子，探讨了秦岭冷杉的生存状况和死亡原因，研究结果对于全面了解小秦岭地区秦岭冷杉生存现状和死亡原因具有重要意义，同时对秦岭冷杉的扦插繁殖、引种育苗等具有参考价值。

河南小秦岭国家级自然保护区，位于豫陕两省交界的灵宝市西部、小秦岭北麓；地理坐标为北纬34°23′—34°31′，东经110°23′—110°44′；保护区东西长31 km，南北宽12 km，总面积15160 hm²，其主峰老鸦岔高达2413.8m，为河南省境内最高峰。本区地处暖温带向亚热带的过渡区，为暖温带大陆性季风型半干旱气候，四季分明。年平均气温依不同的海拔，变化在5.5—14.2℃之间，极端最低气温为-17.0℃，极端最高气温为42.7℃；年均降水量在506—719.2 mm。

通过设置样方的方法只能调查到部分秦岭冷杉。为实现对小秦岭地区全部秦岭冷杉的监测，2014年6月，由保护区工作人员带领，在对保护区全面踏查的基础上，对保护区所有的秦岭冷杉进行实地调查，包括秦岭冷杉活树个体和枯死树个体。采用GPS对每棵秦岭冷杉定位，长期监测，建立冷杉档案。对每棵树挂牌标记，记录树干1.3 m处胸径(DBH)，树高低于1.3 m的秦岭冷杉测量其树干DBH最大处，并采用激光测高仪测量树高。调查过程中为降低人为观测影响造成的差异，参考树木健康性的评价定级^[23]，由同一人记录秦岭冷杉生长生境，分良好、一般和枯死评价秦岭冷杉个体存活状态。同时采用GPS记录每棵秦岭冷杉经纬度和海拔，地质罗盘来测量坡向和坡度。在树基不同方向1 m范围内选取3个点，采用土壤三参数速测仪测定了土壤表层土(0—10 cm)的土壤温度、土壤含水量和土壤电导率。

同时，为了解秦岭冷杉群落的概况，从海拔2000 m开始(秦岭冷杉主要生长于海拔2000 m以上)，选择秦岭冷杉群落设置样方，共 16个20 m ×20 m样地。样方采用GPS 定位，记录其经纬度、海拔、坡度、坡向，及胸径≥1 cm 乔木种类、数量、胸径、树高和冠幅。

本研究采用Spearman 相关分析分别探讨活树和死树与环境因子之间的相关性。活树和死树可能同时受到多个环境因素的影响，因此同时采用偏相关分析，在剔除其他环境因子影响的条件下，研究了活树和死树与某一环境因子的相关程度。环境因子由海拔、坡度、坡向、土壤温度、土壤含水量和土壤电导率组成。土壤温度、土壤含水量和土壤电导率采用的是3个点的平均值。坡向用0—360°表示(0°表示正北方向、90°为正东方向、180°为正南方向、270°表示正西方向)，分析时经过正弦转化。在分析过程中对环境数据进行了标准化处理。

结构等式模型(SEM)所能解决问题的范围和能力比主成分分析、聚类分析等传统多变量统计方法增大和提高了许多，其侧重于确证性检验，更加灵活，能够同时分析多个变量间的因果关系及强度大小^[16]。SEM主要可以处理观察变量、隐变量、综合变量和误差变量4类变量，观察变量指在观察或试验过程中能够直接测量的变量，在图示中用方框来代表；隐变量指无法从观察或实验中直接测量到的变量，用圆形表示；综合变量指多个显变量或隐变量的总体影响，以圆形表示；误差变量指数据分析过程中无法通过预测变量所解释的那部分效应，即残差项，误差变量不带任何边框。在SEM图形框架中，因果关系用(→)表示，箭头起点表示预测变量，箭头所指方向为响应变量。SEM方法主要分为三步：首先基于现有的知识对因果关系进行假定，建立“先验模型”；其次根据实验获得的数据对先验模型进行拟合；最后，对路径进行修改，得到最终模型^{[24, 25]}。

模型的基本假设有：(1)秦岭冷杉的生长受土壤环境的影响，并且活树主要为正相关关系，而枯死树以负相关关系为主；(2)空间位置和地形不但可以影响到土壤含水率、土壤温度和土壤电导率，还可以影响到秦岭冷杉的生长。基于SEM模型，分别分析了土壤(土壤含水率、土壤温度和土壤电导率)、坡度、海拔与活树和死树生长状况的关系。其中树体个体大小的观察变量为DBH和树高。分析前为使变量符合正态分布，所有变量进行了log转换。同时使用卡方统计、拟合优度指数和本特勒-波内特指数来评价模型拟合的优劣。SEM分析过程中，3个评价指标可以通过参数“fit.indices=c("GFI","RMSEA","NFI")得到。

本研究所有分析和绘图均在统计学软件R 2.15.3 (R Development Core Team,http://www.Rproject.org)中进行，其中偏相关分析通过‘ggm’数据包实现，SEM分析通过sem包进行，SEM图形由内置的Graphviz生成^[26]。

对秦岭冷杉普查结果如表 1所示，表 2为秦岭冷杉分布范围。在小秦岭自然保护区，秦岭冷杉有232棵，其中活树仅剩158棵，枯死树有74棵，枯死树占秦岭冷杉总个体数的31.9%。从多年监测和死亡状态上判断枯死树的死亡时间段为近3—4a。枯死树的针叶已全部凋落，生长状态一般的秦岭冷杉生长势也有所减弱，针叶不同程度的出现枯黄，并且结实量较少。图 1为秦岭冷杉活树和枯死树分布图，从图中可看出，枯死树并没有大片呈聚集分布，而是与活树分布基本相一致的，这说明秦岭冷杉的死亡原因不是局部原因造成的。对秦岭冷杉全面普查过程中，很难发现幼苗，个体数极其少。从图 2可看出，活树的胸径主要集中在10 cm到35 cm，树高主要集中在5 m到10 m，也就是说，小树和大树相对较少。枯死树的个体主要集中在胸径0 cm 到30 cm和树高为0 m到10 m。

图 1 秦岭冷杉空间分布图 Fig.1 Abies chensiensis space distribution 图中经纬度为高斯坐标，圆代表胸径大小，灰色代表活树，黑色代表死树

图选项

图 2 秦岭冷杉活树、死树胸径和树高直方图 Fig.2 The distribution structure of DBH and height ofAbies chensiensis of living and dead trees

图选项

从秦岭冷杉群落调查中可以看出，秦岭冷杉多数个体立地条件优越，群落的森林环境主要由秦岭冷杉决定，表 3为秦岭冷杉样地概况。近几年的气象条件变化可能影响秦岭冷杉的生长，小秦岭自然保护区近3a的大气温度、大气湿度、太阳辐射和降雨量资料结果如图 3所示。

图 3 小秦岭自然保护区近3a气象条件变化 Fig.3 The variation of meteorological condition in recent years in Xiaoqinling national nature reserve

图选项

Spearman相关分析结果如表 4所示，活树与海拔、土壤含水量和土壤温度表现出显著正相关或负相关，死树与海拔、坡度和土壤温度表现出显著负相关。由于6个环境因子之间可能存在一定的相互影响，因此偏相关分析更能说明6种因子的真实影响。结果表明(表 5)，活树、死树与海拔和土壤温度都表现出显著相关，死树与坡度和土壤电导率都表现出显著的负相关作用。

根据基本假设建立初始模型(图 4)，依据各变量之间的路径系数和修正指数，对初始模型进行修正，2个最终模型见图 4，图中路径系数为标准化过后的。拟合优度指数、近似误差均方根和本特勒-波内特指数均显示各模型是可以接受的(表 6)。各变量之间的路径系数和显著度见表 7。活树模型中，模型可以解释秦岭冷杉个体大小变异的31.68%，死树模型中，模型可以解释个体大小变异的64.99%。SEM分析结果与相关和偏相关分析结果是相一致的，并且SEM模型可以更好地拟合出各变量间直接、间接的相互作用，有助于更好地认识秦岭冷杉与土壤、地形之间复杂关系。图 4中2个最终模型结果显示，秦岭冷杉活树与土壤环境(土壤含水量和土壤温度)呈明显的正相关系，而死树却表现出负相关关系，这与文章的假设一致。枯死树和活树与空间位置和地形都表现出负相关作用，枯死树与空间位置和地形的路径系数分别为-0.9479和-0.7510，而活树分别为-0.0980和-0.712。

图 4 基于SEM建立的各变量间相互关系修正后的最终模型 Fig.4 The final models of the relationship among variables based on SEM 地形 Topography; 位置 Position; 坡度 Slope; 土壤环境 Soil Environment; 海拔 Elevation; 土壤含水量 Soil Water Content; 土壤温度 Soil Temperature; 空间分布 Spatial Distribution; 经度 Longitude; 纬度 Latitude; 土壤电导率 Soil Electric Conductivity; 个体大小 Individual Size; 胸径 DBH; 树高 Tree Hight

图选项

从秦岭冷杉普查结果可以发现，秦岭冷杉不仅死亡个体数量大，而且活树中很大一部分长势减弱，生长状况一般。从活树和死树分布来看，枯死树分布并没有呈聚集或成片分布，而是与活树分布是相一致的，这说明秦岭冷杉的死亡不是个别现象和局部环境变化所造成的。从秦岭冷杉活树个体的胸径和树高可以看出，幼苗和小树较少，更新能力差，这和许多有关秦岭冷杉种子和结实的研究结果是相一致，和结实率、种子的繁殖特征和鼠类啮食等因素存在一定关系^{[3, 8, 9]}。秦岭冷杉的生长状况和发展趋势不容乐观。

通过秦岭冷杉群落调查发现，秦岭冷杉多数个体立地条件优越，因此不存在因物种竞争而导致物种群体死亡的可能。并且近3年的大气温度、大气湿度、太阳辐射和降水量没有异常变化，也没有发生大规模的自然灾害，因此，气象条件并不是造成秦岭冷杉生长势减弱和死亡的限制性因子。这就很可能与环境因素变化或者人类干扰有关。

从以上分析可以看出，小秦岭地区秦岭冷杉死亡很大可能与环境因素的变化存在关系。张文辉等在研究中也表明环境因素在秦岭冷杉种群生殖过程中发挥了作用^[4]。本研究相关和偏相关分析结果表明秦岭冷杉的生长与海拔、坡度、土壤含水量和土壤温度关系密切。秦岭冷杉的垂直分布跨越很大的海拔范围，不同海拔高度温度和降水等环境因素存在很大差异^[20]。孙玉玲等在陕西秦岭地区研究表明秦岭冷杉球果和种子在不同海拔高度特征存在显著差异^[8]。因此秦岭冷杉的生长受海拔影响比较大。不同坡度的土壤养分和含水量存在很大差异，缓坡养分含量相对于陡坡更丰富^{[21, 22]}，而土壤含水量可以直接影响到秦岭冷杉的生长。并且秦岭冷杉自身对环境因素变化敏感，因此与海拔、坡度、土壤含水量和土壤温度等环境因素关系明显。同时本实验中测定的环境因子可以作为秦岭冷杉就地保护和引种育苗参考指标。

环境因素显著影响着秦岭冷杉的生长，那么，活树和枯死树与周围环境因素的关系有何差异呢？是哪种环境因子造成了秦岭冷杉的死亡？SEM分析结果进一步表明，活树与土壤含水量和温度呈强烈的正相关系，说明土壤含水量和温度对秦岭冷杉生长具有促进作用，而枯死树却表现出负相关作用。王婷等在同区域的伏牛山嵩县研究表明该区域随着5、6月份气温升高，蒸发作用增强，容易造成土壤水分亏缺，从而限制了华山松的生长^[27]，而秦岭冷杉和华山松都属于松科植物。同时王铭等在研究中表明土壤温度是影响土壤呼吸作用的主要因素^[28]。结合秦岭冷杉活树生长势减弱以及枯死树的死亡状态，秦岭冷杉的死亡可能与土壤含水量有着密切的关系，因为保护区建立之前的采矿活动可能造成地下水位下降，加之秦岭冷杉为浅根性物种，从而导致水分成为秦岭冷杉生长的限制性因素。因此，土壤含水量和温度是影响秦岭冷杉的重要环境因素，很有可能是造成秦岭冷杉死亡的一个关键因素。

同时，SEM分析显示枯死树和活树都与海拔和坡度表现出的负相关作用，并且枯死树表现出的负相关作用要明显大于活树。海拔和坡度一定程度上反映了该地点的湿度、温度等微环境条件^[29]，可能是人类活动范围、活动强度和大范围内气候变化逐渐的影响了小秦岭地区降水、温度等环境的变化，从而造成秦岭冷杉逐渐不适应该海拔地段、坡度范围的微环境，从而导致死亡以及活树中部分个体长势减弱。因此秦岭冷杉枯死树和活树都不程度的与海拔和坡度表现出负相关作用，这也进一步证明了两个最终模型的合理性。在今后秦岭冷杉经营保护过程中要注意微环境的变化，减少人为因素对环境的干扰。

影响秦岭冷杉死亡的原因是复杂的而且是一个逐渐的过程。从全国大范围来看，秦岭冷杉死亡的现象在分布地陕西秦岭，神农架等都有报道，特别是河南省鲁山县石人山。这可能与秦岭冷杉分布的地区存在一定关系。秦岭冷杉主要分布于亚热带和暖温带过渡区^[4]，生态过渡区植物群落对全球气候变化反应敏感^[30]，加之人类干扰，从而造成了在秦岭冷杉分布地区都不同程度的出现死亡和生长势减弱现象，这与SEM分析结果也是相一致的。

本研究在小秦岭自然保护区秦岭冷杉普查的基础上，通过对活树、死树分布、胸径和树高的对比分析，并采用相关分析和现代统计分析方法(SEM)结合地形和土壤等环境因子，探讨了秦岭冷杉的生存状况和死亡原因，研究结果对于全面了解小秦岭地区秦岭冷杉生存现状和死亡原因具有重要意义，同时对秦岭冷杉的扦插繁殖、引种育苗等具有参考价值。当然，本研究还有很大的不足之处，所调查的数据没能考虑秦岭冷杉的生理指标等，并且环境因子只能解释秦岭冷杉生长的部分原因，今后研究中应考虑更多的环境因子，比如土壤理化性质和腐殖质等。同时有必要结合石人山、伏牛山、秦岭等地秦岭冷杉的生长状况，在大尺度内研究秦岭冷杉的死亡原因，为秦岭冷杉的保护和种群恢复提供理论和技术上的支持。