中国谷子潜在地理分布的多模型比较

谷子是中国干旱和半干旱区主要的粮食作物之一.它耐旱、耐瘠薄、抗逆性强、适应性广,是未来应对干旱形势的重要战略储备作物.本文基于谷子的157个地理分布点数据,利用中国谷子产量与环境指标的相关性分析,选出10个气候指标、7个土壤指标和3个地形指标,采用MaxEnt、EMFA、RF和GAM共4个物种分布模型,分析中国谷子的潜在适宜性分布.结果表明: 4种模型均可成功模拟谷子的潜在地理分布,其中,MaxEnt模型的模拟效果最好.选用的环境指标中,水热条件对谷子生长最敏感.模型结果结合ArcGIS空间分析模块的结果表明,中国谷子的潜在适宜生长区(最适宜区和适宜区)总面积为55.68×10⁴ km²,远远大于当前谷子的实际种植面积,主要集中在东北地区的东北平原、长白山以南与牡丹江流域,华北地区的淮河以北,华中地区汉江以东与大别山以北,西北地区的黄土高原、鄂尔多斯高原南部、祁连山脉东部、天山山脉东部与阿尔泰山脉,西南地区的重庆以北和贵州西部局地区域.     

利用最大熵模型和规则集遗传算法模型预测孑遗植物裸果木的潜在地理分布及格局

利用野外调查的16个居群分布点和7个环境因子图层, 选择最大熵模型(MAXENT)和规则集遗传算法模型(GARP), 在地理和环境空间上模拟了第三纪孑遗植物裸果木(Gymnocarpos przewalskii)在中国西北地区的潜在分布。结果表明: (1)裸果木的潜在适生区全部集中在西北荒漠区, 其中最佳适生区主要集中在3个区域, 一是河西走廊中部和玉门以西、宁夏北部及内蒙古乌拉特后旗; 二是塔里木盆地西北缘; 三是柴达木盆地西北缘两片极小的高度适生区。裸果木的生态位被确定在一个较广的干旱环境空间: 适生区极端最高气温基本上在29.2-36.8 ℃之间, 极端最低气温在-18.3至-13.4 ℃之间; 年平均降水量40-200 mm; 潜在蒸发率在3-15之间。(2) MAXENT和GARP模型都较好地预测了裸果木的潜在分布, 但GARP产生了相对较大、较连续的潜在分布区, 部分过预测了破碎化生境; 而MAXENT预测到的潜在分布区, 在不同区域具有不同的环境适生性指数, 而且成功地排除了不合理的破碎化分布, 从而更直观地展示了裸果木的潜在分布格局和生态位要求。     

基于MaxEnt模型不同气候变化情景下的豆梨潜在地理分布

针对豆梨的原生境保护和资源利用问题,本研究基于豆梨全球236个分布点和19个环境因子,利用最大熵模型(MaxEnt)和地理信息系统(GIS)预测了豆梨在不同气候条件下的全球生态适宜区.结果表明: 豆梨的生态适宜区主要集中在北美洲、亚洲等地区,面积共约1.6×10⁷ km².其中,中国生态适宜度较高的地区主要分布在湖南省、湖北省、安徽省、江西省、江苏省、浙江省、福建省等地.影响豆梨地理分布的主要气候因子是年平均气温和年降水量,气温季节性变化次之.由模型预测可知,在不同的气候背景下,豆梨适宜生境和低适宜生境的面积有所不同.在空间分布上,豆梨适宜生境和低适宜生境的范围和几何中心都由东部向西部地区扩散,北美洲的适宜生境增长较快,而欧洲地区的低适宜生境增长较快.     

四川省岷江冷杉对气候变化的响应及其潜在分布格局

