Behavioural state‐dependent habitat selection and implications for animal translocations

Associations between glucocorticoids and habitat selection reflect daily and seasonal energy requirements

BACKGROUND

Glucocorticoids are often associated with stressful environments, but they are also thought to drive the best strategies to improve fitness in stressful environments. Glucocorticoids improve fitness in part by regulating foraging behaviours in response to daily and seasonal energy requirements. However, many studies demonstrating relationships between foraging behaviour and glucocorticoids are experimental, and few observational studies conducted under natural conditions have tested whether changing glucocorticoid levels are related to daily and seasonal changes in energy requirements.

METHODS

We integrated glucocorticoids into habitat selection models to test for relationships between foraging behaviour and glucocorticoid levels in elk (Cervus canadensis) as their daily and seasonal energy requirements changed. Using integrated step selection analysis, we tested whether elevated glucocorticoid levels were related to foraging habitat selection on a daily scale and whether that relationship became stronger during lactation, one of the greatest seasonal periods of energy requirement for female mammals.

RESULTS

We found stronger selection of foraging habitat by female elk with elevated glucocorticoids (eß = 1.44 95% CI 1.01, 2.04). We found no difference in overall glucocorticoid levels after calving, nor a significant change in the relationship between glucocorticoids and foraging habitat selection at the time of calving. However, we found a gradual increase in the relationship between glucocorticoids and habitat selection by female elk as their calves grew over the next few months (eß = 1.01, 95% CI 1.00, 1.02), suggesting a potentially stronger physiological effect of glucocorticoids for elk with increasing energy requirements.

CONCLUSIONS

We suggest glucocorticoid-integrated habitat selection models demonstrate the role of glucocorticoids in regulating foraging responses to daily and seasonal energy requirements. Ultimately, this integration will help elucidate the implications of elevated glucocorticoids under natural conditions.

The potential influence of genome-wide adaptive divergence on conservation translocation outcome in an isolated greater sage-grouse population

Conservation translocations are an important conservation tool commonly employed to augment declining or reestablish extirpated populations. One goal of augmentation is to increase genetic diversity and reduce the risk of inbreeding depression (i.e., genetic rescue). However, introducing individuals from significantly diverged populations risks disrupting coadapted traits and reducing local fitness (i.e., outbreeding depression). Genetic data are increasingly more accessible for wildlife species and can provide unique insight regarding the presence and retention of introduced genetic variation from augmentation as an indicator of effectiveness and adaptive similarity as an indicator of source and recipient population suitability. We used 2 genetic data sets to evaluate augmentation of isolated populations of greater sage-grouse (Centrocercus urophasianus) in the northwestern region of the species range (Washington, USA) and to retrospectively evaluate adaptive divergence among source and recipient populations. We developed 2 statistical models for microsatellite data to evaluate augmentation outcomes. We used one model to predict genetic diversity after augmentation and compared these predictions with observations of genetic change. We used the second model to quantify the amount of observed reproduction attributed to transplants (proof of population integration). We also characterized genome-wide adaptive divergence among source and recipient populations. Observed genetic diversity (HO = 0.65) was higher in the recipient population than predicted had no augmentation occurred (HO = 0.58) but less than what was predicted by our model (HO = 0.75). The amount of shared genetic variation between the 2 geographically isolated resident populations increased, which is evidence of periodic gene flow previously assumed to be rare. Among candidate adaptive genes associated with elevated fixation index (FST) (143 genes) or local environmental variables (97 and 157 genes for each genotype-environment association method, respectively), we found clusters of genes with related functions that may influence the ability of transplants to use local resources and navigate unfamiliar environments and their reproductive potential, all possible reasons for low genetic retention from augmentation.

How to account for behavioral states in step-selection analysis: a model comparison

Step-selection models are widely used to study animals' fine-scale habitat selection based on movement data. Resource preferences and movement patterns, however, often depend on the animal's unobserved behavioral states, such as resting or foraging. As this is ignored in standard (integrated) step-selection analyses (SSA, iSSA), different approaches have emerged to account for such states in the analysis. The performance of these approaches and the consequences of ignoring the states in step-selection analysis, however, have rarely been quantified. We evaluate the recent idea of combining iSSAs with hidden Markov models (HMMs), which allows for a joint estimation of the unobserved behavioral states and the associated state-dependent habitat selection. Besides theoretical considerations, we use an extensive simulation study and a case study on fine-scale interactions of simultaneously tracked bank voles (Myodes glareolus) to compare this HMM-iSSA empirically to both the standard and a widely used classification-based iSSA (i.e., a two-step approach based on a separate prior state classification). Moreover, to facilitate its use, we implemented the basic HMM-iSSA approach in the R package HMMiSSA available on GitHub.

