Epigenetic clocks, sex markers and age-class diagnostics in three harvested large mammals.

The development of epigenetic clocks, or the DNA methylation-based inference of age, is an emerging tool for ageing in free ranging populations. In this study, we developed epigenetic clocks for three species of large mammals that are the focus of extensive management throughout their range in North America: white-tailed deer, black bear and mountain goat. We quantified differential DNA methylation patterns at over 30,000 cytosine-guanine sites (CpGs) from tissue samples of all three species (black bear n = 49; white-tailed deer n = 47; mountain goat n = 45). We used a penalized regression model (elastic net) to build explanatory (black bear r = .95; white-tailed deer r = .99; mountain goat r = .97) and robust (black bear Median Absolute Error or MAE = 1.33; white-tailed deer MAE = 0.29; mountain goat MAE = 0.61) models of age or clocks. We also characterized individual CpG sites within each species that demonstrated clear differences in methylation levels between age classes and sex, which can be used to develop a suite of accessible diagnostic markers. This tool has the potential to contribute to wildlife monitoring by providing easily obtainable representations of age structure in managed populations.

Fragmentation effects on an endangered species across a gradient from the interior to edge of its range.

Understanding how habitat fragmentation affects individual species is complicated by challenges associated with quantifying species-specific habitat and spatial variability in fragmentation effects within a species' range. We aggregated a 29-year breeding survey data set for the endangered marbled murrelet (Brachyramphus marmoratus) from >42,000 forest sites throughout the Pacific Northwest (Oregon, Washington, and northern California) of the United States. We built a species distribution model (SDM) in which occupied sites were linked with Landsat imagery to quantify murrelet-specific habitat and then used occupancy models to test the hypotheses that fragmentation negatively affects murrelet breeding distribution and that these effects are amplified with distance from the marine foraging habitat toward the edge of the species' nesting range. Murrelet habitat declined in the Pacific Northwest by 20% since 1988, whereas the proportion of habitat comprising edges increased by 17%, indicating increased fragmentation. Furthermore, fragmentation of murrelet habitat at landscape scales (within 2 km of survey stations) negatively affected occupancy of potential breeding sites, and these effects were amplified near the range edge. On the coast, the odds of occupancy decreased by 37% (95% confidence interval [CI] -54 to 12) for each 10% increase in edge habitat (i.e., fragmentation), but at the range edge (88 km inland) these odds decreased by 99% (95% CI 98 to 99). Conversely, odds of murrelet occupancy increased by 31% (95% CI 14 to 52) for each 10% increase in local edge habitat (within 100 m of survey stations). Avoidance of fragmentation at broad scales but use of locally fragmented habitat with reduced quality may help explain the lack of murrelet population recovery. Further, our results emphasize that fragmentation effects can be nuanced, scale dependent, and geographically variable. Awareness of these nuances is critical for developing landscape-level conservation strategies for species experiencing broad-scale habitat loss and fragmentation.

Efectos de la fragmentación sobre las especies en peligro a lo largo de un gradiente desde el interior hasta el borde de su distribución Resumen Es complicado entender el efecto de la fragmentación del hábitat sobre las especies individuales debido a los retos asociados con la cuantificación de hábitats específicos por especie y la variabilidad espacial de los efectos de la fragmentación dentro de la distribución de la especie. Combinamos los datos de un censo reproductivo realizado durante 29 años para el mérgulo jaspeado (Brachyramphus marmoratus) de >42,000 sitios boscosos a lo largo del noroeste del Pacífico (Oregón, Washington, y el norte de California, EE. UU.). Construimos un modelo de distribución de especie (MDE) en el cual los sitios ocupados estuvieron vinculados con imágenes de Landsat para cuantificar el hábitat específico del mérgulo y después usamos los modelos de ocupación para comprobar la hipótesis de que la fragmentación afecta negativamente la distribución reproductiva de la especie y que estos efectos se amplifican con la distancia entre el hábitat de forrajeo marino y el borde de la distribución de anidación de la especie. El hábitat del mérgulo declinó en la zona en un 20% a partir de 1988, mientras que la proporción de hábitat que comprende bordes incrementó en un 17%, lo que indica un aumento en la fragmentación. Además, la fragmentación del hábitat del mérgulo a escala de paisaje (a de 2 km de las estaciones de censo) afectó negativamente a la ocupación de sitios potenciales de reproducción y estos efectos se amplificaron cerca del borde de la distribución. La probabilidad de ocupación disminuyó en un 37% (95% IC -54 a 12) por cada 10% de incremento en el hábitat de borde (es decir, fragmentación) en la costa, pero en el borde de la distribución (88 km tierra adentro), esta probabilidad disminuyó en un 99% (95% IC 98 a 99). De forma contraria, la probabilidad de ocupación incrementó en un 31% (95% IC 14 a 52) por cada 10% de incremento en el hábitat de borde local (a 100 m de las estaciones de censo). La evasión de la fragmentación a gran escala y el uso de hábitats con calidad reducida y fragmentados a nivel local podría explicar la falta de recuperación poblacional del mérgulo. Más allá, nuestros resultados resaltan que los efectos de la fragmentación pueden estar matizados, depender de la escala y tener variación geográfica. Es importante tener conciencia de estos matices para desarrollar estrategias de conservación a nivel paisaje para las especies que experimentan fragmentación y pérdida del hábitat a gran escala.

