- Methodology
- Open access
- Published:
Unveiling the roles of temporal periodicity, the spatial environment and behavioural modes in terrestrial animal movement
Movement Ecology volume 12, Article number: 57 (2024)
Abstract
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.
Background
Movement is ubiquitous to animal life and central to emergent ecological phenomena such as territoriality and home range formation [1, 2], predator-prey dynamics [3, 4], distribution patterns [5] and population dynamics [6, 7]. To understand such phenomena it is thus key to study what drives the movement of animals through their environment [7, 8]. However, this understanding is hampered by the complex relationships between animal movement and a multitude of behavioural and spatial-temporal environmental drivers, since a specific state of a driver can give rise to different movement patterns depending on the state of other drivers. For example, animals might move slowly through thick vegetation when foraging, or fast if they are translocating from one preferred high-visibility habitat to another. Temporal variation in the response of animals to their environment in diel, lunar or seasonal patterns [e.g. 9], complicates the analysis of animal movement drivers even further. In order to structure the study of animal movement and advance the interpretation of movement patterns, a movement ecology framework has been put forward that divides the movement process into multiple components: an internal motive regarding why to move and the temporally and spatially varying environment influencing when and where to move respectively [10].
First, ‘why move’ relates to the internal motives of movement (e.g. reproduction, resource acquisition and threat avoidance), which are often manifested in different behaviours (e.g. mating, foraging and fleeing) that produce distinct movement modes (i.e., movement trajectory geometries that can be quantified by metrics like speed and turning angle distributions) [10,11,12]. Distinguishing different movement modes in animal trajectories has yielded insights in different foraging strategies within individuals [13], predator-prey interactions [14] and effects of anthropogenic disturbance on behaviour [15]. Internal motives and their emergent movement modes should thus be acknowledged to further the interpretation of animal movement, specifically regarding its iteraction with the environment [11,12,13, 16].
Second, ‘when to move’ deals with temporal variation in animal movement, which often exhibits cyclic patterns. For example, seasonal migratory movements [e.g. 17, 18] have long been recognized, as well as the influence of lunar rhythms [e.g. 19, 20]. Moreover, most animals show strong diel activity patterns [21,22,23,24,25], which can arise directly due to variation in light and thermal conditions [26] or indirectly due to predation risk and human activity [22, 27,28,29]. Temporal periodicity is thus a key characteristic of animal movement [9]. Several methods of varying complexity have been developed to address temporal variation in movement [24, 25, 30,31,32,33]. Basis functions can be used to accommodate temporal heterogeneity in continuous-time movement models [32], for example by modelling processes related to energy discharge and recharge with a temporally explicit recharge function [33], or by augmenting spatial movement models with a generalised additive model that describes temporal movement heterogeneity (such as seasonal and daily patterns) [25]. In discrete-time descriptions of movement, resource selection can be weighted with time-specific distributions to reveal time-varying effects of resource covariates on selection [31], or transition probabilities between movement modes can be modelled with temporal heterogeneity [30].
Third, ‘where to move’ relates the movement of animals to the characteristics of their spatial surroundings. For example, topography influences movement through energetic costs associated with moving through sloping or rugged terrain [34,35,36,37]. Daily cycles can interact strongly with the spatial surrounding, for example in animals that select areas with high vegetation cover during daytime and more open areas during the night to decrease predation risk [38, 39]. These examples illustrate that relationships between spatial environmental characteristics and movement patterns are conditional on the context of internal motive and time, as this context determines the costs and benefits that the spatial environmental characteristics represent [e.g. 40].
Internal motives, temporal periodicity and spatial environmental heterogeneity are thus three intricately linked components of the animal movement process. However, animal movement research has often been hampered by the focus on univariate analyses (Fig. 1) [41], while ideally the relationships between movement and each of these three components are considered within each other’s context [10]. Although the multi-component movement ecology framework has proven helpful in elucidating drivers of animal movement [42] and a plethora of methods exists to study the various components, challenges remain regarding how to (1) include the components in a single coherent modelling workflow, (2) quantify the relative contributions of the components in shaping animal movement, and (3) assess the influence of interactions between the components on the quantitative and qualitative inferences made about the movement process.
Here, we demonstrate an integrative analytic workflow, drawing from and combining existing methods, to evaluate the role of movement modes, time (temporal rhythms) and space (environmental heterogeneity) in shaping animal movement. We apply this workflow to a case study of large savanna ungulates. We quantify the relative contribution of these three components to the movement process in isolation as well as in interaction, and we assess how ecological interpretations of movement change when omitting components and interactions from the analysis. We use fine-scale (10-minute resolution) GPS tracking data of eland antelope (Taurotragus oryx), blue wildebeest (Connochaetes taurinus) and plains zebra (Equus quagga) in a South-African game reserve, and fit Hidden Markov Models (HMMs) with input from tree cover, terrain slope and a constructed daily cyclic covariate (which captures the diel effect on the movement behaviour of these crepuscular species). Since HMMs can describe the movement process as a multi-modal random walk with input from spatial-temporal environmental covariates [43], they provide an ideal method to quantify the effects of modes (as a proxy for internal motives) and temporal and spatial environmental heterogeneity on movement. We include effects of the time of day and the spatial environment in both the mode-specific movement properties speed and turning angle, as well as the switching probabilities between modes. Besides demonstrating our analytic workflow, we use our case study of savanna ungulates to demonstrate that inferences about the behavioural ecology of animals are made best when explicitly modelling the movement process with these three components in interaction.
Methods
Animals and study area
From June to August 2017, 34 eland antelope, 34 blue wildebeest and 35 plains zebra were captured in Welgevonden Game Reserve, Limpopo, South Africa (24°13’S, 27°54’E), and equipped with neck collars containing GPS and accelerometer sensors. For more details about the collaring process and sensors, see [44]. The animals were released in a fenced study area of 1200 hectares at the northern edge of the reserve (Fig. 2). As a part of the Waterberg mountain massif and Bushveld ecoregion [45], the study area is characterised by nutrient poor sandy soils in the lower northern part and rocky, flat hilltop plateaus dissected by ravines extending towards the south. Woody vegetation varies from mixed to open dry deciduous woodland, with the rocky plateaus being the most open.
Animal movement data
The GPS sensors collected positions with an interval ranging between 15 and 2 min, depending on the amount of body activity as determined by the accelerometer [44]. The position data were corrected, smoothed and modelled to regular 10 min resolution trajectories as described in [44]. Due to sensor failures, we selected from each species only the ten individuals with the most recorded data in the period from September to December 2017 (the first four months of the study period, during which the sensors collected most data). Our subset contained a total of 357,815 animal positions across 30 individuals of the three species, comprising 481 bouts (i.e. an uninterrupted and regular sequence of positions of a single individual) with a median bout length of 496 positions (3.4 days).
Environmental data
A digital elevation model of the study area with resolution 2 × 2 m was available, from which terrain slope (hereafter referred to as slope) was calculated (Fig. 2a). We used aerial images of the study area from August 2013 to create an orthomosaic with a 0.25 × 0.25 m resolution, from which pixels were classified into trees, grass, bare ground and other/built up area using the Semi-Automatic Classification Plugin for QGIS [46]. After having applied an isotropic Gaussian kernel smoother with a 10 m standard deviation to the classified pixels, we derived fractional tree cover (hereafter referred to as tree cover) and fractional grass cover maps of the study area at a 2 × 2 m resolution (Fig. 2b). Tree cover was strongly negatively correlated with grass cover (r = 0.95), so we used only tree cover in our analyses and considered it a measure for the inverse of grass cover as well. We extracted tree cover and slope values for each animal position.
Movement model structure
We constructed a Hidden Markov Model (HMM) for each species separately to examine how their movement was shaped by movement modes, diel patterns and spatial environmental heterogeneity. We described the movement process as a random walk consisting of steps between subsequent animal positions with lengths s and turning angles φ. Since intervals between steps were constant at ten minutes, we refer to s as speed throughout. We fitted gamma (Γ) and von Mises (vM) distributions to s and φ, respectively, which were unique to each discrete movement mode z. Thus, z was defined as a unique combination of three movement parameters: the mean μ and standard deviation σ of the speed (converted from the shape and scale parameter of the Γ distribution for easier interpretation) and the vM concentration parameter κ of the turning angle (the vM mean was fixed at 0 since there was no reason to assume a movement bias to either the right or left direction). μ, σ and κ were in turn modelled as functions of several covariates (see section Movement model covariates) using a log-link. Mode switching between the steps was governed by a switching matrix M, where the off-diagonal elements are logit-linear functions of covariates. We considered a three-mode, two-mode and a ‘single-mode’ HMM (i.e. 3 × 3, 2 × 2 and 1 × 1 dimensions for M respectively), where the single-mode model is analogous to a movement model that does not include the mode component. Thus, the HMM had the general form:
Where γ1→2 is the probability of switching from mode z = i to z = j between consecutive steps, fn are functions for mode switching probabilities and Γ and vM parameters, θn are sets of to-be-fitted coefficients of function fn (and for mode z = i in case of θn, i), and ξn are sets of covariates (see next section).
Movement model covariates
We modelled the parameters μ, σ and κ and the mode switching probabilities γ1→2 each as a function of three covariates. Two of those were the spatial environmental covariates tree cover and slope. The third covariate was a constructed daily cyclic covariate (DCC) that describes the diel patterns in μ, σ and κ (driven by e.g. ambient temperature, solar position, perceived predation risk and internal diel rhythm). To determine those diel patterns, we first divided all the animal positions in 10-minute time bins based on the time of day that they were taken (i.e., there were a total of 144 time bins of 10 min, totalling 24 h). We then fitted Γ distributions (parameters μ and σ) to the movement speeds and vM distributions (parameter κ) to the turning angles of the positions within a centred moving window of three time bins wide. This gave us empirical diel patterns in parameters μ, σ and κ (Fig. S1). We rescaled and averaged those empirical diel patterns to one aggregate diel pattern per species (Fig. 3), because the separate patterns in μ, σ and κ were strongly correlated (Fig. S2). Finally, we fitted to that aggregate diel pattern a mixture of a cosine and four wrapped normal distributions with 24 h periods, separately per species, with the wrapped normal distributions representing two peaks and two lows per day that were fixed at the peaks and lows of the aggregate diel patterns to capture the crepuscular movement activity of these species (Fig. 3). We defined this mixture distribution as our daily cyclic covariate (DCC). Thus, DCC had the following form:
Where tday is the time of day (range 0–2π), fWN is the probability density function of the wrapped normal distribution with a 24 h period (describing 1 out of 4 peaks or lows per day in the aggregate diel pattern), and δn are scalars for those distributions. Note that μn and σn here denote the mean and standard deviation of the peaks and lows in the aggregate diel pattern, and are not related to the mean and standard deviation of movement steps. As a result, each animal position had a corresponding DCC value, which depended on the time of day of that position and the species (Fig. 3).
By its design, DCC captures both the purely temporal (space-invariant) effect of time of day on animal movement and the combined spatial-temporal effect of diel variation in habitat selection, but it excludes the purely spatial (time-invariant) environmental effects. Using DCC in our analytic workflow therefore enables quantifying which fraction of the animals’ movement is tied directly or indirectly to time of day. If we were to quantify the influence of all three components in one HMM modelling step, without using the constructed DCC, that would not allow us to distinguish between variation caused by purely temporal (space-invariant) effects, purely spatial (time-invariant) effects and their combined spatial-temporal effects.
We standardized tree cover, slope and DCC to zero-mean and unit variance, and included all three variables and their two-way interactions as covariates of μ, σ and κ in the HMMs. We included tree cover, slope and DCC without interactions as covariates of the mode switching probabilities γ1→2.
Aggregate diel patterns averaged from μ, σ and κ for eland, wildebeest and zebra (solid line). We fitted a daily cyclic covariate (DCC, dashed line) to those aggregate diel patterns, which summarizes the diel activity patterns per species as expressed in the movement properties μ, σ and κ. DCC is a mixture distribution of a cosine and four wrapped normal distributions with 24 h periods
Model fitting
We fitted a series of HMMs per species separately. The HMM series consisted of models containing any combination of the three components of the movement process (modes, time and space), ranging from a null model that contained none, to a full model that contained all three in interaction (Table 1). The models were fit by numerical maximisation of the log-likelihood function, using the fitHMM routine from the R package momentuHMM [47]. We fitted the models in order of increasing complexity, every step adding covariates and/or interactions between covariates to the distribution parameters and switching probabilities. The full model included all three movement components modes, time and space, namely: DCC, tree cover and slope as covariates to the mode switching probabilities γ1→2, and these three variables and their two-way interactions as covariates to μ, σ and κ. This way, we could determine whether the covariates influenced movement within the separate movement modes, or through affecting the probability of switching to another mode, or both. The fitted parameter values from each fit were used as starting values for the next; starting values for other parameters were set to 0 (note that all covariates had been scaled to mean 0 and standard deviation 1). For each iteration, four sets of perturbations of those starting values β were sampled from \(\sim {\rm{N}}(\beta ,\sigma = 0.25 \cdot \beta )\). The model fit with the highest likelihood out of all starting sets was the optimal fit of that iteration, and it was used to initiate the model fit of the next iteration. However, in practice all iteration sets of fitted models converged to almost identical coefficient values and model likelihoods. We used the Akaike information criterion [AIC; 48] to assess model parsimony throughout the workflow. The model structure and fitting sequence were identical for all three species.
Variance partitioning of movement process components
Based on the fitted HMM with all predictors included (Model 3 in Table 1), we partitioned the explained variance of the modelled parameters μ, σ and κ for each species into contributions from the three components of the movement process: movement mode, time (the DCC covariate) and space (the covariates tree cover and slope). To quantify what proportion of the total variance in the predictions of μ, σ and κ was accounted for by each of the three components, we used linear regressions with the predicted values for μ, σ and κ as response variables (on the log-scale). The predicted values for μ, σ and κ were conditional on the most probable movement mode for each step along the movement path, which was determined using Viterbi decoding [49]. As predictors, we alternately omitted or included the movement mode component (i.e. the Viterbi path with the most likely hidden state sequence, included as predictor via a dummy variable), the time component (the DCC covariate), the space component (the tree cover and slope covariates) or their interactions as predictors and summarised the regression fit using R2 values. By doing so, we partitioned the variance of the fitted model parameters μ, σ and κ into the independent model components (i.e. mode, time and space), as well as their combined effects. All analyses were performed in R 3.6.1 [50].
Results
Hidden Markov models
We fitted Hidden Markov Models (HMMs) that included combinations of the three components of the movement process: modes, time and space (Table 1). When including the modes component (Models 1a, 2a, 2b and 3), we explored models with two and three movement modes. The full three-mode models achieved a higher model fit than the full two-mode models, albeit with a relatively small difference (ΔAIC = 32,093; 29,170 and 27,381 for eland, wildebeest and zebra respectively). However, AIC as a model selection criterium generally favours the inclusion of more modes, so the number of included modes should rather be informed by qualitative inspection of interpretability and biological meaningfulness of the modes’ properties [51]. Both the two- and three-mode models detected one faster, more directed mode and one slower, more tortuous mode (Fig. S3). We labelled these as ‘transit’ and ‘encamped’ modes, considering the labels to be reasonable abstractions while retaining spatially descriptive meaning and avoiding to project too specific behavioural annotations on the modes. The three-mode HMMs detected an additional slow and directional mode (Fig. S3). This third mode responded relatively weak to the temporal and spatial environmental covariates, while the response displayed by the other two modes was qualitatively similar for both the two- and three-mode HMMs (Figs. S4-S6). Therefore, we do not further consider the three-mode HMMs, and hereafter refer to the two-mode HMMs when discussing models that include movement modes.
The AIC values of all HMMs are presented in Table 1. The inclusion of movement modes led to much higher model fits, with the models including modes but not time and space (Model 1a) outperforming the models including time and space but not modes (Model 2c; ΔAIC = 63,185; 39,136 and 41,029 for eland, wildebeest and zebra respectively). Furthermore, the models including time but not space were superior to those including space but not time (Model 1c versus 1b; ΔAIC = 9,797; 14,153 and 5,507 for eland, wildebeest and zebra respectively; Model 2b versus 2a; ΔAIC = 9,335; 7,301 and 3,087 for eland, wildebeest and zebra respectively), indicating that DCC is a dominant covariate accounting for substantial variation, which is not accounted for by tree cover and slope. For each species the full model, which includes modes, time and space (Model 3), was the most parsimonious model according to AIC value.
Variance partitioning of movement process components
Partitioning the variance in the modelled parameters μ, σ and κ for each species between the modes, time and space components showed that modes explained most of the variation (see Fig. S7). The contributions of the model components to the within-mode variation in model parameters for the three species are shown in Fig. 4. For most combinations of species, parameter and modes, the temporal component explained most of the variance. However, the interactive effects of spatial-temporal predictors contributed notably to the variance in κ, especially for zebra (Fig. 4). By themselves, the spatial predictors did not explain much variation in movement parameters except for κ in zebra, and σ in eland (encamped mode) and zebra (transit mode; Fig. 4).
Variance partitioning, separately per movement parameter, species and movement mode (enc.= encamped, trans = transit). For variance partitioning that includes modes as a partition, see Fig. S7
Different model descriptions of the movement process
Incrementally expanding the predictor set according to Table 1 (from Model 0 to 3) yielded various interpretations of the movement process (Figs. S8-S14). Model 0, which included neither space, time nor mode as predictors, showed that wildebeest moved on average slowest and most tortuous, with eland and zebra moving equally directed and zebra fastest (Fig. S8). When we included movement modes (Model 1a), zebra showed the same speed as eland in the encamped mode and wildebeest moved faster than eland in the transit mode, instead of being slower overall (Fig. S9). Model 1b, which included only time as a predictor, showed large fluctuations in speed and tortuosity throughout the day, with both movement parameters peaking around sunrise and sunset and showing lows before sunrise and after sunset (Fig. S10). Model 1c, including only spatial predictors, showed that speed decreased with tree cover for eland, but increased for the other species (Fig. S11). Furthermore, speed decreased with slope irrespective of tree cover. Tortuosity decreased with tree cover but increased very slighty with slope. When we included both mode and time (Model 2a), we found the same temporal patterns in speed and tortuosity as the model including only time, but also indications that daily peaks in speed were much higher for animals in the transit mode than those in the encamped mode (Fig. S12). Model 2b, including both modes and space as predictors, showcased that often only in one of the two modes did animals respond substantially to spatial environmental heterogeneity (Fig. S13). For example, only in the transit modes did animals strongly decrease (eland) or increase (wildebeest and zebra) their speed. Zebra did not respond to slope, whereas the model including only spatial predictors suggested it slowed down on sloping terrain. Including space and time as predictors (Model 2c) reveiled that tree cover, slope and time of day interacted in their effect on speed and tortuosity, with their effects being most pronounced in the evening (Fig. S14). Model 3, including space, time and modes all as predictors, yielded the most complete description of movement, with aspects that were not apparent from all previous models (Fig. 5a, b). The movement properties displayed a strong diel pattern, with the effects of tree cover and slope on speed being largest in the evening when animals generally moved faster. Wildebeest and zebra in their transit mode moved faster with increasing tree cover whereas eland moved slower, but only so on steep terrain. In the encamped mode, tree cover had almost no effect on speed, but a strong effect on tortuosity. The animals also generally moved more tortuously with higher tree cover, even more so on sloping terrain. The model indicated that the spatial or temporal environment did not impact movement in all modes, and that the interaction effect between slope and tree cover on movement was also mode-specific.
Model predictions of movement speed μ as a result of tree cover and slope in the morning (04:00, top row), at noon (12:00, middle row) and in the evening (18:00, bottom row) for eland, wildebeest and zebra. The times were chosen to capture a dip, baseline and peak in DCC (see Fig. 3). Blue and red indicate the transit and encamped mode respectively. Colour mapping within blue and red indicates a progression of slope from 0 (light) till 0.15 (the 90th percentile of slope values in the complete dataset; dark). Small ticks at the inside of tree cover axes mark the 10th, 20th … 90th percentile values of tree cover in the respective data sets. Note that the y-axis scales differ between species, but not between the times of day. Confidence intervals have been omitted for visual clarity, but these were generally small
Idem to Fig. 5a, instead with model predictions of turning angle concentration on the y-axis. Concentration parameter κ (range 0-∞) has been converted to ρ (range 0–1) in this figure for visualisation purposes
Discussion
We developed an integrative analytic workflow to quantify and understand the roles of movement modes (‘why move’), temporal rhythms (‘when to move’) and spatial environmental heterogeneity (‘where to move’) in shaping animal movement. In this workflow, multi-modal (why) Hidden Markov Models (HMMs) are fitted on animal movement data, with spatial environmental covariates (where) and our constructed daily cyclic covariate (DCC; when, which aggregates the empirical diel patterns in the animals’ movement parameters) influencing both the movement parameters and the mode-switching probabilities. Finally, by using a variance partitioning approach on the modelled parameters, we quantified the relative importance of the why-, when- and where-components and their interactions in shaping animal movement.
We demonstrated the usefulness of our workflow by applying it on a case study of African ungulates. The variance partitioning showed that movement modes impacted movement the most, followed by diel rhythms and the spatial environment (i.e., tree cover and terrain slope) being the least important. Furthermore, the interactions between these components were important contributors as well, often explaining more of the movement variation than the marginal effect of the spatial environment did. This contrasts recent movement ecology research, where most studies seem to focus on the spatial environment without considering movement modes and temporal rhythms (Fig. 1). That is not to say that the spatial environment is the least important driver of animal movement across all study systems, or that research should focus solely on the most important movement driver of a system. However, when studying a specific component of interest of the movement process (e.g. a spatial environmental variable), that component should still be considered properly in context of the full (why, when, where) movement process [10], to avoid drawing misleading conclusions.
To showcase that misleading conclusions can follow from models in which one or more components of the movement process are ignored (‘omitted-variable bias’), we applied the same analytic pipeline to all combinations of the movement components in our case study. We identified two types of misleading conclusions that followed from models that only considered a subset of the movement process components: (1) failure to detect relationships and (2) detecting questionable relationships. First, the failure to detect relationships leads to overgeneralisations when drawing conclusions about the movement process (resembling a type II error). For example, wildebeest appeared to move slowest of the three species when not considering any movement process component in the model, but when movement modes were considered, the transit mode of wildebeest turned out to be faster than that of eland (Figs. S8 and S9 respectively). While this does not equate an erroneous conclusion (i.e., on average wildebeest did indeed move slowest), it does inhibit uncovering ecologically relevant complexities in the movement process [52]. Second, various questionable (likely spurious) relationships between input and response variables were revealed in models that included only a subset of all movement components that were present in the most parsimonious model (as determined by model AIC value; resembling a type I error). For example, zebra appeared to slow down substantially with increasing slope in the model with only spatial components, but with the addition of movement modes this effect disappeared completely for both modes (Figs. S11 and S13 respectively). Apparently, mode occupancy was a confounding variable that wrongly implied a causal relationship between slope and speed (i.e. zebras were more likely to be in their encamped mode in steep areas than they were in level areas; Fig. S15). Finally, besides these two types of misleading conclusions, the omission of a movement process component that is uncorrelated to other components does not lead to overgeneralisations or questionable relationships, but can nonetheless still lead to simply ‘incomplete’ conclusions. For example, a model with both movement modes and temporal predictors (Fig. S12) displayed only additive effects when compared to the models that solely contained movement modes (Fig. S9) or temporal predictors (Fig. S10).
An important consideration when using our analytic workflow, is that we developed DCC to capture both the direct, space-invariant effects of time of day on the animals’ movement properties, as well as the indirect effects of time of day through diel variation in habitat selection (given that varying environmental conditions between habitats also influence movement properties). For example, crepuscular animals often move faster during early morning because they forage before the midday heat sets in, but on top of that they also move faster at that time of day since they select for open foraging habitat which provides less movement resistance. Therefore, DCC includes both a purely temporal and a spatial-temporal component, which allows to separate those from purely spatial, time-invariant effects on the movement process (as these are captured by the spatial environmental, time-invariant covariates tree cover and slope). The influence of the temporal and spatial-temporal components together on movement was large compared to that of purely spatial environmental heterogeneity. Future research could focus on further disentangling these temporal and spatial-temporal components. For example, in the construction of the daily cyclic covariate, instead of solely averaging the movement parameters μ, σ and κ within each time bin, these movement parameters could be modelled based on time bin-specific spatial environmental covariates, after which the model predictions for e.g. the mean of the spatial environmental covariates could be used as the DCC value of that time bin. In this way, the purely temporal variation in movement could be separated from the temporally varying spatial influence on movement. Moreover, we see opportunities to build upon our workflow with datasets that are larger in spatial-temporal scale. Namely, seasonality (e.g. through climatic variation) should come into play on top of diel rhythms as a temporal covariate when the movement data spans more months, and extra spatial environmental covariates such as vegetation quality or water availability should be added when animals roam over larger or more heterogenous areas.
Conclusions
Our integrative analytic workflow can be used to analyse the effect and relative importance of all components (why, when and where to move) of the movement process. Compared to analyses that focus on a subset of these components, our workflow prevents the drawing of misleading conclusions through overgeneralisations and spurious correlations. Understanding the drivers of animal movement, and ultimately of ecological phenomena that emerge from it, critically depends on considering the various components of the movement process in concert, and especially the interactions between them.
Data availability
The data and script are provided in Additional file 3, and they are available in the 4TU.ResearchData repository at https://doi.org/10.4121/4ed0e7cd-50e7-4bc4-aae7-289fc0e38dc2.
Abbreviations
- AIC:
-
Aikaike information criterion
- DCC:
-
Daily cyclic covatiate
- HMM:
-
Hidden Markov model
References
Moorcroft PR, Lewis MA, Crabtree RL. Home range analysis using a mechanistic home range model. Ecology. 1999;80:1656–65.
Kie JG, Bowyer RT, Nicholson MC, Boroski BB, Loft ER. Landscape heterogeneity at differing scales: effects on spatial distribution of mule deer. Ecology. 2002;83:530–44.
Abrams PA. Habitat choice in predator-prey systems: spatial instability due to interacting adaptive movements. Am Nat. 2007;169:581–94.
Flaxman SM, Lou Y. Tracking prey or tracking the prey’s resource? Mechanisms of movement and optimal habitat selection by predators. J Theor Biol. 2009;256:187–200.
Bailey DW, Gross JE, Laca EA, Rittenhouse LR, Coughenour MB, Swift DM. et al. Mechanisms that result in large herbivore grazing distribution patterns. J Range Manag. 1996;49:386–400.
Revilla E, Wiegand T. Individual movement behavior, matrix heterogeneity, and the dynamics of spatially structured populations. PNAS. 2008;105:19120–5.
Morales JM, Moorcroft PR, Matthiopoulos J, Frair JL, Kie JG, Powell RA. Building the bridge between animal movement and population dynamics. Philosophical Trans Royal Soc B: Biol Sci. 2010;365:2289–301.
Schick RS, Loarie SR, Colchero F, Best BD, Boustany A, Conde DA. et al. Understanding movement data and movement processes: current and emerging directions. Ecol Lett. 2008;11:1338–50.
Riotte-Lambert L, Benhamou S, Chamaillé-Jammes S. Periodicity analysis of movement recursions. J Theor Biol. 2013;317:238–43.
Nathan R, Getz WM, Revilla E, Holyoak M, Kadmon R, Saltz D. et al. A movement ecology paradigm for unifying organismal movement research. PNAS. 2008;105:19052–9.
Morales JM, Haydon DT, Frair J, Holsinger KE, Fryxell JM. Extracting more out of relocation data: building movement models as mixtures of random walks. Ecology. 2004;85:2436–45.
McClintock BT, King R, Thomas L, Matthiopoulos J, McConnell BJ, Morales JM. A general discrete-time modeling framework for animal movement using multistate random walks. Ecol Monogr. 2012;82:335–49.
McClintock BT, London JM, Cameron MF, Boveng PL. Bridging the gaps in animal movement: hidden behaviors and ecological relationships revealed by integrated data streams. Ecosphere. 2017;8:e01751.
Forester JD, Ives AR, Turner MG, Anderson DP, Fortin D, Beyer HL. et al. State–space models link elk movement patterns to landscape characteristics in Yellowstone national park. Ecol Monogr. 2007;77:285–99.
Mews S, Elkenkamp S, Schuhmann P, Tsolak D, Wobbe T. et al. Classifying grey seal behaviour in relation to environmental variability and commercial fishing activity - a multivariate hidden Markov model. Sci Rep. 2019;9:5642.
Jonsen ID, Myers RA, James MC. Robust hierarchical state–space models reveal diel variation in travel rates of migrating leatherback turtles. J Anim Ecol. 2006;75:1046–57.
Dingle H, Drake VA. What is migration? BioScience. 2007;57:113–21.
Dingle H. Seasonal migration patterns. In: Dingle H, editor. Migration: the biology of life on the move. Oxford University Press; 2014. pp. 163–82.
Polansky L, Wittemyer G, Cross PC, Tambling CJ, Getz WM. From moonlight to movement and synchronized randomness: fourier and wavelet analyses of animal location time series data. Ecology. 2010;91:1506–18.
Riotte-Lambert L, Matthiopoulos J. Environmental predictability as a cause and consequence of animal movement. Trends Ecol Evol. 2020;35:163–74.
Green RA, Bear GD. Seasonal cycles and daily activity patterns of rocky mountain elk. J Wildl Manag. 1990;54:272–9.
Ensing EP, Ciuti S, Wijs FALM de, Lentferink DH, Hoedt A ten, Boyce MS, et al. GPS based daily activity patterns in European red deer and North American elk (Cervus elaphus): indication for a weak Circadian clock in ungulates. PLOS ONE. 2014;9:e106997.
Owen-Smith N, Goodall V. Coping with savanna seasonality: comparative daily activity patterns of African ungulates as revealed by GPS telemetry. J Zool. 2014;293:181–91.
Vazquez C, Rowcliffe JM, Spoelstra K, Jansen PA. Comparing diel activity patterns of wildlife across latitudes and seasons: time transformations using day length. Methods Ecol Evol. 2019;10:2057–66.
Buderman FE, Hooten MB, Alldredge MW, Hanks EM, Ivan JS. Time-varying predatory behavior is primary predictor of fine-scale movement of wildland-urban cougars. Mov Ecol. 2018;6:22.
Dussault C, Ouellet J-P, Courtois R, Huot J, Breton L, Larochelle J. Behavioural responses of moose to thermal conditions in the boreal forest. Écoscience. 2004;11:321–8.
Fischhoff IR, Sundaresan SR, Cordingley J, Rubenstein D. Habitat use and movements of plains zebra (Equus burchelli) in response to predation danger from lions. Behav Ecol. 2007;18:725–9.
Owen-Smith N, Traill LW. Space use patterns of a large mammalian herbivore distinguished by activity state: fear versus food? J Zool. 2017;303:281–90.
Skarin A, Danell Ö, Bergström R, Moen J. Reindeer movement patterns in alpine summer ranges. Polar Biol. 2010;33:1263–75.
Li M, Bolker BM. Incorporating periodic variability in hidden Markov models for animal movement. Mov Ecol. 2017;5:1.
Hooten MB, Hanks EM, Johnson DS, Alldredge MW. Temporal variation and scale in movement-based resource selection functions. Stat Methodol. 2014;17:82–98.
Hooten MB, Johnson DS. Basis function models for animal movement. J Am Stat Assoc. 2017;112:578–89.
Hooten MB, Scharf HR, Morales JM. Running on empty: recharge dynamics from animal movement data. Ecol Lett. 2019;22:377–89.
Dailey TV, Hobbs NT. Travel in alpine terrain: energy expenditures for locomotion by mountain goats and bighorn sheep. Can J Zool. 1989;67:2368–75.
Parker KL, Robbins CT, Hanley TA. Energy expenditures for locomotion by mule deer and elk. J Wildl Manag. 1984;48:474–88.
Wall J, Douglas-Hamilton I, Vollrath F. Elephants avoid costly mountaineering. Curr Biol. 2006;16:R527–9.
White RG, Yousef MK. Energy expenditure in reindeer walking on roads and on tundra. Can J Zool. 1978;56:215–23.
Leblond M, Dussault C, Ouellet J-P. What drives fine-scale movements of large herbivores? A case study using moose. Ecography. 2010;33:1102–12.
Ager AA, Johnson BK, Kern JW, Kie JG. Daily and seasonal movements and habitat use by female rocky mountain elk and mule deer. J Mammal. 2003;84:1076–88.
Lubitz N, Bradley M, Sheaves M, Hammerschlag N, Daly R, Barnett A. The role of context in elucidating drivers of animal movement. Ecol Evol. 2022;12:e9128.
Eikelboom JAJ, de Knegt HJ, Klaver M, van Langevelde F, van der Wal T, Prins HHT. Inferring an animal’s environment through biologging: quantifying the environmental influence on animal movement. Mov Ecol. 2020;8:40.
Joo R, Picardi S, Boone ME, Clay TA, Patrick SC, Romero-Romero VS. et al. Recent trends in movement ecology of animals and human mobility. Mov Ecol. 2022;10:26.
McClintock BT, Langrock R, Gimenez O, Cam E, Borchers DL, Glennie R. et al. Uncovering ecological state dynamics with hidden Markov models. Ecol Lett. 2020;23:1878–903.
de Knegt HJ, Eikelboom JAJ, van Langevelde F, Spruyt WF, Prins HHT. Timely poacher detection and localization using sentinel animal movement. Sci Rep. 2021;11:4596.
Mucina L, Rutherford MC, Powrie L, AG R, KGT C, Lötter MC. Vegetation Atlas of South Africa, Lesotho and Swaziland. The Vegetation of South Africa, Lesotho and Swaziland.Pretoria: South African National Biodiversity Institute; 2006. p.748–90.
Congedo L. Semi-automatic classification plugin: a Python tool for the download and processing of remote sensing images in QGIS. J Open Source Softw. 2021;6:3172.
McClintock BT, Michelot T, momentuHMM. R package for generalized hidden Markov models of animal movement. Methods Ecol Evol. 2018;9:1518–30.
Akaike H. Information theory and an extension of the maximum likelihood principle. In: Petrov BN, Csaki BF, editors. Second International Symposium on Information Theory. Akadémiai Kiadó; 1973. p. 267–81.
Zucchini W, Macdonald I. Hidden Markov models for time series: an introduction using R. 2009.
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2023 [cited 2022 Jan 7]. Available from: https://cran.r-project.org/bin/windows/base/old/3.6.1/
Pohle J, Langrock R, van Beest FM, Schmidt NM. Selecting the number of states in hidden Markov models: pragmatic solutions illustrated using animal movement. JABES. 2017;22:270–93.
Spake R, Bowler DE, Callaghan CT, Blowes SA, Doncaster CP, Antão LH. et al. Understanding ‘it depends’ in ecology: a guide to hypothesising, visualising and interpreting statistical interactions. Biol Rev. 2023;98:983–1002.
Acknowledgements
We thank Welgevonden Game Reserve for the animal monitoring, MTN South Africa for the telecommunication and IBM for the data warehousing. We thank two reviewers for providing constructive feedback on earlier versions of the manuscript.
Funding
This research was funded by the Dutch Research Council (NWO program ‘Advanced Instrumentation for Wildlife Protection’).
Author information
Authors and Affiliations
Contributions
H.L., H.J.d.K. and J.A.J.E. conceived the the ideas and methodology for the study; H.J.d.K. and J.A.J.E. provided the data; H.L. and H.J.d.K. analysed the data; H.L. led the writing of the manuscript. All authors provided critical feedback on the manuscript and approved the final version for publication.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Capturing and tagging animals was approved by Welgevonden Game Reserve as a management action and carried out in accordance with relevant guidelines and regulations under the National Environmental Management Protected Areas Act (NEM: PAA; Act No. 57 of 2003; see also [44] and Additional file 2). All people involved in this study consent to participate in its publication.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Supplementary Material 1: Additional file 1.
Systematic literature search queries for Fig. 1 and supplementary figures.
Supplementary Material 2: Additional file 2.
Welgevonden Game Reserve GPS collaring approval letter
Supplementary Material 3: Additional file 3.
Code and data for reproduction of analysis and results
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Linssen, H., de Knegt, H.J. & Eikelboom, J.A. Unveiling the roles of temporal periodicity, the spatial environment and behavioural modes in terrestrial animal movement. Mov Ecol 12, 57 (2024). https://doi.org/10.1186/s40462-024-00489-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40462-024-00489-3