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

Redefining 'state-of-the-art' for integrated population models with immigration

Research Highlight: Christian, M., Oosthuizen, W. C., Bester, M. N., & de Bruyn, P. N. (2024). Robustly estimating the demographic contribution of immigration: Simulation, sensitivity analysis and seals. Journal of Animal Ecology. https://doi.org/10.1111/1365-2656.14053. Immigration can have profound consequences for local population dynamics and demography, but collecting data to accurately quantifying it is challenging. The recent rise of integrated population models (IPMs) offers an alternative by making it possible to estimate immigration without the need for explicit data, and to quantify its contribution to population dynamics through transient Life Table Response Experiments (tLTREs). Simulation studies have, however, highlighted that this approach can be prone to bias and overestimation. In their new study, Christian et al. address one of the root causes of this issue by improving the estimation of time variation in vital rates and immigration using Gaussian processes in lieu of traditionally used temporal random effects. They demonstrate that IPM-tLTRE frameworks with Gaussian processes produce more accurate and less biased estimates of immigration and its contribution to population dynamics and illustrate the applicability of this approach using a long-term data set on elephant seals (Mirounga leonida). Results are validated with a simulation study and suggest that immigration of breeding females has been central for population recovery of elephant seals despite the species' high female site fidelity. Christian et al. thus present new insights into population regulation of long-lived marine mammals and highlight the potential for using Gaussian process priors in IPMs. They also illustrate a suite of 'best practices' for state-of-the-art IPM-tLTRE analyses and provide an inspirational example for the kind of ecological modelling workflow that can be invaluable not just as a starting point for fellow ecologists picking up or improving their own IPM-tLTRE analyses, but also for teaching and in contexts where model estimates are used for informing management and conservation decision-making.

A flexible framework for spatial capture-recapture with unknown identities

Camera traps or acoustic recorders are often used to sample wildlife populations. When animals can be individually identified, these data can be used with spatial capture-recapture (SCR) methods to assess populations. However, obtaining animal identities is often labor-intensive and not always possible for all detected animals. To address this problem, we formulate SCR, including acoustic SCR, as a marked Poisson process, comprising a single counting process for the detections of all animals and a mark distribution for what is observed (eg, animal identity, detector location). The counting process applies equally when it is animals appearing in front of camera traps and when vocalizations are captured by microphones, although the definition of a mark changes. When animals cannot be uniquely identified, the observed marks arise from a mixture of mark distributions defined by the animal activity centers and additional characteristics. Our method generalizes existing latent identity SCR models and provides an integrated framework that includes acoustic SCR. We apply our method to estimate density from a camera trap study of fisher (Pekania pennanti) and an acoustic survey of Cape Peninsula moss frog (Arthroleptella lightfooti). We also test it through simulation. We find latent identity SCR with additional marks such as sex or time of arrival to be a reliable method for estimating animal density.

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.

Crabs ride the tide: incoming tides promote foraging of Giant Mud Crab (Scylla serrata)

BACKGROUND

Effective fisheries management of mobile species relies on robust knowledge of animal behaviour and habitat-use. Indices of behaviour can be useful for interpreting catch-per-unit-effort data which acts as a proxy for relative abundance. Information about habitat-use can inform stocking release strategies or the design of marine protected areas. The Giant Mud Crab (Scylla serrata; Family: Portunidae) is a swimming estuarine crab that supports significant fisheries harvest throughout the Indo-West Pacific, but little is known about the fine-scale movement and behaviour of this species.

METHODS

We tagged 18 adult Giant Mud Crab with accelerometer-equipped acoustic tags to track their fine-scale movement using a hyperbolic positioning system, alongside high temporal resolution environmental data (e.g., water temperature), in a temperate south-east Australian estuary. A hidden Markov model was used to classify movement (i.e., step length, turning angle) and acceleration data into discrete behaviours, while also considering the possibility of individual variation in behavioural dynamics. We then investigated the influence of environmental covariates on these behaviours based on previously published observations.

RESULTS

We fitted a model with two well-distinguished behavioural states describing periods of inactivity and foraging, and found no evidence of individual variation in behavioural dynamics. Inactive periods were most common (79% of time), and foraging was most likely during low, incoming tides; while inactivity was more likely as the high tide receded. Model selection removed time (hour) of day and water temperature (°C) as covariates, suggesting that they do not influence Giant Mud Crab behavioural dynamics at the temporal scale investigated.

CONCLUSIONS

Our study is the first to quantitatively link fine-scale movement and behaviour of Giant Mud Crab to environmental variation. Our results suggest Giant Mud Crab are a predominantly sessile species, and support their status as an opportunistic scavenger. We demonstrate a relationship between the tidal cycle and foraging that is likely to minimize predation risk while maximizing energetic efficiency. These results may explain why tidal covariates influence catch rates in swimming crabs, and provide a foundation for standardisation and interpretation of catch-per-unit-effort data-a commonly used metric in fisheries science.

Daily activity timing in the Anthropocene

Animals are facing novel 'timescapes' in which the stimuli entraining their daily activity patterns no longer match historical conditions due to anthropogenic disturbance. However, the ecological effects (e.g., altered physiology, species interactions) of novel activity timing are virtually unknown. We reviewed 1328 studies and found relatively few focusing on anthropogenic effects on activity timing. We suggest three hypotheses to stimulate future research: (i) activity-timing mismatches determine ecological effects, (ii) duration and timing of timescape modification influence effects, and (iii) consequences of altered activity timing vary biogeographically due to broad-scale variation in factors compressing timescapes. The continued growth of sampling technologies promises to facilitate the study of the consequences of altered activity timing, with emerging applications for biodiversity conservation.