Socioecological factors influencing intraspecific variation in ranging dynamics of western lowland gorillas (Gorilla gorilla gorilla) in Ndoki Forest

Ranging dynamics are physical and behavioral representations of how different socioecological factors affect an organism's spatial decisions and space use strategies. Western lowland gorillas (Gorilla gorilla gorilla) are a model species to investigate the drivers of spatial dynamics based on both the natural variation in socioecological factors within the species and compared with their mountain gorilla counterparts. In this study, we evaluate the influences of resource seasonality and social dynamics on variation in home range size, utilization, and intergroup overlap among multiple gorilla groups over an 8-year study period in the northern Republic of Congo. This study shows that western lowland gorillas can have small home ranges comparable to mountain gorillas, rather than universally larger home ranges as previously supposed, and that home ranges are stable through time. The largest source of variation in space use was the degree of intergroup home range overlap. The study groups did not demonstrate intraspecific variation in range size nor changes in intergroup overlap with respect to seasonality of fruit resources, but all groups demonstrated expansion of monthly range and core area with group size, matching predictions of intragroup feeding competition. These findings highlight the potential impact of intergroup relationships on space use and prompt further research on the role of social dynamics in ranging strategies. In this study, we reveal a greater degree of variability and flexibility in gorilla ranging behavior than previously realized which is relevant to improving comparative studies and informing conservation strategies on behalf of these endangered primates.

A multistate Langevin diffusion for inferring behavior-specific habitat selection and utilization distributions

The identification of important habitat and the behavior(s) associated with it is critical to conservation and place-based management decisions. Behavior also links life-history requirements and habitat use, which are key to understanding why animals use certain habitats. Animal population studies often use tracking data to quantify space use and habitat selection, but they typically either ignore movement behavior (e.g., foraging, migrating, nesting) or adopt a two-stage approach that can induce bias and fail to propagate uncertainty. We develop a habitat-driven Langevin diffusion for animals that exhibit distinct movement behavior states, thereby providing a novel single-stage statistical method for inferring behavior-specific habitat selection and utilization distributions in continuous time. Practitioners can customize, fit, assess, and simulate our integrated model using the provided R package. Simulation experiments demonstrated that the model worked well under a range of sampling scenarios as long as observations were of sufficient temporal resolution. Our simulations also demonstrated the importance of accounting for different behaviors and the misleading inferences that can result when these are ignored. We provide case studies using plains zebra (Equus quagga) and Steller sea lion (Eumetopias jubatus) telemetry data. In the zebra example, our model identified distinct "encamped" and "exploratory" states, where the encamped state was characterized by strong selection for grassland and avoidance of other vegetation types, which may represent selection for foraging resources. In the sea lion example, our model identified distinct movement behavior modes typically associated with this marine central-place forager and, unlike previous analyses, found foraging-type movements to be associated with steeper offshore slopes characteristic of the continental shelf, submarine canyons, and seamounts that are believed to enhance prey concentrations. This is the first single-stage approach for inferring behavior-specific habitat selection and utilization distributions from tracking data that can be readily implemented with user-friendly software. As certain behaviors are often more relevant to specific conservation or management objectives, practitioners can use our model to help inform the identification and prioritization of important habitats. Moreover, by linking individual-level movement behaviors to population-level spatial processes, the multistate Langevin diffusion can advance inferences at the intersection of population, movement, and landscape ecology.

Behavioral-dependent recursive movements and implications for resource selection

Within home ranges, animals repeatedly visit certain areas. Recursive movement patterns are widespread throughout the animal kingdom, but are rarely considered when developing resource selection models. We examined how behavioral state-dependent recursive movements influenced reource selection of eastern wild turkey (Meleagris gallopavo silvestris) broods as they aged from day 1 to 28. Because broods become more plastic in behaviors once they begin roosting off the ground, we separated data into broods that were ground roosting (1-13 days) and tree roosting (14-28 days). We used Hidden Markov Models to identify 2 behavioral states (restricted and mobile). We extracted state-specific recursive movements based on states and specific step lengths, which we integrated into a step selection analysis to evaluate resource selection. We found that in a restricted state, ground roosting broods spent less time in areas of mixed pine-hardwoods and more time in areas with greater vegetation density. Tree roosting broods revisited areas closer to shrub/scrub landcover types, and areas with greater vegetation density. Tree roosting broods also spent less time near mixed pine-hardwoods, while spending more time in areas with greater vegetation density. We found that in a mobile state, ground roosting broods revisited areas closer to secondary roads and mixed pine-hardwoods, but farther from hardwoods. Tree roosting broods revisited areas farther from secondary roads and with greater vegetation density. Tree roosting broods also spent more time in areas closer to pine. Resource selection varied depending on behavioral state and recursive movements. However, revisitation and residence time impacted selection in both ground and tree roosting broods. Our findings highlight the need to consider how behaviors can influence movement decisions and ultimately resource selection.