四川省是我国气候变化的敏感区之一,区域气候的暖干化趋势严重影响植物物种组成与分布。岷江冷杉(Abies faxoniana)作为我国特有种,其分布的动态变化对气候变化具有十分重要的指示作用。基于现有岷江冷杉分布数据、气候、土壤、地形等环境因子,运用最大熵模型(MaxEnt)预测当代气候条件下岷江冷杉潜在分布区,并分析未来时期(2050s和2070s)不同气候变化情景下(RCP2.6、RCP4.5和RCP8.5)岷江冷杉潜在适生分布区,筛选影响岷江冷杉分布的主导环境因子及阈值,探讨岷江冷杉分布对气候变化的响应机制。结果表明：(1)当前岷江冷杉的高适生区集中分布于岷江流域上游地区,在未来两个时期岷江冷杉潜在中、高适生区的面积较当代气候条件下适生区面积均有所增加,且适生区总体向四川南部扩张,北部适宜生境丧失。(2)岷江冷杉潜在中适生区在低排放浓度下(RCP2.6)面积占比最高,而潜在高适生区在高排放浓度下(RCP8.5)的面积占比最高。(3)影响岷江冷杉分布的主要环境因子分别是：降水季节性变异系数、气温年变化幅度、年降水量和海拔(累计贡献>70%)。适宜岷江冷杉潜在分布的环境条件是气温...     

气候变化对阔叶红松林潜在地理分布区的影响

物种地理分布主要取决于它对气候、地形等环境因子的适应性。基于22个环境因子和阔叶红松林的4类主要建群树种——红松、紫椴、水曲柳和蒙古栎的地理分布数据,采用最大熵模型模拟了阔叶红松林的潜在分布区域,并分析决定阔叶红松林地理分布的主要气候和地形因子,最后利用政府间气候变化专门委员会(IPCC)发布的3种排放场景(SRES-A2、SRES-A1B、SRES-B1)下2020、2050、2080年的气候数据预测阔叶红松林的未来潜在分布区。结果表明:各树种的受试者工作特征曲线下面积(AUC值)都大于0.8,说明模型有很好的预测能力;影响阔叶红松林分布的主导环境因子是年降雨量、季节性降雨量、海拔、年平均温度、最湿季度的平均温度。在基准气候条件下,阔叶红松林的高度适宜分布区主要分布在长白山和小兴安岭地区,占研究区总面积的11.69%,低度适宜区面积、不适宜区面积分别占研究区总面积的23%和65.31%。模型预测结果显示,未来在A2、A1B和B1气候情景下,阔叶红松林高度适宜区的南界与北界都向北移动,其面积有缩减的趋势,而低度适宜区的面积有增加的趋势。     

基于最大熵模型的中药材木贼麻黄潜在适生区

木贼麻黄(Ephedra equisetina)具有重要的药用和经济价值,在水土保持、生态平衡、荒漠化防治等方面也发挥重要作用,被列为国家二级保护植物。全球气候变化和人类活动严重威胁到木贼麻黄资源保护与利用,因而预测其潜在适生区空间变化具有重要意义。基于最大熵(MaxEnt)模型和ArcGIS软件,结合201条木贼麻黄在中国的有效分布记录和21个环境变量,预测在自然环境影响及人类活动干扰下中药材木贼麻黄适生区分布,并将木贼麻黄高适生区与土地利用类型进行叠加分析。结果表明:(1)模型添加人类活动因子后曲线下面积(AUC)值由0.919升至0.948,预测结果精度提高,而添加人类活动因子使木贼麻黄适生区面积严重减少,部分适生区等级退化,其分布愈破碎化稀疏零散;(2)木贼麻黄适生区多集中分布于新疆、甘肃、宁夏、陕西和山西地区,未来时期两种情景下木贼麻黄适生区都出现不同程度的扩张,RCP2.6-2050年相对当前时期高适生区面积增加最多;(3)参与建模的环境因子中贡献度较高的酸碱度(pH)、海拔(Altitude)和年平均温度(Bio1)对木贼麻黄地理分布有关键影响;(4)未来时期下木贼麻黄高适生区的几何质心出现不同程度向东迁移的现象,多分布于甘肃省境内;(5)木贼麻黄高适生区部分土地已经被开发利用,剩余未利用的土地分布于新疆、甘肃、内蒙古和宁夏地区,可在这些地区因地制宜发展木贼麻黄种植业。为木贼麻黄的就地、迁地保护、优先保护区以及人工种植选地提供理论依据,这对木贼麻黄资源保护与利用具有重要的现实意义。     

