基于MaxEnt 模型预测鹅掌楸在中国的潜在分布区

邱浩杰， 孙杰杰， 徐 达， 沈爱华， 江 波， 袁位高， 李 胜

（1. 浙江农林大学 风景园林与建筑学院， 浙江 杭州311300； 2. 南京林业大学 生物与环境学院， 江苏南京210037； 3. 浙江省林业科学研究院， 浙江 杭州310023）

北美鹅掌楸Liriodendron tulipifera和鹅掌楸Liriodendron chinense是鹅掌楸属Liriodendron中残存的2 个种。 受人类活动影响和天然更新不良等因素制约， 鹅掌楸的种群数量和分布范围逐年减少， 基本以1 株或数株的零散分布存在。 在中国， 鹅掌楸被列为国家二级珍稀濒危植物［1］， 主要分为东、 西2 个亚区： 西部分区由北部的大巴山经鄂西、 湘西、 川东和黔北中山（武陵山和大娄山）向西南延伸到云贵高原东部中山（苗岭）， 南端沿云贵高原南坡向西南分布至北回归线以南的金平和麻栗坡， 直至越南北部； 东部分区包括地处鄂东、 皖北交界地带的“大别山岛”， 赣西北、 鄂东、 湘东三省交界地带的 “九幕岛”，浙南、 闽北交界地带的 “武夷山岛” 和浙西北、 皖南交界地带的 “天黄岛” 等4 个岛状分布区［2-4］。 目前， 鹅掌楸的研究主要集中在生理生化、 育苗和观赏价值等方面。 如陆畅等［5］对不同种源鹅掌楸的抗寒性能进行了研究； 刘丹等［6］运用随机扩增多态性DNA（RAPD）分子标记手段揭示了鹅掌楸的遗传多样性和变异规律； 张晓平等［7］研究了淹水胁迫对浙江种源鹅掌楸光合特征的影响， 间接揭示光照与降水对鹅掌楸生长的影响。 对鹅掌楸地理分布上的研究主要以调查和整理为主， 对其分布区域和影响因子的定量分析未见报道。 方炎明［2］和贺善安等［3］研究认为： 影响鹅掌楸生长的主要因素是生理障碍、 气候环境和人类活动。 郝日明等［4］进一步研究发现： 由于不合适的生境引起的鹅掌楸天然更新不良是造成鹅掌楸生理障碍的主要原因。 环境因素成为当前全球气候变暖背景下影响生物分布的决定性因素［8］， 预测受威胁物种如鹅掌楸适宜分布区的空间变化和迁移趋势已成为最新研究热点［9-10］。 MaxEnt 最大熵模型（maximum entropy model, MaxEnt）是以最大熵理论为基础的密度估计和物种分布预测模型， 具有运算结果稳定， 运算时间短等优点［11-13］， 被广泛应用于动植物生长环境分析、 害虫预警、 栖息地保护和潜在分布区预测等。 关心怡等［14］借助该模型预测了未来气候情景下的麻栎Quercus acutissima的地理分布； 李绒等［15］基于MaxEnt 模型分析了影响绿孔雀Pavo muticus生长的主要环境因子； 王茹琳等［16］借助此模型分析了影响西藏飞蝗Locusta migratoria tibetensis潜在分布的环境变量和适宜分布区， 在害虫预警方面得到较好反响。 另外， 该模型在交通［17］、 用地条件［18］等新领域也能发挥较好效果。 本研究以鹅掌楸在中国的标本分布记录以及相应分布地区的生物气候环境、 人类活动强度等数据为基础， 模拟鹅掌楸在当代气候环境下的潜在分布区， 分析影响其分布的主导环境因子和适宜的环境范围； 选取21 世纪50 年代和21 世纪70 年代不同气候情景， 分析未来气候环境对鹅掌楸空间分布的影响， 为今后鹅掌楸的栽培保护提供参考。

1 数据来源和研究方法

