基于Maxent生态位模型的互花米草在我国沿海的潜在分布

以分布在我国沿海滩涂的96个互花米草分布记录点及覆盖东部沿海区域的海洋环境数据和气候环境数据为材料,利用Maxent生态位模型,研究外来物种互花米草在我国沿海的潜在分布情况.结果表明:互花米草适宜分布区占我国沿海区域的85%,其中高度适宜分布区占18%,中度适宜分布区占34%,低度适宜分布区占33%,不适宜分布区仅为15%.互花米草在我国沿海的地理分布主要受到年均最低海水温度、年均海水温度、年平均气温和1月最低气温4个环境因素的影响,而年均降水量、年均日较差、海水盐度、最高海水温度、7月最高气温及海流速度对互花米草地理分布的影响较小.互花米草高度适宜分布区的最低海水温度为0.62~24.81℃,平均海水温度为10.46~27.29℃,年均气温为9~25℃,1月最低气温为-13.5~16.7℃.互花米草地理分布概率在我国北部沿海区域达到20%以上,互花米草存在向我国北部进一步入侵的趋势,特别是在渤海湾地区,互花米草入侵潜力较大.互花米草在我国不适宜分布区主要集中在海南中部和南部海岸,以及台湾省大部分区域,依据当前我国分布记录及气候数据,这些区域互花米草入侵风险较小,但不排除未来入侵的可能性.     

Maxent模型复杂度对物种潜在分布区预测的影响

生态位模型在入侵生物学和保护生物学中具有广泛的应用, 其中Maxent模型最为流行, 被越来越多地应用在预测物种的现实分布和潜在分布的研究中。在Maxent模型中, 多数研究者采用默认参数来构建模型, 这些默认参数源自早期对266个物种的测试, 以预测物种的现实分布为目的。近期研究发现, Maxent模型采用复杂机械学习算法, 对采样偏差敏感, 易产生过度拟合, 模型转移能力仅在低阈值情况下较好。基于默认参数的Maxent模型不仅预测结果不可靠, 而且有时很难解释。在本研究中, 作者以入侵害虫茶翅蝽(Halyomorpha halys)为例, 采用经典模型构建方案(即构建本土模型然后将其转移至入侵地来评估), 利用ENMeval数据包来调整本土Maxent模型调控倍频和特征组合参数, 分析各种参数条件下模型的复杂度, 然后选取最低复杂度的模型参数(即为最优模型), 综合比较默认参数和调整参数后Maxent模型的响应曲线和预测结果, 探讨Maxent模型复杂度对预测结果的影响及Maxent模型构建时所需注意事项, 以期对物种潜在分布进行合理的预测, 促进Maxent模型在我国的合理运用和发展。作者认为, 环境变量的选择至关重要, 需要综合分析其对所模拟物种分布的限制作用和环境变量之间的空间相关性。构建Maxent模型前需对物种分布采样偏差及模型的构建区域进行合理地判断, 模型构建时需要比较不同参数下模型的预测结果和响应曲线, 选取复杂度较低的模型参数来最终建模。在茶翅蝽的分析中, Maxent模型的默认参数和最优模型参数不同, 与Maxent模型默认参数相比, 采用调整参数后所构建的模型预测效果较好, 响应曲线较为平滑, 模型转移能力较高, 能够较为合理反映物种对环境因子的响应和准确地模拟该物种的潜在分布。     

基于Maxent模型的假臭草潜在分布区预测