Deriving spatially explicit direct and indirect interaction networks from animal movement data.

Quantifying spatiotemporally explicit interactions within animal populations facilitates the understanding of social structure and its relationship with ecological processes. Data from animal tracking technologies (Global Positioning Systems ["GPS"]) can circumvent longstanding challenges in the estimation of spatiotemporally explicit interactions, but the discrete nature and coarse temporal resolution of data mean that ephemeral interactions that occur between consecutive GPS locations go undetected. Here, we developed a method to quantify individual and spatial patterns of interaction using continuous-time movement models (CTMMs) fit to GPS tracking data. We first applied CTMMs to infer the full movement trajectories at an arbitrarily fine temporal scale before estimating interactions, thus allowing inference of interactions occurring between observed GPS locations. Our framework then infers indirect interactions-individuals occurring at the same location, but at different times-while allowing the identification of indirect interactions to vary with ecological context based on CTMM outputs. We assessed the performance of our new method using simulations and illustrated its implementation by deriving disease-relevant interaction networks for two behaviorally differentiated species, wild pigs (Sus scrofa) that can host African Swine Fever and mule deer (Odocoileus hemionus) that can host chronic wasting disease. Simulations showed that interactions derived from observed GPS data can be substantially underestimated when temporal resolution of movement data exceeds 30-min intervals. Empirical application suggested that underestimation occurred in both interaction rates and their spatial distributions. CTMM-Interaction method, which can introduce uncertainties, recovered majority of true interactions. Our method leverages advances in movement ecology to quantify fine-scale spatiotemporal interactions between individuals from lower temporal resolution GPS data. It can be leveraged to infer dynamic social networks, transmission potential in disease systems, consumer-resource interactions, information sharing, and beyond. The method also sets the stage for future predictive models linking observed spatiotemporal interaction patterns to environmental drivers.

Genomic correlates for migratory direction in a free-ranging cervid.

Animal migrations are some of the most ubiquitous and one of the most threatened ecological processes globally. A wide range of migratory behaviours occur in nature, and this behaviour is not uniform among and within species, where even individuals in the same population can exhibit differences. While the environment largely drives migratory behaviour, it is necessary to understand the genetic mechanisms influencing migration to elucidate the potential of migratory species to cope with novel conditions and adapt to environmental change. In this study, we identified genes associated with a migratory trait by undertaking pooled genome-wide scans on a natural population of migrating mule deer. We identified genomic regions associated with variation in migratory direction, including FITM1, a gene linked to the formation of lipids, and DPPA3, a gene linked to epigenetic modifications of the maternal line. Such a genetic basis for a migratory trait contributes to the adaptive potential of the species and might affect the flexibility of individuals to change their behaviour in the face of changes in their environment.

Explaining detection heterogeneity with finite mixture and non-Euclidean movement in spatially explicit capture-recapture models.

