Safe space in the woods: Mechanistic spatial models for predicting risks of human–bear conflicts in India

Human–wildlife interactions can have negative consequences when they involve large carnivores. Spatial risk modelling could serve as a useful management approach for predicting and pre-emptively mitigating negative interactions. We present a mechanistic modelling framework and examine interactions between humans and sloth bears (Melursus ursinus) in a multi-use forest landscape of central India. We first assessed patterns and determinants of bear distribution across the landscape using indirect sign surveys. At the same spatial scale, we then estimated spatial probabilities of bear attacks on people using information from 675 interviews with local residents, incorporating estimates of distribution probabilities from the previous step. We found the average occupancy probability across 128 grid-cells to be 0.77 (SE = 0.03). Bear occupancy was influenced by terrain ruggedness, forest composition and configuration, vegetation productivity and size of human settlements. The average probability of a bear attack in any given grid-cell was 0.61 (SE = 0.03), mostly determined by bear occurrence patterns, forest cover, terrain ruggedness, and size of human settlements. Using spatial information on people's dependence on forest resources, we identified locations with the highest risk of bear attacks. Our study demonstrates that human attacks by bears—generally believed to be random or incidental—in fact showed deterministic patterns. Our framework can be applied to other scenarios involving human–wildlife conflicts. Based on our findings, we propose that a proactive co-management approach which involves collaboration between wildlife managers and local residents could help better manage human–bear conflicts in central India and elsewhere across the species' range.     

Catchment- and site-scale influences of forest cover and longitudinal forest position on the distribution of a diadromous fish

1. The hydrologic connectivity between landscape elements and streams means that fragmentation of terrestrial habitats could affect the distribution of stream faunas at multiple spatial scales. We investigated how catchment‐ and site‐scale influences, including proportion and position of forest cover within a catchment, and presence of riparian forest cover affected the distribution of a diadromous fish. 2. The occurrence of koaro (Galaxias brevipinnis) in 50‐m stream reaches with either forested or non‐forested riparian margins at 172 sites in 24 catchments on Banks Peninsula, South Island, New Zealand was analysed. Proportions of catchments forested and the dominant position (upland or lowland) of forest within catchments were determined using geographical information system spatial analysis tools. 3. Multivariate analysis of variance indicated forest position and proportion forested at the catchment accounted for the majority of the variation in the overall proportion of sites in a catchment with koaro. 4. Where forest was predominantly in the lower part of the catchments, the presence of riparian cover was important in explaining the proportion of sites with koaro. However, where forest was predominantly in the upper part of the catchment, the effect of riparian forest was not as strong. In the absence of riparian forest cover, no patterns of koaro distribution with respect to catchment forest cover or forest position were detected. 5. These results indicate that landscape elements, such as the proportion and position of catchment forest, operating at catchment‐scales, influence the distribution of diadromous fish but their influence depends on the presence of riparian vegetation, a site‐scale factor.     

Spatial genetic structure and landscape connectivity in black bears: Investigating the significance of using different land cover datasets and classifications in landscape genetics analyses

Landscape genetic analyses allow detection of fine‐scale spatial genetic structure (SGS) and quantification of effects of landscape features on gene flow and connectivity. Typically, analyses require generation of resistance surfaces. These surfaces characteristically take the form of a grid with cells that are coded to represent the degree to which landscape or environmental features promote or inhibit animal movement. How accurately resistance surfaces predict association between the landscape and movement is determined in large part by (a) the landscape features used, (b) the resistance values assigned to features, and (c) how accurately resistance surfaces represent landscape permeability. Our objective was to evaluate the performance of resistance surfaces generated using two publicly available land cover datasets that varied in how accurately they represent the actual landscape. We genotyped 365 individuals from a large black bear population (Ursus americanus) in the Northern Lower Peninsula (NLP) of Michigan, USA at 12 microsatellite loci, and evaluated the relationship between gene flow and landscape features using two different land cover datasets. We investigated the relative importance of land cover classification and accuracy on landscape resistance model performance. We detected local spatial genetic structure in Michigan''s NLP black bears and found roads and land cover were significantly correlated with genetic distance. We observed similarities in model performance when different land cover datasets were used despite 21% dissimilarity in classification between the two land cover datasets. However, we did find the performance of land cover models to predict genetic distance was dependent on the way the land cover was defined. Models in which land cover was finely defined (i.e., eight land cover classes) outperformed models where land cover was defined more coarsely (i.e., habitat/non‐habitat or forest/non‐forest). Our results show that landscape genetic researchers should carefully consider how land cover classification changes inference in landscape genetic studies.     