为有效开展假臭草的防除工作,运用最大熵值模型(Maxent)和地理信息技术(GIS)对假臭草(Praxelis clematidea)在中国的潜在分布区和影响其分布的主导环境变量进行预测。运用AUC法对预测精度进行评价,达到"优秀"水平。结果表明,Maxent模型预测结果具有较高的可信度。假臭草在中国的潜在适生面积为785985km2,占总国土面积的8.19%,主要分布在华南及东南区域,其中福建、广东、海南和广西省为假臭草的中高适生区。福建省的假臭草分布最为广泛,61.98%的福建区域适合假臭草的生长,其中闽南沿海为假臭草的高适生区,而闽北区域假臭草的种群密度较低。2月最低温可能是限制假臭草分布的主要气候变量,贡献率为61.70%;最干季度降水量、9月最高温对假臭草的分布具有一定的影响。假臭草的分布预测为早期监测和预警提供理论依据和技术支持,防止假臭草的进一步扩散蔓延。     

基于Maxent和GARP模型的日本双棘长蠹在中国的潜在地理分布分析

【目的】日本双棘长蠹Sinoxylon japonicum是一种重要的林木蛀干害虫。该虫自1981年首次报道于云南昆明以来,先后在中国的10多个省市被发现。近年来该害虫在北京、天津和河北等地对国槐Sophora japonica的危害日趋严重。【方法】在本研究中,作者按时间顺序梳理了日本双棘长蠹在中国的分布记录,根据日本双棘长蠹已有分布记录及其主要寄主植物国槐在我国的种植区域采用了Maxent和GARP两种生态位模型对日本双棘长蠹的潜在地理分布进行分析。【结果】日本双棘长蠹在中国南方地区呈零星分布,而在中国北方地区较为集中。分布记录的报道时间呈现出由南向北和自东向西的局势,推测这种分布格局是由于其寄主植物所导致的：中国北方地区危害严重可能与国槐在中国北方的集中种植有关,寄主树苗在不同地区间的调运是其种群扩散主要原因。基于生态位模拟的结果显示日本双棘长蠹在中国具有较大的适宜生态空间,潜在地理分布范围较广,从北到南其适生区主要有：辽宁西部、北京、天津、宁夏、河北、山西南部、山东、陕西、江苏、安徽、湖北、重庆、浙江、江西、湖南以及四川和贵州西北部。【结论】这些地区间的树苗转运需要做好检验检疫工作,以防止日本双棘长蠹的扩散。     

基于GIS与Maxent模型的金花茶潜在适生区与保护研究

为明确珍稀植物金花茶(Camellia nitidissima)的潜在分布,利用最大熵模型(Maxent)模拟金花茶在中国的潜在适生区,运用受试者工作特征(ROC)曲线分析法进行验证。结果表明,金花茶的适生区域主要集中在我国西南片区,与已知原生地重合但不局限于已有报道的地点。ROC曲线分析法的面积AUC值为0.992,表明预测结果可靠。Jackknife检验表明,影响金花茶分布的最主要环境因子(阈值)是:年均温(24.3℃)、积温(7 515.84℃)、最干季湿度(14.5%)和最湿季均温(26.2℃)。缓冲区分析定量评价了人类活动对金花茶适生生境的干扰。这些有助于制定金花茶的保护规划。     

基于Maxent模型的檵木分布格局模拟

应用最大熵模型(Maxent)模拟檵木(Loropetalum chinense)在末次间冰期(last inter glacial, LIG)、末次盛冰期(last glacial maximum, LGM)和当代(Present)等不同时期的中国潜在分布格局,分析影响其分布的主导生物气候因子。结果显示:(1)历史气候的变迁,檵木由末次间冰期经末次盛冰期进入当代,适生区面积呈现增加趋势。当代适生区面积占比最大,适生程度也较高;(2)训练数据和测试数据的AUC值分别为0.947和0.954,均达到了极高的精度;(3)刀切法检测表明,影响檵木分布的主导环境因子依次为bio14(最干月份平均降雨量)、bio17(极干季降雨量)、为bio19(极冷季平均降雨量)和bio15(湿度变化方差),其适生值范围分别是23—93 mm, 98—300 mm, 110—350 mm和42%—65%;(4)较纬度而言,经度是影响我国檵木分布格局的主要因子;(5)我国檵木当代的潜在地理分布主要在北回归线以北的区域,总面积为162.55万km²,占国土总面积的16.93%,高适生区集中于...     

