Robert A. Montgomery

Implications of Ignoring Telemetry Error on Inference in Wildlife Resource use Models

ABSTRACT Global Positioning System (GPS) and very high frequency (VHF) telemetry data redefined the examination of wildlife resource use. Researchers collar animals, relocate those animals over time, and utilize the estimated locations to infer resource use and build predictive models. Precision of these estimated wildlife locations, however, influences the reliability of point-based models with accuracy depending on the interaction between mean telemetry error and how habitat characteristics are mapped (categorical raster resolution and patch size). Telemetry data often foster the assumption that locational error can be ignored without biasing study results. We evaluated the effects of mean telemetry error and categorical raster resolution on the correct characterization of patch use when locational error is ignored. We found that our ability to accurately attribute patch type to an estimated telemetry location improved nonlinearly as patch size increased and mean telemetry error decreased. Furthermore, the exact shape of these relationships was directly influenced by categorical raster resolution. Accuracy ranged from 100% (200-ha patch size, 1- to 5-m telemetry error) to 46% (0.5-ha patch size, 56- to 60-m telemetry error) for 10 m resolution rasters. Accuracy ranged from 99% (200-ha patch size, 1- to 5-m telemetry error) to 57% (0.5-ha patch size, 56- to 60-m telemetry error) for 30-m resolution rasters. When covariate rasters were less resolute (30 m vs. 10 m) estimates for the ignore technique were more accurate at smaller patch sizes. Hence, both fine resolution (10 m) covariate rasters and small patch sizes increased probability of patch misidentification. Our results help frame the scope of ecological inference made from point-based wildlife resource use models. For instance, to make ecological inferences with 90% accuracy at small patch sizes (≤5 ha) mean telemetry error ≤5 m is required for 10-m resolution categorical rasters. To achieve the same inference on 30-m resolution categorical rasters, mean telemetry error ≤10 m is required. We encourage wildlife professionals creating point-based models to assess whether reasonable estimates of resource use can be expected given their telemetry error, covariate raster resolution, and range of patch sizes.

Importance of visibility when evaluating animal response to roads

Roads increase risk to animals via direct and indirect mechanisms yet, both positive and negative effects of animal space use in relation to roads have been reported. These contrasting reports may not actually represent animal ecology, but could be a product of the primary variable used to test the relationship between animals and roads. Animal-road associations are often evaluated using Euclidean distance. Euclidean, or straight-line, distance fails to account for the screening effects of vegetation and topography and may document spurious relationships. We evaluated the influence of Euclidean distance, visibility from road and forage quality on summer space use for male elk Cervus elaphus and female elk subherds in Custer State Park, South Dakota, USA. Models that included interactions with visibility from road metrics outperformed models that included only Euclidean distance to road as main effects. Elk response to roads varied by sex and road type, which functioned as an index for vehicle use. Male elk selected habitat away from roads with the greatest vehicle use, an effect that was greater if habitat was visible from those roads. Female elk tended to select habitat with high forage quality in areas visible from roads closed to vehicle use. Interestingly, both male and female elk selected habitat away from roads with secondary vehicle use and near to roads devoid of traffic, regardless of visibility. Our analysis highlights the importance of including both Euclidean distance and visibility from road metrics. Road effects research may be incomplete without consideration of visibility from roads, particularly for large mammals in landscapes with intense road use.

Accuracy of Techniques Used to Assign Mourning Dove Age and Gender

Sex ratios of mourning doves (Zenaida macroura) reported in banding studies have differed markedly between states and years suggesting that techniques used to assign gender may be biased. We tested this hypothesis with experienced mourning dove biologists (n = 4) who assigned age, gender, primary molt progression, and certainty of classification, using external characteristics of harvested doves (n = 468) in Missouri during September 1993. The probability of correctly classifying adult mourning doves, verified with gonadal examination, was 0.904-0.985 and 0.957-0.990 for immatures. Assessment of age and gender combined was 0.944-0.990 for adult females and 0.533-0.993 for adult males. Error probabilities of observer ability to correctly assign age and gender was 0.0-0.045 for adult females and 0.007-0.451 for adult males. Among individual mourning doves examined, only 1 adult and 6 immature doves were incorrectly assigned age and gender by all 4 participants. Error in assigning gender to adult doves could have contributed to disparate sex ratios reported in previous mourning dove banding studies. We recommend that standard training procedures be developed and implemented across states and regions in future banding studies of mourning doves.

Can We Accurately Characterize Wildlife Resource Use When Telemetry Data Are Imprecise

