四川铁布梅花鹿自然保护区梅花鹿夏季栖息地选择

赵成， 李艳红， 石琴， 王彬， 3， 徐康宁， 赵思勤， 胡杰*， 熊远清

(1.宜宾学院生命科学与食品工程学院，四川宜宾644000； 2. 西南野生动植物资源保护教育部重点实验室， 西华师范大学珍稀动植物研究所，四川南充637002； 3. 西华师范大学生态研究院，四川南充637002； 4. 四川省若尔盖县林业局，四川若尔盖624500)

梅花鹿Cervusnippon隶属偶蹄目Atiodactyla鹿科Cervidae，国家Ⅰ级重点保护野生动物，历史上几乎遍布我国华东、华北、华中、华南、西南和青藏高原的东部(Whitehead，1993；郭延蜀，郑慧珍，2000)，由于气候环境的改变以及人为影响，目前中国仅存 3个亚种：东北亚种C.n.hortulorm、华南亚种C.n.kopschi和四川亚种C.n.sichuanicus，分布于黑龙江、吉林、辽宁、安徽、浙江、江西、四川和甘肃迭部交界区域(郭延蜀，2000；宋延龄，刘志涛，2005)；其中，四川铁布梅花鹿自然保护区分布有我国现存最大的野生种群(戚文华等，2014)。

栖息地是野生动物赖以生存的场所，了解和掌握濒危物种栖息地的利用特征是制定保护和管理策略的基础(孟秀祥等，2011；Zhangetal.，2011；Khadka & James，2016)。迄今为止，国外学者对梅花鹿栖息地利用方面的研究成果较为丰富(Asada & Ochiai，1996；Borkowski，2000；Yokoyamaetal.，2000； Sakuragietal.，2003；Diazetal.，2017；Endoetal.，2017)；国内关于梅花鹿栖息地利用方面的研究相对较少，且主要集中于华南亚种(付义强等，2006；李佳等，2015；Lietal.，2017)，针对四川亚种栖息地利用方面的研究较少(Zhaoetal.，2014)。

由于环境的差异，不同地理种群的梅花鹿所面临的生存压力也不一样，影响其栖息地选择的主要因素也存在一定差异。例如，生活于俄罗斯西部的梅花鹿主要选择的生境为蒙古栎Quercusmongolica林、落叶阔叶林和红松Pinuskoraiensis针阔混交林(Aramilev，2009)；分布于吉林珲春东北虎国家级自然保护区的东北亚种却不偏好红松针阔混交林和蒙古栎林(杨海涛等，2018)；而分布于江西的华南亚种不偏好森林，而偏好灌草丛(付义强等，2006；Lietal.，2017)。鉴于此，作者对四川铁布梅花鹿自然保护区梅花鹿夏季的栖息地选择进行了研究，以期进一步丰富四川亚种的生态生物学资料，揭示影响其生境选择的关键因素，并为其保护与管理提供参考。

1 研究区域与研究方法

1.1 研究区域概况

四川省若尔盖县铁布梅花鹿自然保护区(102°52′～103°08′E，33°58′～34°11′N)位于青藏高原东沿，总面积260 km2。区内地貌属中切割山原，谷底海拔2 400 m，峰岭海拔3 000 m左右，最高可达3 850 m。受西风环流及东南季风的影响，夏季温凉、冬季寒冷、干湿季明显，年均降水量625.5 mm。区内植被垂直分布明显，且受坡向影响强烈。河岸及河漫滩地带以柳Salixspp.为主，其间分布有沙棘Hippophaerhamnoides、小檗Berberisspp.、高丛珍珠梅Sorbariaarborea等的河滩灌丛；阴坡以油松Pinustabulaeformis林、粗枝云杉P.asperata林、紫果云杉P.purpurea林、岷江冷杉Abiesfaxoniana林、白桦Betulaplatyphylla林和山杨Populusdavidiana林为主；阳坡是由小檗、沙棘等组成的灌丛草甸和山坡灌丛(郭延蜀，2000)。森林、灌丛、草甸构成的复合生境为梅花鹿生命活动的各个阶段提供了适宜的生存环境(Zhaoetal.，2014)。

1.2 调查方法

