Imputation Approaches for Animal Movement Modeling

Evaluating the Risk Landscape of Hawaiian Monk Seal Exposure to Toxoplasma gondii

Toxoplasmosis is a disease of primary concern for Hawaiian monk seals (Neomonachus schauinslandi), due to its apparently acute lethality and especially heavy impacts on breeding female seals. The disease-causing parasite, Toxoplasma gondii, depends on cats to complete its life cycle; thus, in order to understand how this pathogen infects marine mammals, it is essential to understand aspects of the terrestrial ecosystem and land-to-sea transport. In this study, we constructed a three-tiered model to assess risk of Hawaiian monk seal exposure to T. gondii oocysts: (1) oocyst contamination as a function of cat population characteristics; (2) land-to-sea transport of oocysts as a function of island hydrology, and (3) seal exposure as a function of habitat and space use. We were able to generate risk maps highlighting watersheds contributing the most to oocyst contamination of Hawaiian monk seal habitat. Further, the model showed that free-roaming cats most associated with humans (pets or strays often supplementally fed by people) were able to achieve high densities leading to high levels of oocyst contamination and elevated risk of T. gondii exposure.

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.

Comparison of imputation methods for missing production data of dairy cattle

Nowadays, vast amounts of data representing feed intake, growth, and environmental impact of individual animals are being recorded in on-farm settings. Despite their apparent use, data collected in real-world applications often have missing values in one or several variables, due to reasons including human error, machine error, or sampling frequency misalignment across multiple variables. Since incomplete datasets are less valuable for downstream data analysis, it is important to address the missing value problem properly. One option may be to reduce the dataset to a subset that contains only complete data, but considerable data may be lost via this process. The current study aimed to compare imputation methods for the estimation of missing values in a raw dataset of dairy cattle including 454 553 records collected from 629 cows between 2009 and 2020. The dataset was subjected to a cleaning process that reduced its size to 437 075 observations corresponding to 512 cows. Missing values were present in four variables: concentrate DM intake (CDMI, missing percentage = 2.30%), forage DM intake (FDMI, 8.05%), milk yield (MY, 15.12%), and BW (64.33%). After removing all missing values, the resulting dataset (n = 129 353) was randomly sampled five times to create five independent subsets that exhibit the same missing data percentages as the cleaned dataset. Four univariate and nine multivariate imputation methods (eight machine learning methods and the MissForest method) were applied and evaluated on the five repeats, and average imputation performance was reported for each repeat. The results showed that Random Forest was overall the best imputation method for this type of data and had a lower mean squared prediction error and higher concordance correlation coefficient than the other imputation methods for all imputed variables. Random Forest performed particularly well for imputing CDMI, MY, and BW, compared to imputing FDMI.

Latent trajectory models for spatio-temporal dynamics in Alaskan ecosystems

The Alaskan landscape has undergone substantial changes in recent decades, most notably the expansion of shrubs and trees across the Arctic. We developed a Bayesian hierarchical model to quantify the impact of climate change on the structural transformation of ecosystems using remotely sensed imagery. We used latent trajectory processes to model dynamic state probabilities that evolve annually, from which we derived transition probabilities between ecotypes. Our latent trajectory model accommodates temporal irregularity in survey intervals and uses spatio-temporally heterogeneous climate drivers to infer rates of land cover transitions. We characterized multi-scale spatial correlation induced by plot and subplot arrangements in our study system. We also developed a Pólya-Gamma sampling strategy to improve computation. Our model facilitates inference on the response of ecosystems to shifts in the climate and can be used to predict future land cover transitions under various climate scenarios.

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.

Intrapopulation differences in polar bear movement and step selection patterns

BACKGROUND

