生态位模型的基本原理及其在生物多样性保护中的应用

生态位模型是利用物种已知的分布数据和相关环境变量,根据一定的算法来推算物种的生态需求,然后将运算结果投射至不同的空间和时间中来预测物种的实际分布和潜在分布.近年来,该类模型被越来越多地应用在入侵生物学、保护生物学、全球气候变化对物种分布影响以及传染病空间传播的研究中.然而,由于生态位模型的理论基础未被深入理解,导致得出入侵物种生态位迁移等不符合实际的结论.作者从生态位与物种分布的关系、生态位模型构建的基本原理以及生态位模型和生态位的关系等方面探讨了生态位模型的理论基础.非生物的气候因素、物种间的相互作用和物种的迁移能力是影响物种分布的3个主要因素,它们在不同的空间尺度下作用于物种的分布.生态位模型是利用物种分布点所关联的环境变量来模拟物种的分布,这些分布点本身关联着该物种和其他物种间的相互作用,因此生态位模型所模拟的是现实生态位(realized niche)或潜在生态位(potential niche),而不是基础生态位(fundamental niche).Grinnell生态位和Elton生态位均在生态位模型中得到反映,这取决于环境变量类型的选择、所采用环境变量的分辨率以及物种自身的迁移能力.生态位模型在生物多样性保护中的应用主要包括物种的生态需求分析、未知物种或种群的探索和发现、自然保护区的选择和设计、物种入侵风险评价、气候变化对物种分布的影响、近缘物种生态位保守性及基于生态位分化的物种界定等方面.     

基于本土最优模型模拟入侵物种水盾草在中国的潜在分布

生态位模型预测存在不确定性, 不同模型预测结果差别较大。在生态位保守的前提下, 在本土区域构建经典生态位模型, 利用入侵地独立样本数据检验并选择最优模型, 具有独特优势, 可为入侵物种风险分析提供可靠参考。水盾草(Cabomba caroliniana)是一种恶性水生入侵杂草, 原产于南美洲, 已在我国多个省市建立种群, 本文基于本土最优模型预测其在我国的潜在分布, 以期为其风险分析和综合治理提供依据, 并通过水盾草案例探讨如何提高生态位模型预测准确性的方法。本文按时间顺序梳理了水盾草在我国的分布记录, 然后根据水盾草已有分布记录和其所关联的环境因子比较了不同地理种群所占有的气候生态空间, 测试水盾草在世界入侵过程中的现实生态位保守性。采用两组环境变量和5种算法在南美洲本土地区构建10种生态位模型, 并将其转移至我国, 基于最小遗漏率和记账错率, 利用我国(入侵地)的样本数据选择最优模型预测水盾草在我国的适宜生态空间和潜在分布。研究发现当前水盾草在我国的分布集中在东部水域充沛地区, 沿京杭运河和南水北调工程等向北扩散。生态空间比对中发现水盾草在亚洲与其他大洲所占有的生态空间具有一定的重叠, 其在我国的入侵过程中生态位是保守的。与本土空间相比, 水盾草在我国所占有的生态空间存在较大的生态位空缺, 表明水盾草在我国的潜在分布范围较大。生态位模型预测显示水盾草的适生区主要分布于我国的北京、上海、山东、浙江、江苏、安徽、湖北和湖南等省(市)。水盾草的潜在分布区多聚集在我国东南部, 该地区河流、湖泊、运河和渠道较为密集, 人类活动及自然天敌的缺乏容易助长其入侵趋势, 应在这些适宜地区开展调查, 及时发现疫情并采取相应措施。     

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

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

Realized niche shift associated with Galinsoga quadriradiata (Asteraceae) invasion in China

