Losing your edge: climate change and the conservation value of range‐edge populations

Populations occurring at species' range edges can be locally adapted to unique environmental conditions. From a species' perspective, range‐edge environments generally have higher severity and frequency of extreme climatic events relative to the range core. Under future climates, extreme climatic events are predicted to become increasingly important in defining species' distributions. Therefore, range‐edge genotypes that are better adapted to extreme climates relative to core populations may be essential to species' persistence during periods of rapid climate change. We use relatively simple conceptual models to highlight the importance of locally adapted range‐edge populations (leading and trailing edges) for determining the ability of species to persist under future climates. Using trees as an example, we show how locally adapted populations at species' range edges may expand under future climate change and become more common relative to range‐core populations. We also highlight how large‐scale habitat destruction occurring in some geographic areas where many species range edge converge, such as biome boundaries and ecotones (e.g., the arc of deforestation along the rainforest‐cerrado ecotone in the southern Amazonia), can have major implications for global biodiversity. As climate changes, range‐edge populations will play key roles in helping species to maintain or expand their geographic distributions. The loss of these locally adapted range‐edge populations through anthropogenic disturbance is therefore hypothesized to reduce the ability of species to persist in the face of rapid future climate change.     

Intraspecific trait variation across elevation predicts a widespread tree species' climate niche and range limits

Global change is widely altering environmental conditions which makes accurately predicting species range limits across natural landscapes critical for conservation and management decisions. If climate pressures along elevation gradients influence the distribution of phenotypic and genetic variation of plant functional traits, then such trait variation may be informative of the selective mechanisms and adaptations that help define climatic niche limits. Using extensive field surveys along 16 elevation transects and a large common garden experiment, we tested whether functional trait variation could predict the climatic niche of a widespread tree species (Populus angustifolia) with a double quantile regression approach. We show that intraspecific variation in plant size, growth, and leaf morphology corresponds with the species' total climate range and certain climatic limits related to temperature and moisture extremes. Moreover, we find evidence of genetic clines and phenotypic plasticity at environmental boundaries, which we use to create geographic predictions of trait variation and maximum values due to climatic constraints across the western US. Overall, our findings show the utility of double quantile regressions for connecting species distributions and climate gradients through trait‐based mechanisms. We highlight how new approaches like ours that incorporate genetic variation in functional traits and their response to climate gradients will lead to a better understanding of plant distributions as well as identifying populations anticipated to be maladapted to future environments.     

Dispersal evolution in temporally variable environments: implications for plant range dynamics

Dispersal—the movement of an individual from the site of birth to a different site for reproduction—is an ecological and evolutionary driver of species ranges that shapes patterns of colonization, connectivity, gene flow, and adaptation. In plants, the traits that influence dispersal often vary within and among species, are heritable, and evolve in response to the fitness consequences of moving through heterogeneous landscapes. Spatial and temporal variation in the quality and quantity of habitat are important sources of selection on dispersal strategies across species ranges. While recent reviews have evaluated the interactions between spatial variation in habitat and dispersal dynamics, the extent to which geographic variation in temporal variability can also shape range-wide patterns in dispersal traits has not been synthesized. In this paper, we summarize key predictions from metapopulation models that evaluate how dispersal evolves in response to spatial and temporal habitat variability. Next, we compile empirical data that quantify temporal variability in plant demography and patterns of dispersal trait variation across species ranges to evaluate the hypothesis that higher temporal variability favors increased dispersal at plant range limits. We found some suggestive evidence supporting this hypothesis while more generally identifying a major gap in empirical work evaluating plant metapopulation dynamics across species ranges and geographic variation in dispersal traits. To address this gap, we propose several future research directions that would advance our understanding of the interplay between spatiotemporal variability and dispersal trait variation in shaping the dynamics of current and future species ranges.     

How are tree species distributed in climatic space? A simple and general pattern