Spatial Distribution of Black Bear Incident Reports in Michigan

Interactions between humans and carnivores have existed for centuries due to competition for food and space. American black bears are increasing in abundance and populations are expanding geographically in many portions of its range, including areas that are also increasing in human density, often resulting in associated increases in human-bear conflict (hereafter, bear incidents). We used public reports of bear incidents in Michigan, USA, from 2003–2011 to assess the relative contributions of ecological and anthropogenic variables in explaining the spatial distribution of bear incidents and estimated the potential risk of bear incidents. We used weighted Normalized Difference Vegetation Index mean as an index of primary productivity, region (i.e., Upper Peninsula or Lower Peninsula), primary and secondary road densities, and percentage land cover type within 6.5-km² circular buffers around bear incidents and random points. We developed 22 a priori models and used generalized linear models and Akaike’s Information Criterion (AIC) to rank models. The global model was the best compromise between model complexity and model fit (w = 0.99), with a ΔAIC 8.99 units from the second best performing model. We found that as deciduous forest cover increased, the probability of bear incident occurrence increased. Among the measured anthropogenic variables, cultivated crops and primary roads were the most important in our AIC-best model and were both positively related to the probability of bear incident occurrence. The spatial distribution of relative bear incident risk varied markedly throughout Michigan. Forest cover fragmented with agriculture and other anthropogenic activities presents an environment that likely facilitates bear incidents. Our map can help wildlife managers identify areas of bear incident occurrence, which in turn can be used to help develop strategies aimed at reducing incidents. Researchers and wildlife managers can use similar mapping techniques to assess locations of specific conflict types or to address human impacts on endangered species.     

Identifying Habitat Linkages for American Black Bears in North Carolina,USA

ABSTRACT Understanding landscape structure and the role of habitat linkages is important to managing wildlife populations in fragmented landscapes. We present a data-based method for identifying local- and regional-scale habitat linkages for American black bears (Ursus americanus) on the Albemarle-Pamlico Peninsula of North Carolina, USA. We used weights-of-evidence, a discrete multivariate technique for combining spatial data, to make predictions about bear habitat use from 1,771 telemetry locations on 2 study areas (n = 35 bears). The model included 3 variables measured at a 0.2-km² scale: forest cohesion, forest diversity, and forest-agriculture edge density, adequately describing important habitat characteristics for bears on our study area. We used 2 categories of unique habitat conditions to delineate favorable bear habitat, which correctly classified 79.5% of the bear locations in a 10-fold model validation. Forest cohesion and forest-agriculture edge density were the most powerful predictors of black bear habitat use. We used predicted probabilities of bear occurrence from the model to delineate habitat linkages among local and regional areas where bear densities were relatively high. Our models clearly identified 2 of the 3 sites previously recommended for wildlife underpasses on a new, 4-lane highway in the study area. Our approach yielded insights into how landscape metrics can be integrated to identify linkages suitable as habitat and dispersal routes.     

Landscape partitioning and spatial inferences of competition between black and grizzly bears

Population effects of competition between large carnivore species may be evident by contrasting actual distributions of putative competitors against predictions of inherent landscape quality for each species. Such comparison can be insightful if covariation with external factors known to influence the occurrence, density, or persistence of each species over space and time can be controlled. We used systematically‐distributed DNA hair‐trap stations to sample the occurrence of black bears (Ursus americanus) and grizzly bears (U. arctos) across 5496 km² in southeastern British Columbia, Canada. We describe interspecific landscape partitioning according to terrain, vegetation and land‐cover variables at 2 spatial scales. We developed multivariate models to predict the potential distribution of each species. At sampling site‐session combinations that detected either species, we then investigated whether the expected or actual occurrence of each influenced the likelihood of detecting the other while controlling for human influence and inherent landscape quality. Black bears were more likely than grizzly bears to occur in gentle, valley bottom terrain with lower proportions of open habitats. Each species also was detected less frequently with the other species than predicted by their respective models; however, the strength of this relationship decreased as landscapes became more characteristic of black bear habitat. As landscapes showed higher inherent potential to support grizzly bears, black bears occurred more than model prediction in areas with higher human access and proximity to major highways but less in national parks. As potential to support black bears increased, grizzly bears occurred more than model prediction only in national parks and less with increasing human access and proximity to major highways. Results suggest that competition is occurring between the species, and that the differential response of each species to human disturbance or excessive mortality may influence the outcome and hence landscape partitioning. Moreover, black bears are more likely to benefit from human encroachment into landscapes of high inherent value for grizzly bears than vice versa. Conservation implications relate to potential mediating effects of habitat and human influence on competitive interactions between the species.     

