Analyses of least cost paths for determining effects of habitat types on landscape permeability: wolves in Poland

Determining ecological corridors is crucial for conservation efforts in fragmented habitats. Commonly employed least cost path (LCP) analysis relies on the underlying cost matrix. By using Ecological Niche Factor Analysis, we minimized the problems connected with subjective cost assessment or the use of presence/absence data. We used data on the wolf presence/absence in Poland to identify LCPs connecting patches of suitable wolf habitat, factors that influence patch occupancy, and compare LCPs between different genetic subpopulations. We found that a lower proportion of cities and roads surrounds the most densely populated patches. Least cost paths between areas where little dispersal takes place (i.e., leading to unpopulated patches or between different genetic subpopulations) ran through a higher proportion of roads and human settlements. They also crossed larger maximal distances over deforested areas. We propose that, apart from supplying the basis for direct conservation efforts, LCPs can be used to determine what factors might facilitate or hinder dispersal by comparing different subsets of LCPs. The methods employed can be widely applicable to gain more in-depth information on potential dispersal barriers for large carnivores.     

Microhabitat selection by Eurasian lynx and its implications for species conservation

We studied microhabitat selection of the Eurasian lynxLynx lynx (Linnaeus, 1758) at 116 hunting and 88 resting sites in Białowieża Primeval Forest (Poland) to describe its characteristics and determine the importance of habitat structure for stalking prey and for security during resting. We identified lynx-used sites by radio-tracking 3 male and 3 female lynx. When hunting, the lynx did not select for any type or age class of forest. During both summer and winter, the lynx selected sites characterized by high complexity (number of structures useful for stalking: fallen logs and branches, root plates, patches of dense bushes) and low visibility. In summer, hunting sites were often located in the vicinity of small forest glades that provided good stalking opportunities for lynx and rich foraging resources for roe deer — the main prey of lynx. The habitat at kill sites was more open than at sites where the prey was cached, with higher visibility, lower density of trees and poorer undergrowth. The most important characteristic of resting sites was very low visibility that resulted mainly from using young pine or spruce thickets in the winter and dense undergrowth of oak-lime-hornbeam and ash-alder forests in the summer. The information provided by this study could have direct implications for Eurasian lynx conservation by guiding forest restructuring to better suit the species’ biological requirements.     

Habitat selection by wolvesCanis lupus in the uplands and mountains of southern Poland

Using data from the National Wolf Census, carried out in Poland in 2000–2001, and GIS techniques we analysed habitat selection by wolvesCanis lupus Linnaeus, 1758 in uplands and mountains of southern Poland. We compared ten habitat variables and two parameters related to wolf abundance in 52 circular plots (154 km² each) with recorded wolves and 97 randomly selected plots with no signs of wolf presence. Wolf plots were characterized by higher elevation and closer location to the state border than wolf-free plots. Furthermore, wolf plots had higher forest cover, but smaller number of villages and towns and shorter railways and roads than plots without wolves. The best model explaining wolf distribution included forest cover, number of villages, length of roads and railway lines. Compared to northern Poland, the southern part of the country offers worse habitats for wolves due to significantly denser network of settlements and transportation routes.     

Genetic structure of the Eurasian lynx population in north-eastern Poland and the Baltic states

We analyzed the genotypes of Eurasian lynx (Lynx lynx) from three populations in the westernmost part of the species main range. One population was situated at the distribution edge (NE Poland) and the two other (Latvia and Estonia) were located within the main, contiguous range of the species. The aim was to determine if the genetic composition varied among these populations and if there was evidence of isolation among them. Based on microsatellite allele frequencies, we found the allelic richness in Polish lynx to be lower than that in lynx from Latvia and Estonia. We also found significant differentiation among the lynx populations, with the NE Poland population forming a distinct genetic group relative to the two other populations (R _ST = 0.15 and 0.22, P < 0.0001). We suggest that genetic differentiation among lynx populations is the result of habitat insularisation that limits gene flow. This finding emphasizes the necessity to consider the lynx genetic differentiation in conservation planning of this species in Poland.     