Aim Although many factors undoubtedly affect species geographic distributions, can a single, simple model nonetheless capture most of the spatial variation in the probability of presence/absence in a large set of species? For 482 North American tree species that occur east of the Rocky Mountains, we investigated the shape(s) of the relationship between the probability of occupancy of a given location and macroclimate, and its consistency among species and regions. Location North America. Methods Using Little's tree range maps, we tested four hypothetical shapes of response relating occupancy to climate: (1) high occupancy of all suitable climates; (2) threshold response (i.e. unsuitable climates exclude species, but within the thresholds, species presence is independent of climate); (3) occupancy is a bivariate normal function of annual temperature and precipitation; and (4) asymmetric limitation (i.e. abiotic factors set abrupt range limits in stressful climates only). Finally, we compared observed climatic niches with the occupancy of similar climates on off‐shore islands as well as west of the Rockies. Results (a) Species' distributions in climatic space do not have strong thresholds, nor are they systematically skewed towards less stressful climates. (b) Occupancy can generally be described by a bivariate normal function of temperature and precipitation, with little or no interaction between the two variables. This model, averaged over all species, accounts for 82% of the spatial variation in the probability of occupancy of a given area. (c) Occupied geographic ranges are typically ringed by unoccupied, but climatically suitable areas. (d) Observed climatic niche positions are largely conserved between regions. Main conclusions We conclude that, despite the complexities of species histories and biologies, to a first approximation most of the variation in their geographic distributions relates to climate, in similar ways for nearly all species.     

Applying landscape metrics to species distribution model predictions to characterize internal range structure and associated changes

Distributional shifts in species ranges provide critical evidence of ecological responses to climate change. Assessments of climate-driven changes typically focus on broad-scale range shifts (e.g. poleward or upward), with ecological consequences at regional and local scales commonly overlooked. While these changes are informative for species presenting continuous geographic ranges, many species have discontinuous distributions—both natural (e.g. mountain or coastal species) or human-induced (e.g. species inhabiting fragmented landscapes)—where within-range changes can be significant. Here, we use an ecosystem engineer species (Sabellaria alveolata) with a naturally fragmented distribution as a case study to assess climate-driven changes in within-range occupancy across its entire global distribution. To this end, we applied landscape ecology metrics to outputs from species distribution modelling (SDM) in a novel unified framework. SDM predicted a 27.5% overall increase in the area of potentially suitable habitat under RCP 4.5 by 2050, which taken in isolation would have led to the classification of the species as a climate change winner. SDM further revealed that the latitudinal range is predicted to shrink because of decreased habitat suitability in the equatorward part of the range, not compensated by a poleward expansion. The use of landscape ecology metrics provided additional insights by identifying regions that are predicted to become increasingly fragmented in the future, potentially increasing extirpation risk by jeopardising metapopulation dynamics. This increased range fragmentation could have dramatic consequences for ecosystem structure and functioning. Importantly, the proposed framework—which brings together SDM and landscape metrics—can be widely used to study currently overlooked climate-driven changes in species internal range structure, without requiring detailed empirical knowledge of the modelled species. This approach represents an important advancement beyond predictive envelope approaches and could reveal itself as paramount for managers whose spatial scale of action usually ranges from local to regional.     

Both life‐history plasticity and local adaptation will shape range‐wide responses to climate warming in the tundra plant Silene acaulis

Many predictions of how climate change will impact biodiversity have focused on range shifts using species‐wide climate tolerances, an approach that ignores the demographic mechanisms that enable species to attain broad geographic distributions. But these mechanisms matter, as responses to climate change could fundamentally differ depending on the contributions of life‐history plasticity vs. local adaptation to species‐wide climate tolerances. In particular, if local adaptation to climate is strong, populations across a species’ range—not only those at the trailing range edge—could decline sharply with global climate change. Indeed, faster rates of climate change in many high latitude regions could combine with local adaptation to generate sharper declines well away from trailing edges. Combining 15 years of demographic data from field populations across North America with growth chamber warming experiments, we show that growth and survival in a widespread tundra plant show compensatory responses to warming throughout the species’ latitudinal range, buffering overall performance across a range of temperatures. However, populations also differ in their temperature responses, consistent with adaptation to local climate, especially growing season temperature. In particular, warming begins to negatively impact plant growth at cooler temperatures for plants from colder, northern populations than for those from warmer, southern populations, both in the field and in growth chambers. Furthermore, the individuals and maternal families with the fastest growth also have the lowest water use efficiency at all temperatures, suggesting that a trade‐off between growth and water use efficiency could further constrain responses to forecasted warming and drying. Taken together, these results suggest that populations throughout species’ ranges could be at risk of decline with continued climate change, and that the focus on trailing edge populations risks overlooking the largest potential impacts of climate change on species’ abundance and distribution.     

Integrating mechanistic and empirical model projections to assess climate impacts on tree species distributions in northwestern North America

