Running on empty: recharge dynamics from animal movement data

Unveiling the roles of temporal periodicity, the spatial environment and behavioural modes in terrestrial animal movement

BACKGROUND

Animal movement arises from complex interactions between animals and their heterogeneous environment. To better understand the movement process, it can be divided into behavioural, temporal and spatial components. Although methods exist to address those various components, it remains challenging to integrate them in a single movement analysis.

METHODS

We present an analytic workflow that integrates the behavioural, temporal and spatial components of the movement process and their interactions, which also allows for the assessment of the relative importance of those components. We construct a daily cyclic covariate to represent temporally cyclic movement patterns, such as diel variation in activity, and combine the three components in a multi-modal Hidden Markov Model framework using existing methods and R functions. We compare the trends and statistical fits of models that include or exclude any of the behavioural, spatial and temporal components, and perform variance partitioning on the model predictions that included all components to assess their relative importance to the movement process, both in isolation and in interaction.

RESULTS

We apply our workflow to a case study on the movements of plains zebra, blue wildebeest and eland antelope in a South African reserve. Behavioural modes impacted movement the most, followed by diel rhythms and then the spatial environment (viz. tree cover and terrain slope). Interactions between the components often explained more of the movement variation than the marginal effect of the spatial environment did on its own. Omitting components from the analysis led either to the inability to detect relationships between input and response variables, resulting in overgeneralisations when drawing conclusions about the movement process, or to detections of questionable relationships that appeared to be spurious.

CONCLUSIONS

Our analytic workflow can be used to integrate the behavioural, temporal and spatial components of the movement process and quantify their relative contributions, thereby preventing incomplete or overly generic ecological interpretations. We demonstrate that understanding the drivers of animal movement, and ultimately the ecological phenomena that emerge from it, critically depends on considering the various components of the movement process, and especially the interactions between them.

Bridging the gap between movement data and connectivity analysis using the Time-Explicit Habitat Selection (TEHS) model

BACKGROUND

Understanding how to connect habitat remnants to facilitate the movement of species is a critical task in an increasingly fragmented world impacted by human activities. The identification of dispersal routes and corridors through connectivity analysis requires measures of landscape resistance but there has been no consensus on how to calculate resistance from habitat characteristics, potentially leading to very different connectivity outcomes.

METHODS

We propose a new model, called the Time-Explicit Habitat Selection (TEHS) model, that can be directly used for connectivity analysis. The TEHS model decomposes the movement process in a principled approach into a time and a selection component, providing complementary information regarding space use by separately assessing the drivers of time to traverse the landscape and the drivers of habitat selection. These models are illustrated using GPS-tracking data from giant anteaters (Myrmecophaga tridactyla) in the Pantanal wetlands of Brazil.

RESULTS

The time model revealed that the fastest movements tended to occur between 8 p.m. and 5 a.m., suggesting a crepuscular/nocturnal behavior. Giant anteaters moved faster over wetlands while moving much slower over forests and savannas, in comparison to grasslands. We also found that wetlands were consistently avoided whereas forest and savannas tended to be selected. Importantly, this model revealed that selection for forest increased with temperature, suggesting that forests may act as important thermal shelters when temperatures are high. Finally, using the spatial absorbing Markov chain framework, we show that the TEHS model results can be used to simulate movement and connectivity within a fragmented landscape, revealing that giant anteaters will often not use the shortest-distance path to the destination patch due to avoidance of certain habitats.

CONCLUSIONS

The proposed approach can be used to characterize how landscape features are perceived by individuals through the decomposition of movement patterns into a time and a habitat selection component. Additionally, this framework can help bridge the gap between movement-based models and connectivity analysis, enabling the generation of time-explicit connectivity results.

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.

An integrated path for spatial capture-recapture and animal movement modeling