Natural,not urban,barriers define population structure for a coastal endemic butterfly

Habitat loss and fragmentation are the major causes of biodiversity loss, and, increasingly, habitat is fragmented by urbanization. Yet, the degree to which urbanization creates barriers to animal dispersal remains poorly understood. We used population genetic techniques to determine whether urbanization and/or natural landscape features are dispersal barriers to a butterfly, Atrytonopsis new species 1, throughout its range on coastal sand dunes that are increasingly threatened by development. Using AFLP markers that produced 89 polymorphic loci, we found significant population structure across the range of Atrytonopsis sp1. We found no indication that existing levels of urbanization were barriers to Atrytonopsis sp1 dispersal. Rather, two natural barriers, an ocean inlet and maritime forest, explained the genetic structure. Even in areas with long histories of urbanization, we found no significant isolation-by-distance relationship, and there was very low genetic differentiation between sampling locations. Consequently, conservation strategies for Atrytonopsis sp1, and potentially for other mobile insects that use open-structured habitats, should not focus explicitly on habitat corridors through urban areas, but rather should seek to preserve and restore as much habitat as possible across the butterfly’s range.     

How do low dispersal species establish large range sizes? The case of the water beetle Graphoderus bilineatus

Species’ dispersal abilities have been considered a major driving force in establishment and maintenance of large range sizes. However, recent studies question the general validity of this relationship because the relationship between dispersal ability and range size might in some cases be less important than species phylogeny or local spatial attributes. In this study we used the water beetle Graphoderus bilineatus a philopatric species of conservation concern in Europe as a model to explain large range size and to support effective conservation measures for such species that also have limited dispersal. We recorded the presence/absence of G. bilineatus and measured 14 habitat and 20 landscape variables at 228 localities in Estonia, Poland and Sweden within the core range of the species. Using information theory and average multivariate logistic regression models we determined that presence of G. bilineatus depended on landscape connectivity, distance to a possible source habitat, and stability of the site; however, specificity of habitat characteristics was not vital for the species. We reason that the large range of G. bilineatus is best explained by the historical combination of lakes, river systems and wetlands which used to be highly connected throughout the central plains of Europe. Our data suggest that a broad habitat niche can prevent landscape elements from becoming barriers for species like G. bilineatus. Therefore, we question the usefulness of site protection as conservation measures for G. bilineatus and similar philopatric species. Instead, conservation actions should be focused at the landscape level to ensure a long‐term viability of such species across their range.     

Anthropogenic Habitats Facilitate Dispersal of an Early Successional Obligate: Implications for Restoration of an Endangered Ecosystem

Landscape modification and habitat fragmentation disrupt the connectivity of natural landscapes, with major consequences for biodiversity. Species that require patchily distributed habitats, such as those that specialize on early successional ecosystems, must disperse through a landscape matrix with unsuitable habitat types. We evaluated landscape effects on dispersal of an early successional obligate, the New England cottontail (Sylvilagus transitionalis). Using a landscape genetics approach, we identified barriers and facilitators of gene flow and connectivity corridors for a population of cottontails in the northeastern United States. We modeled dispersal in relation to landscape structure and composition and tested hypotheses about the influence of habitat fragmentation on gene flow. Anthropogenic and natural shrubland habitats facilitated gene flow, while the remainder of the matrix, particularly development and forest, impeded gene flow. The relative influence of matrix habitats differed between study areas in relation to a fragmentation gradient. Barrier features had higher explanatory power in the more fragmented site, while facilitating features were important in the less fragmented site. Landscape models that included a simultaneous barrier and facilitating effect of roads had higher explanatory power than models that considered either effect separately, supporting the hypothesis that roads act as both barriers and facilitators at all spatial scales. The inclusion of LiDAR-identified shrubland habitat improved the fit of our facilitator models. Corridor analyses using circuit and least cost path approaches revealed the importance of anthropogenic, linear features for restoring connectivity between the study areas. In fragmented landscapes, human-modified habitats may enhance functional connectivity by providing suitable dispersal conduits for early successional specialists.     