1.1 标本数据来源

以中国国家标本资源平台（NSII， http：//www.nsii.org.cn/）、 中国数字植物标本馆（CVH， http：//www.cvh.ac.cn/）和教学标本资源共享平台（http：//mnh.scu.edu.cn）等标本信息、 平台信息和文献检索信息为依据， 剔除位置信息未精确到镇的标本数据和重复数据。 根据气候精度， 在每个2.5′×2.5′网格处选取1 个分布点， 共搜集到鹅掌楸标本信息154 条， 按照MaxEnt 软件要求录入并保存信息。

1.2 环境数据来源

从世界气候数据库（http：//www.worldclim.org/）下载1970-2000 年坐标系为WGS84、 栅格大小为25 km2的19 种气候环境数据（表1）， 数据空间分辨率为2.5 min。 其中气候数据采用中国气象局北京气候中心的BCC-CSM1-1 大气环流模型； 未来温室气体排放量（GHG）采用联合国政府间气候变化专门委员会（IPCC）第5 次报告中的温室气体排放场景（representative concentration pathways， RCP）， 选择RCP 2.6 和RCP 8.5 等2 种气候变化情景， 分别代表温室气体增加对未来气候的最低和最高影响［14］； 中国行政区划图来源于国家基础地理信息中心网站（http：//www.ngcc.cn/）。 ArcGIS 软件版本为10.2， MaxEnt 软件版本为3.3.3k。

人类活动对鹅掌楸影响较大。 本研究将人类活动强度数据（hf-v2geo）作为干扰变量进行研究。 该数据来自于国际地球科学信息网络中心（Center for Earth Science Information Network, CIESIN）， 能全面和客观体现人类活动的程度及空间分布［19］。 将下载数据进行格式与坐标转换， 导入到GIS 软件中裁剪得到中国地区的人类活动强度。

表1 模拟鹅掌楸适生区的环境变量及其贡献率Table 1 Environmental factors used for simulating suitable area of L. chinense and percent contribution

1.3 模型的设置与选择

将19 个变量和人类活动强度数据导入MaxEnt 软件。 设分布数据的25%作为测试集（test data）， 剩余的75%作为训练集（trainning data）。 为保证结果的准确性， 设置软件重复运算10 次， 输出分布值为逻辑斯蒂值（logistic）， 表示分布概率（简称为分布值）， 勾选Do Jackknife to measure variable importance 和Creat response Curves 2 项， 其余选项保持软件默认设置。 以省级行政区划图为底图， 以MaxEnt 10 次运算的平均值为依据， 绘制鹅掌楸在当代、 21 世纪50 年代和21 世纪70 年代等3 个时期不同气候变化情景下的适宜性分布图。

2 结果与分析

2.1 模型准确性检验

接受者操作特性曲线（receiver operating characteristic curve， ROC）分析法在物种潜在分布预测模型评价中得到了广泛的应用， 是目前认可度较高的诊断试验评价指标［20］。 ROC 曲线与横坐标围成的区域面积（the area under the ROC curve， AUC）， 取值范围为0～1， 值越接近于1， 表明模型模拟的精确度越高。 10 组训练集AUC 均值为0.973， 测试集的均值为0.953， 表示软件的模拟效果达到了极高的水平［21］，实验数据可信度高。 本研究10 次运算结果中数值最大的1 组训练集AUC 值为0.971， 测试集值为0.968， 可作为分析对象（图1）。

2.2 环境因子的贡献率

MaxEnt 模型通过刀切法（Jackknife）得出各环境因子对模拟鹅掌楸分布的贡献率（contribution rate），由表1 可知： Bio12（年均降水量）、 Bio16（最湿季度降水量）、 Bio6（最冷月最低温）、 Bio15（降水量变异系数）、 和Bio2（昼夜温差月均值）是对鹅掌楸的潜在分布影响最大的5 个环境变量， 累计贡献率达80.4%； 与此同时， 人类活动强度对鹅掌楸分布的贡献率达2.3%， 可见人类活动（如人口压力、 土地利用、 基础设施建设、 交通运输等）对鹅掌楸繁衍和生长也存在着一定程度的影响［19］。

