Delineating important killer whale foraging areas using a spatiotemporal logistic model

Conservation management planning for highly mobile species requires an understanding of the distribution of areas that are biologically important to the species of concern. Collecting data on the locations of animal behaviors linked to biological characteristics, such as foraging, can be used to spatially describe biological important areas. However, spatial modeling of free-ranging animal behavior can be challenging, as behavioral observations of animals are often clustered, and sampling is commonly conducted at a higher frequency than changes in behavioral states, resulting in data that are usually highly autocorrelated in space and time. Here, we fit latent Gaussian process models to observational behavioral data to generate spatially-explicit pre-dictions of foraging behavior within the critical habitat of an endangered population of fish-eating killer whales (Orcinus orca) in southern British Columbia, Canada, and northern Washington State, USA. We compare spatial models treating temporal autocorrelation in behavior in three ways: (1) ignoring temporal autocorrelation entirely; (2) traditional data-thinning to remove temporal autocorrelation, and; (3) using a temporal Gaussian process to account for temporal autocorrelation. Comparisons of autocorrelative structures for each model and visual comparison of broad spatial patterns demonstrate that our third approach yields more accurate results than when ignoring temporal autocorrelation entirely and higher precision results than when applying data-thinning methods. Within the identified areas of critical habitat, our models indicate two primary regions of intense killer whale foraging activity, and we delineate areas wherein the probability of foraging was particularly high as candidate locations for conservation management actions. This study underscores the value of refining our understanding of high-use areas for at-risk species by incorporating animal behavior data to inform area-based conservation measures.