Landscape structure affects animal movement. Differences between landscapes may induce heterogeneity in home range size and movement rates among individuals within a population. These types of heterogeneity can cause bias when estimating population size or density and are seldom considered during analyses. Individual heterogeneity, attributable to unknown or unobserved covariates, is often modelled using latent mixture distributions, but these are demanding of data, and abundance estimates are sensitive to the parameters of the mixture distribution. A recent extension of spatially explicit capture-recapture models allows landscape structure to be modelled explicitly by incorporating landscape connectivity using non-Euclidean least-cost paths, improving inference, especially in highly structured (riparian & mountainous) landscapes. Our objective was to investigate whether these novel models could improve inference about black bear (Ursus americanus) density. We fit spatially explicit capture-recapture models with standard and complex structures to black bear data from 51 separate study areas. We found that non-Euclidean models were supported in over half of our study areas. Associated density estimates were higher and less precise than those from simple models and only slightly more precise than those from finite mixture models. Estimates were sensitive to the scale (pixel resolution) at which least-cost paths were calculated, but there was no consistent pattern across covariates or resolutions. Our results indicate that negative bias associated with ignoring heterogeneity is potentially severe. However, the most popular method for dealing with this heterogeneity (finite mixtures) yielded potentially unreliable point estimates of abundance that may not be comparable across surveys, even in data sets with 136-350 total detections, 3-5 detections per individual, 97-283 recaptures, and 80-254 spatial recaptures. In these same study areas with high sample sizes, we expected that landscape features would not severely constrain animal movements and modelling non-Euclidian distance would not consistently improve inference. Our results suggest caution in applying non-Euclidean SCR models when there is no clear landscape covariate that is known to strongly influence the movement of the focal species, and in applying finite mixture models except when abundant data are available.

Estimating animal abundance at multiple scales by spatially explicit capture-recapture.

Information about how animal abundance varies across landscapes is needed to inform management action but is costly and time-consuming to obtain; surveys of a single population distributed over a large area can take years to complete. Surveys employing small, spatially replicated sampling units improve efficiency, but statistical estimators rely on assumptions that constrain survey design or become less reasonable as larger areas are sampled. Efficient methods that avoid assumptions about similarity of detectability or density among replicates are therefore appealing. Using simulations and data from >3500 black bears sampled on 73 independent study areas in Ontario, Canada, we (1) quantified bias induced by unmodeled spatial heterogeneity in detectability and density; (2) evaluated novel, design-based estimators of average density across replicate study areas; and (3) evaluated two estimators of the variance of average density across study areas: an analytic estimator that assumed an underlying homogeneous spatial Poisson point process for the distribution of animals' activity centers, and an empirical estimator of variance across study areas. In simulations where detectability varied in space, assuming spatially constant detectability yielded density estimates that were negatively biased by 20% to 30%; estimating local detectability and density from local data and treating study areas as independent, equal replicates when estimating average density across study areas using the design-based estimator yielded unbiased estimates at local and landscape scales. Similarly, detectability of black bears varied among study areas and estimates of bear density at landscape scales were higher when no information was shared across study areas when estimating detectability. This approach also maximized precision (relative SEs of estimates of average black bear density ranged from 7% to 18%) and computational efficiency. In simulations, the analytic variance estimator was robust to threefold variation in local densities but the empirical estimator performed poorly. Conducting multiple, similar SECR surveys and treating them as independent replicates during analyses allowed us to efficiently estimate density at multiple scales and extents while avoiding biases caused by pooling spatially heterogeneous data. This approach enables researchers to address a wide range of ecological or management-related questions and is applicable with most types of SECR data.

Harvest and density-dependent predation drive long-term population decline in a northern ungulate.

The relative effect of top-down versus bottom-up forces in regulating and limiting wildlife populations is an important theme in ecology. Untangling these effects is critical for a basic understanding of trophic dynamics and effective management. We examined the drivers of moose (Alces alces) population growth by integrating two independent sources of observations within a hierarchical Bayesian population model. We used one of the largest existing spatiotemporal data sets on ungulate population dynamics globally. We documented a 20% population decline over the period examined. There was negative density-dependent population growth of moose. Although we could not determine the mechanisms producing density-dependent suppression of population growth, the relatively low densities at which we documented moose populations suggested it could be due to density-dependent predation. Predation primarily limited population growth, except at low density, where it was regulating. After we simulated several harvest scenarios, it appeared that harvest was largely additive and likely contributed to population declines. Our results highlight how population dynamics are context dependent and vary strongly across gradients in climate, forest type, and predator abundance. These results help clarify long-standing questions in population ecology and highlight the complex relationships between natural and human-caused mortality in driving ungulate population dynamics.

