Resistance surfaces are often used to fill gaps in our knowledge surrounding animal movement and are frequently the basis for modeling connectivity associated with conservation initiatives. However, ...the methods for quantifying resistance surfaces are varied and there is no general consensus on the appropriate choice of environmental data or analytical approaches. We provide a comprehensive review of the literature on this topic to highlight methods used and identify knowledge gaps. Our review includes 96 papers that parameterized resistance surfaces (sometimes using multiple approaches) for a variety of taxa. Data types used included expert opinion (
n
= 76), detection (
n
= 23), relocation (
n
= 8), pathway (
n
= 2), and genetic (
n
= 28). We organized the papers into three main analytical approaches; one-stage expert opinion, one-stage empirical, and two-stage empirical, each of which was represented by 43, 22, and 36 papers, respectively. We further organized the empirical approaches into five main resource selection functions; point (
n
= 16), matrix (
n
= 38), home range (
n
= 3), step (
n
= 1), and pathway (
n
= 1). We found a general lack of justification for choice of environmental variables and their thematic and spatial representation, a heavy reliance on expert opinion and detection data, and a tendency to confound movement behavior and resource use. Future research needs include comparative analyses on the choice of environmental variables and their spatial and thematic scales, and on the various biological data types used to estimate resistance. Comparative analyses amongst analytical processes is also needed, as well as transparency in reporting on uncertainty in parameter estimates and sensitivity of final resistance surfaces, especially if the resistance surfaces are to be used for conservation and planning purposes.
The importance of examining multiple hierarchical levels when modeling resource use for wildlife has been acknowledged for decades. Multi-level resource selection functions have recently been ...promoted as a method to synthesize resource use across nested organizational levels into a single predictive surface. Analyzing multiple scales of selection within each hierarchical level further strengthens multi-level resource selection functions. We extend this multi-level, multi-scale framework to modeling resistance for wildlife by combining multi-scale resistance surfaces from two data types, genetic and movement. Resistance estimation has typically been conducted with one of these data types, or compared between the two. However, we contend it is not an either/or issue and that resistance may be better-modeled using a combination of resistance surfaces that represent processes at different hierarchical levels. Resistance surfaces estimated from genetic data characterize temporally broad-scale dispersal and successful breeding over generations, whereas resistance surfaces estimated from movement data represent fine-scale travel and contextualized movement decisions. We used telemetry and genetic data from a long-term study on pumas (Puma concolor) in a highly developed landscape in southern California to develop a multi-level, multi-scale resource selection function and a multi-level, multi-scale resistance surface. We used these multi-level, multi-scale surfaces to identify resource use patches and resistant kernel corridors. Across levels, we found puma avoided urban, agricultural areas, and roads and preferred riparian areas and more rugged terrain. For other landscape features, selection differed among levels, as did the scales of selection for each feature. With these results, we developed a conservation plan for one of the most isolated puma populations in the U.S. Our approach captured a wide spectrum of ecological relationships for a population, resulted in effective conservation planning, and can be readily applied to other wildlife species.
Species across the planet are shifting their ranges to track suitable climate conditions in response to climate change. Given that protected areas have higher quality habitat and often harbor higher ...levels of biodiversity compared to unprotected lands, it is often assumed that protected areas can serve as steppingstones for species undergoing climate‐induced range shifts. However, there are several factors that may impede successful range shifts among protected areas, including the distance that must be traveled, unfavorable human land uses and climate conditions along potential movement routes, and lack of analogous climates. Through a species‐agnostic lens, we evaluate these factors across the global terrestrial protected area network as measures of climate connectivity, which is defined as the ability of a landscape to facilitate or impede climate‐induced movement. We found that over half of protected land area and two‐thirds of the number of protected units across the globe are at risk of climate connectivity failure, casting doubt on whether many species can successfully undergo climate‐induced range shifts among protected areas. Consequently, protected areas are unlikely to serve as steppingstones for a large number of species under a warming climate. As species disappear from protected areas without commensurate immigration of species suited to the emerging climate (due to climate connectivity failure), many protected areas may be left with a depauperate suite of species under climate change. Our findings are highly relevant given recent pledges to conserve 30% of the planet by 2030 (30 × 30), underscore the need for innovative land management strategies that allow for species range shifts, and suggest that assisted colonization may be necessary to promote species that are adapted to the emerging climate.
There are several factors that may impede successful climate‐induced range shifts among protected areas, including the distance that must be traveled, unfavorable human land uses and climate conditions along potential movement routes, and lack of analogous climates. We found that over half of protected land area and two‐thirds of the number of protected units across the globe are at risk of climate connectivity failure, casting doubt on whether many species can successfully undergo climate‐induced range shifts among protected areas. Consequently, protected areas are unlikely to serve as steppingstones for a large number of species under a warming climate.
Aim: There is enormous interest in applying connectivity modelling to resistance surfaces for identifying corridors for conservation action. However, the multiple analytical approaches used to ...estimate resistance surfaces and predict connectivity across resistance surfaces have not been rigorously compared, and it is unclear what methods provide the best inferences about population connectivity. Using a large empirical data set on puma (Puma concolor), we are the first to compare several of the most common approaches for estimating resistance and modelling connectivity and validate them with dispersal data. Location: Southern California, USA. Methods: We estimate resistance using presence-only data, GPS telemetry data from puma home ranges and genetic data using a variety of analytical methods. We model connectivity with cost distance and circuit theory algorithms. We then measure the ability of each data type and connectivity algorithm to capture GPS telemetry points of dispersing pumas. Results: We found that resource selection functions based on GPS telemetry points and paths outperformed species distribution models when applied using cost distance connectivity algorithms. Point and path selection functions were not statistically different in their performance, but point selection functions were more sensitive to the transformation used to convert relative probability of use to resistance. Point and path selection functions and landscape genetics outperformed other methods when applied with cost distance; no methods outperformed one another with circuit theory. Main conclusions: We conclude that path or point selection functions, or landscape genetic models, should be used to estimate landscape resistance for wildlife. In cases where resource limitations prohibit the collection of GPS collar or genetic data, our results suggest that species distribution models, while weaker, may still be sufficient for resistance estimation. We recommend the use of cost distance-based approaches, such as least-cost corridors and resistant kernels, for estimating connectivity and identifying functional corridors for terrestrial wildlife.
Context
Roads have several negative effects on large mammals including restricting movements, isolating populations, and mortality due to vehicle collisions. Where large mammals regularly cross ...roads, driver safety is also a concern. Wildlife road crossing structures are often proposed to mitigate the negative effects on wildlife and human safety. However, few studies have compared the efficacy of different methods in identifying wildlife road crossing and mitigation locations.
Objectives
We examined five commonly used methods to identify road crossing locations for wildlife to determine which method best captured empirical crossing locations.
Methods
We used GPS collar data on black bears in Massachusetts, USA, to estimate road crossing locations with a road crossing resource selection function, factorial least-cost paths, resistant kernels, Circuitscape, and an individual based movement model. We evaluated model performance on each road class and all road classes combined with a hold out road crossing data set.
Results
All methods performed well, but Circuitscape consistently outperformed the other methods, both for each road class and all road classes combined.
Conclusions
Road mitigation efforts for wildlife are often costly and have a long-term footprint on the landscape, therefore, it is important to select locations for these efforts that provide functional connectivity for wildlife. We are the first to compare different road crossing models for large mammals across a large study area and on different road types. Our results provide information to assist researchers and managers in selecting an analytical method for identifying potential road mitigation locations for large mammals.
Landscape connectivity is increasingly promoted as a conservation tool to combat the negative effects of habitat loss, fragmentation, and climate change. Given its importance as a key conservation ...strategy, connectivity science is a rapidly growing discipline. However, most landscape connectivity models consider connectivity for only a single snapshot in time, despite the widespread recognition that landscapes and ecological processes are dynamic. In this paper, we discuss the emergence of dynamic connectivity and the importance of including dynamism in connectivity models and assessments. We outline dynamic processes for both structural and functional connectivity at multiple spatiotemporal scales and provide examples of modeling approaches at each of these scales. We highlight the unique challenges that accompany the adoption of dynamic connectivity for conservation management and planning in the context of traditional conservation prioritization approaches. With the increased availability of time series and species movement data, computational capacity, and an expanding number of empirical examples in the literature, incorporating dynamic processes into connectivity models is more feasible than ever. Here, we articulate how dynamism is an intrinsic component of connectivity and integral to the future of connectivity science.
Broad scale population estimates of declining species are desired for conservation efforts. However, for many secretive species including large carnivores, such estimates are often difficult. Based ...on published density estimates obtained through camera trapping, presence/absence data, and globally available predictive variables derived from satellite imagery, we modelled density and occurrence of a large carnivore, the jaguar, across the species' entire range. We then combined these models in a hierarchical framework to estimate the total population. Our models indicate that potential jaguar density is best predicted by measures of primary productivity, with the highest densities in the most productive tropical habitats and a clear declining gradient with distance from the equator. Jaguar distribution, in contrast, is determined by the combined effects of human impacts and environmental factors: probability of jaguar occurrence increased with forest cover, mean temperature, and annual precipitation and declined with increases in human foot print index and human density. Probability of occurrence was also significantly higher for protected areas than outside of them. We estimated the world's jaguar population at 173,000 (95% CI: 138,000-208,000) individuals, mostly concentrated in the Amazon Basin; elsewhere, populations tend to be small and fragmented. The high number of jaguars results from the large total area still occupied (almost 9 million km2) and low human densities (< 1 person/km2) coinciding with high primary productivity in the core area of jaguar range. Our results show the importance of protected areas for jaguar persistence. We conclude that combining modelling of density and distribution can reveal ecological patterns and processes at global scales, can provide robust estimates for use in species assessments, and can guide broad-scale conservation actions.
Conservation and management of wide-ranging carnivores like cougars (
Puma concolor
), which occur across human-altered landscapes can benefit from an in-depth understanding of their genetic status. ...Here, we apply the largest collection of multi-locus genotypes currently available for cougars (
n
= 1,903) to provide a comprehensive assessment of genetic diversity, gene flow, and source-sink dynamics for cougars occurring across Washington, United States and south-central British Columbia, Canada. We found that cougars in the Olympic, Cascade, Kettle, Selkirk, and Blue Mountains ecosystems are genetically differentiated into two clusters with varying degrees of admixture, indicating moderate levels of gene flow across the area with the exception of the Olympic Peninsula and the Blue Mountains which form more distinct genetic groups. We detected several first-generation migrants confirming long-distance movements within our study system, but also observed that migration rates between areas were asymmetrical, which is an indication of genetic source-sink dynamics. Genetic diversity and inbreeding followed a clinal east-to-west pattern with Olympic Peninsula cougars having the lowest genetic diversity and highest inbreeding coefficients among all sites. Spatial autocorrelation results for cougars did not follow sex-specific patterns suggesting that anthropogenic pressures such as habitat fragmentation and/or mortality sources may have an impact on their spatial dynamics. As cougar habitat in the northwestern United States continues to be affected by rising levels of urbanization and anthropogenic activities, long-term regional genetic monitoring represents a critical decision-support tool for formulating effective cougar conservation and management actions to prevent further genetic decline and promote long-term persistence of cougar populations.
With the growth and expansion of human development, large mammals will increasingly encounter humans, elevating the likelihood of human-wildlife conflicts. Understanding the behavior and movement of ...large mammals, particularly around human development, is important for crafting effective conservation and management plans for these species.
We used GPS collar data from American black bears (
) to determine how seasonal food resources and human development affected bear movement patterns and resource use across the Commonwealth of Massachusetts.
We found that though bears moved more and avoided human development during crepuscular and daylight hours than at night, bears preferentially moved through human dominated areas at night. This indicates bears were mitigating the risk of human development by altering their behavior to exploit these areas when human activity is low. This behavioral shift was most prominent in the spring, when natural foods are scarce, and fall, when energetic demands are high. We also observed a high degree of inter-individual variability among our sample of bears. Bears with a higher density of houses in their home ranges (~ 75 houses/km2) displayed less avoidance of human development than more rural bears. Furthermore, bear movement models had different explanatory variables, with preference or avoidance of a variable being dependent on the individual bear. To account for this individuality in our predictive surfaces, we projected the probability of movement for each season and time of day using a spatially weighted surface centered on each bear's home range.
We found that black bears in Massachusetts are operating in a landscape of fear and are altering their movement patterns to use developed areas when human activity is low. We also found seasonal and diel differences among individual bears in resource selection during movement. Accounting for these individual, seasonal, and diel differences when assessing movement for large mammals is especially important if predictive surfaces are to be used in identifying areas for conservation and management.
Estimating landscape resistance to animal movement is the foundation for connectivity modeling, and resource selection functions based on point data are commonly used to empirically estimate ...resistance. In this study, we used GPS data points acquired at 5-min intervals from radiocollared pumas in southern California to model context-dependent point selection functions. We used mixed-effects conditional logistic regression models that incorporate a paired used/available design to examine the sensitivity of point selection functions to the scale of available habitat and to the behavioral state of individual animals. We compared parameter estimates, model performance, and resistance estimates across 37 scales of available habitat, from 250 to 10,000 m, and two behavioral states, resource use and movement. Point selection functions and resistance estimates were sensitive to the chosen scale of the analysis. Multiple characteristic scales were found across our predictor variables, indicating that pumas in the study area are responding at different scales to different landscape features and that multi-scale models may be more appropriate. Additionally, point selection functions and resistance estimates were sensitive to behavioral state; specifically, pumas engaged in resource use behavior had an opposite selection response to some land cover types than pumas engaged in movement behavior. We recommend examining a continuum of scales and behavioral states when using point selection functions to estimate resistance.