图1 MaxEnt 模拟的AUC 曲线分析Figure 1 ROC curve analysis of MaxEnt simulation

2.3 鹅掌楸分布地区的气候特征

由表1 模拟结果可知： 降水因子对鹅掌楸分布影响最大。 选取年均降水量、 最湿季度降水量、 最冷月最低温和降水量变异系数4 个影响最大的变量进行具体分析， 通过MaxEnt 自动输出各变量因子的响应曲线（图2）， 分析最适合鹅掌楸生长的生物气候变化范围。 一般认为， 当存在概率大于0.50 时， 其对应的环境变量数值更有利于鹅掌楸生长［22］。 由图2A 可知： 年均降水量小于1 000 mm 时， 鹅掌楸的存在概率很低（＜0.05）； 随年均降水量增大， 存在概率大幅增长。 当存在概率大于0.50 时可认为该变量处于适宜范围。 鹅掌楸适宜生长的年均降水量范围为1 200～1 800 mm， 1 350 mm 为其最适值。 图2B 显示： 最湿季度降水量小于450 mm 时， 鹅掌楸存在概率极低； 随着降水量的增大， 存在概率急速上升，并在500～900 mm 时达到并大于0.50， 表明500～900 mm 是鹅掌楸适宜的最湿季度降水量， 600～850 mm是鹅掌楸生长最适宜的最湿季度降水量。 最冷月最低温≤-4.0 ℃时， 鹅掌楸存在概率较低； 随着温度的升高， 鹅掌楸存在概率迅速增大， 在-2.0～-3.0 ℃时存在概率超过0.50； 最适温度点为-1.0 ℃。 之后又快速回落， 最冷月最低温超过5.0 ℃后存在概率接近于0（图2C）。 降水量变异系数＜45 时， 鹅掌楸存在概率很低， 适宜范围为45～65； 当降水量变异系数为55 时， 鹅掌楸达最适生长， 之后不断下降（图2D）。

图2 鹅掌楸存在概率对主要气候环境因子的响应Figure 2 Response curve of the existence probability of L. chinense to major climatic environmental factors

2.4 当代环境条件下鹅掌楸在中国的适生区分布

利用ArcGIS 将MaxEnt 模拟结果与中国省级行政区划图叠加， 对鹅掌楸潜在分布区进行适生等级评价。 采用人工分级方法对鹅掌楸潜在分布区进行适宜性等级划分， 适生指数0～0.1 为不适宜分布区，0.1～0.3 为低适宜分布区， 0.3～0.5 为中适宜分布区， 0.5～1.0 为高适宜分布区。 统计每个等级内的栅格数并计算每个等级的面积（表2）。 当代生物气候环境下， 鹅掌楸分布范围为22°～32°N， 103°～123°E， 覆盖了中国19 个（省、 市自治区）， 总适宜面积约129.8 万km2。 湖南、 贵州和江西的适宜区最大， 分别达20.0 万km2， 15.9 万km2和15.6 万km2。 按适宜等级划分， 高适宜分布区面积约20.7 万km2， 主要为中国西南和华东地区。 相比而言， 西南地区面积较大， 分布也较集中； 该区域北起大巴山南的重庆巫溪县和湖北恩施州， 经方斗山、 大娄山和武陵山， 包括重庆石柱县、 武隆县、 彭水县、 酉阳县， 湖南湘西州等大部分地区， 沿西南延伸到遵义等贵州中北部地区， 往南沿苗岭、 凤凰山呈带状， 分布较为稀疏。 华东地区北起天目山， 沿仙霞岭、 武夷山脉延伸至浙南和闽北交界的丘陵山区， 自西向东包括江西上饶， 安徽黄山， 浙江丽水、 台州和宁波等地。 中适宜分布区面积约51.9 万km2， 从高适宜分布区向北延伸到长江以北、 秦岭以南， 陕西安康， 四川达州和巴中都有分布， 向西至华蓥山一带， 且在大凉山附近的四川内江、 乐江、 宜宾和资阳呈现小范围集中分布。 低适宜分布区的范围则进一步扩大， 面积约57.2 万km2， 覆盖了贵州、 重庆、 湖南、 江西及浙江全境， 且在龙门山以东的四川中东部、 武夷山至戴云山脉间的福建中北部、 大别山以南的湖北南部、 南岭山脉的广西北部和广东北部等地出现大面积适宜区。 另外， 在花莲、 新竹等中国台湾中北部地区， 苏州、 无锡和南通等江苏南部地区和上海部分地区也有零星分布。