Individual variation in the habitat selection of upstream migrating fish near a barrier

BACKGROUND

Migration is a vital element of the life cycle of many freshwater fish species but is increasingly hampered globally by riverine barriers. Fish passes are a common approach to enable migration past barriers but are often ineffective. More knowledge is required on fish behaviour as they approach barriers such as habitat preferences.

METHODS

We evaluate the habitat selection of two upstream migrating fish species, barbel Barbus barbus and grayling Thymallus thymallus, at a hydropower plant in southern Germany, considering individual variation and population trends. Fish were tracked via fine-scale 2D acoustic telemetry in 2018 during their spawning migration. Step selection functions were used to evaluate selection of hydraulic parameters by the fish for a time step of 20 s. Exploratory models were built via model selection for each individual fish, to evaluate the extent of individual variation in model structure. A population model was developed for each species by averaging coefficients from individual models to describe general trends. The extent of individual variation was determined and confidence intervals for the population model coefficients were calculated.

RESULTS

Fish varied greatly in individual model structure though common terms were apparent in both species, such as depth, flow velocity, the angular difference between fish and velocity, and the logarithm of the step length. Final population models for barbel included several parameters describing habitat selection and displacement. Barbel selected for faster flows, deeper water, and higher spatial velocity gradients. In addition, they selected to move more with the flow than against. Interactions were also present between habitat parameters, suggesting selection is context dependent. Barbel movement speed also changed with depth, flow velocity and spatial velocity gradient. With grayling, terms often had contrasting effects among individuals and thus general trends could not be distinguished for most terms.

CONCLUSION

Our findings demonstrate habitat selection by upstream migrating fish approaching a fish pass and differences in individual selection which may have an impact on barrier management. Step selection functions are a promising approach and can provide useful insight into habitat selection and movement by migrating freshwater fish in an altered river system.

Flexible hidden Markov models for behaviour-dependent habitat selection

BACKGROUND

There is strong incentive to model behaviour-dependent habitat selection, as this can help delineate critical habitats for important life processes and reduce bias in model parameters. For this purpose, a two-stage modelling approach is often taken: (i) classify behaviours with a hidden Markov model (HMM), and (ii) fit a step selection function (SSF) to each subset of data. However, this approach does not properly account for the uncertainty in behavioural classification, nor does it allow states to depend on habitat selection. An alternative approach is to estimate both state switching and habitat selection in a single, integrated model called an HMM-SSF.

METHODS

We build on this recent methodological work to make the HMM-SSF approach more efficient and general. We focus on writing the model as an HMM where the observation process is defined by an SSF, such that well-known inferential methods for HMMs can be used directly for parameter estimation and state classification. We extend the model to include covariates on the HMM transition probabilities, allowing for inferences into the temporal and individual-specific drivers of state switching. We demonstrate the method through an illustrative example of plains zebra (Equus quagga), including state estimation, and simulations to estimate a utilisation distribution.

RESULTS

In the zebra analysis, we identified two behavioural states, with clearly distinct patterns of movement and habitat selection ("encamped" and "exploratory"). In particular, although the zebra tended to prefer areas higher in grassland across both behavioural states, this selection was much stronger in the fast, directed exploratory state. We also found a clear diel cycle in behaviour, which indicated that zebras were more likely to be exploring in the morning and encamped in the evening.

CONCLUSIONS

This method can be used to analyse behaviour-specific habitat selection in a wide range of species and systems. A large suite of statistical extensions and tools developed for HMMs and SSFs can be applied directly to this integrated model, making it a very versatile framework to jointly learn about animal behaviour, habitat selection, and space use.

Biologging as an important tool to uncover behaviors of cryptic species: an analysis of giant armadillos (Priodontes maximus)