Abstract Telemetry data have been widely used to quantify wildlife habitat relationships despite the fact that these data are inherently imprecise. All telemetry data have positional error, and failure to account for that error can lead to incorrect predictions of wildlife resource use. Several techniques have been used to account for positional error in wildlife studies. These techniques have been described in the literature, but their ability to accurately characterize wildlife resource use has never been tested. We evaluated the performance of techniques commonly used for incorporating telemetry error into studies of wildlife resource use. Our evaluation was based on imprecise telemetry data (mean telemetry error  =  174 m, SD  =  130 m) typical of field-based studies. We tested 5 techniques in 10 virtual environments and in one real-world environment for categorical (i.e., habitat types) and continuous (i.e., distances or elevations) rasters. Technique accuracy varied by patch size for the categorical rasters, with higher accuracy as patch size increased. At the smallest patch size (1 ha), the technique that ignores error performed best on categorical data (0.31 and 0.30 accuracy for virtual and real data, respectively); however, as patch size increased the bivariate-weighted technique performed better (0.56 accuracy at patch sizes >31 ha) and achieved complete accuracy (i.e., 1.00 accuracy) at smaller patch sizes (472 ha and 1,522 ha for virtual and real data, respectively) than any other technique. We quantified the accuracy of the continuous covariates using the mean absolute difference (MAD) in covariate value between true and estimated locations. We found that average MAD varied between 104 m (ignore telemetry error) and 140 m (rescale the covariate data) for our continuous covariate surfaces across virtual and real data sets. Techniques that rescale continuous covariate data or use a zonal mean on values within a telemetry error polygon were significantly less accurate than other techniques. Although the technique that ignored telemetry error performed best on categorical rasters with smaller average patch sizes (i.e., ≤31 ha) and on continuous rasters in our study, accuracy was so low that the utility of using point-based approaches for quantifying resource use is questionable when telemetry data are imprecise, particularly for small-patch habitat relationships.

The influence of winter severity, predation and senescence on moose habitat use

Habitat use is widely known to be influenced by abiotic and biotic factors, such as climate, population density, foraging opportunity and predation risk. The influence of the life-history state of an individual organism on habitat use is less well understood, especially for terrestrial mammals. There is good reason to expect that life-history state would affect habitat use. For example, organisms exhibiting poor condition associated with senescence have an increased vulnerability to predation and that vulnerability is known to alter habitat use strategies. We assessed the influence of life-history stage on habitat use for 732 moose (Alces alces) killed by wolves (Canis lupus) over a 50-year period in Isle Royale National Park, an island ecosystem in Lake Superior, USA. We developed regression models to assess how location of death was associated with a mooses life-history stage (prime-aged or senescent), presence or absence of senescent-associated pathology (osteoarthritis and jaw necrosis), and annual variation in winter severity, moose density and ratio of moose to wolves, which is an index of predation risk. Compared to senescent moose, prime-aged moose tend to make greater use of habitat farther from the shoreline of Isle Royale. That result is ecologically relevant because shoreline habitat on Isle Royale tends to provide better foraging opportunities for moose but is also associated with increased predation risk. During severe winters prime-aged moose tend to make greater use of habitat that is closer to shore in relation to senescent-aged moose. Furthermore, moose of both age classes were more likely to die in riskier, shoreline habitat during years when predation risk was lower in the preceding year. Our results highlight a complicated connection between life history, age-structured population dynamics and habitat-related behaviour. Our analysis also illustrates why intraspecific competition should not be the presumed mechanism underlying density-dependent habitat use, if predation risk is related to density, as it is expected to be in many systems.

Where wolves kill moose: the influence of prey life history dynamics on the landscape ecology of predation.

The landscape ecology of predation is well studied and known to be influenced by habitat heterogeneity. Little attention has been given to how the influence of habitat heterogeneity on the landscape ecology of predation might be modulated by life history dynamics of prey in mammalian systems. We demonstrate how life history dynamics of moose (Alces alces) contribute to landscape patterns in predation by wolves (Canis lupus) in Isle Royale National Park, Lake Superior, USA. We use pattern analysis and kernel density estimates of moose kill sites to demonstrate that moose in senescent condition and moose in prime condition tend to be wolf-killed in different regions of Isle Royale in winter. Predation on senescent moose was clustered in one kill zone in the northeast portion of the island, whereas predation on prime moose was clustered in 13 separate kill zones distributed throughout the full extent of the island. Moreover, the probability of kill occurrence for senescent moose, in comparison to prime moose, increased in high elevation habitat with patches of dense coniferous trees. These differences can be attributed, at least in part, to senescent moose being more vulnerable to predation and making different risk-sensitive habitat decisions than prime moose. Landscape patterns emerging from prey life history dynamics and habitat heterogeneity have been observed in the predation ecology of fish and insects, but this is the first mammalian system for which such observations have been made.

Moving GIS Research Indoors: Spatiotemporal Analysis of Agricultural Animals

