Multivariate Bayesian clustering using covariate-informed components with application to boreal vegetation sensitivity.

Climate change is impacting both the distribution and abundance of vegetation, especially in far northern latitudes. The effects of climate change are different for every plant assemblage and vary heterogeneously in both space and time. Small changes in climate could result in large vegetation responses in sensitive assemblages but weak responses in robust assemblages. But, patterns and mechanisms of sensitivity and robustness are not yet well understood, largely due to a lack of long-term measurements of climate and vegetation. Fortunately, observations are sometimes available across a broad spatial extent. We develop a novel statistical model for a multivariate response based on unknown cluster-specific effects and covariances, where cluster labels correspond to sensitivity and robustness. Our approach utilizes a prototype model for cluster membership that offers flexibility while enforcing smoothness in cluster probabilities across sites with similar characteristics. We demonstrate our approach with an application to vegetation abundance in Alaska, USA, in which we leverage the broad spatial extent of the study area as a proxy for unrecorded historical observations. In the context of the application, our approach yields interpretable site-level cluster labels associated with assemblage-level sensitivity and robustness without requiring strong a priori assumptions about the drivers of climate sensitivity.

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.

Running on empty: recharge dynamics from animal movement data.

Vital rates such as survival and recruitment have always been important in the study of population and community ecology. At the individual level, physiological processes such as energetics are critical in understanding biomechanics and movement ecology and also scale up to influence food webs and trophic cascades. Although vital rates and population-level characteristics are tied with individual-level animal movement, most statistical models for telemetry data are not equipped to provide inference about these relationships because they lack the explicit, mechanistic connection to physiological dynamics. We present a framework for modelling telemetry data that explicitly includes an aggregated physiological process associated with decision making and movement in heterogeneous environments. Our framework accommodates a wide range of movement and physiological process specifications. We illustrate a specific model formulation in continuous-time to provide direct inference about gains and losses associated with physiological processes based on movement. Our approach can also be extended to accommodate auxiliary data when available. We demonstrate our model to infer mountain lion (Puma concolor; in Colorado, USA) and African buffalo (Syncerus caffer; in Kruger National Park, South Africa) recharge dynamics.

The basis function approach for modeling autocorrelation in ecological data.

Analyzing ecological data often requires modeling the autocorrelation created by spatial and temporal processes. Many seemingly disparate statistical methods used to account for autocorrelation can be expressed as regression models that include basis functions. Basis functions also enable ecologists to modify a wide range of existing ecological models in order to account for autocorrelation, which can improve inference and predictive accuracy. Furthermore, understanding the properties of basis functions is essential for evaluating the fit of spatial or time-series models, detecting a hidden form of collinearity, and analyzing large data sets. We present important concepts and properties related to basis functions and illustrate several tools and techniques ecologists can use when modeling autocorrelation in ecological data.