Landscape genetics of the Alpine newt (<Emphasis Type="Italic">Mesotriton alpestris</Emphasis>) inferred from a strip-based approach

Habitat destruction and fragmentation are known to strongly affect dispersal by altering the quality of the environment between populations. As a consequence, lower landscape connectivity is expected to enhance extinction risks through a decrease in gene flow and the resulting negative effects of genetic drift, accumulation of deleterious mutations and inbreeding depression. Such phenomena are particularly harmful for amphibian species, characterized by disjunct breeding habitats. The dispersal behaviour of amphibians being poorly understood, it is crucial to develop new tools, allowing us to determine the influence of landscape connectivity on the persistence of populations. In this study, we developed a new landscape genetics approach that aims at identifying land-uses affecting genetic differentiation, without a priori assumptions about associated ecological costs. We surveyed genetic variation at seven microsatellite loci for 19 Alpine newt (Mesotriton alpestris) populations in western Switzerland. Using strips of varying widths that define a dispersal corridor between pairs of populations, we were able to identify land-uses that act as dispersal barriers (i.e. urban areas) and corridors (i.e. forests). Our results suggest that habitat destruction and landscape fragmentation might in the near future affect common species such as M. alpestris. In addition, by identifying relevant landscape variables influencing population structure without unrealistic assumptions about dispersal, our method offers a simple and flexible tool of investigation as an alternative to least-cost models and other approaches.     

Evaluating Landscape Options for Corridor Restoration between Giant Panda Reserves

The establishment of corridors can offset the negative effects of habitat fragmentation by connecting isolated habitat patches. However, the practical value of corridor planning is minimal if corridor identification is not based on reliable quantitative information about species-environment relationships. An example of this need for quantitative information is planning for giant panda conservation. Although the species has been the focus of intense conservation efforts for decades, most corridor projects remain hypothetical due to the lack of reliable quantitative researches at an appropriate spatial scale. In this paper, we evaluated a framework for giant panda forest corridor planning. We linked our field survey data with satellite imagery, and conducted species occupancy modelling to examine the habitat use of giant panda within the potential corridor area. We then conducted least-cost and circuit models to identify potential paths of dispersal across the landscape, and compared the predicted cost under current conditions and alternative conservation management options considered during corridor planning. We found that due to giant panda''s association with areas of low elevation and flat terrain, human infrastructures in the same area have resulted in corridor fragmentation. We then identified areas with high potential to function as movement corridors, and our analysis of alternative conservation scenarios showed that both forest/bamboo restoration and automobile tunnel construction would significantly improve the effectiveness of corridor, while residence relocation would not significantly improve corridor effectiveness in comparison with the current condition. The framework has general value in any conservation activities that anticipate improving habitat connectivity in human modified landscapes. Specifically, our study suggested that, in this landscape, automobile tunnels are the best means to remove current barriers to giant panda movements caused by anthropogenic interferences.     

Identifying refugia and corridors under climate change conditions for the Sichuan snub‐nosed monkey (Rhinopithecus roxellana) in Hubei Province,China

Using a case study of an isolated management unit of Sichuan snub‐nosed monkey (Rhinopithecus roxellana), we assess the extent that climate change will impact the species’ habitat distribution in the current period and projected into the 2050s. We identify refugia that could maintain the population under climate change and determine dispersal paths for movement of the population to future suitable habitats. Hubei Province, China. We identified climate refugia and potential movements by integrating bioclimatic models with circuit theory and least‐cost model for the current period (1960–1990) and the 2050s (2041–2060). We coupled a maximum entropy algorithm to predict suitable habitat for the current and projected future periods. Suitable habitat areas that were identified during both time periods and that also satisfied home range and dispersal distance conditions were delineated as refugia. We mapped potential movements measured as current flow and linked current and future habitats using least‐cost corridors. Our results indicate up to 1,119 km² of currently suitable habitat within the study range. Based on our projections, a habitat loss of 67.2% due to climate change may occur by the 2050s, resulting in a reduced suitable habitat area of 406 km² and very little new habitat. The refugia areas amounted to 286 km² and were located in Shennongjia National Park and Badong Natural Reserve. Several connecting corridors between the current and future habitats, which are important for potential movements, were identified. Our assessment of the species predicted a trajectory of habitat loss following anticipated future climate change. We believe conservation efforts should focus on refugia and corridors when planning for future species management. This study will assist conservationists in determining high‐priority regions for effective maintenance of the endangered population under climate change and will encourage increased habitat connectivity.     