The spatial ecology of individuals often varies within a population or species. Identifying how individuals in different classes interact with their environment can lead to a better understanding of population responses to human activities and environmental change and improve population estimates. Most inferences about polar bear (Ursus maritimus) spatial ecology are based on data from adult females due to morphological constraints on applying satellite radio collars to other classes of bears. Recent studies, however, have provided limited movement data for adult males and sub-adults of both sexes using ear-mounted and glue-on tags. We evaluated class-specific movements and step selection patterns for polar bears in the Chukchi Sea subpopulation during spring.

METHODS

We developed hierarchical Bayesian models to evaluate polar bear movement (i.e., step length and directional persistence) and step selection at the scale of 4-day step lengths. We assessed differences in movement and step selection parameters among the three classes of polar bears (i.e., adult males, sub-adults, and adult females without cubs-of-the-year).

RESULTS

Adult males had larger step lengths and less directed movements than adult females. Sub-adult movement parameters did not differ from the other classes but point estimates were most similar to adult females. We did not detect differences among polar bear classes in step selection parameters and parameter estimates were consistent with previous studies.

CONCLUSIONS

Our findings support the use of estimated step selection patterns from adult females as a proxy for other classes of polar bears during spring. Conversely, movement analyses indicated that using data from adult females as a proxy for the movements of adult males is likely inappropriate. We recommend that researchers consider whether it is valid to extend inference derived from adult female movements to other classes, based on the questions being asked and the spatial and temporal scope of the data. Because our data were specific to spring, these findings highlight the need to evaluate differences in movement and step selection during other periods of the year, for which data from ear-mounted and glue-on tags are currently lacking.

Discrete-space continuous-time models of marine mammal exposure to Navy sonar

Assessing the patterns of wildlife attendance to specific areas is relevant across many fundamental and applied ecological studies, particularly when animals are at risk of being exposed to stressors within or outside the boundaries of those areas. Marine mammals are increasingly being exposed to human activities that may cause behavioral and physiological changes, including military exercises using active sonars. Assessment of the population-level consequences of anthropogenic disturbance requires robust and efficient tools to quantify the levels of aggregate exposure for individuals in a population over biologically relevant time frames. We propose a discrete-space, continuous-time approach to estimate individual transition rates across the boundaries of an area of interest, informed by telemetry data collected with uncertainty. The approach allows inferring the effect of stressors on transition rates, the progressive return to baseline movement patterns, and any difference among individuals. We apply the modeling framework to telemetry data from Blainville's beaked whale (Mesoplodon densirostris) tagged in the Bahamas at the Atlantic Undersea Test and Evaluation Center (AUTEC), an area used by the U.S. Navy for fleet readiness training. We show that transition rates changed as a result of exposure to sonar exercises in the area, reflecting an avoidance response. Our approach supports the assessment of the aggregate exposure of individuals to sonar and the resulting population-level consequences. The approach has potential applications across many applied and fundamental problems where telemetry data are used to characterize animal occurrence within specific areas.

Recursive Bayesian computation facilitates adaptive optimal design in ecological studies

Optimal design procedures provide a framework to leverage the learning generated by ecological models to flexibly and efficiently deploy future monitoring efforts. At the same time, Bayesian hierarchical models have become widespread in ecology and offer a rich set of tools for ecological learning and inference. However, coupling these methods with an optimal design framework can become computationally intractable. Recursive Bayesian computation offers a way to substantially reduce this computational burden, making optimal design accessible for modern Bayesian ecological models. We demonstrate the application of so-called prior-proposal recursive Bayes to optimal design using a simulated data binary regression and the real-world example of monitoring and modeling sea otters in Glacier Bay, Alaska. These examples highlight the computational gains offered by recursive Bayesian methods and the tighter fusion of monitoring and science that those computational gains enable.

Accounting for phenology in the analysis of animal movement