2011年6—9月在保护区设置20条垂直于等高线的样线，每条长2 km，间距大于2 km，且尽可能覆盖区内各种植被类型。在满足以下条件时设置样方：1)样线的起始点与终点，2)发现梅花鹿实体或新鲜粪便、尸体等痕迹处，3)植被类型发生明显变化，4)若无上述情况，则海拔每上升100 m或行走1 h设置1个样方(张广文等，2007)。

样方设置参考Wei等(2000)的方法：1个1 m×1 m的正方形，1个20 m×20 m的正方形和2个相互垂直的2 m×10 m的长方形，另在20 m×20 m正方形样方中每1/4小样方(10 m×10 m)的中心分别设1个1 m×1 m的正方形样方。每个样方记录如下信息：

坐标：用GPS记录20 m×20 m的正方形样方中心的地理坐标。

样方类型：分为1和0，1表示发现梅花鹿实体(状态：觅食、卧息等)或其粪便、食痕等的样方，0表示没发现梅花鹿及其痕迹的样方(对照样方)。

坡向：东、南、西、北和无。

坡位：上坡、中坡和下坡。

坡度：用坡度仪直接测量。

海拔：在大样方的中心点用GPS测量。

距水源距离：样方中心距最近水源的距离，根据GPS定位，后期结合地形图和卫星影像图测量。

距林缘距离：样方中心距林缘的垂直距离，从林缘到林内取“+”值，从林缘到林外取“-”值。

隐蔽度：在每个大样方的对角线上，相距30 m处，各执一花杆(高1.5 m)，对视对方花杆，记录所见比例，取2次的平均数。分为4级：<25%(所见杆比例>75%)、25%～50%(所见杆比例为50%～75%)、50%～75%(所见杆比例为25%～50%)、>75%(所见杆比例<25%)。

植被类型：针叶林、阔叶林、灌丛、草甸、耕地和其他。

乔木密度：大样方内乔木的数量。

乔木高度：大样方内所有乔木的平均高度。

乔木胸径：测量4个10 m×10 m正方形样方中距中心点最近的乔木胸径，并计算其平均值。

灌木密度：2个2 m×10 m长方形样方的灌木平均数量。

灌木高度：2个2 m×10 m长方形样方的灌木平均高度。

灌木盖度：2个2 m×10 m长方形样方的灌木平均盖度。

草本高度：5个1 m×1 m正方形样方的草本平均高度。

草本盖度：5个1 m×1 m正方形样方的草本平均盖度。

距居民区距离：样方中心距最近居民区的距离，根据GPS定位，后期结合地形图和卫星影像图测量。

距大路距离：距机动车行驶道路的距离，根据GPS定位，后期结合地形图和卫星影像图测量。

1.3 数据分析

针对植被类型、隐蔽度、坡度、坡向4个离散型变量，采用卡方检验分析利用样方与对照样方的差异性。

通过单样本K-S检验分析数值型变量是否符合正态分布，如果符合正态分布，则使用独立样本t检验分析利用样方和对照样方的差异性，否则，使用Mann-WhitneyU检验。

通过信息理论方法(Burnham & Anderson，2002)预测梅花鹿栖息地选择的关键因子。为保证变量的独立性，对所有变量进行相关分析，对相关系数绝对值>0.5的变量进行筛选，保留其中更具生物学意义的变量进入下一步分析。使用赤池信息量准则(Akaike information criterion，AICC)评估所有可能的变量组合，并通过ΔAICC值进行模型筛选；对ΔAICC<2的模型进行模型平均，然后比较每个变量的相对重要性和95%置信区间，以此确定关键的影响因子(Wangetal.，2018)。数据分析在SPSS 20和R中完成。

2 结果

2.1 利用样方和对照样方各因子的差异显著性分析

对保护区梅花鹿夏季利用样方(61个)和对照样方(48个)中4个离散型变量分析的结果表明，植被类型和隐蔽度之间的差异有高度统计学意义(P<0.01)，而坡向和坡位之间的差异无统计学意义(P>0.05)(表1)。梅花鹿偏好灌丛、草甸、针叶林和隐蔽度较高的生境。