Empirical and mechanistic models have both been used to assess the potential impacts of climate change on species distributions, and each modeling approach has its strengths and weaknesses. Here, we demonstrate an approach to projecting climate‐driven changes in species distributions that draws on both empirical and mechanistic models. We combined projections from a dynamic global vegetation model (DGVM) that simulates the distributions of biomes based on basic plant functional types with projections from empirical climatic niche models for six tree species in northwestern North America. These integrated model outputs incorporate important biological processes, such as competition, physiological responses of plants to changes in atmospheric CO₂ concentrations, and fire, as well as what are likely to be species‐specific climatic constraints. We compared the integrated projections to projections from the empirical climatic niche models alone. Overall, our integrated model outputs projected a greater climate‐driven loss of potentially suitable environmental space than did the empirical climatic niche model outputs alone for the majority of modeled species. Our results also show that refining species distributions with DGVM outputs had large effects on the geographic locations of suitable habitat. We demonstrate one approach to integrating the outputs of mechanistic and empirical niche models to produce bioclimatic projections. But perhaps more importantly, our study reveals the potential for empirical climatic niche models to over‐predict suitable environmental space under future climatic conditions.     

Life on the edge: carnivore body size variation is all over the place

Evolutionary biologists have long been fascinated by both the ways in which species respond to ecological conditions at the edges of their geographic ranges and the way that species'' body sizes evolve across their ranges. Surprisingly, though, the relationship between these two phenomena is rarely studied. Here, we examine whether carnivore body size changes from the interior of their geographic range towards the range edges. We find that within species, body size often varies strongly with distance from the range edge. However, there is no general tendency across species for size to be either larger or smaller towards the edge. There is some evidence that the smallest guild members increase in size towards their range edges, but results for the largest guild members are equivocal. Whether individuals vary in relation to the distance from the range edges often depends on the way edge and interior are defined. Neither geographic range size nor absolute body size influences the tendency of size to vary with distance from the range edge. Therefore, we suggest that the frequent significant association between body size and the position of individuals along the edge-core continuum reflects the prevalence of geographic size variation and that the distance to range edge per se does not influence size evolution in a consistent way.     

Topographic,latitudinal and climatic distribution of Pinus coulteri: geographic range limits are not at the edge of the climate envelope

With changing climate, many species are projected to move poleward or to higher elevations to track suitable climates. The prediction that species will move poleward assumes that geographically marginal populations are at the edge of the species' climatic range. We studied Pinus coulteri from the center to the northern (poleward) edge of its range, and examined three scenarios regarding the relationship between the geographic and climatic margins of a species' range. We used herbarium and iNaturalist.org records to identify P. coulteri sites, generated a species distribution model based on temperature, precipitation, climatic water deficit, and actual evapotranspiration, and projected suitability under future climate scenarios. In fourteen populations from the central to northern portions of the range, we conducted field studies and recorded elevation, slope and aspect (to estimate solar insolation) to examine relationships between local and regional distributions. We found that northern populations of P. coulteri do not occupy the cold or wet edge of the species' climatic range; mid‐latitude, high elevation populations occupy the cold margin. Aspect and insolation of P. coulteri populations changed significantly across latitudes and elevations. Unexpectedly, northern, low‐elevation stands occupy north‐facing aspects and receive low insolation, while central, high‐elevation stands grow on more south‐facing aspects that receive higher insolation. Modeled future climate suitability is projected to be highest in the central, high elevation portion of the species range, and in low‐lying coastal regions under some scenarios, with declining suitability in northern areas under most future scenarios. For P. coulteri, the lack of high elevation habitat combined with a major dispersal barrier may limit northward movement in response to a warming climate. Our analyses demonstrate the importance of distinguishing geographically vs. climatically marginal populations, and the importance of quantitative analysis of the realized climate space to understand species range limits.     

Reduced variability in range‐edge butterfly populations over three decades of climate warming

Populations at the high latitude edge of species’ geographical ranges are thought to show larger interannual population fluctuations, with subsequent higher local extinction risk, than those within the ‘core’ climatic range. As climate envelopes shift northward under climate warming, however, we would expect populations to show dampened variability. We test this hypothesis using annual abundance indices from 19 butterfly species across 79 British monitoring sites between 1976 and 2009, a period of climatic warming. We found that populations in the latter (warmer) half of the recording period show reduced interannual population variability. Species with more southerly European distributions showed the greatest dampening in population variability over time. Our results suggest that increases in population variability occur towards climatic range boundaries. British sites, previously existing at the margins of suitable climate space, now appear to fall closer to the core climatic range for many butterfly species.     