Forest degradation drives widespread avian habitat and population declines.

In many regions of the world, forest management has reduced old forest and simplified forest structure and composition. We hypothesized that such forest degradation has resulted in long-term habitat loss for forest-associated bird species of eastern Canada (130,017 km2) which, in turn, has caused bird-population declines. Despite little change in overall forest cover, we found substantial reductions in old forest as a result of frequent clear-cutting and a broad-scale transformation to intensified forestry. Back-cast species distribution models revealed that breeding habitat loss occurred for 66% of the 54 most common species from 1985 to 2020 and was strongly associated with reduction in old age classes. Using a long-term, independent dataset, we found that habitat amount predicted population size for 94% of species, and habitat loss was associated with population declines for old-forest species. Forest degradation may therefore be a primary cause of biodiversity decline in managed forest landscapes.

Conceptual and methodological advances in habitat-selection modeling: guidelines for ecology and evolution.

Habitat selection is a fundamental animal behavior that shapes a wide range of ecological processes, including animal movement, nutrient transfer, trophic dynamics and population distribution. Although habitat selection has been a focus of ecological studies for decades, technological, conceptual and methodological advances over the last 20 yr have led to a surge in studies addressing this process. Despite the substantial literature focused on quantifying the habitat-selection patterns of animals, there is a marked lack of guidance on best analytical practices. The conceptual foundations of the most commonly applied modeling frameworks can be confusing even to those well versed in their application. Furthermore, there has yet to be a synthesis of the advances made over the last 20 yr. Therefore, there is a need for both synthesis of the current state of knowledge on habitat selection, and guidance for those seeking to study this process. Here, we provide an approachable overview and synthesis of the literature on habitat-selection analyses (HSAs) conducted using selection functions, which are by far the most applied modeling framework for understanding the habitat-selection process. This review is purposefully non-technical and focused on understanding without heavy mathematical and statistical notation, which can confuse many practitioners. We offer an overview and history of HSAs, describing the tortuous conceptual path to our current understanding. Through this overview, we also aim to address the areas of greatest confusion in the literature. We synthesize the literature outlining the most exciting conceptual advances in the field of habitat-selection modeling, discussing the substantial ecological and evolutionary inference that can be made using contemporary techniques. We aim for this paper to provide clarity for those navigating the complex literature on HSAs while acting as a reference and best practices guide for practitioners.

Characterizing the demographic history and prion protein variation to infer susceptibility to chronic wasting disease in a naïve population of white-tailed deer (Odocoileus virginianus).

Assessments of the adaptive potential in natural populations are essential for understanding and predicting responses to environmental stressors like climate change and infectious disease. Species face a range of stressors in human-dominated landscapes, often with contrasting effects. White-tailed deer (Odocoileus virginianus; deer) are expanding in the northern part of their range following decreasing winter severity and increasing forage availability. Chronic wasting disease (CWD), a prion disease affecting deer, is likewise expanding and represents a major threat to deer and other cervids. We obtained tissue samples from free-ranging deer across their native range in Ontario, Canada, which has yet to detect CWD in wild populations. We used high-throughput sequencing to assess neutral genomic variation and variation in the prion protein gene (PRNP) that is partly responsible for the protein misfolding when deer contract CWD. Neutral variation revealed a high number of rare alleles and no population structure, and demographic models suggested a rapid historical population expansion. Allele frequencies of PRNP variants associated with CWD susceptibility and disease progression were evenly distributed across the landscape and consistent with deer populations not infected with CWD. We estimated the selection coefficient of CWD, with simulations showing an observable and rapid shift in PRNP allele frequencies that coincides with the start of a novel CWD outbreak. Sustained surveillance of genomic and PRNP variation can be a useful tool for guiding management practices, which is especially important for CWD-free regions where deer are managed for ecological and economic benefits.

When are hypotheses useful in ecology and evolution?

