海南热带雨林特有兰科植物华石斛潜在适宜生境预测

宁瑶, 昌秋霞, 张哲, 赵莹*

(1. 热带特色林木花卉遗传与种质创新教育部重点实验室，海南省热带特色花木资源生物学重点实验室，海南大学林学院，海口 570228；2. 广西壮族自治区中国科学院广西植物研究所，广西 桂林 541006)

自18 世纪工业革命以来，社会经济快速发展、人类无节制地开发自然资源等，大气中的CO2、CH4和N2O 等温室气体含量迅速增长，全球气候呈现变暖趋势。全球气候变暖是21 世纪生物多样性保护的重大挑战之一[1]。气候变暖可能带来物种的分布格局[2]和物候期[3]改变、物种丰富度减少[4]、增加外来植物入侵风险[5]、物种灭绝或加速灭绝[6]等不利影响。研究在气候变化背景下的物种潜在地理分布对揭示物种分布格局的形成、迁移以及制定生物多样性保护策略具有非常重要的意义[7]。

物种分布模型(species distribution models, SDMs)是利用某一特定物种的地理分布数据、气候数据及地形数据等相关数据，通过回归、分类、机器学习等多种算法对物种的生态位进行度量[8–11], 再将模型投影至当前或未来的环境情景中，以分布概率的形式反映特定生境对于该物种的适宜程度, 可以解读为物种出现的概率或生境对物种的适宜度[6]。SDMs 在保护生物多样性[12–14]、入侵生物预警机制构建[15]、探究环境变化对生物的影响[16–17]等方面具有重要的应用价值。在众多模型中，最大熵(MaxEnt)模型具有操作简易，样本量对模拟结果影响小，且在大多数情况下模拟精度高于其他模型的特点，成为物种分布预测研究中应用最广泛的模型之一[18–19]。MaxEnt 模型已成功应用于珍稀濒危物种，如岛杨梅(Canacomyricamonticola)[17]、双色万代兰(Vanda bicolor)[20]、青钱柳(Cyclocaryapaliurus)[21]、蒙古扁桃(Amygdalusmongolica)[22]等生境分布预测研究中,并取得较好效果。当研究物种的分布点数据很少(<10)时，MaxEnt 仍然能够产生较高精度的结果[18,20]。目前，虽然物种分布模型已广泛用于动植物的潜在生境预测和分布格局变化研究中，但大多集中于广布种或分布点数据较多的物种，而对于具有明显狭域性和地域性的特有种，尤其是单一生境或分布点较少的特有种鲜有报道[23]。因此，开展对于狭域种、特有种的适生区预测及影响其分布格局的主要环境因子的研究，不仅有利于探明物种适应单一生境的机制，而且对于物种分布模型的发展具有推动作用。

华石斛(Dendrobiumsinense)隶属兰科(Orchidaceae)石斛属，是多年生草本花卉，热带雨林特有兰科植物，目前发现其仅分布于海南岛热带雨林国家公园内海拔1 000 m 以上的山地常绿林或高山矮林的乔灌木树干或树梢上[24–25]，分布范围极其狭窄。华石斛具有较高的药用价值和观赏价值，是海南重要的传统黎药，此外它还是优良的培育长花期品种的育种材料，常作盆栽观赏[24,26]。作为典型的附生植物，华石斛通常附生在离地面2.5～10.0 m 的树干或一级树杈上，生长所需水分及矿物质主要来源于空气中，表明华石斛对于气候与环境具有较高的依赖性[24]。此外，持续高温会抑制华石斛正常的生理代谢，导致叶片萎蔫、失绿、掉落等，推测其对持续高温或极端高温抵御能力差[26]。目前华石斛的研究主要集中在组织培养[27]、菌根生物学[28–29]、种群动态[30–31]、遗传多样性[32]、种子活力[33]、化学成分[34]、高温胁迫[26]等方面，而华石斛种群响应气候变化的研究鲜见报道[35]。明确华石斛的潜在适宜生境分布，以及限制其分布的主导环境因子成为华石斛保育工作和资源可持续利用的关键。为此，本研究利用野外实地调查分布数据，选取气候数据及地形数据作为预测因子，基于MaxEnt 模型和GIS 技术，估测当前、2050 和2070 年3 个时期气候情景下华石斛在海南岛的潜在分布范围，探究华石斛对气候变化的响应，同时探讨影响该物种分布格局的主要环境因素，旨在为华石斛种质资源的保护与合理利用、人工引种栽培及重引入工作提供理论依据。

1 材料和方法

1.1 数据来源与预处理