A test of the central–marginal hypothesis using population genetics and ecological niche modelling in an endemic salamander (Ambystoma barbouri)

The central–marginal hypothesis (CMH) predicts that population size, genetic diversity and genetic connectivity are highest at the core and decrease near the edges of species' geographic distributions. We provide a test of the CMH using three replicated core‐to‐edge transects that encompass nearly the entire geographic range of the endemic streamside salamander (Ambystoma barbouri). We confirmed that the mapped core of the distribution was the most suitable habitat using ecological niche modelling (ENM) and via genetic estimates of effective population sizes. As predicted by the CMH, we found statistical support for decreased genetic diversity, effective population size and genetic connectivity from core to edge in western and northern transects, yet not along a southern transect. Based on our niche model, habitat suitability is lower towards the southern range edge, presumably leading to conflicting core‐to‐edge genetic patterns. These results suggest that multiple processes may influence a species' distribution based on the heterogeneity of habitat across a species' range and that replicated sampling may be needed to accurately test the CMH. Our work also emphasizes the importance of identifying the geographic range core with methods other than using the Euclidean centre on a map, which may help to explain discrepancies among other empirical tests of the CMH. Assessing core‐to‐edge population genetic patterns across an entire species' range accompanied with ENM can inform our general understanding of the mechanisms leading to species' geographic range limits.     

Cold range edges of marine fishes track climate change better than warm edges

Species around the world are shifting their ranges in response to climate change. To make robust predictions about climate‐related colonizations and extinctions, it is vital to understand the dynamics of range edges. This study is among the first to examine annual dynamics of cold and warm range edges, as most global change studies average observational data over space or over time. We analyzed annual range edge dynamics of marine fishes—both at the individual species level and pooled into cold‐ and warm‐edge assemblages—in a multi‐decade time‐series of trawl surveys conducted on the Northeast US Shelf during a period of rapid warming. We tested whether cold edges show stronger evidence of climate tracking than warm edges (due to non‐climate processes or time lags at the warm edge; the biogeography hypothesis or extinction debt hypothesis), or whether they tracked temperature change equally (due to the influence of habitat suitability; the ecophysiology hypothesis). In addition to exploring correlations with regional temperature change, we calculated species‐ and assemblage‐specific sea bottom and sea surface temperature isotherms and used them to predict range edge position. Cold edges shifted further and tracked sea surface and bottom temperature isotherms to a greater degree than warm edges. Mixed‐effects models revealed that for a one‐degree latitude shift in isotherm position, cold edges shifted 0.47 degrees of latitude, and warm edges shifted only 0.28 degrees. Our results suggest that cold range edges are tracking climate change better than warm range edges, invalidating the ecophysiology hypothesis. We also found that even among highly mobile marine ectotherms in a global warming hotspot, few species are fully keeping pace with climate.     

A resurrection study reveals limited evolution of thermal performance in response to recent climate change across the geographic range of the scarlet monkeyflower

Evolutionary rescue can prevent populations from declining under climate change, and should be more likely at high-latitude, “leading” edges of species’ ranges due to greater temperature anomalies and gene flow from warm-adapted populations. Using a resurrection study with seeds collected before and after a 7-year period of record warming, we tested for thermal adaptation in the scarlet monkeyflower Mimulus cardinalis. We grew ancestors and descendants from northern-edge, central, and southern-edge populations across eight temperatures. Despite recent climate anomalies, populations showed limited evolution of thermal performance curves. However, one southern population evolved a narrower thermal performance breadth by 1.31°C, which matches the direction and magnitude of the average decrease in seasonality experienced. Consistent with the climate variability hypothesis, thermal performance breadth increased with temperature seasonality across the species’ geographic range. Inconsistent with performance trade-offs between low and high temperatures across populations, we did not detect a positive relationship between thermal optimum and mean temperature. These findings fail to support the hypothesis that evolutionary response to climate change is greatest at the leading edge, and suggest that the evolution of thermal performance is unlikely to rescue most populations from the detrimental effects of rapidly changing climate.     

Mutualist‐mediated effects on species' range limits across large geographic scales