Research hypotheses have been a cornerstone of science since before Galileo. Many have argued that hypotheses (1) encourage discovery of mechanisms, and (2) reduce bias-both features that should increase transferability and reproducibility. However, we are entering a new era of big data and highly predictive models where some argue the hypothesis is outmoded. We hypothesized that hypothesis use has declined in ecology and evolution since the 1990s, given the substantial advancement of tools further facilitating descriptive, correlative research. Alternatively, hypothesis use may have become more frequent due to the strong recommendation by some journals and funding agencies that submissions have hypothesis statements. Using a detailed literature analysis (N = 268 articles), we found prevalence of hypotheses in eco-evo research is very low (6.7%-26%) and static from 1990-2015, a pattern mirrored in an extensive literature search (N = 302,558 articles). Our literature review also indicates that neither grant success nor citation rates were related to the inclusion of hypotheses, which may provide disincentive for hypothesis formulation. Here, we review common justifications for avoiding hypotheses and present new arguments based on benefits to the individual researcher. We argue that stating multiple alternative hypotheses increases research clarity and precision, and is more likely to address the mechanisms for observed patterns in nature. Although hypotheses are not always necessary, we expect their continued and increased use will help our fields move toward greater understanding, reproducibility, prediction, and effective conservation of nature.

Heritability Estimates of Antler and Body Traits in White-Tailed Deer (Odocoileus virginianus) From Genomic-Relatedness Matrices.

Estimating heritability (h2) is required to predict the response to selection and is useful in species that are managed or farmed using trait information. Estimating h2 in free-ranging populations is challenging due to the need for pedigrees; genomic-relatedness matrices (GRMs) circumvent this need and can be implemented in nearly any system where phenotypic and genome-wide single-nucleotide polymorphism (SNP) data are available. We estimated the heritability of 5 body and 3 antler traits in a free-ranging population of white-tailed deer (Odocoileus virginianus) on Anticosti Island, Quebec, Canada. We generated classic and robust GRMs from >10,000 SNPs: hind foot length, dressed body mass, and peroneus muscle mass had high h2 values of 0.62, 0.44, and 0.55, respectively. Heritability in male-only antler features ranged from 0.07 to 0.33. We explored the influence of filtering by minor allele frequency and data completion on h2: GRMs derived from fewer SNPs had reduced h2 estimates and the relatedness coefficients significantly deviated from those generated with more SNPs. As a corollary, we discussed limitations to the application of GRMs in the wild, notably how skewed GRMs, specifically many unrelated individuals, can increase variance around h2 estimates. This is the first study to estimate h2 on a free-ranging population of white-tailed deer and should be informative for breeding designs and management as these traits could respond to selection.

Improving spatial predictions of animal resource selection to guide conservation decision making.

Resource selection is often studied by ecologists interested in the environmental drivers of animal space use and movement. These studies commonly produce spatial predictions, which are of considerable utility to resource managers making habitat and population management decisions. It is thus paramount that predictions from resource selection studies are accurate. We evaluated model building and fitting strategies for optimizing resource selection function predictions in a use-availability framework. We did so by simulating low- and high-intensity spatial sampling data that respectively predicted study area and movement-based resource selection. We compared one of the most commonly used forms of statistical regularization, Akaike's Information Criterion (AIC), with the lesser used least absolute shrinkage and selection operator (LASSO). LASSO predictions were less variable and more accurate than AIC and were often best when considering additive and interacting variables. We explicitly demonstrate the predictive equivalence using the logistic and Poisson likelihoods and how it is lost when the available sample is too small. Regardless of modeling approach, interpreting the sign of coefficients as a measure of selection can be misleading when optimizing for prediction.

Behavioural valuation of landscapes using movement data.

