珍贵用材树种红椿4个变种栽培的潜在气候适生区预测*

张春华，雷晨雨，王 储，冯德枫，孙永玉

(1.中国林业科学研究院 高原林业研究所，云南 昆明 650233；2.云南元谋干热河谷生态系统国家定位观测研究站，云南 昆明 650233；3.西南林业大学 林学院，云南 昆明 650224)

最大熵(MaxEnt)物种分布模型源于统计力学，已发展成为保护生态学、恢复生态学和植物地理学研究中广泛应用的物种分布模型[1-3]。相较于其他生态位模型，MaxEnt 物种分布模型在模拟物种适生区分布过程中具有数据集不完整(仅存分布数据即可)、模型运行时间短、操作简单、样本小且模拟精度高等优点，同时能揭示物种适生区分布与特定环境变量之间的关系[4-7]。应用MaxEnt 物种分布模型的研究方向主要涉及物种应对全球气候变化的适生区变迁[8]、入侵物种潜在扩散区预测及控制[9]、寻找高经济价值物种和特有濒危物种潜在适生区及引种区划[10-11]以及各类自然保护区的核心区、缓冲区和边缘区界定[12-13]等方面。诸多国内外研究证实：MaxEnt 物种分布模型能较好地模拟物种(动植物)的潜在适生区。

红椿(Toona ciliata)在中国有“中国桃花心木”之美称，属楝科(Meliaceae)香椿属(Toona)落叶树种，主干通直，木材纹理优美，是中国南方丘陵山地零星分布的乡土珍贵用材树种之一；包括红椿原变种(T.ciliatavar.ciliata)、毛红椿(T.ciliatavar.pubescens)、思茅红椿(T.ciliatavar.henryi)、滇红椿(T.ciliatavar.yunnanensis)和疏花红椿(T.ciliatavar.sublaxiflora)共5个变种，且在云南都有分布[14]，但因过渡利用、生境破碎及自然更新能力低下而被纳入了国家二级重点保护野生植物名录[15-16]。虽然在20 世纪初期就已对红椿引种、繁育和栽培等进行了探索与研究，但由于未考虑红椿变种对环境需求的异质性，除少量科学试验林[17]外，尚未有规模化栽培的成功案例，主要表现为保存率低、生长缓慢和病虫害严重。红椿为热带、亚热带树种，自然分布于云南海拔1600 m 以下地区[14]，而该区域的土壤类型较为单一(红壤)且立体气候显著[18]，因此，区域内气候差异及水热配比是影响红椿变种人工栽培成败的关键因素。本研究以红椿原变种、毛红椿、思茅红椿和滇红椿为研究对象，利用Max-Ent 模型模拟当前气候下红椿变种在云南的潜在气候适生区分布，并确定影响潜在适生区分布的主要环境因素，为红椿变种的气候适宜性评价、资源保护和人工林种植提供理论支持。

1 材料与方法

1.1 红椿变种地理分布数据