表2 当代鹅掌楸在中国各省（市）的潜在分布区面积Table 2 Potential distribution area of L. chinense in various provinces of China in current

2.5 鹅掌楸未来的适生分布区

选取未来21 世纪50 年代和21 世纪70 年代的不同气候环境情景， 利用MaxEnt 模型对鹅掌楸的潜在地理分布进行模拟， 得到未来不同气候情景下鹅掌楸适宜分布区分布面积（表3）。 结果表明： 随着年代的推移， 鹅掌楸在中国的适宜分布面积先稳定后减少， 并呈现向高纬度地区和东部沿海地区轻微扩展趋势。 当温室气体排放量达到典型目标值（RCP 8.5）时， 该趋势更加明显。 相比当代， RCP 8.5 情景下，21 世纪50 年代适宜分布区将减少0.5 万km2， 21 世纪70 年代将减少6.9 万km2， 下降幅度达到了5.3%。RCP2.6 情景下， 面积变化趋势同样为先稳定后下降， 但变化幅度较小。 具体来看， 当代、 21 世纪50 年代和21 世纪70 年代的总适宜分布区面积分别为129.8、 131.3 和126.9 万km2； 低适宜分布区面积和高适宜分布区面积变化趋势相似， 都是先增加后减少； 中适宜分布区面积在当代最高， 之后不断下降并保持稳定。 从分布地区来看， 与当代分布区相比， 未来低适宜分布区较当代有比较明显的北抬趋势， 江苏中北部和陕西南部都会出现一定范围的分布， 中国台湾中北部的低适宜分布区也出现了明显的增加； 中高适宜分布区开始呈现东西分布区连接中断的趋势， 且分布区域出现破碎化， 并出现向东部迁移的倾向。 四川资阳、 内江和宜宾的大凉山附近分布区域逐渐萎缩； 丰城、 新余、 抚州等江西中北部， 岳阳、常德等湖南北部以及鄂州、 黄石等湖北东部的中高适宜分布区消失， 自西向东的中高适宜分布区在此出现不同程度的断裂。 集中分布区仅存在于大巴山以南的湖北西南部和重庆东南部以及浙皖交界处的天目山一带。

3 讨论

3.1 主导气候因子对鹅掌楸的分布影响

本研究的20 个环境因子中， 降水相关因子对鹅掌楸分布影响最大， 累计贡献率达到71.4%，其中， 年均降水量贡献率最大。 MaxEnt 模型预测1 200～1 800 mm 为鹅掌楸生长的最适年降水量范围， 该结果与王章荣［23］提出的鹅掌楸分布于雨量充沛、 湿度大、 季节分明的亚热带、 暖温带气候区， 年均降水量为800～2 300 mm 的结论相符。 事实上， 从中国气象局所查询到的各个站点的气候数据显示， 已存在的154 个标本分布点年均降水量都在这一范围， 特别是高适宜分布区（贵州、 四川等）的年均降水量都接近最适降水量（1 350 mm）。 最冷月最低温对鹅掌楸的分布也有一定的影响， 低于-4.0 ℃和高于5.0 ℃以上都不适宜鹅掌楸生长。 正因如此， 在当代， 气温和降水两大影响因子将鹅掌楸的适生范围限定在了秦岭—淮河以南； 而福建南部、 广东等低纬度地区虽然冬季温度较高， 但鹅掌楸也较少分布。 此外， 季节性降雨（Bio15）， 极端气温（Bio6、 Bio10、 Bio11）和降水（Bio16）也对其生长存在影响。