Wildlife tracking is one of the most frequently employed approaches to monitor and study wildlife populations. To date, the application of tracking data to applied objectives has focused largely on the intensity of use by an animal in a location or the type of habitat. While this has provided valuable insights and advanced spatial wildlife management, such interpretation of tracking data does not capture the complexity of spatio-temporal processes inherent to animal behaviour and represented in the movement path. Here, we discuss current and emerging approaches to estimate the behavioural value of spatial locations using movement data, focusing on the nexus of conservation behaviour and movement ecology that can amplify the application of animal tracking research to contemporary conservation challenges. We highlight the importance of applying behavioural ecological approaches to the analysis of tracking data and discuss the utility of comparative approaches, optimization theory and economic valuation to gain understanding of movement strategies and gauge population-level processes. First, we discuss innovations in the most fundamental movement-based valuation of landscapes, the intensity of use of a location, namely dissecting temporal dynamics in and means by which to weight the intensity of use. We then expand our discussion to three less common currencies for behavioural valuation of landscapes, namely the assessment of the functional (i.e. what an individual is doing at a location), structural (i.e. how a location relates to use of the broader landscape) and fitness (i.e. the return from using a location) value of a location. Strengthening the behavioural theoretical underpinnings of movement ecology research promises to provide a deeper, mechanistic understanding of animal movement that can lead to unprecedented insights into the interaction between landscapes and animal behaviour and advance the application of movement research to conservation challenges. This article is part of the theme issue 'Linking behaviour to dynamics of populations and communities: application of novel approaches in behavioural ecology to conservation'.

Impacts of the Northwest Forest Plan on forest composition and bird populations.

The Northwest Forest Plan (NWFP) initiated one of the most sweeping changes to forest management in the world, affecting 10 million hectares of federal land. The NWFP is a science-based plan incorporating monitoring and adaptive management and provides a unique opportunity to evaluate the influence of policy. We used >25 years of region-wide bird surveys, forest data, and land-ownership maps to test this policy's effect on biodiversity. Clearcutting decreased rapidly, and we expected populations of older-forest-associated birds to stabilize on federal land, but to continue declining on private industrial lands where clearcutting continued. In contrast, we expected declines in early-seral-associated species on federal land because of reduced anthropogenic disturbance since the NWFP. Bayesian hierarchical models revealed that bird species' population trends tracked changes in forest composition. However, against our expectations, declines of birds associated with older forests accelerated. These declines are partly explained by losses of older forests due to fire on federal land and continued clearcutting elsewhere. Indeed, the NWFP anticipated that reversing declines of older forests would take time. Overall, the early-seral ecosystem area was stable, but declined in two ecoregions-the Coast Range and Cascades-along with early-seral bird populations. Although the NWFP halted clearcutting on federal land, this has so far been insufficient to reverse declines in older-forest-associated bird populations. These findings underscore the importance of continuing to prioritize older forests under the NWFP and ensuring that the recently proposed creation of early-seral ecosystems does not impede the conservation and development of older-forest structure.

Correction: A comment on priors for Bayesian occupancy models.

[This corrects the article DOI: 10.1371/journal.pone.0192819.].

Synergistic effects of climate and land-use change influence broad-scale avian population declines.

Climate and land-use changes are expected to be the primary drivers of future global biodiversity loss. Although theory suggests that these factors impact species synergistically, past studies have either focused on only one in isolation or have substituted space for time, which often results in confounding between drivers. Tests of synergistic effects require congruent time series on animal populations, climate change and land-use change replicated across landscapes that span the gradient of correlations between the drivers of change. Using a unique time series of high-resolution climate (measured as temperature and precipitation) and land-use change (measured as forest change) data, we show that these drivers of global change act synergistically to influence forest bird population declines over 29 years in the Pacific Northwest of the United States. Nearly half of the species examined had declined over this time. Populations declined most in response to loss of early seral and mature forest, with responses to loss of early seral forest amplified in landscapes that had warmed over time. In addition, birds declined more in response to loss of mature forest in areas that had dried over time. Climate change did not appear to impact populations in landscapes with limited habitat loss, except when those landscapes were initially warmer than the average landscape. Our results provide some of the first empirical evidence of synergistic effects of climate and land-use change on animal population dynamics, suggesting accelerated loss of biodiversity in areas under pressure from multiple global change drivers. Furthermore, our findings suggest strong spatial variability in the impacts of climate change and highlight the need for future studies to evaluate multiple drivers simultaneously to avoid potential misattribution of effects.

A comment on priors for Bayesian occupancy models.