基于Maxent生态位模型预测球果角胫象在云南的潜在地理分布

掌握球果角胫象Shirahoshizo conifera Chao在云南潜在分布区,对了解其危害和防控具有重要的意义。本文利用球果角胫象的分布点数据和环境因子数据,通过Maxent生态位模型预测其在云南的潜在分布区。结果表明：球果角胫象的潜在分布区主要集中滇中地区;在ArcGIS中进行显示与风险等级划分,按栅格数值的大小将球果角胫象的适生范围分为4级,分级标准为：高适生区、中适生区、低适生区、非适生区,经ROC曲线分析法验证,Maxent生态位模型的AUC值为0998,表明预测获得了较好的效果。     

基于生态位模型的番荔枝实蝇潜在适生性分布预测

番荔枝实蝇Ceratitis artortoe(Graham)是一种重要的外来人侵性检疫害虫.在广东口岸,其幼虫连续从入境旅客所携带的水果中被榆出.目前关于番荔枝实蝇潜在适生性分布的研究进行得很少,但对于我国的生物生态安全却有重要意义.本研究中,我们使用3种牛态位模型(ENFA模型,马氏典型性模型和Maxent模型)对番荔枝实蝇在中斟以及全球范尉内的潜在适生性分布区域进行了预测分析.结果显示:Maxent模型拥有最好的预测精确度,马氏典型件模型次之,而ENFA模型的预测精确度最差;Maxent模型和马氏典型性模型的预测精确度无显著性差异;根据Maxent模型的预测结果,番荔枝实蝇在中国的潜在适生区主要是广西、广东、海南以及云南的少部分地区.分析结果显示,番荔枝实蝇从境外传人中国南部地区并最终在上述地区定殖的风险可能性存在,但风险较小.另外,折刀法(Jackknife)分析显示,6种环境因子,例如地面霜冻频率、年平均降雨量、十月降雨量、四月降雨量、年最低温度以及蒸气压,对于番荔枝实蝇在全球和局部地区的分布模式有显著的影响.     

芒和五节芒在中国的潜在分布

利用最大熵模型, 将我国现有芒(Miscanthus sinensis)和五节芒(M. floridulus)的地理分布信息与19个降水及温度等气候因子相拟合, 预测了芒和五节芒在我国的潜在分布区域, 并推测出芒和五节芒的基本生态位。结果显示: 芒的潜在适生区包括四川西部、陕西北部、宁夏、内蒙古中部、黑龙江、吉林西部、辽宁西部、青海东南部等地区, 其基本生态位参数为: 最暖季节降水量为400-1 000 mm, 平均8月降水量为100-350 mm, 7月平均最低气温为15 ℃, 平均7月降水量为100-350 mm, 11月平均最高气温为-10-22 ℃, 最干月平均气温为-15-20 ℃, 平均12月降水量为100 mm以下; 五节芒的潜在适生区为云南、陕西、山西、宁夏、河南、山东、吉林、辽宁以及四川西部、甘肃南部和内蒙古东部等地区, 其基本生态位参数为: 最暖季节降水量至少在400 mm以上, 平均6月降水量为150-550 mm, 7月平均最低气温在15-30 ℃之间, 6月平均最低气温为10 ℃, 平均4月降水量为50-100 mm。结果表明, 在进行遗传改良的前提下, 我国有丰富的适合栽植芒和五节芒的土地资源。     

入侵害虫菊方翅网蝽在中国的潜在分布预测