表3 鹅掌楸不同气候情景下的潜在分布区的面积Table 3 Areas of potential distribution under different climate scenarios of L. chinense

3.2 气候变化下鹅掌楸地理空间分布的变化

目前， 学者大都认为未来物种的适宜分布面积将不断减少， 且有向高纬度和高海拔地区迁移的趋势［24］。 如MCKENNEY 等［25］发现： 在未来北美130 种树木分布范围都将出现一定程度的缩减。 吴建国［26］对国内7 种乔木的研究发现， 气候变化下这些乔木未来分布面积将不断萎缩， 且有向高海拔区域扩散的趋势。 近年来， 借助于MaxEnt 等模型分析认为， 未来3 种气候情景下黄山花楸Sorbus amabilis适生区面积较当代缩减［27］； 白栎Quercus fabri可能向北扩张， 且面积也有一定程度的缩小［28］。 但也有学者得出不同的结论。 关心怡等［14］在研究麻栎的适宜性分布时认为： 气候变暖对部分植物生长有促进作用， 在未来气候条件下麻栎适生区域面积总体上是增加的。 本研究预测了21 世纪50 年代和21 世纪70 年代不同RCP 环境情景下鹅掌楸的潜在分布区， 结果表明： 鹅掌楸适宜分布区面积呈现先稳定后不断减少的趋势。 作为喜湿植物， 鹅掌楸对年均降水量的适宜范围较大， 有一定的耐寒性， 可以推测这是未来鹅掌楸适宜面积保持稳定的一个原因。 但21 世纪50 年代鹅掌楸低适宜分布范围有明显的北抬趋势， 在陕西南部、 江苏中部等地出现较大面积的分布区。 这是温度和降水量双重影响的结果。 一方面降水作为最大的影响因子， 在全球变暖的背景下， 中国区域内的降水量会不断增加， 尤以北方地区最为明显［29］； 另一方面不断增大的降水势必会破坏低海拔生境［27］。 到21 世纪70 年代， 气温和降水增大趋势更加明显， 鹅掌楸原生长区域的适宜度进一步降低。 中适宜分布区和高适宜分布区呈现岛状分布， 由当前西南、 华东2个高适生区缩减为以天目山到黄山为中心的东部片区和以湖北、 湖南和贵州交界处为中心的西部片区；剩余中高适宜分布区较为零散， 分布面积也因此缩小。 数据显示， 在RCP 8.5 的情景下， 总适宜面积相比当代将下降5.3%。

3.3 鹅掌楸的栽培与保护

就预测结果而言， 鹅掌楸未来适宜分布区不容乐观。 当前及今后， 鹅掌楸的保护工作需更加深入。一是扩大种植范围， 在进一步认识鹅掌楸适生气候环境前提下， 在未知的适生区域和未来新增区域上划定保护区， 减少适生区尤其是新增适生区域的人为干扰； 加强已有适生区域， 特别是高适宜分布区的保护， 防止乱砍乱伐， 使之成为鹅掌楸应对气候变化的安全地和避难所［30］。 另一方面， 鹅掌楸是优美的园林树种， 可在相关区域作为园林绿化树种进行推广， 通过多种途径提高种群数量。 本研究建模时未将土壤因子， 地形数据（包括坡度和坡向）等作为影响因子， 栽培时建议在本研究基础上， 选择排水良好的酸性土（pH 4.5～6.5）且具有适宜光照的区域种植［31］， 为其提供优越的生长环境。