The analysis of animal tracking data provides important scientific understanding and discovery in ecology. Observations of animal trajectories using telemetry devices provide researchers with information about the way animals interact with their environment and each other. For many species, specific geographical features in the landscape can have a strong effect on behavior. Such features may correspond to a single point (eg, dens or kill sites), or to higher dimensional subspaces (eg, rivers or lakes). Features may be relatively static in time (eg, coastlines or home-range centers), or may be dynamic (eg, sea ice extent or areas of high-quality forage for herbivores). We introduce a novel model for animal movement that incorporates active selection for dynamic features in a landscape. Our approach is motivated by the study of polar bear (Ursus maritimus) movement. During the sea ice melt season, polar bears spend much of their time on sea ice above shallow, biologically productive water where they hunt seals. The changing distribution and characteristics of sea ice throughout the year mean that the location of valuable habitat is constantly shifting. We develop a model for the movement of polar bears that accounts for the effect of this important landscape feature. We introduce a two-stage procedure for approximate Bayesian inference that allows us to analyze over 300 000 observed locations of 186 polar bears from 2012 to 2016. We use our model to estimate a spatial boundary of interest to wildlife managers that separates two subpopulations of polar bears from the Beaufort and Chukchi seas.

Time-varying predatory behavior is primary predictor of fine-scale movement of wildland-urban cougars

BACKGROUND

While many species have suffered from the detrimental impacts of increasing human population growth, some species, such as cougars (Puma concolor), have been observed using human-modified landscapes. However, human-modified habitat can be a source of both increased risk and increased food availability, particularly for large carnivores. Assessing preferential use of the landscape is important for managing wildlife and can be particularly useful in transitional habitats, such as at the wildland-urban interface. Preferential use is often evaluated using resource selection functions (RSFs), which are focused on quantifying habitat preference using either a temporally static framework or researcher-defined temporal delineations. Many applications of RSFs do not incorporate time-varying landscape availability or temporally-varying behavior, which may mask conflict and avoidance behavior.

METHODS

Contemporary approaches to incorporate landscape availability into the assessment of habitat selection include spatio-temporal point process models, step selection functions, and continuous-time Markov chain (CTMC) models; in contrast with the other methods, the CTMC model allows for explicit inference on animal movement in continuous-time. We used a hierarchical version of the CTMC framework to model speed and directionality of fine-scale movement by a population of cougars inhabiting the Front Range of Colorado, U.S.A., an area exhibiting rapid population growth and increased recreational use, as a function of individual variation and time-varying responses to landscape covariates.

RESULTS

We found evidence for individual- and daily temporal-variability in cougar response to landscape characteristics. Distance to nearest kill site emerged as the most important driver of movement at a population-level. We also detected seasonal differences in average response to elevation, heat loading, and distance to roads. Motility was also a function of amount of development, with cougars moving faster in developed areas than in undeveloped areas.

CONCLUSIONS

The time-varying framework allowed us to detect temporal variability that would be masked in a generalized linear model, and improved the within-sample predictive ability of the model. The high degree of individual variation suggests that, if agencies want to minimize human-wildlife conflict management options should be varied and flexible. However, due to the effect of recursive behavior on cougar movement, likely related to the location and timing of potential kill-sites, kill-site identification tools may be useful for identifying areas of potential conflict.

Monitoring dynamic spatio-temporal ecological processes optimally

Population dynamics vary in space and time. Survey designs that ignore these dynamics may be inefficient and fail to capture essential spatio-temporal variability of a process. Alternatively, dynamic survey designs explicitly incorporate knowledge of ecological processes, the associated uncertainty in those processes, and can be optimized with respect to monitoring objectives. We describe a cohesive framework for monitoring a spreading population that explicitly links animal movement models with survey design and monitoring objectives. We apply the framework to develop an optimal survey design for sea otters in Glacier Bay. Sea otters were first detected in Glacier Bay in 1988 and have since increased in both abundance and distribution; abundance estimates increased from 5 otters to >5,000 otters, and they have spread faster than 2.7 km/yr. By explicitly linking animal movement models and survey design, we are able to reduce uncertainty associated with forecasting occupancy, abundance, and distribution compared to other potential random designs. The framework we describe is general, and we outline steps to applying it to novel systems and taxa.