Advances in biologging have increased the understanding of how animals interact with their environment, especially for cryptic species. For example, giant armadillos (Priodontes maximus) are the largest extant species of armadillo but are rarely encountered due to their fossorial and nocturnal behavior. Through the analysis of speed, turning angles, and accelerometer activity counts, we estimated behavioral states, characterized activity budgets, and investigated the state-habitat associations exhibited by individuals monitored with GPS telemetry in the Brazilian Pantanal from 2019 to 2020. This methodology is proposed as a useful framework for the identification of priority habitat. Using the non-parametric Bayesian mixture model for movement (M3), we estimated four latent behavioral states that were named 'vigilance-excavation', 'local search', 'exploratory', and 'transit'. These states appeared to correspond with behavior near burrows or termite mounds, foraging, ranging, and rapid movements, respectively. The first and last hours of activity presented relatively high proportions of the vigilance-excavation state, while most of the activity period was dominated by local search and exploratory states. The vigilance-excavation state occurred more frequently in regions between forest and closed savannas, whereas local search was more likely in high proportions of closed savanna. Exploratory behavior probability increased in areas with high proportions of both forest and closed savanna. Our results establish a baseline for behavioral complexity, activity budgets, and habitat associations in a relatively pristine environment that can be used for future work to investigate anthropogenic impacts on giant armadillo behavior and fitness. The integration of accelerometer and GPS-derived movement data through our mixture model has the potential to become a powerful methodological approach for the conservation of other cryptic species.

Drivers of polar bear behavior and the possible effects of prey availability on foraging strategy

BACKGROUND

Change in behavior is one of the earliest responses to variation in habitat suitability. It is therefore important to understand the conditions that promote different behaviors, particularly in areas undergoing environmental change. Animal movement is tightly linked to behavior and remote tracking can be used to study ethology when direct observation is not possible.

METHODS

We used movement data from 14 polar bears (Ursus maritimus) in Hudson Bay, Canada, during the foraging season (January-June), when bears inhabit the sea ice. We developed an error-tolerant method to correct for sea ice drift in tracking data. Next, we used hidden Markov models with movement and orientation relative to wind to study three behaviors (stationary, area-restricted search, and olfactory search) and examine effects of 11 covariates on behavior.

RESULTS

Polar bears spent approximately 47% of their time in the stationary drift state, 29% in olfactory search, and 24% in area-restricted search. High energy behaviors occurred later in the day (around 20:00) compared to other populations. Second, olfactory search increased as the season progressed, which may reflect a shift in foraging strategy from still-hunting to active search linked to a shift in seal availability (i.e., increase in haul-outs from winter to the spring pupping and molting seasons). Last, we found spatial patterns of distribution linked to season, ice concentration, and bear age that may be tied to habitat quality and competitive exclusion.

CONCLUSIONS

Our observations were generally consistent with predictions of the marginal value theorem, and differences between our findings and other populations could be explained by regional or temporal variation in resource availability. Our novel movement analyses and finding can help identify periods, regions, and conditions of critical habitat.

Sex Differences Dictate the Movement Patterns of Striped Hyenas, Hyaena hyaena, in a Human-Dominated Landscape

Large-carnivore populations have experienced significant declines in the past centuries in extended parts of the world. Habitat loss, fragmentation, and depletion of natural resources are some of the main causes of this decline. Consequently, behavioral flexibility, enabling the exploitation of anthropogenic food resources in highly disturbed human-dominated landscapes, is becoming critical for the survival of large carnivores. These behavioral changes increase the potential for human-large carnivore conflict and can further intensify carnivore persecution. Here, we examine how land cover types (representing a gradient of anthropogenic disturbance) alter the behavior of striped hyenas (Hyaena hyaena) in a human-dominated landscape in Israel, and whether differences in life history between males and females affect their reaction to such disturbances and consequently their level of exposure to humans. We used a Hidden Markov Model on GPS-tracking data for seven striped hyenas to segment individual-night trajectories into behavioral states (resting, searching, and traveling). We then used multinomial logistic regression to model hyenas’ behavioral state as a function of the interaction between land cover and sex. Females traveled less than males both in terms of average distance traveled per hour, per night, and nightly net displacement. Most steps were classified as “searching” for females and as “traveling” for males. Female hyenas spent a higher proportion of time in human-dominated areas and a lower proportion in natural areas compared to males, and they were also more likely to be found close to settlements than males. Females changed their time allocation between natural and human-dominated areas, spending more time resting than traveling in natural areas but not in human-dominated ones; males spent more time searching than resting in human-dominated areas but were equally likely to rest or search in natural ones. The differences in life history between male and female hyenas may reflect different motivations for space use as a means to optimize fitness, which affects their exposure to humans and therefore their potential involvement in human-hyenas conflict. Understanding the mechanisms that lead to behavioral change in response to human disturbance is important for adaptive management and promoting human large-carnivores co-existence in general.