A human-induced landscape of fear influences foraging behavior of brown bears

Animals adapt their foraging behavior to variations in food availability and predation risk. In Sweden, brown bears (Ursus arctos) depend on a nearly continuous intake of berries, especially bilberries (Vaccinium myrtillus) during late summer and early autumn to fatten up prior to hibernation. This overlaps with the bear hunting season that starts on 21 August. Bilberry occurrence varies across space, as does human-induced mortality risk. Here, we hypothesize that brown bears select for areas with a high probability of bilberry occurrence across a boreal forest ecosystem in Sweden (H1), and that human-induced mortality risk reduces bear selection for bilberries (H2). In addition, we hypothesized that bears that survived the hunting season avoided bilberry areas associated with high risk, whereas bears that were later killed selected more strongly for berries and less against risk prior to the hunting season (H3). To evaluate our hypotheses, we used resource selection functions to contrast bear GPS relocation data (N = 35, 2012–2015) and random positions within the bearś home range with generalized linear mixed effect models against two focal variables: a map predicting bilberry occurrence and a map predicting human-induced mortality risk. We found that bears selected for areas with a high probability of bilberry occurrence (supporting H1), but avoided these areas if they were associated with and high risk of hunting mortality (supporting H2). The killed and surviving bears did not differ in their selection for bilberries, but they did differ in their selection against risk (partially supporting H3). Surviving bears strongly avoided high risk areas, whereas killed bears responded less to risk and selected for high-risk areas with a low probability of bilberry occurrence. This suggests that killed bears selected for other food sources than berries in high risk areas, which exposed them to human hunters. We conclude that bears respond to a landscape of fear during the berry season and that different foraging strategies may have a direct impact on individual mortality during the hunting season.     

Resource Selection by Recolonizing American Black Bears

American black bears (Ursus americanus) were extirpated from Oklahoma, USA, in the early twentieth century but have since recolonized eastern portions of the state after immigrating from Arkansas, where they were successfully translocated. Within the last 2 decades, a population of black bears was detected in the Oklahoma Ozark region, prompting studies to determine population size, growth rate, and genetic makeup. To understand how black bears were recolonizing the human-dominated landscape, we investigated resource selection at 2 scales. Between 2011 and 2016, we collected global positioning system collar spatial data for 10 males and 13 females. We calculated average kernel density home ranges on a seasonal scale for all collared bears. We used generalized linear mixed models to calculate resource selection functions at the study area, defined by locations of all radio-collared black bears (second order) and the scale of individual black bear home ranges (third order). Resource selection did not differ significantly by sex. Black bears across seasons and scales selected riparian forest and moist oak (Quercus spp.) forest land cover types and mostly selected against indicators of human activity (e.g., pasture-prairie, anthropogenic land cover types, roads, and areas of high human population density). Black bears also selected areas with rugged terrain at high elevations, although not consistently across seasons and scales. Black bear recolonization appeared to be negatively affected by areas and features characterized as human-altered. Further expansion of the range of black bears may be limited by anthropogenic disturbance in the region. © 2021 The Wildlife Society.     

陕西汉中地区黑熊的现状,分布及保护措施

１９９１－１９９３年,作者在陕西省汉中林区采用三阶抽样设计调查并询访群众,求算每个所选择的最小样本单位（三级区）的方差,推算黑熊整体数量。根据黑熊数量将该地林区划分为３个密度区级：高密度区（０．０３头／ｋｍ￣２以上）,秦岭中山区;中密度区（０．０１－０．０２９头／ｋｍ￣２）,巴山中山区;低密度区（０．０１头／ｋｍ￣２以下）,秦巴低山丘陵区。全区有黑熊３１２头,平均密度０．０３５／ｋｍ￣２,遍布１０县１６６个乡。分析了各地黑熊数量多寡的原因,提出了相应的保护措施。     

Distribution and abundance of trees in floodplain forests of the Wisconsin River: Environmental influences at different scales