Understanding the processes determining species range limits is central to predicting species distributions under climate change. Projected future ranges are extrapolated from distribution models based on climate layers, and few models incorporate the effects of biotic interactions on species' distributions. Here, we show that a positive species interaction ameliorates abiotic stress, and has a profound effect on a species' range limits. Combining field surveys of 92 populations, 10 common garden experiments throughout the range, species distribution models and greenhouse experiments, we show that mutualistic fungal endophytes ameliorate drought stress and broaden the geographic range of their native grass host Bromus laevipes by thousands of square kilometres (~ 20% larger) into drier habitats. Range differentiation between fungal‐associated and fungal‐free grasses was comparable to species‐level range divergence of congeners, indicating large impacts on range limits. Positive biotic interactions may be underappreciated in determining species' ranges and species' responses to future climates across large geographic scales.     

Climate structures genetic variation across a species' elevation range: a test of range limits hypotheses

Gene flow may influence the formation of species range limits, and yet little is known about the patterns of gene flow with respect to environmental gradients or proximity to range limits. With rapid environmental change, it is especially important to understand patterns of gene flow to inform conservation efforts. Here we investigate the species range of the selfing, annual plant, Mimulus laciniatus, in the California Sierra Nevada. We assessed genetic variation, gene flow, and population abundance across the entire elevation‐based climate range. Contrary to expectations, within‐population plant density increased towards both climate limits. Mean genetic diversity of edge populations was equivalent to central populations; however, all edge populations exhibited less genetic diversity than neighbouring interior populations. Genetic differentiation was fairly consistent and moderate among all populations, and no directional signals of contemporary gene flow were detected between central and peripheral elevations. Elevation‐driven gene flow (isolation by environment), but not isolation by distance, was found across the species range. These findings were the same towards high‐ and low‐elevation range limits and were inconsistent with two common centre‐edge hypotheses invoked for the formation of species range limits: (i) decreasing habitat quality and population size; (ii) swamping gene flow from large, central populations. This pattern demonstrates that climate, but not centre‐edge dynamics, is an important range‐wide factor structuring M. laciniatus populations. To our knowledge, this is the first empirical study to relate environmental patterns of gene flow to range limits hypotheses. Similar investigations across a wide variety of taxa and life histories are needed.     

Habitat associations of thermophilous butterflies are reduced despite climatic warming

Climate warming threatens the survival of species at their warm, trailing‐edge range boundaries but also provides opportunities for the ecological release of populations at the cool, leading edges of their distributions. Thus, as the climate warms, leading‐edge populations are expected to utilize an increased range of habitat types, leading to larger population sizes and range expansion. Here, we test the hypothesis that the habitat associations of British butterflies have expanded over three decades of climate warming. We characterize the habitat breadth of 27 southerly distributed species from 77 monitoring transects between 1977 and 2007 by considering changes in densities of butterflies across 11 habitat types. Contrary to expectation, we find that 20 of 27 (74%) butterfly species showed long‐term contractions in their habitat associations, despite some short‐term expansions in habitat breadth in warmer‐than‐usual years. Thus, we conclude that climatic warming has ameliorated habitat contractions caused by other environmental drivers to some extent, but that habitat degradation continues to be a major driver of reductions in habitat breadth and population density of butterflies.     

Is the spatial distribution of European beech (Fagus sylvatica L.) limited by its potential height growth?

Aim To improve our understanding of species range limits by studying how height growth, a trait related to plant survival, varies throughout the geographic range of Fagus sylvatica L. in France. Location The geographic range of beech in France, representing the western area of its European distribution, within which this species exhibits range distribution limits in both plains and mountainous areas. Methods A generalized linear regression model was used to link beech growth performance to environmental variables using data from 819 plots of the French National Forest Inventory (IFN) database. This model was applied to predict potential growth on 97,281 IFN plots covering the geographic range of beech in France. A kriging technique was used to interpolate estimated growth potential. Finally, the performance of plot‐based predictions of potential growth from the map (i.e. map quality) was evaluated against an independent data set. Results The beech growth performance model highlighted the major impact of climate on potential tree growth at a broad spatial scale. The relevant climatic factors were related mainly to spring cold, summer heat, and winter temperatures and rainfall. The study also revealed the predictive power of soil parameters, which explained a large proportion of the variation in potential beech growth (c. 30%). Analyses of height growth patterns near the boundary of the species range in France showed that the limit only partly coincides with the growth decline caused by climatic and soil factors. Along parts of the range limit, the predicted potential for growth was high, suggesting that in these areas the limit of the range could be explained by other factors, such as competition or constraints on reproduction. Main conclusions The spatial variation in the potential height growth of Fagus sylvatica can be explained by environmental factors and is partly correlated with its regional range limits. By identifying areas where growth potential constrains the geographic range of species, environmental growth models can help to improve our knowledge of the spatial drivers of species geographic range limits and shed light on their response to future environmental changes.     