Physical barriers and corridors in urban habitats affect colonisation and parasitism rates of a specialist leaf miner

1. Urban environments are fragmented habitats characterised by the presence of physical barriers, which may negatively affect dispersal and colonisation by insect herbivores and their natural enemies. Conversely, plants growing along pavements may function as dispersal corridors, helping to moderate the harmful effects of resource patch isolation on organism movement and population persistence. 2. We experimentally tested the effects of walls as physical barriers to the dispersal of the leaf miner Liriomyza commelinae Frost and colonisation of its host plant, Commelina erecta L., in urban habitats. We also evaluated whether plants along pavements could act as corridors for this species. 3. We exposed experimental host plants to the leaf miner in houses with front gardens and back yards, the latter being completely surrounded by walls. The front gardens had walls but none separating them from the pavement. Previously mined plants were also exposed to parasitoids in the yards to determine parasitoid attack. 4. Liriomyza commelinae took longer to colonise back yards with higher walls, and the abundance of mined plants along pavements reduced the colonisation time. Leaf‐miner abundance was marginally affected by the yard type, and was lower in back yards. Cumulative parasitism rates decreased with increasing distance at which mined plants were placed from pavements. 5. Constructions act as physical barriers, having a negative impact on colonisation of host plants by leaf miners. The function of pavements as corridors seems to depend on the abundance of mined plants. Parasitism may be affected by distance from the corridor rather than physical barriers or other potential hosts.     

Dispersal through stream networks: modelling climate‐driven range expansions of fishes

Aim To incorporate dispersal through stream networks into models predicting the future distribution of a native, freshwater fish given climate change scenarios. Location Sweden. Methods We used logistic regression to fit climate and habitat data to observed pike (Esox lucius Linnaeus) distributions in 13,476 lakes. We used GIS to map dispersal pathways through streams. Lakes either (1) contained pike or were downstream from pike lakes, (2) were upstream from pike lakes, but downstream from natural dispersal barriers, or (3) were isolated from streams or were upstream from natural dispersal barriers. We then used climate projections to model future distributions of pike and compared our results with and without including dispersal. Results Given climate and habitat, pike were predicted present in all of 99,249 Swedish lakes by 2100. After accounting for dispersal barriers, we only predicted pike presence in 31,538 lakes. Dispersal barriers most strongly limited pike invasion in mountainous regions, but low connectivity also characterized some relatively flat regions. Main conclusions The dendritic network structure of streams and interconnected lakes makes a two‐dimensional representation of the landscape unsuitable for predicting range shifts of many freshwater organisms. If dispersal through stream networks is not accounted for, predictions of future fish distributions in a warmer climate might grossly overestimate range expansions of warm and cool‐water fishes and underestimate range contractions of cold‐water fishes. Dispersal through stream networks can be modelled in any region for which a digital elevation model and species occurrence data are available.     

Behavioural and spatial adaptation of the Eurasian lynx to a decline in prey availability