菊方翅网蝽Corythucha marmorata(Uhler,1878)是我国新近发现的外来入侵害虫,研究明确菊方翅网蝽在我国的潜在分布范围对其监测预警及科学防控具有重要意义。本研究根据菊方翅网蝽的地理分布数据及相关环境变量,运用Maxent生态位模型与ArcGIS预测了菊方翅网蝽在中国的潜在地理分布范围。预测结果表明:菊方翅网蝽在我国的适生区主要分布于100°～125°E,20°～40°N的亚热带、暖温带区域,其中高适生区主要集中在长江中下游地区,包括浙江、江苏、湖南、上海大部分地区、安徽南部、湖北南部、江西西部及南部、贵州东部、福建东部、广西北部、山东中部、河南南部以及重庆、台湾局部;此外,极端气温、平均气温、最干月份降雨量对菊方翅网蝽的潜在分布影响较大。菊方翅网蝽已在我国成功入侵并迅速蔓延成灾,应在疫区边缘地带加强监测,并采取措施防止其进一步扩散。     

气候变化情境下木荷潜在地理分布及生态适宜性分析

木荷（Schima superba）是中国南方常见的速生阔叶耐火树种之一,在生态修复和维持环境稳定方面具有重要作用,因而预测其潜在适生区变化具有重要意义。基于Maxent模型,结合158条木荷在中国的有效分布记录和筛选后的8个生物气候变量,预测木荷现代和2041-2060年3种气候排放浓度情景下木荷的潜在分布区。Maxent模型的可靠性以受试者工作特征曲线来判断,综合气候因子的贡献率、置换重要值和刀切法检验,来分析并探讨制约木荷地理分布的重要气候因子。结果表明：（1） Maxent模型训练和测试数据的AUC值分别为0.936±0.001和0.943±0.008,表明模型的预测精度极高;（2）现代木荷潜在地理分布的总适生区面积198.87×10⁴ km²,占中国国土面积的20.7%,主要位于在中国长江以南的大部分地区,包括福建、浙江、台湾、江西、湖南、广东、海南、广西、贵州、云南、重庆,以及四川、河南和安徽部分区域;（3）制约木荷潜在地理分布的重要生物气候因子为降水因子（最干月降水量26-214 mm、年降水量800-2500 mm和最湿季降水量590-1030 mm）和温度因子（气温年变化8.5-28℃）;（4）使用Maxent模型模拟的木荷适生区面积大于其实际分布面积,暗示着木荷人工林在未来有较大的发展潜力。本研究结果可为木荷种质资源收集和人工林种植奠定理论基础,对维持区域性生态系统安全和稳定、缓解气候变化和2060年实现碳中和目标具有战略性意义。     

气候变化情景下基于最大熵模型的青海云杉潜在分布格局模拟

青海云杉(Picea crassifolia)是我国青藏高原东北缘特有树种,在维系我国西北地区生态平衡、水土保持、水源涵养和生物多样性等方面发挥着重要作用。基于其分布范围内的69个地理分布样点,利用最大熵(Maxent)模型对现实气候条件下青海云杉的潜在分布及其分布的主导气候因子进行分析,同时结合3种大气环流模型模拟青海云杉在3种气候变化情景(温室气候排放量不同)下未来2050s和2080s潜在分布区的变化。结果表明:Maxent模型对青海云杉潜在分布区的预测具有极高的准确度,所有模型的平均受试者工作特征曲线下面积(AUC测试值)均高于0.99;Jackknife检验和气候因子响应曲线表明年最低降雨量是限制青海云杉分布的主导因子;当前青海云杉的潜在分布区主要集中于青海东部、甘肃东南部、宁夏大部分地区、西藏东部、四川西部山区以及陕西、新疆和内蒙古部分地区。在未来3种增温情景下,青海云杉在2050s和2080s的潜在分布总面积与当前相比变化不明显,但不同适生等级的潜在分布面积变化较大,其中,中度适生区和低度适生区受气候增温影响显著,中度增温下这些区域在2080s的面积明显增大,而高度适生区(核心分布)则在所有增温情景下均呈缩小趋势。同时,在未来3种增温情景下,青海云杉在2050s和2080s的潜在分布区有向北移动趋势,但其心分布区域(高度适生区)仍然以青海东部、甘肃北部为主,无明显变迁趋势。从气候因素角度考虑,本研究表明未来气候变化情景下,青海云杉依然在西部高山地区,特别是作为我国重要生态屏障的祁连山、贺兰山等山区具有重要的经济价值并将持续其生态服务功能。     