分布数据来源华石斛的25个分布点数据全部来自海南岛(18°10′～20°10′ N, 108°37′～111°03′ E)野外实地调查，利用GPS (global positioning system)全球定位系统获取经纬度，主要分布于海南热带雨林国家公园霸王岭、尖峰岭、鹦哥岭、五指山、吊罗山、黎母山及佳西分局内(图1)。将所有分布点的经度、纬度、海拔、编号及地名输入Excel 表格并保存为“.CSV”格式，以便MaxEnt 模型读取数据。

图1 华石斛野外分布Fig. 1 Distribution of Dendrobium sinense

环境数据来源参与模型模拟的环境变量包括地形因子和气候因子，地形数据来源于中国科学院的地理空间数据云平台(http://www.gscloud.cn)。地形因子包含海拔、坡度和坡向3 个变量，空间分辨率为90 m 的DEM 数据。因获取的DEM 数据是分幅的，通过Arcgis 数据管理工具中的“融合”工具(Masaic To New Raster)拼接为一整幅栅格图层。其次，通过Arcgis 空间分析模块工具可将海南岛行政区划矢量图做腌膜，裁剪得到海南岛海拔、坡度、坡向3 个栅格图层。气候数据来源于世界气候数据库(http://www.worldclim.org/)，版本为2.0，分辨率为30″ (约1 km)。气候数据包中包含了19 个生物气候变量，如环境变量间出现多重共线性的情况，会引起模型的过度拟合，最终导致模型无法正常运行或预测结果不准确。因此，需对所有气候变量进行Pearson 相关性分析，对于相关性较高(|r|>0.8)的变量，仅保留1 个，结果得到bio2、bio3、bio5、bio7、bio12、bio15、bio16 及海拔、坡度、坡向参与模型模拟(表1)。最后，对所有环境图层进行标准化设置，即统一坐标系为WGS-84，统一研究区域边界，统一分辨率为100 m。

表1 用于模拟物种潜在适宜生境分布的10 个环境变量Table 1 Ten Environmental variables used for modeling species potential suitable habitat distribution

1.2 方法

1.2.1 模型构建

预测物种潜在适宜生境，需要针对不同的建模目标和物种的数据特征如物种分布样本量、空间分布特征进行选取预测性能较优的物种分布模型[36–37]。前人将MaxEnt 模型成功应用于岛杨梅[17]、双色万代兰[20]、青钱柳[21]、蒙古扁桃[22]等珍稀濒危植物及海南凤仙花(Impatienshainanensis)[23]等狭域植物的研究中，并取得较好效果，表明该模型输出结果稳定且可信度较高，因此，本研究选择MaxEnt 模型用于华石斛的潜在适宜生境预测。

将华石斛的分布点数据和环境变量数据导入软件MaxEnt 3.3.3k 中，选中“创建反应曲线”和“刀切法检验”，选择交叉验证法，随机抽取60%的样点作为训练数据用于建模，40%的样点作为检验数据用于验证，模拟过程重复10 次，其他参数选择默认设定，最终得到10 套华石斛适宜生境预测模拟数据。模型将会生成0～1 的估测值，0 表示物种出现的最低概率，即完全不适宜物种生长，1 代表物种出现的最高概率，即非常适宜物种生长。最后，在Arcgis 10.2平台中进行模型输出结果的可视化绘图与处理分析。

1.2.2 适宜生境重分类

将模型模拟结果导入Arcgis 10.2中转为栅格数据图层，再运用空间分析模块的“重分类工具”，采取自然间断点分级方法，将生境之于物种的适宜程度划分为5 个等级：非常不适宜、低度适宜、中度适宜、高度适宜和非常适宜。在预测结果中剔除不可能适宜物种生存的土地类型区域，如农田、水域、城镇及农村居住区等。

所生成的适宜生境分布图坐标系为WGS 1984,无法对各适宜区范围进行面积计算，可利用Arcgis中的数据管理工具的“投影变换”功能，打开适宜生境分布栅格图层的属性表，添加“面积”字段，利用“字段计算器”，一个栅格单元的面积大小与各等级生境所占栅格数量的乘积即为各适宜等级的生境总面积。

1.2.3 模型精度评估

受试者工作特征曲线下面积(area under the receiver operating characteristic curve, AUC)、Kappa系数、真实技巧统计值(true skill statistic, TSS)及总体精度是4 种常用于评估物种分布模型精度的检验指标[38]。AUC 是概率值，其值域为0～1，通常AUC<0.7，模型预测效果很差；0.70.8 时，模型结果较为理想；TSS<0.4 时，模型结果较差。本研究选取这3 个指标用于评估模型的精度。

2 结果和分析

2.1 当前潜在适宜生境预测

MaxEnt 模型基于25 个实际分布点对华石斛进行了预测模拟，运用刀切法检验, 3 个评估模型精度的指标分别为：AUC=0.997±0.003，Kappa=0.990±0.03，TSS=0.990±0.03，表明模型具有较高精度, 预测结果的准确度较高。

从图2 可见，华石斛在海南岛的总适宜生境面积为403.4 km2，占海南岛总面积的1.1%，其中非常适宜的生境面积为50. 9 km2，高度适宜为46.5 km2,中度适宜为75.9 km2，低度适宜为230.1 km2和非常不适宜为34 996.6 km2。华石斛的适宜生境呈现严重的生境斑块化，各斑块间地理隔离严重。华石斛的潜在适生区主要包括：横跨海南热带雨林国家公园霸王岭分局、鹦哥岭分局、尖峰岭分局、吊罗山分局以及五指山分局管辖区域。此外，黎母山分局管辖区、五指山市南部、琼中县、万宁市、保亭县和三亚市等地有零星分布。其中猕猴岭东北部及位于猕猴岭南部的毫肉岭、豪雷岭和鹦哥岭南部的蛙岭、白石岭及石板岭存在大片未被实地调查的非常适宜和高度适宜生境。

图2 华石斛在海南岛的潜在适宜生境分布图Fig. 2 Potential suitable habitat of Dendrobium sinense in Hainan Island

2.2 未来气候情景下华石斛分布变化分析

相比当前，未来2050 和2070 年在RCP4.5 和RCP8.5 气候情景下，华石斛的潜在适宜生境明显减少(表2, 图3)，适宜生境分布区块呈现破碎化，连续性不高，分布格局呈现向分布区块中心缩减的趋势。非常适宜和高度适宜生境全都分布在海南热带雨林国家公园边界内，中度适宜生境和低度适宜生境除少部分零星散落外，也几乎都分布在海南热带雨林国家公园边界内。这说明目前华石斛的适宜生境在保护区内能够得到有效保护，应当重视已发现种群及其潜在适宜生境的就地保护。从时间上来看，从当前到2050 年再到2070 年，当前华石斛的适宜生境面积最大，2070 年适宜生境相比2050 年基本保持不变。

表2 华石斛潜在适宜生境面积(km2)的变化Table 2 Changes in potential suitable habitat area (km2) of Dendrobium sinense

图3 华石斛2050s 和2070s 在2 种气候情景(RCP4.5 和RCP8.5)下的适宜生境分布图Fig. 3 Suitable habitat distribution of Dendrobium sinense under two climate scenarios (RCP4.5 and RCP8.5) in 2050s and 2070s

2.3 环境因子的影响

从华石斛模型预测的环境因子贡献率(表3, 图4)来看，海拔(48.0%)、降水量季节性变异系数(39.7%)、最湿季降水量(8.3%)、最暖月最高温(1.2%)的贡献率最大，累积贡献率达97.2%，因此影响华石斛分布的主要环境变量是海拔、降水量季节性变异系数、最湿季降水量和最暖月最高温4 个环境因子。为了阐述华石斛当前潜在适宜生境分布区的环境因子特征，确定其适宜范围，本研究选择这4 个主导环境因子，根据模型输出结果中的单因子响应曲线，确定存在概率>0.3 的范围为该环境因子对华石斛的相应适宜区间，得到适宜海拔为1 070～1 600 m,降水量季节性变异系数为0.8～0.88，最湿季降水量为1 700～2 266 mm，最暖月最高温为18 ℃～25 ℃。因此，华石斛适宜生长在季节性降水变化大、较为凉爽、暖季气温不高且海拔在1 000 m 以上的高山湿润环境。

表3 MaxEnt 模型预测的华石斛环境变量贡献率Table 3 Contribution rate of environment variables for Dendrobium sinense by MaxEnt

图4 华石斛的预测变量相对重要性的刀切法检验。ALL: 包含所有变量; slo、ele、bio7、bio5、bio3、bio2、bio16、bio15、bio12、asp 见表1。Fig. 4 Results of jackknife evaluations of relative importance of predictor variables for Dendrobium sinense Maxent model. ALL: With all variables; slo, ele,bio7, bio5, bio3, bio2, bio16, bio15, bio12 and asp see Table 1.

3 结论和讨论

3.1 华石斛的适宜生境分布

华石斛的适宜生境分布在海南10 个市县，其当前的潜在适宜生境主要位于海南热带雨林国家公园霸王岭、鹦哥岭、尖峰岭、吊罗山、五指山、黎母山等分局范围内，五指山市南部、琼中县、万宁市、保亭县和三亚市等非海南热带雨林国家公园区域也有零星分布，预测结果与杨琦[39]的研究结果基本一致。华石斛的潜在适宜分布区范围极其狭窄，并呈现破碎化状态，存在严重的地理隔离。未来2050 年和2070 年在RCP4.5 和RCP8.5 气候情景下，相比当前，华石斛的潜在适宜生境明显减少。造成华石斛野外种群减少及其适宜生境破碎的可能原因有3 个：一是人为因素，华石斛是海南传统黎药，能起到抗菌、抗肿瘤活性的作用[24]，同时具有较高的观赏价值，因此遭到大肆采摘，导致野生种群越来越少，此外，进入工业时代后，矿山开采、经济林的生产活动也严重破坏了华石斛的适宜生境[39]；二是环境因素，华石斛所需生境条件较为苛刻，对高温极为敏感，对其附生树种也较为挑剔[26]，预测其在海南岛的适宜生境面积大致为403.4 km2，仅占海南岛总面积的1.1%，其适宜生境面积狭窄, 各适宜分布区间存在严重的地理隔离，使得华石斛的繁衍与扩散存在一定的困难；三是自身因素，华石斛受传粉者的影响，结实率较低[30]，幼苗在自然环境下萌发率不高，繁殖能力弱，生长缓慢[40]。

华石斛仅分布在海南岛中部和西部海拔1 000 m以上的高山上，是兰科石斛属的典型附生植物, 对周围环境的变化有着较为敏感的反应[41]。若将华石斛引种栽培到低海拔高温地区，持续高温会引起生长严重不适应，植株易出现叶片枯萎、凋落、褐化甚至死亡的现象。此外，华石斛在30 ℃以上的持续高温下幼苗和植株的生理发生改变甚至死亡[26]。本研究结果也表明影响华石斛分布的4 个主导环境因子是海拔、降水量季节性变异系数、最湿季降水量和最暖月最高温。这表明华石斛已经适应了高海拔地区的冷凉湿润生境, 而对高温环境较为敏感,这也可能是华石斛未分布于低海拔地区或其他温度较高地区的主要原因。

3.2 华石斛的保护建议

华石斛目前已发现的种群主要分布于海南热带雨林国家公园内，因此对华石斛的保护主要加强就地保护，同时兼顾其附生宿主树的保护。华石斛具有珍贵的药用价值，因此遭到人为挖采，建议有关部门加强普及华石斛的相关保护知识，提高民众保护意识，并维护种群的生境稳定。此外，华石斛引种栽培时，环境因子的控制可以参考主导环境因子适宜范围，开展华石斛菌根一体化成苗技术体系[42]，制定合理的重引入计划，扩大种群数量。在我们预测的适宜分布区内，都发现了华石斛新增种群分布点，在各个保护区内存在的非常适宜、高度适宜和中度适宜生境区，尤其是猕猴岭东北部及位于猕猴岭南部的毫肉岭、豪雷岭，和鹦哥岭南部的蛙岭、白石岭及石板岭存在大片非常适宜和高度适宜生境，这些区域至今还未有人去深入调查过，可以作为今后进行野外调查以及华石斛再引入工作的重点目标，不仅有利于指导开展野外调研工作,同时也可以为海南热带雨林国家公园各分局管理工作提供一些科学的指导。

3.3 模型因子选择

本研究选取了来自世界气候数据库19 个生物气候变量，同时还将地形变量纳入模型中，丰富了预测模型的模拟变量，增加了预测结果的可靠性。有研究表明，宿主树组成可能会对附生植物组合有显著影响[43–45]，而华石斛对宿主树选择存在广泛的多样性，其宿主树种类多达25 科46 属71 种，无宿主专一性[24]。研究宿主专一性附生植物的潜在适宜生境分布格局不仅需要了解区域气候，还需要考虑其宿主树的分布以及附生植物所处空间内生物的相互作用。本研究选取的环境变量考虑了温度、水分等气候因子及地形因子，而对于微生境、生物竞争、人为因素等其他因素未进行探讨，对于结果的准确性存在一定影响。如何将人为因素、生物因素量化成预测变量纳入到物种分布模型中参与建模可能是未来的研究热点之一。

致谢感谢海南热带雨林国家公园护林员王进强在野外调查过程中的协助，郑希龙与我们共享华石斛的分布数据，海南大学林学院邱梓轩老师、郝馨老师对本论文提出的修改建议。