Ecologists and conservation biologists increasingly rely on spatial capture-recapture (SCR) and movement modeling to study animal populations. Historically, SCR has focused on population-level processes (e.g., vital rates, abundance, density, and distribution), whereas animal movement modeling has focused on the behavior of individuals (e.g., activity budgets, resource selection, migration). Even though animal movement is clearly a driver of population-level patterns and dynamics, technical and conceptual developments to date have not forged a firm link between the two fields. Instead, movement modeling has typically focused on the individual level without providing a coherent scaling from individual- to population-level processes, whereas SCR has typically focused on the population level while greatly simplifying the movement processes that give rise to the observations underlying these models. In our view, the integration of SCR and animal movement modeling has tremendous potential for allowing ecologists to scale up from individuals to populations and advancing the types of inferences that can be made at the intersection of population, movement, and landscape ecology. Properly accounting for complex animal movement processes can also potentially reduce bias in estimators of population-level parameters, thereby improving inferences that are critical for species conservation and management. This introductory article to the Special Feature reviews recent advances in SCR and animal movement modeling, establishes a common notation, highlights potential advantages of linking individual-level (Lagrangian) movements to population-level (Eulerian) processes, and outlines a general conceptual framework for the integration of movement and SCR models. We then identify important avenues for future research, including key challenges and potential pitfalls in the developments and applications that lie ahead.

Modeling abundance, distribution, movement and space use with camera and telemetry data

Studies of animal abundance and distribution are often conducted independently of research on movement, despite the important links between processes. Movement can cause rapid changes in spatial variation in density, and movement influences detection probability and therefore estimates of abundance from inferential methods such as spatial capture-recapture (SCR). Technological developments including camera traps and GPS telemetry have opened new opportunities for studying animal demography and movement, yet statistical models for these two data types have largely developed along parallel tracks. We present a hierarchical model in which both datasets are conditioned on a movement process for a clearly defined population. We fitted the model to data from 60 camera traps and 23,572 GPS telemetry locations collected on 17 male white-tailed deer in the Big Cypress National Preserve, Florida, USA during July 2015. Telemetry data were collected on a 3-4 h acquisition schedule, and we modeled the movement paths of all individuals in the region with a Ornstein-Uhlenbeck process that included individual-specific random effects. Two of the 17 deer with GPS collars were detected on cameras. An additional 20 male deer without collars were detected on cameras and individually identified based on their unique antler characteristics. Abundance was 126 (95% CI: 88-177) in the 228 km² region, only slightly higher than estimated using a standard SCR model: 119 (84-168). The standard SCR model, however, was unable to describe individual heterogeneity in movement rates and space use as revealed by the joint model. Joint modeling allowed the telemetry data to inform the movement model and the SCR encounter model, while leveraging information in the camera data to inform abundance, distribution and movement. Unlike most existing methods for population-level inference on movement, the joint SCR-movement model can yield unbiased inferences even if non-uniform sampling is used to deploy transmitters. Potential extensions of the model include the addition of resource selection parameters, and relaxation of the closure assumption when interest lies in survival and recruitment. These developments would contribute to the emerging holistic framework for the study of animal ecology, one that uses modern technology and spatio-temporal statistics to learn about interactions between behavior and demography.

Simple metrics to characterize inter-individual and temporal variation in habitat selection behaviour

Individual variation in habitat selection and movement behavior is receiving growing attention, but primarily with respect to characterizing behaviors in different contexts as opposed to decomposing structure in behavior within populations. This focus may be limiting advances in understanding the diversity of individual behavior and its influence on population organization. We propose a framework for characterizing variation in space-use behavior with the aim of advancing interpretation of its form and function. Using outputs from integrated Step Selection Analyses of 20 years of telemetry data from African elephants (Loxodonta Africana), we developed four metrics characterizing differentiation in resource selection behavior within a population [specialization (magnitude of the response independent of direction), heterogeneity (inter-individual variation), consistency (temporal shift in response) and reversal (frequency of directional changes in the response)]. We contrast insight from the developed metrics relative to the mean population response using an example focused on two covariates. We then expanded this contrast by evaluating if the metrics identify structurally important information on seasonal shifts in resource selection behaviors in addition to that provided by mean selection coefficients through Principal Component Analyses (PCAs) and a random forest classification. The simplified example highlighted that for some covariates focusing on the population average failed to capture complex individual variation in behaviors. The PCAs revealed that the developed metrics provided additional information in explaining the patterns in elephant selection beyond that offered by population average covariate values. For elephants, specialization and heterogeneity were informative, with specialization often being a better descriptor of differences in seasonal resource selection behavior than population average responses. Summarizing these metrics spatially and temporally, we illustrate how these metrics can provide insights on overlooked aspects of animal behavior. Our work offers a new approach in how we conceptualize variation in space-use behavior (i.e., habitat selection and movement) by providing ways of encapsulating variation that enables diagnoses of the drivers of individual level variability in a population. The developed metrics explicitly distill how variation in a behavior is structured among individuals and over time which could facilitate comparative work across time, populations, or strata within populations.