气候变化情景下大沙鼠潜在地理分布

大沙鼠（Rhombomys opimus）是中亚地区典型的荒漠啮齿动物,其采食和掘洞行为造成了荒漠林和荒漠草原退化加剧,生态环境恶化。基于大沙鼠分布数据、气候、土壤和地形因子数据,采用MaxEnt模型预测大沙鼠在当前气候和温室气体低、中、高3种浓度排放情景下2050年和2070年的潜在适生区,分析亚洲大陆未来气候条件下大沙鼠适生面积和分布格局的变化趋势,探讨影响大沙鼠分布的主要环境因子。结果表明：模型AUC（Area Under Curve）值达到0.9以上,预测的准确性达到"极好"。经刀切法分析（Jackknife）表明,影响大沙鼠在适生区分布最主要的环境变量为温度季节性变化的标准差、土壤基本饱和度、最干季度降水量、最暖季度降水量和土壤可交换钠盐。Rcp2.6、Rcp4.5和Rcp8.5三种气候场景下2050年高适生区面积较当前分别增长15.78%、15.10%和13.44%;Rcp2.6、Rcp4.5和Rcp8.5三种气候场景下2070年高适生区面积较当前增长8.32%、13.18%和18.18%。中国大沙鼠适生区范围内,新疆所分布的大沙鼠适生区分布范围变化较大,3种情景模式下大沙鼠的适生区位置向新疆北部扩张;甘肃适生区位置向西北部扩张;内蒙西北部和阿拉善地区大沙鼠的适生区位置向四周扩张。研究揭示了未来气候下大沙鼠高适生区范围和空间变化,并得到影响其分布的主要环境变量,对其防控具有重要意义。     

气候变化下栓皮栎潜在地理分布格局及其主导气候因子

栓皮栎(Quercus variabilis)是东亚天然分布最广泛的树种之一。利用最大熵(Maxent)模型对现实气候条件下栓皮栎在东亚地区的潜在分布及其分布的主导气候因子进行分析,同时结合全新世中期(6000年前)和未来气候(2050年)来模拟和预测气候变化背景下栓皮栎潜在分布格局的变化。结果表明:现实气候条件下东亚栓皮栎适生区(适宜生境和低适宜生境)面积占总研究区面积的21.88%,主要集中在东亚南部区域,在我国北起陕西中部、山西和河北南部边缘、山东,西起甘肃东部边缘、四川中东部、云南、西藏东部边缘,一直到东部沿海区域,同时在朝鲜半岛南部和日本中南部也有分布,其中适宜生境面积占研究区总面积的5.69%,主要集中在秦岭山脉、大巴山脉、伏牛山、云南的云贵高原、罗霄山脉、南岭山脉、武夷山和台湾岛;气候变化情景下,栓皮栎的适生区分布面积变化较小,但其适宜生境的分布范围却发生了较大的变化,随着全球气候的波动性变化,适宜生境分布范围逐渐向西部秦岭山脉、大巴山脉、四川、重庆和云贵高原等区域集中,并使该分布中心的适宜生境面积逐渐扩大;影响栓皮栎分布的主要气候因子为最冷月的最低气温(Bio6)、最冷季平均气温(Bio11)和年降水量(Bio12),三者的贡献率分别为48.6%、21.4%和14.2%。     

基于Maxent模型的未来气候变化情景下胡杨在中国的潜在地理分布