Understanding patterns of species occurrence and the processes underlying these patterns is fundamental to the study of ecology. One of the more commonly used approaches to investigate species occurrence patterns is occupancy modeling, which can account for imperfect detection of a species during surveys. In recent years, there has been a proliferation of Bayesian modeling in ecology, which includes fitting Bayesian occupancy models. The Bayesian framework is appealing to ecologists for many reasons, including the ability to incorporate prior information through the specification of prior distributions on parameters. While ecologists almost exclusively intend to choose priors so that they are "uninformative" or "vague", such priors can easily be unintentionally highly informative. Here we report on how the specification of a "vague" normally distributed (i.e., Gaussian) prior on coefficients in Bayesian occupancy models can unintentionally influence parameter estimation. Using both simulated data and empirical examples, we illustrate how this issue likely compromises inference about species-habitat relationships. While the extent to which these informative priors influence inference depends on the data set, researchers fitting Bayesian occupancy models should conduct sensitivity analyses to ensure intended inference, or employ less commonly used priors that are less informative (e.g., logistic or t prior distributions). We provide suggestions for addressing this issue in occupancy studies, and an online tool for exploring this issue under different contexts.

Applying network theory to animal movements to identify properties of landscape space use.

Network (graph) theory is a popular analytical framework to characterize the structure and dynamics among discrete objects and is particularly effective at identifying critical hubs and patterns of connectivity. The identification of such attributes is a fundamental objective of animal movement research, yet network theory has rarely been applied directly to animal relocation data. We develop an approach that allows the analysis of movement data using network theory by defining occupied pixels as nodes and connection among these pixels as edges. We first quantify node-level (local) metrics and graph-level (system) metrics on simulated movement trajectories to assess the ability of these metrics to pull out known properties in movement paths. We then apply our framework to empirical data from African elephants (Loxodonta africana), giant Galapagos tortoises (Chelonoidis spp.), and mule deer (Odocoileous hemionus). Our results indicate that certain node-level metrics, namely degree, weight, and betweenness, perform well in capturing local patterns of space use, such as the definition of core areas and paths used for inter-patch movement. These metrics were generally applicable across data sets, indicating their robustness to assumptions structuring analysis or strategies of movement. Other metrics capture local patterns effectively, but were sensitive to specified graph properties, indicating case specific applications. Our analysis indicates that graph-level metrics are unlikely to outperform other approaches for the categorization of general movement strategies (central place foraging, migration, nomadism). By identifying critical nodes, our approach provides a robust quantitative framework to identify local properties of space use that can be used to evaluate the effect of the loss of specific nodes on range wide connectivity. Our network approach is intuitive, and can be implemented across imperfectly sampled or large-scale data sets efficiently, providing a framework for conservationists to analyze movement data. Functions created for the analyses are available within the R package moveNT.

Movement reveals scale dependence in habitat selection of a large ungulate.

Ecological processes operate across temporal and spatial scales. Anthropogenic disturbances impact these processes, but examinations of scale dependence in impacts are infrequent. Such examinations can provide important insight to wildlife-human interactions and guide management efforts to reduce impacts. We assessed spatiotemporal scale dependence in habitat selection of mule deer (Odocoileus hemionus) in the Piceance Basin of Colorado, USA, an area of ongoing natural gas development. We employed a newly developed animal movement method to assess habitat selection across scales defined using animal-centric spatiotemporal definitions ranging from the local (defined from five hour movements) to the broad (defined from weekly movements). We extended our analysis to examine variation in scale dependence between night and day and assess functional responses in habitat selection patterns relative to the density of anthropogenic features. Mule deer displayed scale invariance in the direction of their response to energy development features, avoiding well pads and the areas closest to roads at all scales, though with increasing strength of avoidance at coarser scales. Deer displayed scale-dependent responses to most other habitat features, including land cover type and habitat edges. Selection differed between night and day at the finest scales, but homogenized as scale increased. Deer displayed functional responses to development, with deer inhabiting the least developed ranges more strongly avoiding development relative to those with more development in their ranges. Energy development was a primary driver of habitat selection patterns in mule deer, structuring their behaviors across all scales examined. Stronger avoidance at coarser scales suggests that deer behaviorally mediated their interaction with development, but only to a degree. At higher development densities than seen in this area, such mediation may not be possible and thus maintenance of sufficient habitat with lower development densities will be a critical best management practice as development expands globally.