基于最大熵模型预测气候变化下河蚬在中国的潜在分布

河蚬是一类在我国广泛分布的底栖贝类,具有重要的经济价值及生态价值。近年来,河蚬野生资源量锐减,了解河蚬在国内的潜在分布能为河蚬的保护和合理利用提供重要参考。基于河蚬在中国的136个分布点和8个环境因子,采用ENMeval包和biasfile优化后的最大熵模型(MaxEnt)预测分别河蚬现代和未来(2041—2060年和2081—2100年) 6个气候情景下的潜在分布。综合Jackknife检验、置换重要值和环境因子贡献率评估影响现代河蚬潜在分布的主要因子,比较未来气候情景下潜在适生区差异从而分析预测河蚬适宜分布的变化。结果表明：(1)优化后的MaxEnt模型预测准确度极高,平均AUC值为0.900±0.037,平均AUC_DIFF值为0.019,现代河蚬潜在分布区域总面积为188.33×10⁴ km²,主要集中在长江流域、海河流域、淮河流域、珠江流域、东南沿海区域以及黄河流域下游和渤海湾沿岸区域。(2)影响河蚬潜在分布的主要环境因子为海拔、温度(年均温和温度年较差)和降水(年降水量)。(3)在未来6种气候情景下,河蚬主要潜...     

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

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

基于MaxEnt模型的大豆蚜全球潜在地理分布分析

【目的】为预测和分析大豆蚜Aphis glycines的全球潜在地理分布,研究大豆蚜分布与环境变量之间的联系。【方法】利用最大熵法生态位模型(maximum entropy niche-based modeling, MaxEnt)和地理信息系统软件ArcGIS,根据收集的大豆蚜已知分布点和环境变量,预测大豆蚜的全球潜在地理分布区,推测环境变量对大豆蚜分布的影响。【结果】结果表明,大豆蚜适生区主要分布在低海拔地区,高 度适生区集中在25°~50°N的中国、日本、韩国、朝鲜、加拿大、美国、意大利和格鲁吉亚。决定大豆蚜分布地点的关键环境变量为最暖季度降水量、最暖季度平均温度、最湿季度平均温度、最干月降水量、月平均昼夜温差和温度季节性变化标准差。【结论】大豆蚜潜在地理分布区域广泛,应在各国大豆农产品贸易时做好检验检疫工作,以防止大豆蚜的扩散。     

利用CLIMEX预测纳塔尔实蝇在中国的潜在地理分布

纳塔尔实蝇Ceratitis rosa Karsch属双翅目实蝇科腊实蝇属,为害30余种经济植物,被我国列为进境植物检疫性有害生物。本研究运用CLIMEX 3.0及ArcGIS 9.3对纳塔尔实蝇在我国目前及未来的潜在地理分布进行了预测。结果显示：在目前的气候条件下,纳塔尔实蝇在我国的潜在地理分布区为18.250°N-30.250°N,其中,华南和东南地区为高度潜在地理分布区,包括四川、重庆、云南、福建、广东、广西、海南以及台湾等地。在未来的气候条件下,2020和2050年潜在地理分布区的北界分别移至31.250°N和32.250°N,且中度潜在地理分布区北移明显;2100年潜在地理分布区的北界移至33.750°N,且高度适生区北移明显。因此建议目前应加强纳塔尔实蝇检疫措施,完善监测体系,监测网点主要设在我国的南方地区,尤其应对云南、广西、广东、福建和海南等地进行长期监测;同时,应对湖南、贵州、江西、湖北、江苏、安徽等省进行定期监测,严防该虫入侵。     

秦岭冷杉球果与种子的形态变异