胡杨(Populus euphratica)是全世界干旱和半干旱区急需优先保护的林木基因资源,预测未来气候变化情景下胡杨在中国的潜在地理分布将为胡杨种群资源的保护和管理提供科学依据,并为绿洲恢复过程中胡杨的合理种植和配置提供有价值的理论指导。本研究基于胡杨在中国地区的92条有效分布记录和10个环境因子变量,利用Maxent模型和ArcGIS软件预测了未来气候变化情景下胡杨在中国的潜在地理分布,综合环境因子变量贡献率及置换重要值、刀切法检验评估制约现代胡杨潜在地理分布的重要因子,采用响应曲线确定环境因子变量的适宜区间,定量确定胡杨未来受威胁的潜在地理分布区域和面积。结果表明：（1）Maxent模型的预测准确度极高,受试者工作曲线面积（AUC值）达0.932,现代胡杨潜在地理分布的总适生区面积为289.94×104km2,主要位于内蒙古中西部地区（额济纳旗和阿拉善地区）、新疆大部分地区、甘肃北部和西北部地区、青海中西部地区和宁夏北部地区;（2）影响胡杨的潜在地理分布的主要环境因子变量为气温因子变量（年均温和最冷月最低温）和降水因子变量（最湿月降水量和最干季降水量）,最湿月降水量是影响胡杨潜在地理分布的关键因素;（3）在未来4种气候变化情景下,胡杨不同等级潜在地理分布区的面积较现代潜在地理分布区面积均有不同程度的缩小,且整体上看胡杨的潜在地理分布区有向高海拔区域迁移的趋势。     

基于MaxEnt模型的气候变化下杏树鬃球蚧在中国的潜在地理分布预测

【目的】杏树鬃球蚧在天山野果林暴发成灾,对野杏的生存以及野果林的生态环境造成严重影响,同时危害栽培果园。预测未来气候情景下杏树鬃球蚧在中国的潜在适生区,可为农林部门开展防控工作提供科学依据。【方法】收集关于杏树鬃球蚧的分布点与气候数据,采用MaxEnt模型预测其在中国的适生区。【结果】MaxEnt模型中的AUC值为0.949,表明预测结果具有较高可靠性。最冷月最低温度、最冷季度降水量、温度季节性变化标准差、等温性是影响杏树鬃球蚧潜在分布的主要环境变量。杏树鬃球蚧在我国分布区较为广阔。该模型预测杏树鬃球蚧的高适生区主要分布在辽宁省环渤海区域、河北、北京、天津、山东、河南、陕西省中部、甘肃和陕西的临界处,以及新疆维吾尔自治区的伊犁哈萨克自治州、巴音郭楞蒙古自治州和阿勒泰地区,中适生区以高适生区为核心向外蔓延。其高、中、低适生区分别占中国陆地面积的4.62%、6.69%、18.77%。未来气候情景下,高适生区呈进一步扩散趋势,尤其在新疆地区,高适生区大面积扩散。【结论】杏树鬃球蚧适生区面积占陆地面积的30.08%,有继续扩张和大面积暴发风险,在我国西北地区、华北地区面临的形势更为严峻,应引起相关部门关注。     

气候变化下中国八种梧桐属树种潜在适生区模拟

