扩展功能
文章信息
- 李文娟, 赵小飞, 饶才红, 黄志旁, 任国鹏, 陈尧, 何晓东, 王新文, 白洪兵, 张淑霞, 肖文
- LI Wenjuan, ZHAO Xiaofei, RAO Caihong, HUANG Zhipang, REN Guopeng, CHEN Yao, HE Xiaodong, WANG Xinwen, BAI Hongbing, ZHANG Shuxia, XIAO Wen
- 高黎贡山雉类生境预测
- Prediction on the Suitable Habitat of Four Pheasants Species Using Ecological Niche Modeling in the Gaoligong Mountains
- 四川动物, 2019, 38(3): 328-337
- Sichuan Journal of Zoology, 2019, 38(3): 328-337
- 10.11984/j.issn.1000-7083.20180332
-
文章历史
- 收稿日期: 2018-11-09
- 接受日期: 2019-04-08
2. 普洱市林业局, 云南普洱 665000;
3. 云南高黎贡山国家级自然保护区福贡管护分局, 云南怒江 673400;
4. 云岭省级自然保护区管护局, 云南怒江 671400;
5. 云南高黎贡山国家级自然保护区怒江管护局, 云南怒江 673200;
6. 云南高黎贡山国家级自然保护区泸水管护分局, 云南怒江 673200;
7. 云南高黎贡山国家级自然保护区贡山管护分局, 云南怒江 673500
2. Forestry Bureau of Pu'er City, Pu'er, Yunnan Province 665000, China;
3. Fugong Administrative Sub-Bureau of Gaoligongshan National Nature Reserve, Nujiang, Yunnan Province 673400, China;
4. Yunling Provincial Nature Reserve Administrative Bureau, Nujiang, Yunnan Province 671400, China;
5. Nujiang Administrative Bureau of Gaoligongshan National Nature Reserve, Nujiang, Yunnan Province 673200, China;
6. Lushui Administrative Sub-Bureau of Gaoligongshan National Nature Reserve, Nujiang, Yunnan Province 673200, China;
7. Gongshan Administrative Sub-Bureau of Gaoligongshan National Nature Reserve, Nujiang, Yunnan Province 673500, China
物种的分布信息对野生动物的保护和管理至关重要(高惠等,2017;崔绍朋等,2018)。但调查野生动植物的分布,尤其是稀有物种的分布,通常需要耗费大量的人力、物力和财力(宓春荣等,2017)。而通过模型方法对物种分布进行预测,能快速有效地了解物种分布,明确生境现状,为受胁物种的有效保护提供指导(黄勇杰等,2017;Ren et al., 2017)。生态位模型能够根据物种的出现数据和环境变量在大空间尺度上预测物种生境的分布,而不需要对物种的生理和行为特征进行广泛的调查和详细的描述,已经发展成为预测物种生境分布的有效工具(Wang et al., 2017)。
雉类是全球最重要的鸟类类群之一,被人类广泛驯养、狩猎获取食物和羽毛。然而,猎杀、生境丧失、人为干扰和城市化正威胁着很多雉类物种的生存。自1990年以来,雉类研究取得了巨大进展,但在宏观生态学、分子遗传学和保护上仍然有研究空缺,保护方面亦需要清晰、科学的规划来改善(Tian et al., 2018)。雉类集中分布区的滇西北高黎贡山(Cai et al., 2018),在区域尺度上开展雉类分布预测研究,可为该区域雉类整体性保护提供必要信息。
本研究以高黎贡山4种雉类(白尾梢虹雉Lophophorus sclateri、血雉Ithaginis cruentus、白鹇Lophura nycthemera和白腹锦鸡Chrysolophus amherstiae)的访查分布记录为基础,基于该区域气候变量数据,应用生境偏好、利比希最小因子定律、谢尔福德耐受性定律构建生态位模型,预测其生境分布,并评估生境质量。通过现有保护区的网络覆盖情况,分析4种雉类保护的空缺状况,以期为高黎贡山雉类的生境保护提供整体性规划依据,提高保护效率。
1 材料与方法 1.1 研究区域高黎贡山位于滇西北中缅边境地区,是怒江、独龙江两大水系的分水岭,地理坐标97°30′~99°00′E,24°40′~28°30′N,面积约2.8万km2,是南北动物迁徙扩散的天然通道和东西生物交汇的过渡纽带,有“世界物种基因库”“雉类和鹛类的乐园”“哺乳类动物祖先的发源地”“东亚植物区系的摇篮”之称(熊清华,朱明育,2006),是三大生物多样性热点地区(印缅、中国西南山地和喜马拉雅生物多样性热点地区)的交汇区域(Marchese,2015)。
1.2 分布数据2012年,大理大学东喜玛拉雅研究院在整个高黎贡山怒江段,南起泸水县上江乡大练地村,北至贡山县丙中洛乡秋那桶村,对78个自然村(图 1)进行访查,根据受访者的辨识,现场在Google Earth上勾绘出各雉类的分布范围,共获得白尾梢虹雉49条、血雉106条、白鹇90条、白腹锦鸡115条分布记录。
1.3 气候数据与处理气候数据来自WorldClim网站(http://www.worldclim.org/),空间分辨率为30 s,包括1970—2000年全球每月的最高气温(tmax1: 12)、最低气温(tmin1: 12)、平均气温(tavg1: 12)以及降水量(prec1: 12),共48个变量,计算多年的月平均值(Hijmans et al., 2005)。按研究区域对上述变量进行裁剪,为了去除环境变量间的相关性,突出主要的环境信息,对这些变量进行了主成分分析,筛选出贡献率大于1%的主成分用于物种分布预测模型构建。
1.4 物种分布预测模型构建生境偏好是生境选择的结果,指物种对某些资源高频率的利用,可以用物种出现点和背景点的环境条件对比来反映(Neu et al., 1974;李欣海等,2001;Zhang et al., 2018)。本研究从每个物种的分布范围中,系统抽取2 000个点作为物种出现点,从整个高黎贡山怒江段系统抽取10 000个点作为背景点。从物种出现点和背景点中各随机抽取75%用于物种分布预测模型构建,其余25%用于模型验证。
参照Zhang等(2018)的方法,将建模所用的物种出现点和背景点在环境因子(e)的取值的概率密度函数分别用核密度函数拟合(图 2)。对于任意的环境因子,物种的生境选择系数(habitat selectivity,w)=KDEpresence(e)/KDEbackground(e)。
其中,KDE表示用核密度函数拟合的概率密度函数。一般说来,如果物种对环境因子有选择,则w与e的关系曲线可能有3种形式:S型(图 3:a)、钟型(图 3:b)和Z型(图 3:c)。当w>1时,表示该物种偏好此环境条件;当w<1时,表示物种回避此环境条件;当w=0时,表示环境条件超过了该物种的耐受范围。
本研究以w与e的关系曲线为基础,用分段线性函数来拟合生境质量(H)与环境因子的关系(图 3)。以图 3-e说明图中点A(e1,0)、B(e2,0.5)、C(e3,1)、D(e4,1)、E(e5,0.5)、F(e6,0)的各个坐标值是如何确定的。分别用emin、emax、Q1、Q2、Q3表示用于建模的物种出现点环境因子的最小值、最大值、下四分位数、中位数和上四分位数。按以下步骤(1~4)构建一个物种在环境因子的生境质量模型H=f(e)。
(1) 预先规定,当w=1时,H=0.5,对应图 3-e中的点B(e2,0.5)和E(e5,0.5)。当w=1时,表示物种既不偏好也不回避此环境条件,因而,规定这2个点生境质量为0.5是合理的。注意,因为点B和点E是确定的,所以,e2和e5为已知值。
(2) 点A为生境质量从0到大于0的临界点,令e1=max[emin,Q2-1.5×(Q3-Q1)]。
(3) 点F为生境质量从大于0到等于0的临界点,令e6=min[emax,Q2+1.5×(Q3-Q1)]。
(4) 点A和点B确定了1条直线,点C在这条直线上,可知点C的横坐标e3=2e2-e1;同理,点D的横坐标e4=2e5-e6。
对于互相独立的全部环境因子,物种的生境质量取决于生境质量最低的生态因子(利比希最小因子定律)。因此,本研究以所有生境因子中的最低H值来反映物种在此环境条件下的整体生境质量。当w>1时,H>0.5,定义为高质量生境;当1>w>0时,0.5>H>0,定义为边缘生境;当w=0时,H=0,定义为非生境。
通过计算研究区域每个位置(像元)的生境质量,得到研究区整体的潜在生境分布预测图。
1.5 模型评估本研究采用受试者工作特征(receiver operating characteristic,ROC)曲线下面积(area under curve,AUC)评估模型的准确度(王运生等,2007)。
1.6 保护优先区的确定及级别划分最理想的情况下,高黎贡山的这4种濒危雉类的潜在栖息地都能得到保护。而当资源有限时,应该优先保护多个物种均有分布的区域,尤其是高质量生境分布区。将4种雉类的生境质量评估结果加权后(高质量生境赋值为2,边缘生境赋值为1,非生境赋值为0)进行叠加分析。生境质量累加值大于0的像元(30 m×30 m)从大到小进行排序,按像元数量进行三等分,分别定义为一级、二级和三级保护优先区。
本研究所有数据处理都在R Statistics 3.4.3平台下使用raster(Hijmans,2017)和dismo(Hijmans et al., 2017)等实现。
2 结果 2.1 气候变量的主成分对48个气候变量做主成分分析,前4个主成分的累积贡献率达到了99%。其中,第一主成分(贡献率86%)与每月最高气温、平均气温和降水量正相关;第二主成分(贡献率7%)与8—12月的最低气温负相关;第三主成分(贡献率3%)与4—8月的最低气温负相关;第四主成分(贡献率2%)受月最低气温影响(表 1)。
气候变量 | 第一主成分 | 第二主成分 | 第三主成分 | 第四主成分 |
prec1 | 0.155 | -0.013 | 0.009 | 0.042 |
prec2 | 0.155 | -0.010 | 0.010 | 0.048 |
prec3 | 0.155 | -0.014 | 0.026 | 0.042 |
prec4 | 0.155 | -0.001 | 0.019 | 0.052 |
prec5 | 0.155 | -0.016 | 0.009 | 0.036 |
prec6 | 0.155 | -0.017 | 0.026 | 0.035 |
prec7 | 0.155 | -0.014 | 0.006 | 0.066 |
prec8 | 0.155 | -0.012 | -0.011 | 0.063 |
prec9 | 0.155 | -0.016 | -0.018 | 0.058 |
prec10 | 0.155 | -0.018 | -0.018 | 0.052 |
prec11 | 0.155 | -0.011 | 0.026 | 0.048 |
prec12 | 0.155 | 0.001 | 0.002 | 0.048 |
tmax1 | 0.155 | -0.015 | 0.032 | 0.053 |
tmax2 | 0.155 | -0.025 | -0.001 | 0.068 |
tmax3 | 0.155 | -0.003 | 0.036 | 0.055 |
tmax4 | 0.155 | -0.022 | 0.053 | 0.034 |
tmax5 | 0.155 | -0.033 | 0.029 | 0.049 |
tmax6 | 0.155 | -0.022 | -0.020 | 0.064 |
tmax7 | 0.155 | -0.021 | -0.019 | 0.045 |
tmax8 | 0.155 | -0.025 | 0.053 | 0.027 |
tmax9 | 0.155 | -0.017 | 0.052 | 0.057 |
tmax10 | 0.154 | -0.028 | -0.003 | 0.078 |
tmax11 | 0.154 | -0.019 | 0.053 | 0.042 |
tmax12 | 0.154 | -0.010 | -0.028 | 0.059 |
tavg1 | 0.154 | -0.026 | -0.021 | 0.075 |
tavg2 | 0.154 | -0.035 | 0.052 | 0.053 |
tavg3 | 0.154 | -0.040 | 0.057 | 0.043 |
tavg4 | 0.154 | -0.044 | 0.057 | 0.035 |
tavg5 | 0.154 | -0.037 | 0.057 | 0.051 |
tavg6 | 0.154 | -0.020 | -0.038 | 0.060 |
tavg7 | 0.154 | -0.051 | 0.058 | 0.053 |
tavg8 | 0.154 | -0.055 | 0.057 | 0.045 |
tavg9 | 0.154 | -0.024 | -0.039 | 0.071 |
tavg10 | 0.154 | -0.019 | 0.071 | 0.061 |
tavg11 | 0.154 | -0.048 | 0.058 | 0.061 |
tavg12 | 0.153 | -0.038 | 0.074 | 0.057 |
tmin1 | 0.141 | -0.090 | -0.014 | -0.321 |
tmin2 | 0.140 | -0.057 | 0.061 | -0.393 |
tmin3 | 0.139 | -0.096 | -0.031 | -0.400 |
tmin4 | 0.136 | 0.069 | -0.221 | -0.366 |
tmin5 | 0.123 | -0.204 | -0.135 | -0.436 |
tmin6 | 0.119 | 0.034 | -0.492 | -0.139 |
tmin7 | 0.065 | 0.172 | -0.650 | 0.151 |
tmin8 | -0.054 | -0.325 | -0.418 | 0.306 |
tmin9 | -0.072 | -0.460 | 0.152 | -0.023 |
tmin10 | -0.073 | -0.467 | -0.061 | 0.018 |
tmin11 | -0.086 | -0.425 | -0.062 | -0.041 |
tmin12 | -0.099 | -0.402 | -0.051 | 0.130 |
注:prec1: 12为1—12月每月降水量,tmax1: 12为1—12月每月最高气温,tavg1: 12为1—12月每月平均气温,tmin1: 12为1—12月每月最低气温 Notes:prec1: 12 is the monthly precipitation from January to December,tmax1: 12 is the maximum temperature in each month from January to December,tavg1: 12 is the average monthly temperature from January to December,tmin1: 12 is the minimum temperature in each month from January to December |
白尾梢虹雉、血雉、白鹇、白腹锦鸡预测的AUC值分别为0.92、0.89、0.83、0.80,表明该模型对白尾梢虹雉的预测达优秀水平,对其余3种预测效果较好。
2.3 模型预测结果在高黎贡山区域,白鹇和白腹锦鸡比白尾梢虹雉和血雉具有更大的潜在生境;但后两者具有更多的高质量生境(表 2,图 4)。
生境质量等级 | 生境面积/km2 | |||
白尾梢虹雉 | 血雉 | 白鹇 | 白腹锦鸡 | |
高质量生境 | 2 357 | 4 097 | 772 | 1 021 |
边缘生境 | 4 075 | 4 367 | 8 801 | 12 671 |
总计 | 6 432 | 8 464 | 9 573 | 13 691 |
2.4 生态位分化
主成分分析发现,白尾梢虹雉和血雉的生境有一定的相似性,表明存在一定程度的生态位重叠;与白腹锦鸡和白鹇的生境存在一定的差异性(图 5),说明它们之间存在一定的生态位分化,其生境预测结果也很好地说明了这点。
2.5 保护优先区分级一级保护优先区集中分布于高黎贡山北段贡山县高黎贡山东西坡、丙中洛区域,在福贡县有零星分布。二级保护优先区主要分布于高黎贡山中段泸水县和福贡县。高黎贡山北段是4种雉类共存的潜在生境分布区,从高黎贡山雉类保护整体规划角度,认为高黎贡山北段是雉类保护的优先区域(图 6)。
3 讨论 3.1 本研究生态位模型特色生态系统内输入太阳辐射的多少和气温的高低、降水的多少和潜在蒸腾作用的大小造成了物种丰富度空间分布格局的差异(Hawkins et al., 2005;McCain, 2007, 2009)。在较大的空间尺度上,如区域或山系尺度,气候条件限制了物种基础生态位的分布,直接影响着物种的分布和生境质量,因此气候变量在物种生境预测中具有重要作用(Ren et al., 2017)。通过主成分分析,本研究在最大化的保留原有多个环境因子信息的同时,又使进入模型的环境因子(主成分)之间互相独立,从而可以根据利比希最小因子定律获得任何地点的生态质量评估体系。本研究构建模型的过程和张桂铭等(2013)、Zhang等(2018)比较类似,但本研究还根据物种/环境因子上的响应曲线形态(钟型、Z型、S型)构建了不同的分段线性模型。本研究将生境选择系数大于1界定为高质量生境,小于1界定为边缘生境,这与生境偏好的定义是相符合的。与张桂铭等(2013)不同的是,我们认为,对于特定的环境因子,每个物种都有一个最适的区间(图 3-e中的平台部分),而不是仅有一个最适的值(点)。使用独立的测试数据评估模型,基于本研究构建的生态位模型的白尾梢虹雉、血雉、白鹇、白腹锦鸡的生境预测效果均达到较好效果,说明本研究构建模型的方法是合理可行的。在未来的研究中,我们将尝试将此物种分布预测模型用于更多的物种以测试该方法的适用性和稳定性。
3.2 高海拔分布的物种在高黎贡山具有较大的高质量生境面积本研究预测结果发现血雉和白尾梢虹雉比白鹇和白腹锦鸡具有更大的高质量生境面积。血雉和白尾梢虹雉分布的海拔较高,主要在高山针叶林、高山灌木丛,白鹇和白腹锦鸡则主要分布于中、低海拔的常绿阔叶林、落叶阔叶林中(杨岚,1995;Cai et al., 2018)。在云南,人类对于中、低海拔动物生境的破坏远大于高海拔人迹罕至的区域(Wu et al., 2017),因此在高黎贡山存留的高海拔生境面积实际亦可能大于中、低海拔的。这也从另一个角度证明了本研究预测的可靠性。
主成分分析发现,白尾梢虹雉和血雉的生境选择有一定的相似性。在第一主成分上,偏向于每月降水量、最高气温、平均气温的最小值;在第二主成分上,偏向于8—12月每月最低气温的最小值,而这2种雉类在高黎贡山高海拔地带确实同域分布,有混群取食现象(罗旭等,2016)。
3.3 保护空缺分析白尾梢虹雉为世界自然保护联盟(IUCN)濒危物种红皮书中的易危(VU)物种,受生境丧失、偷猎等威胁,亟需开展保护研究。本研究发现在高黎贡山北段和中段有该物种大面积成片分布的潜在生境,韩联宪等(2004)对云南白尾梢虹雉的分布调查中亦在高黎贡山北段贡山县实地观察到白尾梢虹雉;但目前针对该物种的生态调查研究集中在高黎贡山南段的腾冲县(施晓春,1999;罗旭等,2004)。未来的白尾梢虹雉种群数量调查应重点在高黎贡山北段和中段优先开展,以确保所有野外种群均被发现和保护。
在雉类的整体保护规划上,高黎贡山北部4种雉类共存的可能性较大,明显具有保护的优先性。而北部地区由于交通原因不易实施保护项目,提示未来的高黎贡山的雉类保护项目应以高黎贡山北段为优先区域。
4 结论在高黎贡山区域尺度上,本研究使用基于生境选择和谢尔福德耐受性定律构建的生态位模型预测4种雉类的潜在分布区。结果发现,高黎贡山北段是4种雉类共存的潜在分布区。从高黎贡山雉类保护整体规划角度,认为高黎贡山北段是雉类保护,尤其是白尾梢虹雉保护的优先区域,建议加大该区域的雉类调查和保护力度。
崔绍朋, 罗晓, 李春旺. 2018. 基于Maxent模型预测白唇鹿的潜在分布区[J]. 生物多样性, 26(2): 171–176. |
高惠, 滕丽微, 汪洋, 等. 2017. 阿拉善马鹿(Cervus alashanicus)生境适宜性评价[J]. 生态学报, 37(11): 3926–3931. |
韩联宪, 黄石林, 罗旭, 等. 2004. 云南白尾梢虹雉的分布与保护[J]. 生物多样性, 12(5): 523–527. DOI:10.3321/j.issn:1005-0094.2004.05.008 |
黄勇杰, 卢佳斌, 王锋堂, 等. 2017. 基于MaxEnt模型预测海南岛海南臭蛙的潜在地理分布[J]. 动物学杂志, 52(1): 30–41. |
李欣海, 马志军, 李典谟, 等. 2001. 应用资源选择函数研究朱鹮的巢址选择[J]. 生物多样性, 9(4): 352–358. DOI:10.3321/j.issn:1005-0094.2001.04.005 |
罗旭, 吴太平, 黄安琪. 2016. 云南高黎贡山两种同域分布雉类冬季的食物分析及取食策略[J]. 生态学杂志, 35(4): 1003–1008. |
罗旭, 韩联宪, 艾怀森. 2004. 高黎贡山冬季白尾梢虹雉运动方式和生境偏好的初步观察[J]. 动物学研究, 25(1): 48–52. DOI:10.3321/j.issn:0254-5853.2004.01.009 |
宓春荣, 郭玉民, Huettmann F, 等. 2017. 基于物种分布模型的精确采样提高目标物种发现率——以黑颈鹤(Grus nigricollis), 白头鹤(Grus monacha)为例[J]. 生态学报, 37(13): 4476–4482. |
施晓春. 1999. 高黎贡山白尾梢虹雉生境利用的初步观察[J]. 动物学研究, 20(1): 50–54. DOI:10.3321/j.issn:0254-5853.1999.01.011 |
王运生, 谢丙炎, 万方浩, 等. 2007. ROC曲线分析在评价入侵物种分布模型中的应用[J]. 生物多样性, 15(4): 365–372. DOI:10.3321/j.issn:1005-0094.2007.04.005 |
熊清华, 朱明育. 2006. 高黎贡山周边社区研究(精)[M]. 北京: 科学出版社. |
杨岚. 1995. 云南鸟类志上卷非雀形目[M]. 昆明: 云南科技出版社. |
张桂铭, 朱阿兴, 杨胜天, 等. 2013. 基于核密度估计的动物生境适宜度制图方法[J]. 生态学报, 33(23): 7590–7600. |
Cai T, Fjeldså J, Wu Y, et al. 2018. What makes the Sino-Himalayan mountains the major diversity hotspots for pheasants?[J]. Journal of Biogeography, 45(3): 640–651. DOI:10.1111/jbi.2018.45.issue-3 |
Hawkins BA, Diniz-Filho JA, Soeller SA. 2005. Water links the historical and contemporary components of the Australian bird diversity gradient[J]. Journal of Biogeography, 32(6): 1035–1042. DOI:10.1111/jbi.2005.32.issue-6 |
Hijmans RJ. 2017. Raster: geographic data analysis and modeling. R package version 2.6-7[EB/OL].[2018-12-20]. https://CRAN.R-project.org/package=raster. |
Hijmans RJ, Phillips S, Leathwick J, et al. 2017. dismo: species distribution modeling. R package version 1.1-4[EB/OL].[2018-12-10]. https://CRAN.R-project.org/package=dismo. |
Hijmans RJ, Cameron SE, Parra JL, et al. 2005. Very high resolution interpolated climate surfaces for global land areas[J]. International Journal of Climatology, 25(15): 1965–1978. DOI:10.1002/(ISSN)1097-0088 |
Marchese C. 2015. Biodiversity hotspots:a shortcut for a more complicated concept[J]. Global Ecology and Conservation, 3: 297–309. DOI:10.1016/j.gecco.2014.12.008 |
McCain CM. 2009. Global analysis of bird elevational diversity[J]. Global Ecology and Biogeography, 18(3): 346–360. DOI:10.1111/geb.2009.18.issue-3 |
McCain CM. 2007. Could temperature and water availability drive elevational species richness patterns? A global case study for bats[J]. Global Ecology and Biogeography, 16(1): 1–13. DOI:10.1111/geb.2007.16.issue-1 |
Neu C, Byers C, Peek J. 1974. A technique for analysis of utilization-availability data[J]. The Journal of Wildlife Management, 38(3): 541–545. DOI:10.2307/3800887 |
Ren G, Yang Y, He X, et al. 2017. Habitat evaluation and conservation framework of the newly discovered and critically endangered black snub-nosed monkey[J]. Biological Conservation, 209: 273–279. DOI:10.1016/j.biocon.2017.02.029 |
Tian S, Xu JL, Li JQ, et al. 2018. Research advances of Galliformes since 1990 and future prospects[J]. Avian Research, 9: 32. DOI:10.1186/s40657-018-0124-7 |
Wang B, Xu Y, Ran JH. 2017. Predicting suitable habitat of the Chinese monal (Lophophorus lhuysii) using ecological niche modeling in the Qionglai Mountains, China[J]. Peerj, 5(7): e3477. |
Wu F, Liu L, Fang J, et al. 2017. Conservation value of human-modified forests for birds in mountainous regions of south-west China[J]. Bird Conservation International, 27(2): 187–203. DOI:10.1017/S0959270916000058 |
Zhang G, Zhu AX, Windels SK, et al. 2018. Modelling species habitat suitability from presence-only data using kernel density estimation[J]. Ecological Indicators, 93: 387–396. DOI:10.1016/j.ecolind.2018.04.002 |