Questions: 1. How do physiography, flooding regime, landscape pattern, land‐cover history, and local soil conditions influence the presence, community structure and abundance of overstorey trees? 2. Can broad‐scale factors explain variation in the floodplain forest community, or are locally measured soil conditions necessary? Location: Floodplain of the lower 370 km of the Wisconsin River, Wisconsin, USA. Methods: Floodplain forest was sampled in 10 m × 20 m plots [n= 405) during summers of 1999 and 2000 in six 12‐ to 15‐km reaches. Results: Species observed most frequently were Fraxinus pennsylvanica, Acer saccharinum and Ulmus americana. Physiography (e.g. geographic province) and indicators of flooding regime (e.g. relative elevation and distance from main channel) were consistently important in predicting occurrence, community composition, and abundance of trees. Correspondence analysis revealed that flood‐tolerant and intolerant species segregated along the primary axis, and late‐successional species segregated from flood‐tolerant species along the secondary axis. Current landscape configuration only influenced species presence or abundance in forests that developed during recent decades. Land‐cover history was important for tree species presence and for the abundance of late‐successional species. Comparison of statistical models developed with and without soils data suggested that broad‐scale factors such as geographic province generally performed well. Conclusions: Physiography and indicators of flood regime are particularly useful for explaining floodplain forest structure and composition in floodplains with a relatively high proportion of natural cover types.     

Presence-absence versus presence-only modelling methods for predicting bird habitat suitability

Habitat suitability models can be generated using methods requiring information on species presence or species presence and absence. Knowledge of the predictive performance of such methods becomes a critical issue to establish their optimal scope of application for mapping current species distributions under different constraints. Here, we use breeding bird atlas data in Catalonia as a working example and attempt to analyse the relative performance of two methods: the Ecological Niche factor Analysis (ENFA) using presence data only and Generalised Linear Models (GLM) using presence/absence data. Models were run on a set of forest species with similar habitat requirements, but with varying occurrence rates (prevalence) and niche positions (marginality). Our results support the idea that GLM predictions are more accurate than those obtained with ENFA. This was particularly true when species were using available habitats proportionally to their suitability, making absence data reliable and useful to enhance model calibration. Species marginality in niche space was also correlated to predictive accuracy, i.e. species with less restricted ecological requirements were modelled less accurately than species with more restricted requirements. This pattern was irrespective of the method employed. Models for wide‐ranging and tolerant species were more sensitive to absence data, suggesting that presence/absence methods may be particularly important for predicting distributions of this type of species. We conclude that modellers should consider that species ecological characteristics are critical in determining the accuracy of models and that it is difficult to predict generalist species distributions accurately and this is independent of the method used. Being based on distinct approaches regarding adjustment to data and data quality, habitat distribution modelling methods cover different application areas, making it difficult to identify one that should be universally applicable. Our results suggest however, that if absence data is available, methods using this information should be preferably used in most situations.     

Unraveling the mystery of the glacier bear: Genetic population structure of black bears (Ursus americanus) within the range of a rare pelage type

Glacier bears are a rare grey color morph of American black bear (Ursus americanus) found only in northern Southeast Alaska and a small portion of western Canada. We examine contemporary genetic population structure of black bears within the geographic extent of glacier bears and explore how this structure relates to pelage color and landscape features of a recently glaciated and highly fragmented landscape. We used existing radiocollar data to quantify black bear home‐range size within the geographic range of glacier bears. The mean home‐range size of female black bears in the study area was 13 km² (n = 11), whereas the home range of a single male was 86.9 km². We genotyped 284 bears using 21 microsatellites extracted from noninvasively collected hair as well as tissue samples from harvested bears. We found ten populations of black bears in the study area, including several new populations not previously identified, divided largely by geographic features such as glaciers and marine fjords. Glacier bears were assigned to four populations found on the north and east side of Lynn Canal and the north and west side of Glacier Bay with a curious absence in the nonglaciated peninsula between. Lack of genetic relatedness and geographic continuity between black bear populations containing glacier bears suggest a possible unsampled population or an association with ice fields. Further investigation is needed to determine the genetic basis and the adaptive and evolutionary significance of the glacier bear color morph to help focus black bear conservation management to maximize and preserve genetic diversity.     

The contributions of topoclimate and land cover to species distributions and abundance: fine‐resolution tests for a mountain butterfly fauna