中国梧桐属（Firmiana）在世界梧桐属中占比较大,且除梧桐外其余种均为中国特有且分布范围狭窄的植物种,灭绝风险大,研究气候变化对中国梧桐属树种的影响对于维护生物多样性具有重要的意义。结合多时期第六次国际气候耦合模式比较计划（CMIP6）气候变量数据和中国八种梧桐属树种的分布数据,基于R语言kuenm程序包优化的最大熵（Maxent）模型模拟分析中国八种梧桐属树种在多尺度下的潜在适生区,得出梧桐属最适宜的模拟尺度、潜在适生区的面积变化和迁移方向、梧桐属多样性保护关键区域及保护空缺。结果表明：（1）梧桐属最适宜的模拟尺度为亚洲;（2） Maxent模型的接收者操作特征曲线下面积（AUC）值均大于0.9,表明模型对梧桐属潜在适生区预测结果具有较高准确度;（3）气候变化影响下除云南梧桐（Firmiana major）外其它树种的潜在适生区都将在未来有所扩大;（4）中国八种梧桐属树种潜在适生区迁移方向主要为东西向,南北向大跨度迁移较少,纬度变化不大;（5）丹霞梧桐（Firmiana danxiaensis）的稳定潜在适生区最小;（6）中国梧桐属多样性保护关键区域主要分布于广西壮族自治区及云南、广东、海南等省区;（7）中国梧桐属多样性保护空缺区域主要分布于广西壮族自治区中部及海南省北部;（8）梧桐属多样性保护关键区域正在为人造地表所侵蚀。研究分析气候变化对中国八种梧桐属树种的影响及其潜在适生区变化、中国梧桐属多样性保护状态,可为中国梧桐属建立多样性保护廊道提供相关建议,为制定多样性保护规划及相应措施提供参考。     

末次盛冰期以来红豆树在不同气候变化情景下的分布动态

红豆树(Ormosia hosiei Hemsl. et Wils.)是中国的特有种,具有极高的经济价值、景观价值和药用价值。由于木材珍贵,人工盗伐严重,其种群数量和分布范围不断减少,被世界自然保护联盟(IUCN)濒危物种红色名录列为近危物种。气候变化会对物种的分布造成严重影响,理解该影响将有助于物种保护策略的制定,尤其是为濒危物种未来的保护提供重要参考。本研究基于红豆树在中国的地理分布数据,借助经相关性分析后筛选出的9个变量因子,利用最大熵模型(MaxEnt)模拟红豆树在末次盛冰期、全新世中期、当代和未来气候情景下的潜在分布区和影响其分布的主导环境因子,并且通过空间分析模拟其在不同气候变化情景下空间分布格局的变化。结果表明MaxEnt在各时期训练集和测试集的AUC(受试者工作特征曲线下的面积)均值均达到0.9以上,表明模型有很好的预测能力。刀切法(Jackknife)表明年均降水量、温度季节性变化标准差和昼夜温差月均值是主导其分布的三大因子,累计贡献率达到91.8%。将模拟结果导入到ArcGIS后,处理得到不同时期红豆树适生区空间分布格局变化。结果表明,自末次盛冰期以来,红豆树的适生...     

不同气候情景下四子柳的亚洲潜在地理分布格局变化预测

四子柳(Salix tetrasperma Roxb.)为杨柳科柳属为数不多的分布区扩展到热带的物种之一,广泛分布于华南和东南亚地区,但其生境破碎化较为严重。该物种具有较高的园林绿化价值和生态价值,预测不同气候情景下该物种的地理分布将为四子柳的资源开发和合理利用乃至柳属的起源和分化研究提供重要的科学依据。利用四子柳全面且精确的分布信息和高分辨率环境数据,基于Maxent模型和ArcGIS空间分析,构建其末次间冰期、末次盛冰期、当代以及未来(2050,2070)的潜在空间分布格局,评价环境因子对分布模型的重要性,定量确定未来受到威胁的适宜生境区域和面积。结果表明,四子柳目前的适宜生境面积为234.65×10~4 km~2,主要位于东亚、南亚和东南亚的热带地区,气温年较差和年均降水量是限制其分布的主要环境因子。总体而言,从冰期至未来,四子柳的分布中心有南北往返迁移的趋势。四子柳在末次盛冰期时向西北方向扩张并开始出现于热带地区,表明该物种开始适应热带气候,进入全新世中期以后种群收缩并退缩到云贵高原的河谷地区和印度尼西亚的平原区域。随着全球气候变暖,在不同二氧化碳浓度路径下2050年和2070年四子柳的潜在适生境有可能破碎化增加,建议对缅甸东部及中部成片减少的边缘群体进行实时监测。