表1 铁布梅花鹿自然保护区梅花鹿夏季利用样方 和对照样方中离散型变量比较Table 1 Discrete variables in used quadrats and control quadrats for Cervus nippon sichuanicus in the Tiebu Sika Deer Nature Reserve during summer

对14个连续性变量分析的结果表明，距水源距离、灌木高度、灌木盖度、灌木密度、草本高度和草本盖度6个变量在利用样方和对照样方之间的差异有统计学意义(F<0.05或U<0.05)，其他变量之间的差异无统计学意义(F>0.05或U>0.05)(表2)。

2.2 逻辑斯蒂回归模型

相关分析结果发现，灌木盖度和灌木密度之间的相关系数较大(r=0.838，P<0.001)。鉴于灌木密度既能反映样方的隐蔽度又能代表样方的食物量情况，仅将灌木密度纳入模型进行分析。各回归模型中ΔAICC<2的有3个，共涉及11个变量(表3)。对ΔAICC<2的模型进行模型平均，结果表明，影响梅花鹿夏季栖息地选择的主要因子为草本盖度、距林缘距离和距水源距离；且梅花鹿偏好草本盖度大、靠近水源和林缘的栖息地，回避远离水源、林缘的生境。

3 讨论

研究表明，食物是影响动物栖息地选择的最重要因素(李伟等，2010；赵成等，2012；郑维超等，2012；王征，张旭晖，2014；刘鹏，张薇薇，2017)。郭延蜀(2001)研究发现，草甸和灌丛为梅花鹿四川亚种的主要食物基地。许多鹿科动物都会选择隐蔽性较好的生境以避开天敌(刘振生等，2004；游章强等，2014)。本研究表明，梅花鹿通常偏好灌木和草本繁茂、隐蔽度较高的生境。

表2 铁布梅花鹿自然保护区梅花鹿利用样方和对照样方中连续型变量比较Table 2 Continuous variables in used quadrats and control quadrats for Cervus nippon sichuanicus in the Tiebu Sika Deer Nature Reserve during summer

表3 最优逻辑斯蒂回归模型集合Table 3 Summary of the best logistic regression analysis

注： 植被类型为分类变量， 包含针叶林、阔叶林、灌丛、草甸、耕地和其他

Note： vegetation types are categorical variables， including coniferous forests， broadleaved forests， shrubs， meadows， cultivated land， and others

水源对野生动物的分布有明显的影响，鹿科动物多偏好在距水源较近的区域活动(胡杰等，2018)，梅花鹿也不例外(付义强等，2006；Zhaoetal.，2014)。本研究结果显示，梅花鹿偏好选择距离水源较近的生境，这进一步印证了水源对梅花鹿分布的重要性，即便是在雨水相对丰富的夏季，梅花鹿依然偏好近水生境。郭延蜀(2000)曾指出，梅花鹿四川亚种多出现在距离水源30～550 m的范围内，本研究中梅花鹿夏季利用样方距离水源324.67 m±247.68 m，与其类似。此外，作为典型的林缘动物，梅花鹿利用的栖息地也相对靠近林缘地带。

表4 模型平均结果Table 4 Model-averaged coefficient

注： 植被类型(针叶林)作为固定参考因子

Note： regard vegetation type (coniferous forest) as a fixed reference factor

研究表明，除了食物、隐蔽地和水源，人为干扰也是影响野生动物栖息地利用的重要因素(刘振生等，2004；冉江洪等，2006；游章强等，2014；Zhaoetal.，2017；胡杰等，2018)。本研究显示，人为干扰在梅花鹿利用样方和对照样方之间的差异无统计学意义，未进入最优模型，可能有两方面因素：一方面，区内道路和居民点多依水而建，梅花鹿对水的需求致使其无法长期回避人为活动频繁的区域；另一方面，区内尚有大量居民从事生活生产活动，人鹿关系密切，梅花鹿已适应一定程度的人为干扰。

本研究结果表明，分布于四川铁布梅花鹿自然保护区的梅花鹿四川亚种在夏季偏好选择隐蔽度高、距离水源和林缘较近、灌木和草本生长繁茂的生境。鉴于此，加强管护靠近水源和林缘的灌草丛，是保护梅花鹿四川亚种栖息地的关键。