Aim Models relating species distributions to climate or habitat are widely used to predict the effects of global change on biodiversity. Most such approaches assume that climate governs coarse‐scale species ranges, whereas habitat limits fine‐scale distributions. We tested the influence of topoclimate and land cover on butterfly distributions and abundance in a mountain range, where climate may vary as markedly at a fine scale as land cover. Location Sierra de Guadarrama (Spain, southern Europe) Methods We sampled the butterfly fauna of 180 locations (89 in 2004, 91 in 2005) in a 10,800 km² region, and derived generalized linear models (GLMs) for species occurrence and abundance based on topoclimatic (elevation and insolation) or habitat (land cover, geology and hydrology) variables sampled at 100‐m resolution using GIS. Models for each year were tested against independent data from the alternate year, using the area under the receiver operating characteristic curve (AUC) (distribution) or Spearman's rank correlation coefficient (r_s) (abundance). Results In independent model tests, 74% of occurrence models achieved AUCs of > 0.7, and 85% of abundance models were significantly related to observed abundance. Topoclimatic models outperformed models based purely on land cover in 72% of occurrence models and 66% of abundance models. Including both types of variables often explained most variation in model calibration, but did not significantly improve model cross‐validation relative to topoclimatic models. Hierarchical partitioning analysis confirmed the overriding effect of topoclimatic factors on species distributions, with the exception of several species for which the importance of land cover was confirmed. Main conclusions Topoclimatic factors may dominate fine‐resolution species distributions in mountain ranges where climate conditions vary markedly over short distances and large areas of natural habitat remain. Climate change is likely to be a key driver of species distributions in such systems and could have important effects on biodiversity. However, continued habitat protection may be vital to facilitate range shifts in response to climate change.     

A Multiscale Assessment of Den Selection by Black Bears in Louisiana

ABSTRACT Understanding den selection for the federally threatened Louisiana black bear (Ursus americanus luteolus) will assist managers in effectively managing bear populations and the landscapes they inhabit. We assessed den selection of 72 female Louisiana black bears at multiple spatial scales in northern and central Louisiana, USA. We used 230 den-years to examine den type preference, microhabitat characteristics of dens, and effects of landscape characteristics on den selection. We also evaluated reuse of tree dens and their availability across several study areas. Bears selected tree dens more frequently (65%) than ground dens. Most (55%) ground dens were associated with standing or downed trees, and most tree dens were in baldcypress (Taxodium distichum; 86%) or surrounded by water (80%). Selection of ground dens was positively associated with landscape metrics related to presence of water, whereas selection of tree dens was positively associated with proximity to edge and greater proportions of swamp and areas containing water. Reuse of tree dens averaged 15% (range = 0–20%) across study areas and our estimates of tree-den availability indicated that tree dens were not a limiting resource for Louisiana black bears in the northern or central portions of Louisiana. Although Louisiana black bears demonstrated plasticity in den selection, we recommend managers use forest management practices that conserve appropriate tree dens and create ground-denning opportunities, because both would promote the recovery of Louisiana black bear populations.     

Highlighting the effects of land-use change on a threatened amphibian in a human-dominated landscape

In Western Europe, habitat loss and landscape fragmentation has led to significant population decline in various animal groups, including amphibians. The extinction of the last natural populations of the yellow-bellied toad in Belgium, Luxembourg and several regions of southern and western France suggests a widespread decline. By using site-occupancy models and adding covariates corresponding to the human-influenced features of the landscape, we tried to identify the relative effects of different land-use types on the species’ distribution pattern in a man-made environment (the Alsatian Rhine floodplain in France). We recorded presence–absence data in 150 forest sample plots (300 × 300 m) and then modeled species distribution while taking into account detection errors in the field. Land-use was recorded on two spatial scales: within the forest sample plots and in a 1500 m radius buffer area around the forest plots. In the forest plots, toad occurrence was negatively correlated with loss of forest cover to agricultural land. In contrast, occurrence is positively correlated with the density of human-made rutted dirt paths and tracks, which provide semi-natural breeding sites. In the 1500 m radius buffer zones around forest plots, toad occurrence was negatively correlated with the density of urbanization and road networks. These results can be used to plan conservation strategies for amphibians in human-dominated landscapes.     

黑龙江省10年来熊类资源的分布及变动趋势

熊类是大型动物,其毛皮、油脂以及胆、掌、肉等都有很高的利用价值,因此,在黑龙江省林区,熊类历来是重要的狩猎对象。然而,熊是大型猛兽,每被击伤或受惊扰,常凶猛袭人,故猎熊并非易事。历史文献中曾记载:“熊矫捷而罴憨猛,皆兽之绝有力者”,又说,“刺虎者,不过五人一排…刺熊罴,非十人不能胜。盖其力倍于虎也”(《满洲源流考》)。现在,黑龙江省的猎人,通常是在冬季待熊在树洞或地穴中冬眠时进行猎杀,即猎取洞中冬眠的熊,俗称“掏仓子”,或在春季捕捉刚出洞的幼熊,故每年都有一定的猎取量。     