The distribution and abundance of food resources is a major factor influencing animal populations. I studied the effect of a roe and red deer population decline on diet composition, home range size and foraging pattern in the Eurasian lynxLynx lynx (Linnaeus, 1758) in the Bia?owie?a Primeval Forest (BPF), eastern Poland. The population of cervids in BPF experienced a nearly two-fold reduction in size from 1991 through 2006 due to severe hunting pressure between 1991 and 1996. Comparison of published data on lynx diet during the high abundance of ungulates with new data obtained for the low abundance period showed that despite a significant decline in their availability, cervids (roe and red deer) continued to form the majority of the diet of lynx, with roe deer being most preferred in both periods. Home range sizes of lynx showed a tendency to increase with declining prey densities, as indicated by relative percentage increases in average yearly home range sizes amongst different sex/age groups. In response to lower availability of their main prey, lynx increased their daily straight-line movement distances by 44% and doubled the ranges covered in 5-day periods. This illustrated that, with declining prey abundance, the lynx increased their hunting efforts by either spending more time actively searching for prey or continuing foraging even after a successful hunt. Spatial analysis of the distribution of ungulates and lynx indicated that deer were evenly distributed throughout lynx ranges in BPF and spatial proximity of the predator to prey sites did not play an important role in the efficiency of hunting. Lynx may adapt to changing prey availability by increasing search effort, but this was not sufficient to prevent the negative influences of the prey decline on the lynx population. Prey depletion has an immediate effect on lynx spatial organization and, in consequence, on their density. This information has to be considered in prioritizing lynx conservation measures and management of ungulates.     

Habitat persistence,habitat availability and the evolution of dispersal

Many organisms live in ephemeral habitats, making dispersal a vital element of life history. Here, we investigate how dispersal rate evolves in response to habitat persistence, mean habitat availability and landscape pattern. We show that dispersal rate is generally lowered by reduced habitat availability and by longer habitat persistence. However, for habitats that persist for an average of ten times the length of a generation, we show a clear non-monotonic relationship between habitat availability and dispersal rate. Some patterns of available habitat result in populations with dispersal polymorphisms. We explain these observations as a metapopulation effect, with the rate of evolution a function of both within-population and between-population selection pressures. Individuals in corridors evolve much lower dispersal rates than those in the mainland populations, especially within long, narrow corridors. We consider the implications of the results for conservation.     

Habitat corridors function as both drift fences and movement conduits for dispersing flies

Corridors connect otherwise isolated habitat patches and can direct movement of animals among such patches. In eight experimental landscapes, we tested two hypotheses of how corridors might affect dispersal behavior. The Traditional Corridor hypothesis posits that animals preferentially leave patches via corridors, following them into adjacent patches. The Drift Fence hypothesis posits that animals dispersing through matrix habitat are diverted into patches with corridors because they follow corridors when encountered. House flies (Musca domestica L.), a species that prefers the habitat of our patches and corridors, were released in a central patch (100×100 m) and recaptured in peripheral patches that were or were not connected by a corridor. Flies were captured more frequently in connected than unconnected patches, thereby supporting the Traditional Corridor hypothesis. The Drift Fence hypothesis was also supported, as flies were captured more frequently in unconnected patches with blind (dead end) corridors than in unconnected patches of equal area without blind corridors. A second experiment tested whether these results might be dependent on the type of patch-matrix boundary encountered by dispersing flies and whether edge-following behavior might be the mechanism underlying the observed corridor effect in the first experiment. We recorded dispersal patterns of flies released along forest edges with dense undergrowth in the forest (“closed” edges) and along edges with little forest understory (“open” edges). Flies were less likely to cross and more likely to follow closed edges than open edges, indicating that when patch and corridor edges are pronounced, edge-following behavior of flies may direct them along corridors into connected patches. Because edges in the first experiment were open, these results also suggest that corridor effects for flies in that experiment would have been even stronger if the edges around the source patches and corridors had been more closed. Taken together, our results suggest that corridors can affect dispersal of organisms in unappreciated ways (i.e., as drift fences) and that edge type can alter dispersal behavior.     

基于野外观测检验动物廊道选址模型的准确性——以白头叶猴为例