在秦岭冷杉的天然分布区 ,分别选取核心区和边缘区进行取样 ,比较球果和种子的形态变异。核心区为秦岭地区 ,边缘区选在湖北神农架地区。根据秦岭冷杉分布特征 ,核心区沿海拔高度连续采样 ,边缘区选择两个不同海拔高度的种群 ,分单株采样 ,测量了球果、种子、种鳞、种翅等器官的 15个形态指标。相关分析表明 ,秦岭冷杉球果和种子形态特征 (变化 )与海拔高度显著相关 ,说明海拔因素影响着球果和种子形态变异。比较发现 ,核心区 (即秦岭地区 )的秦岭冷杉球果和种子的绝大部分指标均值显著大于边缘区 (即神农架地区 ) ,说明秦岭冷杉在核心区的生殖生长好于边缘区。相同海拔条件下 ,秦岭冷杉球果和种子形态在地区间差异显著 ,约 85 %的变异来自地区内个体间和个体内 ,说明秦岭冷杉的球果和种子形态特征变化受遗传控制更显著。在边缘区还比较了不同种群间、人工林和天然林在球果和种子形态特征的差异。球果和种子形态在同一地区两个种群间、在人工林与天然林间都有显著差异 ,大部分人工林球果和种子指标均值和变异幅度大于天然林 ,表明通过人为经营可以改善秦岭冷杉的一些生殖性状     

濒危植物秦岭冷杉种群数量动态

为了对濒危植物秦岭冷杉种群数量动态评价和预测,通过样地调查和数据统计,研究了秦岭冷杉种群的年龄结构、静态生命表及其与环境因子关系,运用时间序列模型预测了种群数量动态.结果表明,多数秦岭冷杉种群幼龄级个体数较少,中老龄个体数量较大,呈衰退趋势.仅处于低海拔地区的秦岭冷杉-木蓝-苔草群丛中的种群(D种群)由于立地条件较好,幼龄级个体数量相对丰富,种群稳定.不同秦岭冷杉种群生命表和存活曲线的分析表明,尽管生境条件差异,但存活曲线基本接近DeeveyⅢ型;不同种群偏离典型存活曲线的程度与幼苗缺乏程度有关,一般Ⅲ～Ⅴ龄级死亡率较高.时间序列分析表明,在未来20、40和80年中,不同秦岭冷杉种群均会呈现老龄级株数先增后减的趋势,种群稳定性长期维持困难.对影响秦岭冷杉种群增长的10个环境因子通过主成分分析(PCA)发现,乔木层盖度、土壤有机质含量和空气湿度对种群发挥有利影响,而人为干扰和光照强度对秦岭冷杉种群增长发挥不利影响.应充分利用秦岭冷杉性喜荫、耐寒、种子活力较强的特点,加强现有林分就地保护,重点是具有结实能力的中老龄个体;在阴坡地带,对林下灌木比较密集的群丛,通过砍灌、清理林下活地被物等抚育措施,为幼苗发育创造良好的环境条件;就地采种育苗,扩大人工种群.     

濒危植物秦岭冷杉地理分布和生物生态学特性研究

秦岭冷杉(Abies chensiensis)为中国特有三级保护植物,仅岛屿化分布于陕西、甘肃、四川、河南的秦巴山地和湖北神农架中高山地区,分布海拔1300—2300m。为了对秦岭冷杉种群进行深入研究和有效保护,通过样地调查和相关资料分析,对其地理分布、群落特性、种群特性进行了系统研究。秦岭冷杉群落包含种子植物40科63属102种,可划分成11个地理区系成分,其中温带分布型属占84．21％,热带分布型属占15．79％,温带分布是群落的基本特性。秦岭山区的秦岭冷杉群落可分为5个群丛,基本特点为：生境阴湿;群落郁闭度大,结构复杂;种群个体年龄差异大;天然更新主要发生在保护较好的阴坡林窗地段。秦岭冷杉苗期生长缓慢,在30—70年龄阶段,高生长旺盛;在40—100年龄阶段,胸径生长较快。秦岭冷杉开始结实年龄在30—40年龄阶段,结实间隔期为3—5年。秦岭山区成熟个体平均产球果8个／株,饱满种子56粒／球果。种群更新的关键问题是种子产量低,林下幼苗稀疏,更新缺乏后继资源。种子向幼苗转化率低是生活史中的脆弱环节。在未来经营保护中,应以就地保护为主,通过抚育管理,促进天然更新;丰年应适时采种,就地育苗,扩大人工种群。     