Climate and land use changes will degrade the configuration of the landscape for titi monkeys in eastern Brazil

Land use changes have profound effects on populations of Neotropical primates, and ongoing climate change is expected to aggravate this scenario. The titi monkeys from eastern Brazil (Callicebus personatus group) have been particularly affected by this process, with four of the five species now allocated to threatened conservation status categories. Here, we estimate the changes in the distribution of these titi monkeys caused by changes in both climate and land use. We also use demographic‐based, functional landscape metrics to assess the magnitude of the change in landscape conditions for the distribution predicted for each species. We built species distribution models (SDMs) based on maximum entropy for current and future conditions (2070), allowing for different global circulation models and contrasting scenarios of glasshouse gas concentrations. We refined the SDMs using a high‐resolution map of habitat remnants. We then calculated habitat availability and connectivity based on home‐range size and the dispersal limitations of the individual, in the context of a predicted loss of 10% of forest cover in the future. The landscape configuration is predicted to be degraded for all species, regardless of the climatic settings. This include reductions in the total cover of forest remnants, patch size and functional connectivity. As the landscape configuration should deteriorate severely in the future for all species, the prevention of further loss of populations will only be achieved through habitat restoration and reconnection to counteract the negative effects for these and several other co‐occurring species.     

Foraging strategies by omnivores: are black bears actively searching for ungulate neonates or are they simply opportunistic predators?

Omnivores feed on animals with dynamic distributions and on plants with static distributions. The search tactics they adopt will not only define the risk for the targeted prey, but also for other prey that may be consumed when encountered. The potential impact of omnivores on the dynamics of multi‐prey systems thus depends on resource selection and on the tactics used to find their prey. We present an approach that can clarify the foraging decisions of omnivores by combining analyses of habitat selection, local residency time, and interpatch movements. We use this framework to evaluate whether predation by omnivorous black bears on ungulate neonates resulted from an active search or from incidental encounters. We monitored 12 bears, 22 forest‐dwelling caribou, and 36 moose during calving seasons. We estimated the spatial patterns in relative occurrence probability of ungulate neonates using Resource Selection Functions (RSFs). We also mapped plant abundance from vegetation surveys. RSF were then built to assess the link between bear distribution and the distribution of these three food types (vegetation, moose calves, caribou fawns). We further evaluated the search tactic used by bears that led to this spatial dependency by exploring patterns of residency times and interpatch movements. Bears did not select areas with a high probability of encounter with neonates, but selected areas with abundant vegetation. Surprisingly, bears displayed shorter residency times in vegetation‐rich areas. The selection for vegetation‐rich areas was therefore achieved by moving preferentially, but frequently, between areas offering abundant vegetation. Such frequent interpatch movements could result in high rates of fortuitous encounters with neonates, even if bears are not actively searching for them. To mitigate the impacts of forest harvesting on threatened caribou populations, vegetation‐rich areas selected by bears (e.g. roadsides) should be segregated from large patches of mature conifer forest suitable for caribou.     

Local coexistence and niche differences between the Lusitanian and Mediterranean pine voles (<Emphasis Type="Italic">Microtus lusitanicus</Emphasis> and <Emphasis Type="Italic">M. duodecimcostatus</Emphasis>)

In the present study, we analyzed the coexistence pattern of the Lusitanian pine vole (Microtus lusitanicus) and the Mediterranean pine vole (Microtus duodecimcostatus) in a potential area of sympatry in a Mediterranean landscape (Portugal). We also determined the relative contribution of local, landscape, and spatial factors explaining the differences in the distribution patterns of the two species in the region. Using a kriging interpolation method, we obtained a map of sympatric and allopatric areas of species occurrence. The estimated sympatry area corresponded to a northwest–southeast belt representing 11.3% of the study area. Habitat niche differences were assessed with binomial GLMs followed by a variance partitioning. At a local scale, higher altitude, higher cover of shrubs, lower clay content in the soil, and lower cover of tree canopy were the most important factors distinguishing M. lusitanicus presence sites from those with M. duodecimcostatus. At a larger scale, the presence of forest landscape units and the low abundance of “montado” units were the most influencing landscape factors in the identification of M. lusitanicus occurrence sites when compared to M. duodecimcostatus. Our results suggested that local coexistence of M. lusitanicus and M. duodecimcostatus in the field is a rare event. The differences in distribution patterns of the two pine vole species were mostly explained by fine-scale environmental factors and by shared spatial effects.