A role for lakes in revealing the nature of animal movement using high dimensional telemetry systems

Movement ecology is increasingly relying on experimental approaches and hypothesis testing to reveal how, when, where, why, and which animals move. Movement of megafauna is inherently interesting but many of the fundamental questions of movement ecology can be efficiently tested in study systems with high degrees of control. Lakes can be seen as microcosms for studying ecological processes and the use of high-resolution positioning systems to triangulate exact coordinates of fish, along with sensors that relay information about depth, temperature, acceleration, predation, and more, can be used to answer some of movement ecology's most pressing questions. We describe how key questions in animal movement have been approached and how experiments can be designed to gather information about movement processes to answer questions about the physiological, genetic, and environmental drivers of movement using lakes. We submit that whole lake telemetry studies have a key role to play not only in movement ecology but more broadly in biology as key scientific arenas for knowledge advancement. New hardware for tracking aquatic animals and statistical tools for understanding the processes underlying detection data will continue to advance the potential for revealing the paradigms that govern movement and biological phenomena not just within lakes but in other realms spanning lands and oceans.

What processes must we understand to forecast regional-scale population dynamics?

An urgent challenge facing biologists is predicting the regional-scale population dynamics of species facing environmental change. Biologists suggest that we must move beyond predictions based on phenomenological models and instead base predictions on underlying processes. For example, population biologists, evolutionary biologists, community ecologists and ecophysiologists all argue that the respective processes they study are essential. Must our models include processes from all of these fields? We argue that answering this critical question is ultimately an empirical exercise requiring a substantial amount of data that have not been integrated for any system to date. To motivate and facilitate the necessary data collection and integration, we first review the potential importance of each mechanism for skilful prediction. We then develop a conceptual framework based on reaction norms, and propose a hierarchical Bayesian statistical framework to integrate processes affecting reaction norms at different scales. The ambitious research programme we advocate is rapidly becoming feasible due to novel collaborations, datasets and analytical tools.

The context dependence of non-consumptive predator effects

Non-consumptive predator effects (NCEs) are now widely recognised for their capacity to shape ecosystem structure and function. Yet, forecasting the propagation of these predator-induced trait changes through particular communities remains a challenge. Accordingly, focusing on plasticity in prey anti-predator behaviours, we conceptualise the multi-stage process by which predators trigger direct and indirect NCEs, review and distil potential drivers of contingencies into three key categories (properties of the prey, predator and setting), and then provide a general framework for predicting both the nature and strength of direct NCEs. Our review underscores the myriad factors that can generate NCE contingencies while guiding how research might better anticipate and account for them. Moreover, our synthesis highlights the value of mapping both habitat domains and prey-specific patterns of evasion success ('evasion landscapes') as the basis for predicting how direct NCEs are likely to manifest in any particular community. Looking ahead, we highlight two key knowledge gaps that continue to impede a comprehensive understanding of non-consumptive predator-prey interactions and their ecosystem consequences; namely, insufficient empirical exploration of (1) context-dependent indirect NCEs and (2) the ways in which direct and indirect NCEs are shaped interactively by multiple drivers of context dependence.

Statistical Implementations of Agent-Based Demographic Models