廊道构建是减少栖息地破碎化负面影响的重要策略之一。目前,已经有许多模型用于动物廊道的选址,而"选址模型是否能准确预测动物迁移的实际发生位置"一直是保护生物学最为关注的问题。最小成本路径模型（LCP）和条件最小成本廊道模型（CMTC）是两种较为常用的廊道选址模型。以白头叶猴（Trachypithecus leucocephalus）为目标物种,分别运用LCP和CMTC模拟生成白头叶猴迁移廊道,将模拟结果与野外观测廊道进行对比,检验两种方法的准确性。结果表明：与野外观测实际廊道相比,LCP模型模拟结果的完全准确率为46.7%,部分准确率为20%,完全不准确率为33.3%;CMTC模型模拟结果的完全准确率为26.7%,其余73.3%为部分准确,无完全不准确的结果;总体上看,CMTC廊道的准确率较LCP高,因而CMTC模型模拟白头叶猴实际迁移廊道位置的准确性优于LCP模型。输入"源"要素类型、阻力面栅格尺度设定、栖息地土地利用类型变化以及动物迁移行为复杂性4个因素是影响该模拟结果准确性的主要原因。     

Diet of Eurasian lynx <Emphasis Type="Italic">Lynx lynx</Emphasis> in the northern Dinaric Mountains (Slovenia and Croatia)

The Eurasian lynx Lynx lynx (Linnaeus 1758) is an opportunistic predator that usually selects the smallest ungulate available. Its diet varies considerably among different regions; therefore it is important to study lynx diet in different parts of the species’ range. We studied the diet of lynx from the endangered Dinaric population in Slovenia and Croatia by analyzing lynx scats, prey remains, and stomach contents. Dinaric lynx mainly killed European roe deer Capreolus capreolus (0.64 frequency of occurrence, 79% of all consumed biomass), which were used more frequently during winter and spring. Ungulates were killed more often by adult males than by lynx of other age and sex groups. In contrast to studies from other regions, lynx in the northern Dinaric Mountains also frequently fed on the edible dormice Glis glis (0.18 frequency of occurrence, 7% of all consumed biomass). This large rodent appears to be an important alternative prey, especially for females and young lynx, and was the reason for the highest use of rodents reported so far for the Eurasian lynx. Edible dormice in Dinaric forests have highly variable numbers of active animals. Seasonal and possibly annual variation in dormouse availability obviously affects lynx diet. This is a rare example where variability in the availability of the alternative prey and not the preferred prey leads to the dietary shift. This study confirms the opportunistic nature of Eurasian lynx and the regional variability of its diet.     

Geoclimatic factors influence the population genetic connectivity of Incarvillea arguta (Bignoniaceae) in the Himalaya–Hengduan Mountains biodiversity hotspot

Geoclimatic factors related to the uplift of the Himalaya and the Quaternary climatic oscillations influence the population genetic connectivity in the Himalaya–Hengduan Mountains (HHM) biodiversity hotspot. Therefore, to explore the relative roles played by these two factors, we examined the population dynamics and dispersal corridors of Incarvillea arguta (Royle) Royle incorporating ensemble species distribution modelling (SDM). Thirty‐seven populations were genotyped using plastid chloroplast DNA and low copy nuclear gene (ncpGS) sequences. Phylogeographic analysis was carried out to reveal the genetic structure and lineage differentiation. Ensemble SDMs were carried out for distributional change in the last glacial maximum, present, and future. Finally, the least cost path method was used to trace out possible dispersal corridors. The haplotypes were divided into four clades with strong geographical structure. The late Miocene origin of I. arguta in the western Himalaya ca. 7.92 Ma indicates lineage diversification related to the uplift of the HHM. The variability in habitat connectivity revealed by SDM is due to change in suitability since the Pleistocene. A putative dispersal corridor was detected along the drainage systems and river valleys, with strong support in the eastern Hengduan Mountains group. Our results support the signature of geoclimatic influence on population genetic connectivity of I. arguta in the HHM. We proposed that the major drainage systems might have assisted the rapid dispersal of isolated riverine plant species I. arguta in the HHM. The population genetic connectivity, using the fine‐tuned ensemble SDMs, enables scientists and policymakers to develop conservation strategies for the species gene pool in the HHM biodiversity hotspots.     