濒危植物秦岭冷杉种群空间分布格局及动态

为阐明秦岭冷杉种群分布格局动态特征,达到有效保护和恢复的目的,通过野外调查,采用相邻格子样方法,离散分布拟合、扩散系数、Morisita指数、均方-曲线等方法,对不同生境和年龄阶段的秦岭冷杉种群的空间分布格局及其动态进行了系统研究。结果表明：（1）秦岭冷杉种群空间分布格局在不同发育阶段和不同生境条件下基本属于聚集型分布,在25m^2、150m^2和300m^2的面积上聚集强度较大。天然条件下,种群聚集强度最大取样尺度,可能是其种群更新的最佳面积。（2）秦岭冷杉种群空间分布格局随海拔上升,聚集强度降低;聚集强度在海拔1300～1700m较小,在海拔1700～2100m地区较大;阴坡种群分布聚集强度大于阳坡。（3）随着种群年龄增加,分布格局由聚集型向随机型过渡,聚集强度减弱;聚集性最强的尺度规模在幼、中龄期为50m^2,老龄期为300m^2。（4）种群的空间分布格局与生物生态学特性、生境条件关系密切。在未来森林经营和种群恢复中,应针对生境和种群年龄,森林抚育管理作业面积应该尽量接近聚集强度最大的尺度。低海拔和阴坡应该适度间伐非目的乔木树种,高海拔和阳坡地区,应适度清理灌木和草本层,为群落内团块状更新创造空间;实行就地保护,减少人为干扰,促进天然更新保护策略。     

濒危植物秦岭冷杉补光育苗技术研究

秦岭冷杉种子生活力差、成苗率低(2.2%)是制约种群恢复的薄弱环节。通过对播种后苗期2年连续补光处理,幼苗表现出良好的生长效应,苗高10.5 cm(对照为7.3 cm),地径4.5 mm(对照为2.7 mm),保存率88.6%(对照66.4%),延长了幼苗生长期,增加了苗木生物量积累,加快了苗木培育进程。     

濒危植物秦岭冷杉群落数量特征及其动态

对石人山自然保护区秦岭冷杉群落特征和动态进行了研究,结果表明:秦岭冷杉群落共有种子植物61种,分属27科40属,群落垂直结构简单,可分为乔木层、灌木层和草本层。秦岭冷杉和华山松种群大、中胸径的个体占绝对优势,幼苗、幼树较少,属于衰退种群;湖北花楸小、中胸径的个体占优势,属于进展种群。通过聚类分析,将秦岭冷杉群落分为5个群丛类型:①秦岭冷杉 湖北花楸-山刺玫-结缕草群丛;②秦岭冷杉 华山松-秀雅杜鹃-藜芦群丛;③秦岭冷杉 血皮槭-黄刺玫-野菊群丛;④秦岭冷杉 华山松-箭竹-牛蒡群丛;⑤秦岭冷杉 锐齿栎-箭竹-水芹群丛。5个群丛类型均处于动态变化中,种群恢复应实行就地保护,减少人为干扰,促进天然更新。     

基于最大熵生态位模型的中华穿山甲潜在适宜生境预测

中华穿山甲（Manis pentadactyla）属于全球极度濒危物种,也是我国一级保护动物。对中华穿山甲的非法捕杀曾导致其种群数量锐减。但是,近年来相关研究报道较少,穿山甲分布状况不明,极大地制约了对该物种的有效保护。搜集了近年来国内中华穿山甲的救护记录和救护新闻,甄别出67个记录分布点,利用最大熵模型软件（MaxEnt）进行因子筛选,结果表明最冷季度降水量、人口密度、年降水量、坡度、坡向、海拔等6个环境变量是与中华穿山甲分布显著相关的影响因子。基于6个主导环境变量构建的MaxEnt模型AUC平均值为0.961±0.014,预测结果达到极好标准。刀切法（Jackknife）表明,其中最冷季度降水量、年降水量、人口密度和海拔是影响中华穿山甲分布的主要因素。中华穿山甲适宜生境（出现概率大于0.498）具有以下特点：最冷季度降水量141.22-439.46 mm,年降水量1471.67-2386.56 mm,人口密度≥390人/km²,海拔<316.98 m。该模型预测中华穿山甲在我国的潜在分布适宜区主要位于我国长江以南地区,总面积约为74.27×10⁴ km²,占国土面积的7.73%,主要集中在江西、广东、湖南和广西省,面积分别占该区域的97.58%,89.65%,76.90%和73.08%;其次是浙江、福建、台湾和安徽省。湖北、江苏、四川、云南、贵州等省份也有中华穿山甲的零星分布。湖北东南部、江苏南部、浙江西南部和福建西北部等与江西接壤的区域也是中华穿山甲的重要潜在分布适宜区。明确中华穿山甲的潜在分布适宜区,可为该物种的种群保护和栖息地管理提供科技支撑。     