依据红椿变种的主要分类性状[19]，于2018 年9 月—2019 年9 月在云南省范围内通过实地调查得到100个红椿分布位置数据，其中思茅红椿20个、毛红椿16个、红椿原变种45个、滇红椿19个，并采用精度为 ±5 m 的多通道全球定位系统(GPS)记录经纬度坐标，在Excel 中整理成只包含有物种名、经度和纬度的CVS 格式文件备用[17]。从国家基础地理信息系统(http://bzdt.ch.mnr.gov.cn/index.jsp)下载云南省行政区划矢量地图(1∶4 000 000)，应用ArcGIS 10.5 绘制红椿4个变种的潜在气候适生区分布图并估算面积[3]。

1.2 环境因子

环境变量是影响物种分布的关键因素，而19个生物气候变量(表1)已被大量研究证实是影响物种潜在分布的最重要变量[5-7]。此外，海拔也是影响目标物种潜在分布的重要环境因子[20]。从WorldClim (www.worldclim.org)下载得到当代气候条件下云南省行政区范围内的19个生物气候变量和海拔，空间分辨率为30 s (约1 km2)[21]。建立红椿各变种适生区预测模型时移除聚集在10 km 范围内的样本，以减少采样偏差的影响[22]。采用PHILLIPS 等[4]的方法删除相关系数大于0.8的环境因子，以避免环境变量间的多重共线性。

1.3 研究方法

MaxEnt 模型以概率论为基础，具有稳定的算法及简明的数学定义。它以目标物种现存分布数据和环境因子数据构建相应的物种潜在适生区分布模型，最终以概率(0～1)的形式表示目标物种在某个生境条件下存在的可能性，且不受样本量大小的影响[4,9-11]。本研究采用MaxEnt (v3.4.1)研究上述4个红椿变种在云南的潜在气候适生区。建模时将红椿各变种对应的物种分布数据和环境因子加载到MaxEnt 中，随机选取75%和25%的物种分布数据分别作为训练集和测试集[23]，最大迭代次数固定为3 000，其他值保持默认水平[24]。通过计算受试者工作特征(receiver operating characteristic，ROC)曲线下的面积(area under ROC curve，AUC)度量模型的准确性，AUC 值低于0.7 为较差，0.7～0.9 为中等，大于0.9 为好[25]。采用刀切法计算环境变量对预测模型的贡献率[20]。最后根据各变种MaxEnt 生态位模拟模型输出的概率(P值)划分高适生区(P≥0.6)、中适生区(0.4≤P＜0.6)、低适生区(0.2≤P＜0.4)和非适生区(P＜0.2)[26]。

2 结果与分析

2.1 红椿变种潜在气候适生区模型预测精度检验

思茅红椿潜在气候适生区分布模型的平均训练AUC=0.999，平均测试AUC=0.999；毛红椿潜在气候适生区分布模型的平均训练AUC=0.992，平均测试AUC=0.983；红椿原变种潜在气候适生区分布模型的平均训练AUC=0.991，平均测试AUC=0.991；滇红椿潜在气候适生区分布模型的平均训练AUC=0.996，平均测试AUC=0.984。据模型预测精度判定标准，红椿变种潜在气候适生区分布模型的预测精度达“好”级别，预测结果的精度和可信度均较高。

2.2 影响红椿变种潜在气候适生区预测模型的主要环境因子判别

经MaxEnt 模型分析(表1)可知：限制红椿原变种潜在气候适生区分布的因子为Bio1、Bio3、Bio4、Bio7 和Alti；限制毛红椿潜在气候适生区分布的因子为Bio1、Bio5、Bio6、Bio12 和Alti；而Bio4、Bio6、Bio12、Bio13、Bio14和Alti 是思茅红椿潜在气候适生区的限制因子；Bio1、Bio2、Bio3、Bio12、Bio18 和Alti 则是滇红椿潜在气候适生区的限制因子。

由表1 还可知：对红椿原变种潜在气候适生区分布模型贡献率最高的前4 位环境变量依次为等温性、温度季节变化标准差、年均温和年均温变化范围，累积贡献率为96.2%；对毛红椿潜在气候适生区分布模型贡献率最高的前4 位环境变量依次为最冷月最低温、年降水量、年均温和海拔，累积贡献率为96.0%；对思茅红椿潜在气候适生区分布模型贡献率最高的前4 位环境变量依次为温度季节变化标准差、最干月降水量、海拔和最湿月降水量，累积贡献率为96.0%；对滇红椿潜在气候适生区分布模型贡献率最高的前4 位环境变量依次为年均温、等温性、海拔和最暖季降水量，累积贡献率为90.5%。

表1 环境变量及其对分布模型的贡献率Tab.1 Environmental variables and their contribution rates in predicting the current distribution

环境因子响应曲线直观体现了环境变量与物种潜在分布的定量关系。根据响应曲线，等温性大于47.5、温度季节变化标准差30～42、年均温17～22.5 ℃和年均温变化范围22～24 ℃是红椿原变种潜在气候适生区分布的最佳环境条件组合(图1a～d)；最冷月最低温-1～8.5 ℃、年降水量600～1 300 mm、年均温大于16 ℃和海拔600～1 700 m 则为毛红椿潜在气候适生区分布的最佳环境条件组合(图1e～h)；限制思茅红椿潜在气候适生区分布的环境变量组合为温度季节变化标准差30～42、最干月降水量11～13 mm、海拔700～1 400 m 和最湿月降水量200～280 mm (图1i～l)；限制滇红椿潜在气候适生区分布的环境变量组合为年均温大于18 ℃、等温性大于48、海拔300～1 100 m 和最暖季降水量170～520 mm (图1m～p)。

图1 红椿各变种前4个最大贡献环境因子响应曲线Fig.1 The response curves for the top four environmental variables of Toona ciliata varietas with the largest contribution

2.3 红椿各变种的潜在气候适生区分布及面积

由图2 可知：红椿原变种的潜在气候适生区除迪庆、昭通、曲靖和文山有少量分布外，其他地州均有大面积分布，其潜在气候低、中、高适生区面积分别约74 353.11、76 696.77 和71 753.55 km2；潜在气候中、高适生区主要分布于怒江、大理、昆明、玉溪和红河一线以南，另在金沙江河谷部分地区存在一定面积高适生区。毛红椿的潜在气候适生区主要分布于滇中高原南缘以北地区，潜在气候低、中、高适生区面积分别约74 326.58、73 370.11 和55 355.32 km2。思茅红椿的潜在气候适生区主要分布于普洱、临沧、大理、楚雄和玉溪的交界处，还包括西双版纳的部分区域，其潜在气候低、中、高适生区面积分别约36 290.01、28 648.56 和9 092.09 km2。滇红椿的潜在气候适生区主要分布于滇西南的德宏、普洱、临沧、西双版纳和红河，其潜在气候低、中、高适生区面积分别约37 640.56、36 243.74 和34 356.04 km2。

图2 红椿变种潜在气候适生区Fig.2 The potentially suitable climate areas for T.ciliata varietas

3 讨论

物种适生区与难以量化的种间关系、居民的生产活动和居群动态密不可分[13,17]，且前人的研究重点是气候变化对目标物种潜在适生区的影响[22-26]，本研究则主要探讨在当前气候条件下红椿变种潜在气候适生区在云南省的分布情况，气候变化对其适生区的影响及复杂的种间关系、居民生产生活和物种居群动态等因子不在模型考虑的参数范围之列。本研究表明：红椿原变种、毛红椿、思茅红椿和滇红椿潜在气候适生区分布预测模型的工作特征曲线下面积均大于0.9，达到Maxent 预测模型精度判别标准中的最优级别[10]，说明对红椿各变种的潜在气候适生区模拟准确、预测结果可靠。

环境因子对物种分布模型的贡献率与前人对众多物种的潜在适生区分布研究结果[27-29]一致，即影响不同物种潜在适生区分布的环境因子各有不同，即使相同环境因子在不同物种预测模型中贡献率也各有不同。红椿为亚热带树种[14]，根据环境因子对红椿变种适生区分布的贡献率可知：红椿原变种对温度变化要求较高，对降水则无特定需求，为温度控制型树种[30]，分布于年均温大于17 ℃的云南大部分区域，即使年平均降雨量最少的金沙江干热河谷[31]也出现了一定面积的潜在气候高适生区；而毛红椿、思茅红椿和滇红椿的潜在气候适生区分布则受温度和降水量共同影响[32-33]。其中，毛红椿的潜在气候适生区则要求年均温大于16 ℃，且有一定的低温条件(-1～8.5 ℃区域)，另外年降水量600～1 300 mm限制了其主要分布于滇中高原以北区域；而温度季节变化标准差、最干月降水量、海拔和最湿月降水量则把思茅红椿的潜在气候适生区限制在以受西南和东南季风交汇地[31]的普洱地区北部为中心的小部分区域；滇红椿潜在气候适生区要求年均温大于18 ℃，且对等温性的要求高于红椿原变种，加之最暖季降水量170～520 mm 的限制，形成了其高适生区与红椿原变种南缘除西双版纳州外高度重合。

国家标本平台(http://www.nsii.org.cn)、数字植物标本馆(http://www.cvh.org.cn)和全球生物多样性数据库(GBIF，https://www.gbif.org)是物种现存分布记录的重要补充，但目前不能从中获取物种分布点的地理坐标信息，因此不能作为有效分布用于物种潜在适生区研究[33]。红椿变种潜在气候适生区与其自然分布[14]高度相似，说明早期获得的红椿变种分布数据可靠。从现有文献记载看，思茅红椿和滇红椿仅分布于云南省行政区域内，而思茅红椿集中分布于普洱地区，导致对其潜在适生区的预测结果看似出现了“牛眼”现象，但与实际分布[14]是高度吻合的。

MaxEnt 模型得到的物种潜在气候高适生区是物种原地保护和资源培育的最佳区域[27,34-35]。从红椿变种的高适生区面积来看，红椿原变种面积最大(71 753.55 km2)，毛红椿(55 355.32 km2)和滇红椿(34 356.04 km2)次之，思茅红椿(仅9 092.09 km2)最小。由此得到红椿变种保护优先顺序应为思茅红椿、滇红椿、毛红椿和红椿原变种。变种对环境的要求不一致，在进行物种保护和种群重建过程中需在遵循“适地适树”原则[36]下做到适地适变种分区域推进，才能实现红椿各变种保护和资源培育双重效益最大化的目标。滇中高原西部至西南部地区，包括德宏、临沧、普洱和西双版纳的大部以及大理、丽江和楚雄的金沙江沿岸区域是红椿原变种原地保护和发展大径材人工用材林的理想区域；滇中高原的北部地区则适合于毛红椿的保护和人工林的发展；普洱北部和临沧东部区域是思茅红椿的重点保护区域，德宏—临沧—普洱一线以南及以西区域则是滇红椿的重点保护区域。

4 结论

MaxEnt 模型能较好地模拟红椿原变种、毛红椿、思茅红椿和滇红椿在云南省的潜在气候适生区范围。等温性、温度季节变化标准差、年均温和年均温变化范围等因子影响红椿原变种潜在气候分布区，潜在气候适生区总面积约22.58×104km2；最冷月最低温、年降水量、年均温和海拔等因子影响毛红椿潜在气候分布区，潜在气候适生区总面积约20.31×104km2；温度季节变化标准差、最湿月降水量、海拔和最干月降水量等因子影响思茅红椿潜在气候分布区，潜在气候适生区总面积约7.40×104km2；年均温、等温性、海拔和最暖季降水量等因子影响滇红椿潜在气候分布区，潜在气候适生区总面积约10.82×104km2。