A proof of concept applying wildlife ecology techniques to animal welfare science in intensive agricultural environments was conducted using non-cage laying hens. Studies of wildlife ecology regularly use Geographic Information Systems (GIS) to assess wild animal movement and behavior within environments with relatively unlimited space and finite resources. However, rather than depicting landscapes, a GIS could be developed in animal production environments to provide insight into animal behavior as an indicator of animal welfare. We developed a GIS-based approach for studying agricultural animal behavior in an environment with finite space and unlimited resources. Concurrent data from wireless body-worn location tracking sensor and video-recording systems, which depicted spatially-explicit behavior of hens (135 hens/room) in two identical indoor enclosures, were collected. The spatial configuration of specific hen behaviors, variation in home range patterns, and variation in home range overlap show that individual hens respond to the same environment differently. Such information could catalyze management practice adjustments (e.g., modifying feeder design and/or location). Genetically-similar hens exhibited diverse behavioral and spatial patterns via a proof of concept approach enabling detailed examinations of individual non-cage laying hen behavior and welfare.

Living amidst a sea of agriculture: predicting the occurrence of Canada lynx within an ecological island

Carnivore populations are increasingly confined to reserves surrounded by anthropogenic development. The boundaries of ecological islands are risky because of habitat loss and human—carnivore conflict. From snow track survey data collected over a 5-winter period we developed a population-level resource selection function for Canada lynx Lynx canadensis in Riding Mountain National Park, Canada. This park has been characterized as an ‘ecological island situated amidst a sea of agricultural land’ and while lynx are protected within the park, they are subject to harvest outside of the park. Winter resource selection of lynx increased with higher elevation and in highly suitable habitat for snowshoe hare Lepus americanus and decreased in habitat with greater proportions of agriculture and grasslands, both of which are common along the edge of the park. Habitat with medium to high relative probabilities of lynx occurrence tended to be distributed in the interior of the park. However, the highest relative probability of lynx occurrence was associated with habitat near the southeast border of the park in close proximity to a human community and a lake where snowmobiling, skiing, and snowshoeing are common recreational pursuits. We attribute this relationship to a 25-year old fire which presently created successional habitat highly suitable for snowshoe hare and correspondingly for lynx in this area. Our results suggest that Canada lynx occurrence tends to be associated with habitat that is highly suitable for their primary prey even if that habitat is located near to sources of human recreational activity.

Mapping the Relative Probability of Common Toad Occurrence in Terrestrial Lowland Farm Habitat in the United Kingdom

Introduction The common toad (Bufo bufo) is of increasing conservation concern in the United Kingdom (UK) due to dramatic population declines occurring in the past century. Many of these population declines coincided with reductions in both terrestrial and aquatic habitat availability and quality and have been primarily attributed to the effect of agricultural land conversion (of natural and semi-natural habitats to arable and pasture fields) and pond drainage. However, there is little evidence available to link habitat availability with common toad population declines, especially when examined at a broad landscape scale. Assessing such patterns of population declines at the landscape scale, for instance, require an understanding of how this species uses terrestrial habitat. Methods We intensively studied the terrestrial resource selection of a large population of common toads in Oxfordshire, England, UK. Adult common toads were fitted with passive integrated transponder (PIT) tags to allow detection in the terrestrial environment using a portable PIT antenna once toads left the pond and before going into hibernation (April/May-October 2012 and 2013). We developed a population-level resource selection function (RSF) to assess the relative probability of toad occurrence in the terrestrial environment by collecting location data for 90 recaptured toads. Results The predicted relative probability of toad occurrence for this population was greatest in wooded habitat near to water bodies; relative probability of occurrence declined dramatically > 50 m from these habitats. Toads also tended to select habitat near to their breeding pond and toad occurrence was negatively related to urban environments.

Examining disease prevalence for species of conservation concern using non‐invasive spatial capture–recapture techniques

Summary Non-invasive techniques have long been used to estimate wildlife population abundance and density. However, recent technological breakthroughs have facilitated non-invasive estimation of the proportion of animal populations with certain diseases. Giraffes Giraffa camelopardalis are increasingly becoming recognized as a species of conservation concern with decreasing population trajectories across their range in Africa. Diseases may be an important component impacting giraffe population declines, and the emerging ‘Giraffe Skin Disease’ (GSD), characterized by the appearance of wrinkled skin and alopecic lesions on the limbs, neck, and chest of infected giraffe, may hinder movement causing increased susceptibility to predation. We examined the prevalence of GSD in Tanzanias Ruaha National Park over a 4-month period in 2015, using photographic capture–recapture surveys via road-based transects. We divided the study area into five circuitous survey units, each approximately 100 km in length (