Effects of matrix heterogeneity on animal dispersal: from individual behavior to metapopulation-level parameters

Mounting theoretical and empirical evidence shows that matrix heterogeneity may have contrasting effects on metapopulation dynamics by contributing to patch isolation in nontrivial ways. We analyze the movement properties during interpatch dispersal in a metapopulation of Iberian lynx (Lynx pardinus). On a daily temporal scale, lynx habitat selection defines two types of matrix habitats where individuals may move: open and dispersal habitats (avoided and used as available, respectively). There was a strong and complex impact of matrix heterogeneity on movement properties at several temporal scales (hourly and daily radiolocations and the entire dispersal event). We use the movement properties on the hourly temporal scale to build a simulation model to reconstruct individual dispersal events. The two most important parameters affecting model predictions at both the individual (daily) and metapopulation scales were related to the movement capacity (number of movement steps per day and autocorrelation in dispersal habitat) followed by the parameters representing the habitat selection in the matrix. The model adequately reproduced field estimates of population-level parameters (e.g., interpatch connectivity, maximum and final dispersal distances), and its performance was clearly improved when including the effect of matrix heterogeneity on movement properties. To assume there is a homogeneous matrix results in large errors in the estimate of interpatch connectivity, especially for close patches separated by open habitat or corridors of dispersal habitat, showing how important it is to consider matrix heterogeneity when it is present. Movement properties affect the interaction of dispersing individuals with the landscape and can be used as a mechanistic representation of dispersal at the metapopulation level. This is so when the effect of matrix heterogeneity on movement properties is evaluated under biologically meaningful spatial and temporal scales.     

Environmental constraints in the colonization sequence of roe deer (Capreolus capreolus Linnaeus, 1758) across the Iberian Mountains, Spain

Aim This paper describes the dispersal process of roe deer (Capreolus capreolus Linnaeus, 1758) with respect to climatic factors, landscape characteristics and human activity. We hypothesized that environmental characteristics constrain the relative abundance and dispersal process of roe deer. Location The study was conducted in the Iberian Mountains, north‐eastern Spain, during 1986–2000. Methods Roe deer colonization dates in the study area were obtained from a survey for large mammals and from direct interviews with the employees of the Fish and Game Agency. We used a 10 × 10 km UTM grid as the sampling unit (n = 91). The relative abundance of roe deer was estimated by counting the number of pellet groups in line transects, which were representative of the habitat availability in the 10 × 10 km UTM grid. Climatic factors were obtained from meteorological stations placed near to the plot. Landscape structure indices, topography and human activity factors were obtained from digital maps using fragstats 3.3 and idrisi 32. We discarded a number of variables with no statistical significance and avoided multicollinearity by using Spearman rank correlation. Then, we used GLMz (with a multinomial error distribution and a logit link function) to analyse the influence of each variable considered in the dispersal process. Finally, GLMz (with a binomial error distribution and a logit link function) were used a posteriori to differentiate between the effects of the explanatory factors on a particular phase of the dispersal process. Results Our results indicate that proximity to a previously colonized grid significantly affected the global process of roe deer colonization. Independently of the proximity of the nearest population in the previous phase of colonization, our results also indicated that the dispersal process was influenced by precipitation and landscape structure, leading the species to colonize even apparently hostile places. Original nuclei of these populations occurred in territories with high and constant precipitation, and a landscape formed by mosaics of agricultural land with a high proportion of forests. During the expansion process, roe deer colonized territories with summer droughts, dominated by large agricultural patches and with few forest patches. Main conclusions These data support the working hypothesis that environmental characteristics constrain roe deer relative abundance and dispersal process through the Iberian Mountains. The location of the original nuclei of these populations probably had environmental conditions that were more favourable for the roe deer. Areas settled in the final phases of the dispersal process had low relative abundances of roe deer, and are likely to act mainly as dispersal corridors rather than being able to sustain viable populations themselves.