粗毛牛膝菊在中国的入侵与生态位漂移有关 在外来物种入侵和扩散过程中,生态位的漂移可能起到了重要作用。粗毛牛膝菊(Galinsoga quadriradiata) 在中国已造成了较为严重的入侵,占据了许多与其原产地不同的气候区。为此,本研究力图揭示粗毛牛膝菊入侵过程中的气候生态位漂移,分析其在该物种入 侵中国过程中可能发挥的作用。本研究结合粗毛牛膝菊原 产地和入侵地的分布点与气候数据, 采用Maxent模型预测了其在中国潜在的分布,并采用主成分分析的方法评估 了在入侵中国过程中粗毛牛膝菊气候生态位的漂移。模型结果显示,该物种原产地种群和入侵地种群之间只 有32.7%的生态位重叠,两个种群的生态位相似性较低(Schoener's D = 0.093, P < 0.005),这暗示了在其入侵过程中发生了生态位漂移。相比于其原产地种群,其在中国的入侵种群气候生态位的整体范围和中心都明 显地漂移向了温度更低、降水更少的区域;中国南方大部分区域属于粗毛牛膝菊的稳定适生区,而位于入侵 前沿的北方地区则存在局域适应和潜在拓殖区域。这些研究结果说明,粗毛牛膝菊在中国的入侵种群仍处于准平衡阶段,未来有可能继续向新的适生区扩散入侵,其生态位的变化有力地解释了为什么该物种在中国的入侵性强、危害范围大。     

MaxEnt模型参数设置对其所模拟物种地理分布和生态位的影响——以茶翅蝽为例

【目的】生态位模型被广泛应用于入侵生物学和保护生物学研究,现有建模工具中,MaxEnt是最流行和运用最广泛的生态位模型。然而最近研究表明,基于MaxEnt模型的默认参数构建模型时,模型倾向于过度拟合,并非一定为最佳模型,尤其是在处理一些分布点较少的物种。【方法】以茶翅蝽为例,通过设置不同的特征参数、调控倍频以及背景拟不存在点数分别构建茶翅蝽的本土模型,然后将其转入入侵地来验证和比较模型,通过检测模型预测的物种对环境因子的响应曲线、潜在分布在生态空间中的生态位映射以及潜在分布的空间差异性,探讨3种参数设置对MaxEnt模型模拟物种分布和生态位的影响。【结果】在茶翅蝽的案例分析中,特征参数的设置对MaxEnt模型所模拟的潜在分布和生态位的影响最大,调控倍频的影响次之,背景拟不存在点数的影响最小。与其他特征相比,基于特征H和T的模型其响应曲线较为曲折;随着调控倍频的增加,响应曲线变得圆滑。【结论】在构建MaxEnt模型时,需要从生态空间中考虑物种的生态需求,分析模型参数对预测物种分布和生态位可能造成的影响。     

洲际入侵植物生态位稳定性研究进展

人类活动引起的大规模洲际物种交换与生物入侵, 改变了当地生态系统结构与功能, 使生物多样性受到日益严重的威胁。本文通过综合分析主要国家和地区入侵植物的地理起源, 发现洲际入侵主要包括东亚—北美、东亚—南美、欧洲—南非、欧洲—北美、欧洲—东亚、北美—大洋洲等, 这些洲际入侵造成的后果往往比陆内入侵更为严重。利用物种分布模型(SDMs)预测入侵物种潜在分布范围是有效管理和提早预防生物入侵的重要依据, 但这些模型的一个关键假定是: 入侵物种的生态位在空间和时间上是保守的、稳定的。然而, 对于远离原产地种群并能快速适应新生境的洲际入侵植物来说, 生态位可能发生显著的变化。入侵种能否在入侵地保持原有的生态位, 取决于制约其生态分布的限制因素和生态过程在不同地区间是否发生变化。本文中作者总结了洲际入侵与陆内入侵的生态与进化过程的异同点, 认为这些限制物种原产地分布的因素如扩散限制、种间互作、适应性进化、生态可塑性和种群遗传特性等均可能导致入侵物种生态位的改变。建议下一步的研究应该重视: (1)对生态位属性进行多尺度的研究, 包括时间、空间、环境或系统发育等几个方面; (2)对比生态位稳定与发生偏移的物种特性, 确定什么样的入侵物种更容易改变原有的生态位; (3)进行生态位时间动态格局研究, 探讨生态位变化的倾向、历史速率和偏移程度, 以便判定生态位变化趋势。这些研究结果将会进一步提高物种分布模型的预测能力, 有助于更为准确地揭示气候变化和物种入侵对生物多样性的影响。     

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

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

加拿大一枝黄花入侵对本土植物群落动态的影响及其机制

了解外来植物入侵对本土植物群落种群动态的影响对于植物入侵的防控极为重要。该文以加拿大一枝黄花(Solidago canadensis)入侵不同阶段的植物群落为研究对象, 对本土植物物种多样性以及常见优势种群的生态位变化进行了定量分析。结果表明: 加拿大一枝黄花氮素积累能力高于其他本土优势种群。随着加拿大一枝黄花入侵的深入, 本土植物群落的物种多样性呈现显著下降趋势; 氮素积累能力高的本土优势种群生态位宽度呈现明显的上升趋势, 而氮素积累能力低的本土优势种群生态位宽度则呈现明显下降的趋势; 本土优势种群的生态位重叠平均值呈现逐步下降的趋势。加拿大一枝黄花的入侵, 显著提高了土壤硝态氮含量, 而土壤铵态氮、有效磷、全磷和全氮含量显著降低。对氮素的积累能力决定了加拿大一枝黄花入侵后, 本土植物种群的动态变化格局。     

紫茎泽兰在中国的气候生态位稳定且其分布范围仍有进一步扩展的趋势

明确入侵物种气候生态位的时空动态格局, 不仅有助于理解入侵物种生态位的变化趋势, 还可以使物种分布模型更真实地预测物种的扩散动态。本文使用恶性入侵植物紫茎泽兰(Ageratina adenophora)在原产地及我国的分布点和环境因子数据, 结合生态位分析及物种分布模型, 分析了紫茎泽兰入侵至我国后其气候生态位的时空变化特征、地理分布格局变化及其限制因子。结果表明: (1)紫茎泽兰入侵至我国后其气候生态位在时间和空间上虽略有扩张, 但仍保持较高的稳定性。尤其是入侵的初始阶段(1980年之前)稳定性最高(NS = 0.995), 此后其气候生态位的稳定性逐渐降低, 具体表现为该植物在中国更倾向于分布在潮湿(高最暖季降水量)、寒冷(低最干季平均温)且平均日较差低的生境。(2)过去50年里紫茎泽兰的适宜分布范围显著增加(151.77%-271.14%), 且其适宜分布区的重心在2010年之前表现出向东北(贵州西南部、四川南部)扩散的趋势, 2010年以后表现出向西南(云南北部、贵州北部)扩散的趋势。在扩散过程中, 该物种的气候生态位稳定性也表现出随着纬度的增加而逐渐降低的趋势。(3)未来气候变化情景下, 紫茎泽兰的适宜分布范围将进一步增加(111.97%), 到2040年该物种基本可以扩散至其整个潜在适宜分布区, 仅在藏东南地区有小范围的适宜生境受扩散限制无法达到。     

应用生态位模型研究外来入侵物种生态位漂移

由于基础生态位和实际生态位的改变,外来入侵物种在入侵地成功定殖、扩散后常会发生生态位漂移,而物种生态位漂移往往很难直接证明。生态位模型在假设入侵物种的生态位需求保守的前提下,以物种在其原产地的生态位需求为基础,预测其在入侵地的潜在分布,通过比较预测分布与实际分布的差异可以从一定程度上得到外来入侵物种的生态位是否发生漂移的间接证据。以我国入侵杂草胜红蓟在原产地的生态位需求为基础,应用生态位模型预测其在其他地区的潜在分布。研究结果表明,生态位模型可以很好地预测胜红蓟在亚太平洋地区和非洲地区的分布,但在我国,其预测分布与实际分布存在较大差别。胜红蓟在我国预测分布主要为云南、海南、台湾部分地区,而胜红蓟入侵我国后现已广泛分布于长江以南地区,其实际分布比预测分布广泛得多,由此推测胜红蓟在入侵我国后其生态位已经产生了漂移。     

An evaluation of transferability of ecological niche models

Ecological niche modeling (ENM) is used widely to study species’ geographic distributions. ENM applications frequently involve transferring models calibrated with environmental data from one region to other regions or times that may include novel environmental conditions. When novel conditions are present, transferability implies extrapolation, whereas, in absence of such conditions, transferability is an interpolation step only. We evaluated transferability of models produced using 11 ENM algorithms from the perspective of interpolation and extrapolation in a virtual species framework. We defined fundamental niches and potential distributions of 16 virtual species distributed across Eurasia. To simulate real situations of incomplete understanding of species’ distribution or existing fundamental niche (environmental conditions suitable for the species contained in the study area; N* _F ), we divided Eurasia into six regions and used 1–5 regions for model calibration and the rest for model evaluation. The models produced with the 11 ENM algorithms were evaluated in environmental space, to complement the traditional geographic evaluation of models. None of the algorithms accurately estimated the existing fundamental niche (N* _F ) given one region in calibration, and model evaluation scores decreased as the novelty of the environments in the evaluation regions increased. Thus, we recommend quantifying environmental similarity between calibration and transfer regions prior to model transfer, providing an avenue for assessing uncertainty of model transferability. Different algorithms had different sensitivity to completeness of knowledge of N* _F , with implications for algorithm selection. If the goal is to reconstruct fundamental niches, users should choose algorithms with limited extrapolation when N* _F is well known, or choose algorithms with increased extrapolation when N* _F is poorly known. Our assessment can inform applications of ecological niche modeling transference to anticipate species invasions into novel areas, disease emergence in new regions, and forecasts of species distributions under future climate conditions.     

Ecological Niche Transferability Using Invasive Species as a Case Study

Species distribution modeling is widely applied to predict invasive species distributions and species range shifts under climate change. Accurate predictions depend upon meeting the assumption that ecological niches are conserved, i.e., spatially or temporally transferable. Here we present a multi-taxon comparative analysis of niche conservatism using biological invasion events well documented in natural history museum collections. Our goal is to assess spatial transferability of the climatic niche of a range of noxious terrestrial invasive species using two complementary approaches. First we compare species’ native versus invasive ranges in environmental space using two distinct methods, Principal Components Analysis and Mahalanobis distance. Second we compare species’ native versus invaded ranges in geographic space as estimated using the species distribution modeling technique Maxent and the comparative index Hellinger’s I. We find that species exhibit a range of responses, from almost complete transferability, in which the invaded niches completely overlap with the native niches, to a complete dissociation between native and invaded ranges. Intermediate responses included expansion of dimension attributable to either temperature or precipitation derived variables, as well as niche expansion in multiple dimensions. We conclude that the ecological niche in the native range is generally a poor predictor of invaded range and, by analogy, the ecological niche may be a poor predictor of range shifts under climate change. We suggest that assessing dimensions of niche transferability prior to standard species distribution modeling may improve the understanding of species’ dynamics in the invaded range.     

Using ecological niche modelling to evaluate niche stability in deep time

Aim To investigate relative niche stability in species responses to various types of environmental pressure (biotic and abiotic) on geological time‐scales using the fossil record. Location The case study focuses on Late Ordovician articulate brachiopods of the Cincinnati Arch in eastern North America. Methods Species niches were modelled for a suite of fossil brachiopod species based on five environmental variables inferred from sedimentary parameters using GARP and Maxent . Niche stability was assessed by comparison of (1) the degree of overlap of species distribution models developed for a time‐slice and those generated by projecting niche models of the previous time‐slice onto environmental layers of a second time‐slice using GARP and Maxent , (2) Schoener’s D statistic, and (3) the similarity of the contribution of each environmental parameter within Maxent niche models between adjacent time‐slices. Results Late Ordovician brachiopod species conserved their niches with high fidelity during intervals of gradual environmental change but responded to inter‐basinal species invasions through niche evolution. Both native and invasive species exhibited similar levels of niche evolution in the invasion and post‐invasion intervals. Niche evolution was related mostly to decreased variance within the former ecological niche parameters rather than to shifts to new ecospace. Main conclusions Although the species examined exhibited morphological stasis during the study interval, high levels of niche conservatism were observed only during intervals of gradual environmental change. Rapid environmental change, notably inter‐basinal species invasions, resulted in high levels of niche evolution among the focal taxa. Both native and invasive species responded with similar levels of niche evolution during the invasion interval and subsequent environmental reorganization. The assumption of complete niche conservatism frequently employed in ecological niche modelling (ENM) analyses to forecast or hindcast species geographical distributions is more likely to be accurate for climate change studies than for invasive species analyses over geological time‐scales.     

Selecting Biological Meaningful Environmental Dimensions of Low Discrepancy among Ranges to Predict Potential Distribution of Bean Plataspid Invasion

Background

The Bean plataspid (Megacopta cribraria) (Hemiptera: Pentatomidae), native to Asia, is becoming an invasive species in North America; its potential spread to soybean producing areas in the US is of great concern. Ecological niche modelling (ENM) has been used increasingly in predicting invasive species'' potential distribution; however, poor niche model transferability was sometimes reported, leading to the artifactual conclusion of niche differentiation during species'' invasion.

Methodology/Principals

We aim to improve the geographical transferability of ENM via environmental variable selection to predict the potential distribution of Bean plataspid invasion. Sixteen environmental dimensions between native and introduced Bean plataspid populations were compared, and classified into two datasets with different degrees of discrepancy by the interquartile range (IQR) overlap in boxplot. Niche models based on these two datasets were compared in native model prediction and invading model projection. Classical niche model approaches (i.e., model calibrated on native range and transferred outside) were used to anticipate the potential distribution of Bean plataspid invasion.

Conclusions/Significance

Niche models based on the two datasets showed little difference in native model predictions; however, when projecting onto the introduced area, models based on the environmental datasets showing low discrepancy among ranges recovered good model transferability in predicting the newly established population of Bean plataspid in the US. Recommendations were made for selecting biological meaningful environmental dimensions of low discrepancy among ranges to improve niche model transferability among these geographically separated areas. Outside of its native range, areas with invasion potential include the southeastern US in North America, southwestern Europe, southeastern South America, southern Africa, and the eastern coastal Australia.     

Effect of geographic background and equilibrium state on niche model transferability: predicting areas of invasion of Leptoglossus occidentalis

Niche conservatism providing support for using ecological niche modeling in biological invasions has been widely noticed, however, the equilibrium state and geographic background effect on niche model transferability has received scant attention. The western conifer seed bug, Leptoglossus occidentalis, native to western North America, has expanded its range eastward and has become an invasive pest in Europe and Asia. Niche models calibrated on the ranges of a small native population and two large expanding populations were compared. We found that the climate niche of L. occidentalis is conserved during its steady expansion in North America and rapid spread in Europe. Models based on the small western native range successfully captured the eastern expanding and introduced European populations, whereas the large area-based models varied with the presumed state of equilibrium. The equilibrium state based model succeeded but the non-equilibrium based model failed to predict the range in Europe. Our study estimates global invasion risk zones for L. occidentalis and suggests that, based on niche conservatism, modeling based on a reasonable geographic distribution at a climatic equilibrium of a species could guarantee the transferability of niche model prediction. Caution is warranted in interpreting low niche model transferability with niche differentiation and forwarding message for management strategy.     

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling–Daba Mountains,a major glacial refugium in Central China

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling–Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling–Daba Mountains of Central China.     

Modeling potential invasion of stored-product pest Cryptamorpha desjardinsii (Guérin-Méneville, 1844) (Coleoptera: Silvanidae) with emphasis on newly recorded areas

Desjardin’s flat bark beetle – Cryptamorpha desjardinsii (Guérin-Méneville, 1844) – is a widespread invasive species that can become a stored products pest and threaten native species. Our objectives were to provide the first detailed records of this species to South America and New Caledonia, and model suitable areas of C. desjardinsii globally based on different Ecological Niche Models (ENM). We combine several techniques to predict the ecological niche of C. desjardinsii using distribution records and a set of environmental predictor variables. Overall, the species appears to be well established in tropical and subtropical regions, particularly in islands, with new and previous records predicted by the consensus model. In South America, the driest and coldest regions were predicted to be less suitable for the occurrence of C. desjardinsii. In the current global spread context, these findings are pivotal to predicting areas of high risk of invasion of C. desjardinsii.     

Ecological niche modeling re‐examined: A case study with the Darwin's fox

Many previous studies have attempted to assess ecological niche modeling performance using receiver operating characteristic (ROC) approaches, even though diverse problems with this metric have been pointed out in the literature. We explored different evaluation metrics based on independent testing data using the Darwin's Fox (Lycalopex fulvipes) as a detailed case in point. Six ecological niche models (ENMs; generalized linear models, boosted regression trees, Maxent, GARP, multivariable kernel density estimation, and NicheA) were explored and tested using six evaluation metrics (partial ROC, Akaike information criterion, omission rate, cumulative binomial probability), including two novel metrics to quantify model extrapolation versus interpolation (E‐space index I) and extent of extrapolation versus Jaccard similarity (E‐space index II). Different ENMs showed diverse and mixed performance, depending on the evaluation metric used. Because ENMs performed differently according to the evaluation metric employed, model selection should be based on the data available, assumptions necessary, and the particular research question. The typical ROC AUC evaluation approach should be discontinued when only presence data are available, and evaluations in environmental dimensions should be adopted as part of the toolkit of ENM researchers. Our results suggest that selecting Maxent ENM based solely on previous reports of its performance is a questionable practice. Instead, model comparisons, including diverse algorithms and parameterizations, should be the sine qua non for every study using ecological niche modeling. ENM evaluations should be developed using metrics that assess desired model characteristics instead of single measurement of fit between model and data. The metrics proposed herein that assess model performance in environmental space (i.e., E‐space indices I and II) may complement current methods for ENM evaluation.