基于MaxEnt模型预测白唇鹿的潜在分布区

研究物种分布格局对制定有效保护措施至关重要。本研究从文献记录和实地调查中获取97个分布点, 筛选出7个重要环境变量, 采用最大熵(MaxEnt)模型分析了中国青藏高原特有种白唇鹿(Przewalskium albirostris)的潜在分布区域。结果显示: 白唇鹿潜在分布区主要位于青藏高原东部, 适宜生境呈现出以西藏、青海和四川三省交界地带为核心向四周扩散的模式。刀切法(Jackknife)分析表明, 地形、气温季节变化和年降水量是决定白唇鹿分布的重要因素, 人类影响变量贡献相对较小。白唇鹿种群现状尚不清晰, 建议对该物种进行全面深入调查。     

Comparative analysis of genetic diversity and population genetic structure in Abies chensiensis and Abies fargesii inferred from microsatellite markers

Abies chensiensis Tieghem and Abies fargesii Franchet are two closely related tree species of Pinaceae endemic to China. A. chensiensis is usually found scattered in small forest fragments, whereas A. fargesii is a dominant member of coniferous forest. To evaluate the genetic effect of fragmentation on A. chensiensis, a total of 24 populations were sampled from the whole distribution of the two species. Seven nuclear microsatellite loci were employed to analyze comparatively the genetic diversity and population genetic differentiation. Both A. chensiensis and A. fargesii have high level within-population genetic diversity and low inter-population genetic differentiation. Low microsatellite differentiation (2.1%) between A. fargesii and A. chensiensis was observed. But microsatellite marker was able to discriminate most populations of these two species. Compared to A. fargesii, A. chensiensi has lower allelic diversity and higher genetic differentiation among populations. It suggested the existence of negative genetic impacts of habitat fragmentation on A. chensiensis.     

Species delimitation and phylogeography of the Abies chensiensis complex inferred from morphological and molecular data

Species delimitation in the Abies chensiensis complex, consisting of A. chensiensis, A. recurvata, A. ernestii and A. ernestii var. salouenensis, has been the subject of a long‐term dispute, as various taxonomic combinations have been proposed. We combined different lines of evidence, including morphological characters and mitochondrial and plastid DNA sequence data, to assess species delimitation. Plastid DNA sequence variation was highly consistent with the morphological result and separated A. recurvata from the other three taxa, but the mitochondrial DNA genealogy was more strongly related to geographical distribution rather than species delimitation, as proposed in previous studies. On the basis of morphological characters and plastid DNA sequence variation, the current species status of A. recurvata was well supported, and we also accepted A. ernestii and A. ernestii var. salouenensis as varieties of A. chensiensis. Furthermore, the phylogeographical history of A. chensiensis was surveyed by molecular data and ecological niche modelling. On the basis of molecular data, the current phylogeographical pattern of the A. chensiensis complex was probably shaped by recent rapid expansions from multiple refugia in the Hengduan Mountains and habitat contraction after the last glaciation. This hypothesis was supported by the results of ecological niche modelling and previous fossil records. Our findings suggest that these fir taxa might have continued to expand their range after the largest glaciation until the end of the last glaciation, but have contracted since then with the increasing temperature and habitat changes. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 177 , 175–188.