Integrated animal movement and spatial capture-recapture models: Simulation, implementation, and inference

Over the last decade, spatial capture-recapture (SCR) models have become widespread for estimating demographic parameters in ecological studies. However, the underlying assumptions about animal movement and space use are often not realistic. This is a missed opportunity because interesting ecological questions related to animal space use, habitat selection, and behavior cannot be addressed with most SCR models, despite the fact that the data collected in SCR studies - individual animals observed at specific locations and times - can provide a rich source of information about these processes and how they relate to demographic rates. We developed SCR models that integrated more complex movement processes that are typically inferred from telemetry data, including a simple random walk, correlated random walk (i.e., short-term directional persistence), and habitat-driven Langevin diffusion. We demonstrated how to formulate, simulate from, and fit these models with standard SCR data using data-augmented Bayesian analysis methods. We evaluated their performance through a simulation study, in which we varied the detection, movement, and resource selection parameters. We also examined different numbers of sampling occasions and assessed performance gains when including auxiliary location data collected from telemetered individuals. Across all scenarios, the integrated SCR movement models performed well in terms of abundance, detection, and movement parameter estimation. We found little difference in bias for the simple random walk model when reducing the number of sampling occasions from T = 25 to T = 15. We found some bias in movement parameter estimates under several of the correlated random walk scenarios, but incorporating auxiliary location data improved parameter estimates and significantly improved mixing during model fitting. The Langevin movement model was able to recover resource selection parameters from standard SCR data, which is particularly appealing because it explicitly links the individual-level movement process with habitat selection and population density. We focused on closed population models, but the movement models developed here can be extended to open SCR models. The movement process models could also be easily extended to accommodate additional "building blocks" of random walks, such as central tendency (e.g., territoriality) or multiple movement behavior states, thereby providing a flexible and coherent framework for linking animal movement behavior to population dynamics, density, and distribution.

Modeling spatiotemporal abundance and movement dynamics using an integrated spatial capture-recapture movement model

Animal movement is a fundamental ecological process affecting the survival and reproduction of individuals, the structure of populations, and the dynamics of communities. Methods to quantify animal movement and spatiotemporal abundances, however, are generally separate and therefore omit linkages between individual-level and population-level processes. We describe an integrated spatial capture-recapture (SCR) movement model to jointly estimate (1) the number and distribution of individuals in a defined spatial region and (2) movement of those individuals through time. We applied our model to a study of polar bears (Ursus maritimus) in a 28,125 km² survey area of the eastern Chukchi Sea, USA in 2015 that incorporated capture-recapture and telemetry data. In simulation studies, the model provided unbiased estimates of movement, abundance, and detection parameters using a bivariate normal random walk and correlated random walk movement process. Our case study provided detailed evidence of directional movement persistence for both male and female bears, where individuals regularly traversed areas larger than the survey area during the 36-day study period. Scaling from individual- to population-level inferences, we found that densities varied from <0.75 bears/625 km² grid cell/day in nearshore cells to 1.6-2.5 bears/grid cell/day for cells surrounded by sea ice. Daily abundance estimates ranged from 53 to 69 bears, with no trend across days. The cumulative number of unique bears that used the survey area increased through time due to movements into and out of the area, resulting in an estimated 171 individuals using the survey area during the study (95% credible interval 124-250). Abundance estimates were similar to a previous multiyear integrated population model using capture-recapture and telemetry data (2008-2016; Regehr et al., Scientific Reports 8:16780, 2018). Overall, the SCR-movement model successfully quantified both individual- and population-level space use, including the effects of landscape characteristics on movement, abundance, and detection, while linking the movement and abundance processes to directly estimate density within a prescribed spatial region and temporal period. Integrated SCR-movement models provide a generalizable approach to incorporate greater movement realism into population dynamics and link movement to emergent properties including spatiotemporal densities and abundances.

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.

Confronting spatial capture-recapture models with realistic animal movement simulations

Spatial capture-recapture (SCR) models have emerged as a robust method to estimate the population density of mobile animals. However, model evaluation has generally been based on data simulated from simplified representations of animal space use. Here, we generated data from animal movement simulated from a mechanistic individual-based model, in which movement emerges from the individual's response to a changing environment (i.e., from the bottom-up), driven by key ecological processes (e.g., resource memory and territoriality). We drew individual detection data from simulated movement trajectories and fitted detection data sets to a basic, resource selection and transience SCR model, as well as their variants accounting for resource-driven heterogeneity in density and detectability. Across all SCR models, abundance estimates were robust to multiple, but low-degree violations of the specified movement processes (e.g., resource selection). SCR models also successfully captured the positive effect of resource quality on density. However, covariate models failed to capture the finer scale effect of resource quality on detectability and space use, which may be a consequence of the low temporal resolution of SCR data sets and/or model misspecification. We show that home-range size is challenging to infer from the scale parameter alone, compounded by reliance on conventional measures of "true" home-range size that are highly sensitive to sampling regime. Additionally, we found the transience model challenging to fit, probably due to data sparsity and violation of the assumption of normally distributed inter-occasion movement of activity centers, suggesting that further development of the model is required for general applicability. Our results showed that further integration of complex movement into SCR models may not be necessary for population estimates of abundance when the level of individual heterogeneity induced by the underlying movement process is low, but appears warranted in terms of accurately revealing finer scale patterns of ecological and movement processes. Further investigation into whether this holds true in populations with other types of realistic movement characteristics is merited. Our study provides a framework to generate realistic SCR data sets to develop and evaluate more complex movement processes in SCR models.