A variety of demographic statistical models exist for studying population dynamics when individuals can be tracked over time. In cases where data are missing due to imperfect detection of individuals, the associated measurement error can be accommodated under certain study designs (e.g. those that involve multiple surveys or replication). However, the interaction of the measurement error and the underlying dynamic process can complicate the implementation of statistical agent-based models (ABMs) for population demography. In a Bayesian setting, traditional computational algorithms for fitting hierarchical demographic models can be prohibitively cumbersome to construct. Thus, we discuss a variety of approaches for fitting statistical ABMs to data and demonstrate how to use multi-stage recursive Bayesian computing and statistical emulators to fit models in such a way that alleviates the need to have analytical knowledge of the ABM likelihood. Using two examples, a demographic model for survival and a compartment model for COVID-19, we illustrate statistical procedures for implementing ABMs. The approaches we describe are intuitive and accessible for practitioners and can be parallelised easily for additional computational efficiency.

Movement reveals reproductive tactics in male elephants

Long-term bio-logging has the potential to reveal how movements, and hence life-history trade-offs, vary over a lifetime. Reproductive tactics in particular may vary as individuals' trade-off current investment versus lifetime fitness. Male African savanna elephants (Loxodona africana) provide a telling example of balancing body growth with reproductive fitness due to the combination of indeterminate growth and strongly delineated periods of sexual activity (musth), which results in reproductive tactics that alter with age. Our study aims to quantify the extent to which male elephants alter their movement patterns, and hence energetic allocation, in relation to (a) reproductive state and (b) age, and (c) to determine whether musth periods can be detected directly from GPS tracking data. We used a combination of GPS tracking data and visual observations of 25 male elephants ranging in age from 20 to 52 years to examine the influence of reproductive state and age on movement. We then used a three-state hidden Markov model (HMM) to detect musth behaviour in a subset of sequential tracking data. Our results demonstrate that male elephants increased their daily mean speed and range size with age and in musth. Furthermore, non-musth speed decreased with age, presumably reflecting a shift towards energy acquisition during non-musth. Thus, despite similar speeds and marginally larger ranges between reproductive states at age 20, by age 50, males were travelling 2.0 times faster in a 3.5 times larger area in musth relative to non-musth. The distinctiveness of musth periods over age 35 meant the three-state HMM could automatically detect musth movement with high sensitivity and specificity, but could not for the younger age class. We show that male elephants increased their energetic allocation into reproduction with age as the probability of reproductive success increases. Given that older male elephants tend to be both the target of legal trophy hunting and illegal poaching, man-made interference could drive fundamental changes in elephant reproductive tactics. Bio-logging, as our study reveals, has the potential both to quantify mature elephant reproductive tactics remotely and to be used to institute proactive management strategies around the reproductive behaviour of this charismatic keystone species.

Optimizing the use of biologgers for movement ecology research

The paradigm-changing opportunities of biologging sensors for ecological research, especially movement ecology, are vast, but the crucial questions of how best to match the most appropriate sensors and sensor combinations to specific biological questions and how to analyse complex biologging data, are mostly ignored. Here, we fill this gap by reviewing how to optimize the use of biologging techniques to answer questions in movement ecology and synthesize this into an Integrated Biologging Framework (IBF). We highlight that multisensor approaches are a new frontier in biologging, while identifying current limitations and avenues for future development in sensor technology. We focus on the importance of efficient data exploration, and more advanced multidimensional visualization methods, combined with appropriate archiving and sharing approaches, to tackle the big data issues presented by biologging. We also discuss the challenges and opportunities in matching the peculiarities of specific sensor data to the statistical models used, highlighting at the same time the large advances which will be required in the latter to properly analyse biologging data. Taking advantage of the biologging revolution will require a large improvement in the theoretical and mathematical foundations of movement ecology, to include the rich set of high-frequency multivariate data, which greatly expand the fundamentally limited and coarse data that could be collected using location-only technology such as GPS. Equally important will be the establishment of multidisciplinary collaborations to catalyse the opportunities offered by current and future biologging technology. If this is achieved, clear potential exists for developing a vastly improved mechanistic understanding of animal movements and their roles in ecological processes and for building realistic predictive models.