Cetacean range and climate in the eastern North Atlantic: future predictions and implications for conservation

There is increasing evidence that the distributions of a large number of species are shifting with global climate change as they track changing surface temperatures that define their thermal niche. Modelling efforts to predict species distributions under future climates have increased with concern about the overall impact of these distribution shifts on species ecology, and especially where barriers to dispersal exist. Here we apply a bio‐climatic envelope modelling technique to investigate the impacts of climate change on the geographic range of ten cetacean species in the eastern North Atlantic and to assess how such modelling can be used to inform conservation and management. The modelling process integrates elements of a species' habitat and thermal niche, and employs “hindcasting” of historical distribution changes in order to verify the accuracy of the modelled relationship between temperature and species range. If this ability is not verified, there is a risk that inappropriate or inaccurate models will be used to make future predictions of species distributions. Of the ten species investigated, we found that while the models for nine could successfully explain current spatial distribution, only four had a good ability to predict distribution changes over time in response to changes in water temperature. Applied to future climate scenarios, the four species‐specific models with good predictive abilities indicated range expansion in one species and range contraction in three others, including the potential loss of up to 80% of suitable white‐beaked dolphin habitat. Model predictions allow identification of affected areas and the likely time‐scales over which impacts will occur. Thus, this work provides important information on both our ability to predict how individual species will respond to future climate change and the applicability of predictive distribution models as a tool to help construct viable conservation and management strategies.     

A changing climate is eroding the geographical range of the Namib Desert tree Aloe through population declines and dispersal lags

While poleward species migration in response to recent climatic warming is widely documented, few studies have examined entire range responses of broadly distributed sessile organisms, including changes on both the trailing (equatorward) and the leading (poleward) range edges. From a detailed population census throughout the entire geographical range of Aloe dichotoma Masson, a long-lived Namib Desert tree, together with data from repeat photographs, we present strong evidence that a developing range shift in this species is a 'fingerprint' of anthropogenic climate change. This is explained at a high level of statistical significance by population level impacts of observed regional warming and resulting water balance constraints. Generalized linear models suggest that greater mortalities and population declines in equatorward populations are virtually certainly the result, due to anthropogenic climate change, of the progressive exceedance of critical climate thresholds that are relatively closer to the species' tolerance limits in equatorward sites. Equatorward population declines are also broadly consistent with bioclimatically modelled projections under anticipated anthropogenic climate change but, as yet, there is no evidence of poleward range expansion into the area predicted to become suitable in future, despite good evidence for positive population growth trends in poleward populations. This study is among the first to show a marked lag between trailing edge population extinction and leading edge range expansion in a species experiencing anthropogenic climate change impacts, a pattern likely to apply to most sessile and poorly dispersed organisms. This provides support for conservative assumptions of species' migration rates when modelling climate change impacts for such species. Aloe dichotoma 's response to climate change suggests that desert ecosystems may be more sensitive to climate change than previously suspected.     

Evaluating the potential causes of range limits of birds of the Colombian Andes

Aim To evaluate how factors acting at different spatial scales influence range limits in bird species of the Colombian Andes. Location Andes Mountains of Colombia. Methods We used Maxent , a climate envelope model (CEM), and environmental and geographic information to study range‐filling (i.e. the extent to which a species occurs in all the areas in which it is predicted to occur) in 70 range‐restricted bird species of the Colombian Andes. Environmental data were taken from the WorldClim database, and species occurrence data were taken from museum data collated by the BioMap project, an observational database, and the literature. We evaluated how climate and geographic barriers may shape range limits at two scales. Results At a broad extent (i.e. across the three main cordilleras within the Colombian Andes), we find that CEMs predict there to be suitable environmental conditions for particular species in regions where the species is absent, possibly as a result of dispersal limitation or biotic interactions. In contrast, at a finer scale (within a given cordillera), species generally occur across the entire area predicted to be suitable by a given CEM. Geographic discontinuities within cordilleras do not generally correspond to range limits; instead, range limits correspond to changes in environmental conditions. Main conclusions Our results suggest that different mechanisms influence the presence of species at different scales. Dispersal limitation, potentially combined with species interactions, may influence range limits at a broad extent (the entire Colombian Andes), while strong environmental gradients correspond to range limits at a finer scale (within a cordillera).