
==== Front
Mov Ecol
Mov Ecol
Movement Ecology
2051-3933
BioMed Central London

39238061
502
10.1186/s40462-024-00502-9
Research
Effects of the social environment on movement-integrated habitat selection
Webber Quinn qwebber@uoguelph.ca

12
Prokopenko Christina 3
Kingdon Katrien 3
Turner Julie 3
Vander Wal Eric 13
1 https://ror.org/04haebc03 grid.25055.37 0000 0000 9130 6822 Cognitive and Behavioural Ecology Interdisciplinary Program, Memorial University of Newfoundland, St. John’s, NL Canada
2 https://ror.org/01r7awg59 grid.34429.38 0000 0004 1936 8198 Present Address: Department of Integrative Biology, University of Guelph, Guelph, ON Canada
3 https://ror.org/04haebc03 grid.25055.37 0000 0000 9130 6822 Department of Biology, Memorial University of Newfoundland, St. John’s, NL Canada
5 9 2024
5 9 2024
2024
12 6110 5 2023
20 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ 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/.
Background

Movement links the distribution of habitats with the social environment of animals using those habitats. Despite the links between movement, habitat selection, and socioecology, their integration remains a challenge due to lack of shared vocabulary across fields, methodological gaps, and the implicit (rather than explicit) historical development of theory in the fields of social and spatial ecology. Given these challenges can be addressed, opportunity for further study will provide insight about the links between social, spatial, and movement ecology. Here, our objective was to disentangle the roles of habitat selection and social association as drivers of movement in caribou (Rangifer tarandus).

Methods

To accomplish our objective, we modelled the relationship between collective movement and selection of foraging habitats using socially informed integrated step selection function (iSSF). Using iSSF, we modelled the effect of social processes, i.e., nearest neighbour distance and social preference, and movement behaviour on patterns of habitat selection.

Results

By unifying social network analysis with iSSF, we identified movement-dependent social association, where individuals took shorter steps in lichen habitat and foraged in close proximity to more familiar individuals.

Conclusions

Our study demonstrates that social preference is context-dependent based on habitat selection and foraging behaviour. We therefore surmise that habitat selection and social association are drivers of collective movement, such that movement is the glue between habitat selection and social association. Here, we put these concepts into practice to demonstrate that movement is the glue connecting individual habitat selection to the social environment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40462-024-00502-9.

Keywords

Caribou
Integrated step selection analysis
Movement ecology
Social preference
Social network analysis
http://dx.doi.org/10.13039/501100000038 Natural Sciences and Engineering Research Council of Canada issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcBackground

Movement is defined by a change in spatial location and is the behavioural link between the physical space an animal occupies and the conditions, risks, and resources that exist in their environment [1, 2]. In the context of the social environment, which is defined as the nature, quality, and patterning of relationships among members of a group or population [3], movement represents the connection between the distribution of resources and the social structure of animals that consume those resources [4]. Disentangling the social and spatial drivers of movement is a formidable challenge within behavioural ecology. In many cases, research omits the social contexts within which animals move to, from, and within the areas that contain foraging resources [5, 6]. Spatially-explicit models of sociality highlight that some gregarious species aggregate at areas associated with profitable foraging resources [7], whereas some territorial species only interact at territory edges [8]. Sharing space, either at foraging sites, territory edges, or elsewhere within an animal’s range is required to form the social environment. For example, animals are predicted to select habitat as a function of the profitability and availability of the habitat [9]. A logical extension can be made to conspecifics; individuals form groups based on their familiarity with conspecifics and the profitability of associating with familiar conspecifics. We aim to quantify the relative importance of habitat and conspecifics by developing a socially informed integrated step selection analysis, a movement-based method that accounts for the relative intensity of selection for habitats and social neighbours.

For social animals, individual movement shapes social encounters and subsequent interactions with conspecifics that can affect consensus decision making associated with collective movement [10, 11]. Further complicating our understanding of collective movement is the idea that the type, quality, and distribution of habitats on the landscape can constrain or promote collective movement [12]. For example, dense vegetation impedes visibility, which could influence movement of individuals and reduce the probability a group remains together [13, 14]. In addition, individual movement and habitat selection are affected by the distribution of resources. For example, patchily distributed foraging resources could facilitate large aggregations, whereas homogenously distributed foraging resources could result in a reduction in social associations [5]. While the physical space an individual, or group, occupies is an important driver of animal movement [4], the familiarity of animals within a social group can also affect how animals move through the landscape [11].

Not all social groups are equal; some groups contain unfamiliar individuals (i.e., anonymous groups) [15], while others contain familiar individuals [16]. For anonymous and familiar groups, social foraging occurs when the costs (e.g. competition) and benefits (e.g. reduced predation risk) of an individual’s foraging behaviour are linked with the foraging behaviour of conspecifics [17]. Social foraging can be most beneficial when social information about foraging resources comes from familiar individuals [18]. For example, when foraging resources are unpredictable, familiar individuals obtain reliable information from conspecifics to increase foraging efficiency [19, 20], such that time searching for forage is reduced in favour of more time spent foraging. In the context of movement and habitat selection, theory on social foraging and the benefits of social familiarity provides a framework through which the costs and benefits of collective movement can be explored [17, 21].

Apparent social familiarity or preference is the long-term repeated social association due to shared space at the same time. Although individuals often associate with many conspecifics, non-random repeated social associations with certain individuals form the basis for social preference [22]. Proximately, long-term social relationships can influence collective movement via the reliability of information transfer about foraging resources or predator risk [23, 24], while ultimately they can enhance fitness [25]. The social environment can be influenced by the availability of foraging resources, but social groups can also be composed of individuals with similar physiological or nutritional requirements that occupy the same locations. Apparent social preference may therefore arise as a function of spatial constraints [26], including physical barriers, such as rivers or mountains. Disentangling social preference from spatial constraint could inform our understanding of collective movement and habitat selection [27, 28].

Here, we use contemporary tools from movement ecology and social network analysis and integrate aspects of socioecology, spatial ecology, and movement ecology. We test the effects of the social environment on patterns of habitat selection behaviour by parameterizing socially informed integrated step selection models (Fig. 1). Using caribou (Rangifer tarandus) as a model system, our over-arching goal was to disentangle the roles of habitat selection and social association as drivers of collective movement, when the availability and distribution of foraging resources are variable. To accomplish our goals, we had an over-arching objective to test how social information influences patterns of selection for foraging resources. Social information was measured using proximity to conspecifics (i.e., nearest neighbour distance) and time-lagged social association as proxies for information flow. Objective 1 was to quantify patterns of habitat selection in the presence of conspecifics. We incorporated proximity to conspecifics (i.e., nearest neighbour distance) as a proxy for social information transfer, such that individuals near one another are more likely to exchange information than individuals that are further away from one another (O1). We therefore expect that individuals will be more likely to select foraging habitat (lichen) when conspecifics are near rather than far. Objective 2 was to quantify patterns of habitat selection in the presence of familiar conspecifics. Here, we incorporated time-lagged social association as a proxy for reliable information flow by decomposing social relationships into familiar and unfamiliar contexts. Specifically, we expect that individuals with stronger short-term social preference should collectively select foraging habitat, which tends to be more open and therefore facilitate grouping (O2). The corollary is that individuals should also take short steps in the presence of conspecifics. For O2 we assessed short-term time-lagged social preference using weekly measures of association and tested whether recency bias of association affected patterns of foraging habitat selection. O2 relies on the assumption that from a movement ecology perspective, shorter steps typically represent foraging behaviour and longer steps represent searching behaviour [29].Fig. 1 Summary of the data pipeline used to generate integrated step selection function (iSSF) models. Primary data sources were landcover data and caribou GPS relocation data, which were combined to determine the physical locations of animals on the landscape. The pairing of animal locations and landcover data was used to generate the comparison of used to available points (panel a), which is the response variable in iSSF models, as well as the habitat type in which a given relocation occurred: lichen (defined in text as open-forage) and forest (panel b). Caribou relocation data were also used to generate two movement parameters (panel c) and aspects of the social environment (panels d and e). Our movement parameter included step length, which is the linear distance between consecutive relocations. The social environment included nearest neighbour distance (panel d) and weekly social networks and the dyadic SRI generated based on a moving-window as a proxy for short-term social preference (panel e). The bottom row represents a graphical formulation of our iSSF models, where habitat selection (1:20 ratio of used to available relocations) was regressed against habitat type (lichen and forest), movement parameters (step length and turn angle), nearest neighbour distance, and weekly dyadic SRI

Methods

Caribou as a model system

We investigated patterns of movement, space use, and social behaviour for caribou on Fogo Island, Newfoundland, Canada. Fogo Island is a small (~ 237km2) island situated ~ 12 km off the northeastern coast of Newfoundland with a humid continental climate. The dominant habitat types on Fogo Island consist of coniferous and mixed forests of balsam fir (Abies balsamea), black spruce (Picea mariana) and white birch (Betula papyrifera), as well as bogs, lakes, and barren rock. Between 1964 and 1967, 26 caribou were introduced to Fogo Island from the Island of Newfoundland [30]. Currently, Fogo Island has a population of approximately 300 caribou and unlike many of the other herds in Newfoundland [31], caribou on Fogo Island have a relatively stable population that has not declined in recent years. In addition, caribou on Fogo Island are sedentary and do not display any migratory or long-distance movements.

Caribou live in fission–fusion societies [32], and throughout much of their range, caribou forage primarily on lichen, grasses, sedges, and other deciduous browse with access to these resources changing between the seasons [33]. During winter (January to March), the landscape is covered by snow, and caribou forage primarily on lichen [34]. Lichen is heterogeneously distributed, and access is impeded by snow and ice cover. Caribou dig holes in the snow, termed craters, to access lichen in the winter, often where snow depth is relatively shallow (~ 30–60 cm deep). Consequently, caribou have limited access to lichen buried under the snow and tend to re-use established craters. To cope with this limitation, caribou use conspecific attraction and social information transfer to gain access to foraging opportunities [35]. In winter, caribou activity budgets suggest that caribou spend ~ 50% of their time foraging, while ~ 40% of their time is spent lying down or ruminating, 7% of their time is spent walking or trotting, and 3% of their time is spent standing [36, 37]. In addition, caribou typically avoid forested habitats due to deep snow in forests and lack of access to forage opportunities [38], whereas most open habitats on Fogo Island are windswept in the winter, facilitating foraging and movement [33].

We used GPS location data collected from Fogo Island caribou (2017–2019) to assess the relationship between social behaviour, habitat selection, and movement (see supplementary information for details on collaring procedures). For all analyses, we restricted locations to only include relocations from the first 75 days of each year (1 January–16 March). Each relocation was assigned to a given habitat classification that was extracted from Landsat images with 30m × 30m pixels [39]. We pooled habitats into three categories. Open habitats that caribou move through consisted of wetland (21.2% of habitats on Fogo Island), rocky outcrops (6.7%), and water/ice (9.2%) habitats. Forest habitats consisted of conifer scrub (27.0%), conifer forest (10.4%), mixed-wood (9.2%), and broadleaf forest (0.2% of habitats) habitats. Finally, we considered lichen barrens (11.7%) as caribou foraging habitat. Overall, the proportion of GPS relocations (see below) that fell within each habitat were: 54.2% of GPS relocations were in open habitats, 30.4% of GPS relocations were in lichen habitat, 10.8% of GPS relocations were in forested habitats, and 4.4% of locations were in undesignated habitats. We then calculated the proportion of each habitat type (i.e., open foraging, open moving, or forest) within 200 m around each used and available point location (see below).

Adult female caribou (n = 26 individual caribou, n = 76 caribou-years) were immobilized and fitted with global positioning system (GPS) collars (Lotek Wireless Inc., Newmarket, ON, Canada, GPS4400M collars, 1250 g). Prior to analyses, we removed all erroneous and outlier GPS locations following Bjørneraas et al. [40]. Collars collected data throughout the year and were programmed to collect locations every two hours. Specifically, we removed implausible GPS fixes, including those that were recorded in the ocean and those that exceeded the maximum diameter of Fogo Island (30 km). We removed individuals with collar failure during the study period (n = 16 caribou-years), individuals that swam to nearby adjacent islands (n = 3 caribou-years), mortality (n = 3), or individuals without data in the winter (n = 16). Given that habitat selection and kernel analyses can be biased by fix success rates [41, 42], we removed individuals with high rates of collar failure (i.e., where > 20% of fixes were missing), resulting in a total of 21 individual caribou (n = 38 caribou-years). Overall, average relocation success rate for the 21 individuals included in subsequent analyses was 98.8%. We did not collar all female caribou in the herd, however, and collared individuals were randomly selected from the population. We therefore assume that our sample of collared animals was randomly distributed. Although associations between collared and uncollared animals were unrecorded, we assumed that our networks (see below) were unbiased representations of the relative degree of social association among all caribou.

Formulating integrated step selection models

To accomplish Objective 2, we formulated integrated step selection functions (iSSF); a type of movement model that simultaneously incorporates movement and habitat selection within a conditional logistic regression framework (Fig. 1) [43–45]. As in other resource and step selection analyses [46], iSSF models habitat selection as a binomial response variable where ‘use’ represents the location an animal was observed and ‘availability’ represents the geographical area an animal could potentially use but was not necessarily observed (Fig. 2). iSSF defines availability based on probability distributions of step lengths and turn angles [43], where a step is the linear connection between consecutive relocations, and turn angle is the angular deviation between the headings of two consecutive steps [47]. We generated available steps and turn angles based on the distributions informed by observed population-level movement behaviour using the amt package in R [48]. First, we sampled step lengths from a gamma distribution of observed step lengths for the study population; values were log-transformed for analysis. The statistical coefficient of log-transformed step length is a modifier of the shape parameter from the gamma distribution originally used to generate available steps [43]. Second, we sampled turn angles (measured in radians) for available steps from observed values between -π and π following a Von Mises distribution. Each observed relocation was paired through a shared start point with 20 available steps generated from step-length and turn-angle distributions and compared in a conditional logistic regression framework. In addition to generating available movement parameters, we also generated an available social environment (see below). To evaluate the predictive performance of our model, we used k-fold (k = 5) cross validation [49] following the methods of Fortin et al. [50]. For details on k-fold cross validation see Appendix 2.Fig. 2 Schematic of integrated step selection functions in the context of social network analysis. Available (random) steps are generated based on the distribution of used step lengths (thin dashed orange and blue lines) and turn angles. We compared used (observed) to available (random) nearest neighbour distance based on fine-scale movement decisions of individuals. Blue lines represent used (dark thick lines) and available (light dashed lines) steps of individual i and orange lines represent used (dark thick lines) and available (light dashed lines) steps of individual j. The dashed grey line (NNi,j) represents the observed nearest neighbour distance between i and j at t3. For each set of available steps, we re-calculated nearest neighbour distance, denoted by a dashed grey line (NNi,j′) between available steps for i and j, which represents the available nearest neighbour distance at a given iteration. Step length is the distance between the used start (e.g. tj,2) point and the step end point (tj,3) and turn angle is calculated as the angular deviation between the previous step heading (grey line and θi and θj) and the subsequent used step

Social network analysis

To accomplish Objective 1, we generated proximity-based social association networks. Nodes in the networks represented individual caribou and edges represented the frequency of association based on proximity between individuals. We generated social networks based on proximity of locations between individual caribou during each week throughout the winter to assess the role of short-term social preference on patterns of habitat selection (O1). For weekly networks (see below for details), we assumed association between two individuals when simultaneous locations (i.e. GPS relocations that occurred within 5 min of each other) were within 50 m of one another [32, 35]. We selected the 50 m threshold based on the standard distance applied to assign individuals to groups in studies of ungulate group size and social behaviour [51] and based on sensitivity analyses in our system. We deem 50 m as an appropriate threshold for caribou as lower thresholds yield relatively sparse networks and higher thresholds yield saturated networks [52]. We applied the ‘chain rule’ to define groups, where each discrete GPS fix was buffered by 50 m and we considered individuals in the same group if 50 m buffers for two or more individuals were contiguous. In some cases, the distance between discrete fixes of two individuals exceeded 50 m because contiguous buffers overlapped with at least one other central individual to form a ‘chain’ of buffers [51]. Groups were defined based on the gambit-of-the-group assumption, which was used to construct social networks [53]. We weighted edges of social networks by the strength of association between dyads of caribou using a dataset specific simplification of the simple ratio index [54], SRI:SRI=xx+yAB

where x is the number of times individuals A and B were within 50 m of each other and yAB is the number of simultaneous fixes from individuals A and B that were separated by > 50 m [55]. To generate networks, we used GPS location data using the R [56] packages spatsoc [57] and igraph [58].

Weekly networks

To accomplish Objective 1, we iteratively generated weekly social networks using a moving window approach and calculated the observed SRI to be included as a covariate in our iSSF model (see formulating iSSF section). The first network was calculated for 1 January to 7 January, the second was 2 January to 8 January, and so on. Weekly networks contained 84 relocations per individual (12 relocations per day). For each of these networks, we used dyadic values of SRI as a proxy for short-term social preference. We used a three-step process. First, to incorporate SRI within the iSSF framework, we determined the identity and distance (m) of each individual’s nearest neighbour at each relocation. Second, for each focal individual and their nearest neighbour at each relocation, we matched the time-lagged dyadic SRI value for the prior week. For example, for individual A at 12:00 on 8 January, we determined the nearest neighbour was individual B and we extracted the dyadic SRI value for these individuals for the previous week. Third, we repeated steps one (i.e., calculate the identity and distance to nearest neighbours) and two (i.e., match the time-lagged dyadic SRI to nearest neighbours) for each observed relocation and the accompanying set of ‘available’ relocations defined by random steps generated in the iSSF. Therefore, each individual at each relocation had one observed weekly dyadic SRI value and 20 ‘available’ weekly dyadic SRI values (see dashed lines on Fig. 2 for example of available steps and how these steps are related to nearest neighbour distance).

Modelling collective movement and habitat selection

To accomplish Objective 2, we fit three iSSF models with a series of fixed and random effects following Muff et al. [62] (see Fig. 1 for iSSF breakdown). When fitting iSSFs, the model is formulated by incorporating a selection-free movement kernel, which considers step length and turn angle distributions, and a selection kernel, which considers habitat types [59]. Parameters in both the movement and selection kernels can be simultaneously estimated and interactions between movement and habitat types can be included [59]. Importantly, an iSSF produces a mechanistic movement model that is empirically parameterized [43, 48].

We took advantage of the fact that the conditional logistic regression model is a likelihood-equivalent to a Poisson model with stratum‐specific fixed intercepts. The approach outlined by Muff et al. [62] uses a mixed modelling approach which allows intercepts and/or slopes to vary by individual, while also incorporating shared information that is present in the data from different individuals [60]. For social species that may move collectively, and therefore have correlated movement trajectories, varying intercepts by individual is recommended to account for correlation within nested groupings of locations [61]. However, we are unable to nest individual within groups in our models because caribou live in fission–fusion societies [32], where group membership changes through space and time.

Following Muff et al. [62], all variables included in the fixed effect structure were also included in the random effect structure. Leveraging the method described by Muff and colleagues allows us to quantify population level effects of covariates (fixed effects) while controlling for individual variation in each covariate (random effects), which ultimately ensures fixed effect coefficients are more robust than previously described two step approaches [62]. For models to yield unbiased results, each fixed effect must also be parameterized as a random slope at the individual animal level. With GPS data, the quantity of data is higher (i.e., ~ 900 observations per animal) than most observational studies over a similar time-period. The quantity of data therefore allows the development of more complex models. Each fixed effect and interaction term has biological relevance for our objectives, indicating that a single large model is more effective for the comparison of effect sizes than many smaller models. From a single model, we test the importance of covariates relative to one another, therefore enabling us to test our objectives related to social foraging.

Our initial model included the proportion of lichen, forest, and open habitat within 200 m of the point location, the natural log-transformed step length, natural log-transformed nearest neighbour distance, and log-transformed time-lagged weekly dyadic SRI (O1, see Fig. 1 for model visualization). Given high collinearity among habitat types, we removed open habitat from our model for all subsequent analyses (Table S2). We natural log-transformed step length, nearest neighbour distance, and weekly dyadic SRI to account for potential decay function in each variable, to scale variables, and to ensure model convergence. We then generated three follow up models. First, we fit a null model that only included log-transformed step length, lichen habitat, and forest habitat within 200 m of the point location as well as interactions between step length and both habitat types (Table S3). Second, we fit a model that include the same terms as model 1, but also included turn angle in an interaction with both habitat types (Table S3). Third, we fit a model that included social variables: nearest neighbour distance (m), which was measured as the distance between a focal individual and the nearest collared conspecific and was calculated for all used and available steps and time-lagged weekly dyadic SRI. In our third model, we also included interactions between: 1) step length and each of the proportion of lichen and forest habitats within 200 m of the point location, 2) nearest neighbour distance and step length, 3) nearest neighbour distance and each of the proportion of lichen and forest habitats within 200 m of the point location as separate interactions, 4) log-transformed SRI and each of the proportion of lichen and forest habitats within 200 m of the point location. For interactions that included nearest neighbour distance, we used either distance at the start of a step or at the end of the step, depending on the other variable in the interaction (Fig. 2). Specifically, for the interaction between step length and nearest neighbour distance (O2), we used distance at the start of the step because the likelihood of taking a shorter or longer step is predicted to vary based on the distance to conspecifics before the step is taken. By contrast, for interactions between habitat variables and nearest neighbour distance, we used distance at the end of the step because the likelihood of selecting a given habitat is predicted to vary based on the distance to conspecifics when that habitat is being selected, i.e., at the end of the step. To assess model fit, we compared the null model to the model with social variables using AIC [63]. We fit our iSSF model using the glmmTMB package in R [64].

Calculating effect sizes

We calculated individual-level relative selection strength (RSS) to demonstrate how habitat features influenced selection [65]. Our model was therefore interpreted using RSS, which is the biological effect size extracted from conditional logistic models [65]. Interpretation of RSS provides information about how a covariate is selected or avoided across an environmental gradient of what is available to the individual. We calculated the strength for selecting one step over another that differed in the habitat value where those steps ended. RSS was calculated for each habitat type (i.e., forest or lichen habitats) as a function of nearest neighbour distance and the shared dyadic SRI between nearest neighbours.

Results

Social networks were generally well connected such that most individuals associated with most other individuals at least once during the winter (Fig.’s S1, S2, and S3). On average, graph density (i.e. the proportion of realized pairwise connections in each network) was 0.83 (2017 = 1.00; 2018 = 0.78; 2019 = 0.71). Moreover, individuals tended to cluster geographically, suggesting a possible spatial driver of social network structure (Fig.’s ,S1, S2, and S3yadic SRI values for weekly networks were skewed to zero with a long tail (median across all networks [± SD]: 0.03 ± 0.05), indicating that even though most individuals were connected at least once throughout the winter, it is clear some dyads consistently demonstrated social preference for one another at shorter timeframes.

We found evidence for an effect of the social environment on patterns of habitat selection (Table 1; Fig. 3). Specifically, our iSSF model which included social variables was considered a better fit than our null models (ΔAIC = 19,000, Table S3). Given the social nature of caribou, the effect of the social environment on selection was nuanced, and we found partial support for our expectation of social foraging. Overall, individuals took longer steps than expected (0.34, 95% CI: 0.18, 0.50; Table 1) and tended to moved more slowly when selecting lichen (– 0.38, 95% CI: – 0.45, − 0.31) and forest habitat (− 0.36, 95% CI: − 0.43, − 0.31) (Fig. 4; Table 1). Meanwhile, individuals tend to be closer to their nearest neighbours when selecting lichen habitat (0.55, 95% CI: 0.46, 0.65; Table 1) and relative selection strength for lichen decayed as nearest neighbours were further away (Fig. 5). Though the interaction between forest habitat and nearest neighbour distance was non-significant (− 0.01, 95% CI: − 0.08, 0.11; Table 1), relative selection strength for forest habitat decayed as nearest neighbours were further away (Fig. 5). Finally, individuals tended to share a higher dyadic SRI value with their nearest neighbours when selecting lichen habitat (− 4.03, 95% CI: − 4.64, 3.43), relative to the availability of lichen habitat (Table 1), while there was no relationship between forest habitat and SRI (− 0.21, 95% CI: − 0.78, 0.36; Table 1). Taken together, relative selection for lichen habitat was stronger than forest habitats as a function of nearest neighbour distance and SRI (Fig. 5). Our k-fold cross-validation had high scores (rho = 0.81 SE ± 0.03), demonstrating our model was better than random at predicting where caribou moved (see Fig. S4 to view coefficients and error for variables in each fold). Table 1 Summary of integrated step selection function. We considered covariates significant if the 95% confidence intervals did not overlap zero

Covariate	ß (95% CI)	z-value	p-value	
Intercept	6.94 (− 3.77, 17.6)	1.27	0.20	
Step length	0.34 (0.18, 0.50)	4.24	 < 0.0001	
Proportion Forest	0.61 (− 0.43, 1.66)	1.14	0.25	
Proportion Lichen	− 6.54 (− 7.67, − 5.42)	− 11.4	 < 0.0001	
Nearest neighbour distance (end)	− 0.75 (− 0.82, − 0.68)	− 20.9	 < 0.0001	
Simple ratio index	4.62 (4.13, 5.11)	18.5	 < 0.0001	
Proportion Forest: Step length	− 0.36 (− 0.43, − 0.31)	− 12.0	 < 0.0001	
Proportion Lichen: Step length	− 0.38 (− 0.45, − 0.31)	− 11.4	 < 0.0001	
Step length: Nearest neighbour distance (start)	− 0.013 (− 0.04, 0.012)	− 1.02	0.31	
Proportion Forest: Nearest neighbour distance (end)	− 0.01 (− 0.08, 0.11)	0.29	0.76	
Proportion Lichen: Nearest neighbour distance (end)	0.55 (0.46, 0.65)	11.9	 < 0.0001	
Step length: Simple ratio index	− 0.03 (− 0.09, 0.03)	− 1.13	0.25	
Proportion Forest: Simple ratio index	− 0.21 (− 0.78, 0.36)	− 0.71	0.47	
Proportion Lichen: Simple ratio index	− 4.03 (− 4.64, − 3.43)	− 13.1	 < 0.0001	

Fig. 3 Habitat selection coefficients and 95% confidence intervals for our integrated step selection model. Coefficients that do not overlap zero are deemed significant (see Table S3). Colours denote which habitat each fixed effect term was associated with and in cases where a term was not associated with habitat, points were coloured grey. Step length, nearest neighbour distances (both start and end), and the simple ratio index were log-transformed. Measures of distance at the end of a step are included when the distance term is a single term in a model or included in an interaction with habitat type, but measures of distance at the start of the step are included when the distance term is in an interaction with step length. These terms model the effect of nearest neighbour distance at the end of a step on the likelihood of selecting a given habitat type because in the case of habitat selection, proximity to conspecifics is more relevant for where an animal is going. By contrast, the effect of nearest neighbour distance at the start of a step on the likelihood of taking a shorter or longer step because in the case of movement, proximity to conspecifics is more relevant for where an animal is starting. Note, the lower panel is a zoom of the upper panel for enhanced visualization of coefficients with relatively small confidence intervals

Fig. 4 Relationship between expected speed (m/hr) as a function of changes in the proportion of lichen and forest habitats within 200 m of a given point location for individual caribou

Fig. 5 Relative selection strength of forest and lichen habitats as a function of nearest neighbour distance (m) (top panels) and shared dyadic simple ratio index between nearest neighbours (bottom panels). The dotted horizontal line represents no response, while values above the line indicate the population is selecting to be closer to that habitat than expected or to have higher shared dyadic SRI than expected, and below the dotted line the population is selecting to be farther from that habitat than expected or to have a lower shared dyadic SRI than expected. Interpretation for RSS values are that individuals generally tend to select to be near to conspecifics when selecting lichen habitat relative to its availability, whereas the response to nearest neighbours in forest habitat relative to its availability is limited. Meanwhile, individuals select for familiar nearest neighbours in all habitats

Discussion

Our study examined apparent social preference in the context of shared space use using socially informed integrated step selection functions. We integrate contemporary social network analysis within a traditional movement ecology and habitat selection framework. Although social networks were well connected at the population level (Fig.’s S1, S2, and S3), we highlight two forms of social preference, including long-term temporal stability of associations among individuals, and an effect of short-term social familiarity on habitat selection. Further, we found that individuals tended to select foraging habitat near familiar individuals, which could be one way that individuals reduce competition for forage at cratering sites. Based on our unification of social network analysis with integrated step selection functions, we highlight the influence of collective movement and preferred associations on habitat selection and foraging.

Our findings reveal that caribou are social in nearly all circumstances, although we observed a social hierarchy of movement-dependent social associations. Specifically, individuals tended to select to be closer to nearest neighbours in lichen habitat relative to forest habitats regardless of the familiarity of nearest neighbours. Within the movement ecology literature for ungulates, there is an assumption that slower movement in a given habitat represents foraging behaviour and faster movement represents searching behaviour [29]. Our results support this assumption; individuals moved slowly in lichen habitat. In conjunction with our observations in the field and past work in our system, we consider lichen habitat as key foraging habitat and dietary item for caribou in winter [33, 34, 52]. Individuals are more likely to trust social information about food sources and predation risk from familiar individuals, but the potential costs are an increase in competition at foraging patches. Despite the potential for competition, individual selection for lichen habitat was strongest in close proximity to conspecifics, suggesting a potential trade-off between the costs (i.e. competition) and benefits (i.e. access to information or anti-predator vigilance) of social foraging. Lichen habitat is typically open, suggesting the possibility that individuals may remain in visual and vocal contact, thereby facilitating and maintaining social cohesion during foraging [66]. Our results are also corroborated by other ungulate systems. In bison (Bison bison), the social environment in combination with recent knowledge of local foraging options dictated whether individuals followed, or left, a group [67]. Moreover, in the bison system, the costs and benefits of foraging in a group are moderated by collective decision making [68] and collective movement [69], both of which are likely involved in the foraging decisions made by caribou. Here, we elucidate potential behavioural mechanisms (i.e., foraging) that influence the frequency and magnitude of social associations.

The emergent geometry of collective movement and spatial arrangement of individuals in a group appears to change as individuals adjust their behaviour based on the availability of resources and the presence of familiar conspecifics [70]. Assamese macaques (Macaca assamensis) distance from one another during foraging, but move collectively between foraging sites [71], while individual giraffes (Giraffa camelopardalis) show social preference for conspecifics during foraging, but not during movement [24]. Interestingly, macaques foraged in closer proximity to individuals of similar dominance rank, but for giraffes it was unclear whether observed social preference was the result of passive or active assortment. For caribou, dominance hierarchies are linear and typically driven by body size [72], suggesting that social preference in caribou could also be related to dominance. Our ability to delineate aspects of the social environment between collective movement and habitat selection is useful for disentangling passive or active assortment, for example dominance rank, conspecific attraction, or the transfer of information about foraging resources.

We assumed that moving with familiar conspecifics is the result of information transfer about the location or quality of cratering sites. Despite clear patterns of social foraging, the question remains how individuals balance potential competition at craters in the winter, which can be substantial [72]. Selection for lichen habitat relative to its availability in groups could reflect the use of social information about the location of foraging sites [32] or predation [73], both of which could counter-act potential for direct competition. Moreover, it is possible that direct competition among familiar individuals is low; an assumption which requires further testing. Craters can vary in size and distribution [33]; however, craters may only be large enough for a single individual to forage at a time [74]. We propose that while caribou generally have larger group sizes in winter [75], groups vary in size based on movement and habitat selection behaviour presumably to balance the trade-off between competition and information acquisition. Furthermore, female caribou often have antlers, which unlike males, persist into winter. Females are hypothesized to use their antlers to defend craters and exert dominance over both males and females without antlers [72, 76]. Our interpretation is corroborated by theory used to explain fission–fusion dynamics, where individuals are expected to split and merge through space and time to reduce conflict and competition during foraging [77, 78].

We wish to highlight an important caveat related to sample size in our study. Despite the small number of individuals each year (n = 21 total individuals; N = 38 caribou-years), the number of GPS relocations per individual was high (n ~ 900 relocations). Street et al. [79] highlight that as few as two individuals (but ideally > 20) with 100 relocations per individual may be used to infer population-level habitat selection using selection function type analyses. Our sample size of individuals and relocations per individual far exceed minimum thresholds, and we are confident our findings are stable for individuals and populations [79].

Within a socially informed integrated step selection framework, we bridge the theoretical and methodological gap between social network analysis, movement ecology, and habitat selection. We demonstrate that although caribou are social in nearly all situations, individuals tend to be the most social when selecting lichen relative to its availability. Our synthesis of integrated step selection functions with social networks to describe the relationships between sociality, movement, and habitat selection is an important step towards identifying the roles of physical space and animal space use as factors influencing the social environment [12]. Moreover, individual variation in phenotypes attributable to movement or habitat selection may affect how individuals experience the social environment [80, 81]. Movement, habitat selection, and social behaviour are clearly linked; as van Moorter et al. [2] described movement as the ‘glue’ connecting habitat selection to the physical location of a given set of habitats, we posit that movement is the glue connecting collective habitat selection to the social environment.

Supplementary Information

Additional file 1.

Acknowledgements

We respectfully acknowledge the territory in which data were collected and analyzed as the ancestral homelands of the Beothuk and the Island of Newfoundland as the ancestral homelands of the Mi’kmaq and Beothuk. We thank M. Laforge, M. Bonar, J. Hendrix, and R. Huang for help in the field, and D. Wilson, K. Lewis, I. Fleming, G. Albery, members of the Wildlife Evolutionary Ecology Lab, and three anonymous reviewers for helpful comments on previous versions of this manuscript. We also thank members of the Newfoundland and Labrador Wildlife Division, including S. Moores, B. Adams, W. Barney, and J. Neville for facilitating animal captures and for logistical support in the field.

Author contributions

QMRW and EVW conceived the ideas and designed method; QMRW conducted fieldwork QMRW and generated social networks; QMRW, CMP, KAK, and JWT conducted movement and spatial analysis; QMRW led the writing of the manuscript; EVW developed the research program. All authors contributed critically to drafts and gave final approval for publication.

Funding

Funding for this study was provided by Natural Sciences and Engineering Research Council (NSERC) Vanier Canada Graduate Scholarships to QMRW and CMP, a NSERC Canada Graduate Scholarship (Masters and Doctoral) to KAK, and a NSERC Discovery Grant to EVW.

Availability of data materials

All code and data used for statistical analysis and figures are archived in Zenodo: 10.5281/zenodo.4549509

Declarations

Ethics approval and consent to participate

All animal captures and handling procedures were consistent with the American Society of Mammologist guidelines and were approved by Memorial University Animal Use Protocol No. 20152067.

Competing interests

The authors declare no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Matthiopoulos J, Fieberg J, Aarts G. Species-habitat associations: spatial data, predictive models, and ecological insights. Species-habitat associations: Spatial data, predictive models, and ecological insights. 2020.
2. Van Moorter B Rolandsen CM Basille M Gaillard J-M Movement is the glue connecting home ranges and habitat selection J Anim Ecol 2016 85 21 31 10.1111/1365-2656.12394 25980987
Van Moorter B, Rolandsen CM, Basille M, Gaillard J-M. Movement is the glue connecting home ranges and habitat selection. J Anim Ecol. 2016;85:21–31.25980987 10.1111/1365-2656.12394
3. Hinde RA Interactions, relationships and social structure Man 1976 11 1 17 10.2307/2800384
Hinde RA. Interactions, relationships and social structure. Man. 1976;11:1–17.10.2307/2800384
4. He P Maldonado-Chaparro AA Farine DR The role of habitat configuration in shaping social structure: a gap in studies of animal social complexity Behav Ecol Sociobiol 2019 73 9 10.1007/s00265-018-2602-7
He P, Maldonado-Chaparro AA, Farine DR. The role of habitat configuration in shaping social structure: a gap in studies of animal social complexity. Behav Ecol Sociobiol. 2019;73:9.10.1007/s00265-018-2602-7
5. Spiegel O Leu ST Bull CM Sih A What’s your move? Movement as a link between personality and spatial dynamics in animal populations Ecol Lett 2017 20 3 18 10.1111/ele.12708 28000433
Spiegel O, Leu ST, Bull CM, Sih A. What’s your move? Movement as a link between personality and spatial dynamics in animal populations. Ecol Lett. 2017;20:3–18.28000433 10.1111/ele.12708
6. Strandburg-Peshkin A Papageorgiou D Crofoot MC Farine DR Inferring influence and leadership in moving animal groups Philos Trans R Soc B 2018 373 20170006 10.1098/rstb.2017.0006
Strandburg-Peshkin A, Papageorgiou D, Crofoot MC, Farine DR. Inferring influence and leadership in moving animal groups. Philos Trans R Soc B. 2018;373:20170006.10.1098/rstb.2017.0006
7. Chamaillé-Jammes S Fritz H Valeix M Murindagomo F Clobert J Resource variability, aggregation and direct density dependence in an open context: The local regulation of an African elephant population J Anim Ecol 2008 77 135 144 10.1111/j.1365-2656.2007.01307.x 17986249
Chamaillé-Jammes S, Fritz H, Valeix M, Murindagomo F, Clobert J. Resource variability, aggregation and direct density dependence in an open context: The local regulation of an African elephant population. J Anim Ecol. 2008;77:135–44.17986249 10.1111/j.1365-2656.2007.01307.x
8. Spiegel O Sih A Leu ST Bull CM Where should we meet? Mapping social network interactions of sleepy lizards shows sex-dependent social network structure Anim Behav 2018 136 207 215 10.1016/j.anbehav.2017.11.001
Spiegel O, Sih A, Leu ST, Bull CM. Where should we meet? Mapping social network interactions of sleepy lizards shows sex-dependent social network structure. Anim Behav. 2018;136:207–15.10.1016/j.anbehav.2017.11.001
9. van Beest FM Uzal A Vander Wal E Laforge MP Contasti AL Colville D Increasing density leads to generalization in both coarse-grained habitat selection and fine-grained resource selection in a large mammal J Anim Ecol 2014 83 147 156 10.1111/1365-2656.12115 23931034
van Beest FM, Uzal A, Vander Wal E, Laforge MP, Contasti AL, Colville D, et al. Increasing density leads to generalization in both coarse-grained habitat selection and fine-grained resource selection in a large mammal. J Anim Ecol. 2014;83:147–56.23931034 10.1111/1365-2656.12115
10. Conradt L Roper TJ Consensus decision making in animals Trends Ecol Evol 2005 20 449 456 10.1016/j.tree.2005.05.008 16701416
Conradt L, Roper TJ. Consensus decision making in animals. Trends Ecol Evol. 2005;20:449–56.16701416 10.1016/j.tree.2005.05.008
11. Jolles JW King AJ Killen SS The role of individual heterogeneity in collective animal behaviour Trends Ecol Evol 2020 35 278 291 10.1016/j.tree.2019.11.001 31879039
Jolles JW, King AJ, Killen SS. The role of individual heterogeneity in collective animal behaviour. Trends Ecol Evol. 2020;35:278–91.31879039 10.1016/j.tree.2019.11.001
12. Strandburg-Peshkin A Farine DR Crofoot MC Couzin ID Habitat structure shapes individual decisions and emergent group structure in collectively moving wild baboons Elife 2017 6 e19505 10.7554/eLife.19505 28139196
Strandburg-Peshkin A, Farine DR, Crofoot MC, Couzin ID. Habitat structure shapes individual decisions and emergent group structure in collectively moving wild baboons. Elife. 2017;6:e19505.28139196 10.7554/eLife.19505
13. Aben J Pellikka P Travis JMJ A call for viewshed ecology: Advancing our understanding of the ecology of information through viewshed analysis Methods Ecol Evol 2018 9 624 633 10.1111/2041-210X.12902
Aben J, Pellikka P, Travis JMJ. A call for viewshed ecology: Advancing our understanding of the ecology of information through viewshed analysis. Methods Ecol Evol. 2018;9:624–33.10.1111/2041-210X.12902
14. Aben J Signer J Heiskanen J Pellikka P Travis JMJ What you see is where you go: visibility influences movement decisions of a forest bird navigating a three-dimensional-structured matrix: What you see is where you go Biol Lett 2021 17 20200478 10.1098/rsbl.2020.0478 33497591
Aben J, Signer J, Heiskanen J, Pellikka P, Travis JMJ. What you see is where you go: visibility influences movement decisions of a forest bird navigating a three-dimensional-structured matrix: What you see is where you go. Biol Lett. 2021;17:20200478.33497591 10.1098/rsbl.2020.0478
15. Harel R Spiegel O Getz WM Nathan R Social foraging and individual consistency in following behaviour: testing the information centre hypothesis in free-ranging vultures Proc R Soc B 2017 284 20162654 10.1098/rspb.2016.2654 28404771
Harel R, Spiegel O, Getz WM, Nathan R. Social foraging and individual consistency in following behaviour: testing the information centre hypothesis in free-ranging vultures. Proc R Soc B. 2017;284:20162654.28404771 10.1098/rspb.2016.2654
16. Lachlan RF Crooks L Laland KN Who follows whom? Shoaling preferences and social learning of foraging information in guppies Anim Behav 1998 56 181 190 10.1006/anbe.1998.0760 9710476
Lachlan RF, Crooks L, Laland KN. Who follows whom? Shoaling preferences and social learning of foraging information in guppies. Anim Behav. 1998;56:181–90.9710476 10.1006/anbe.1998.0760
17. Giraldeau LA Dubois F Social foraging and the study of exploitative behavior Advances in the Study of Behavior 2008 38 59 104 10.1016/S0065-3454(08)00002-8
Giraldeau LA, Dubois F. Social foraging and the study of exploitative behavior. Advances in the Study of Behavior. 2008;38:59–104.10.1016/S0065-3454(08)00002-8
18. Patin R Fortin D Sueur C Chamaillé-Jammes S Space use and leadership modify dilution effects on optimal vigilance under food-safety trade-offs Am Nat 2019 193 E15 28 10.1086/700566 30624106
Patin R, Fortin D, Sueur C, Chamaillé-Jammes S. Space use and leadership modify dilution effects on optimal vigilance under food-safety trade-offs. Am Nat. 2019;193:E15-28.30624106 10.1086/700566
19. Spiegel O Crofoot MC The feedback between where we go and what we know—information shapes movement, but movement also impacts information acquisition Curr Opin Behav Sci 2016 12 90 96 10.1016/j.cobeha.2016.09.009
Spiegel O, Crofoot MC. The feedback between where we go and what we know—information shapes movement, but movement also impacts information acquisition. Curr Opin Behav Sci. 2016;12:90–6.10.1016/j.cobeha.2016.09.009
20. Jones TB Patrick SC Evans JC Wells MR Consistent sociality but flexible social associations across temporal and spatial foraging contexts in a colonial breeder Ecol Lett 2020 23 1085 1096 10.1111/ele.13507 32314533
Jones TB, Patrick SC, Evans JC, Wells MR. Consistent sociality but flexible social associations across temporal and spatial foraging contexts in a colonial breeder. Ecol Lett. 2020;23:1085–96.32314533 10.1111/ele.13507
21. Giraldeau L-A Caraco T Social foraging theory 2018 Princeton Princeton University Press
Giraldeau L-A, Caraco T. Social foraging theory. Princeton: Princeton University Press; 2018.
22. Mourier J Vercelloni J Planes S Evidence of social communities in a spatially structured network of a free-ranging shark species Anim Behav 2012 83 389 401 10.1016/j.anbehav.2011.11.008
Mourier J, Vercelloni J, Planes S. Evidence of social communities in a spatially structured network of a free-ranging shark species. Anim Behav. 2012;83:389–401.10.1016/j.anbehav.2011.11.008
23. Best EC Seddon JM Dwyer RG Goldizen AW Social preference influences female community structure in a population of wild eastern grey kangaroos Anim Behav 2013 86 1031 1040 10.1016/j.anbehav.2013.09.008
Best EC, Seddon JM, Dwyer RG, Goldizen AW. Social preference influences female community structure in a population of wild eastern grey kangaroos. Anim Behav. 2013;86:1031–40.10.1016/j.anbehav.2013.09.008
24. Muller Z Cantor M Cuthill IC Harris S Giraffe social preferences are context dependent Anim Behav 2018 146 37 49 10.1016/j.anbehav.2018.10.006
Muller Z, Cantor M, Cuthill IC, Harris S. Giraffe social preferences are context dependent. Anim Behav. 2018;146:37–49.10.1016/j.anbehav.2018.10.006
25. Silk JB Social components of fitness in primate groups Science 2007 317 1347 1351 10.1126/science.1140734 17823344
Silk JB. Social components of fitness in primate groups. Science. 2007;317:1347–51.17823344 10.1126/science.1140734
26. Spiegel O Leu ST Sih A Bull CM Socially interacting or indifferent neighbours? Randomization of movement paths to tease apart social preference and spatial constraints Methods Ecol Evol 2016 7 971 979 10.1111/2041-210X.12553
Spiegel O, Leu ST, Sih A, Bull CM. Socially interacting or indifferent neighbours? Randomization of movement paths to tease apart social preference and spatial constraints. Methods Ecol Evol. 2016;7:971–9.10.1111/2041-210X.12553
27. Pinter-Wollman N Hobson EA Smith JE Edelman AJ Shizuka D De Silva S The dynamics of animal social networks: Analytical, conceptual, and theoretical advances Behav Ecol 2013 25 242 255 10.1093/beheco/art047
Pinter-Wollman N, Hobson EA, Smith JE, Edelman AJ, Shizuka D, De Silva S, et al. The dynamics of animal social networks: Analytical, conceptual, and theoretical advances. Behav Ecol. 2013;25:242–55.10.1093/beheco/art047
28. Croft DP Darden SK Wey TW Current directions in animal social networks Curr Opin Behav Sci 2016 12 52 58 10.1016/j.cobeha.2016.09.001
Croft DP, Darden SK, Wey TW. Current directions in animal social networks. Curr Opin Behav Sci. 2016;12:52–8.10.1016/j.cobeha.2016.09.001
29. Owen-Smith N Fryxell JM Merrill EH Foraging theory upscaled: the behavioural ecology of herbivore movement Philos Trans R Soc B 2010 365 2267 2278 10.1098/rstb.2010.0095
Owen-Smith N, Fryxell JM, Merrill EH. Foraging theory upscaled: the behavioural ecology of herbivore movement. Philos Trans R Soc B. 2010;365:2267–78.10.1098/rstb.2010.0095
30. Bergerud AT Mercer WE Caribou introductions in eastern North America Wildl Soc Bull 1989 17 111 120
Bergerud AT, Mercer WE. Caribou introductions in eastern North America. Wildl Soc Bull. 1989;17:111–20.
31. Bastille-Rousseau G Schaefer JA Mahoney SP Murray DL Population decline in semi-migratory caribou (Rangifer tarandus): intrinsic or extrinsic drivers? Can J Zool 2013 91 820 828 10.1139/cjz-2013-0154
Bastille-Rousseau G, Schaefer JA, Mahoney SP, Murray DL. Population decline in semi-migratory caribou (Rangifer tarandus): intrinsic or extrinsic drivers? Can J Zool. 2013;91:820–8.10.1139/cjz-2013-0154
32. Lesmerises F Johnson CJ St-Laurent M-H Landscape knowledge is an important driver of the fission dynamics of an alpine ungulate Anim Behav 2018 140 39 47 10.1016/j.anbehav.2018.03.014
Lesmerises F, Johnson CJ, St-Laurent M-H. Landscape knowledge is an important driver of the fission dynamics of an alpine ungulate. Anim Behav. 2018;140:39–47.10.1016/j.anbehav.2018.03.014
33. Bergerud AT Relative abundance of food in winter for Newfoundland caribou Oikos 1974 25 379 387 10.2307/3543960
Bergerud AT. Relative abundance of food in winter for Newfoundland caribou. Oikos. 1974;25:379–87.10.2307/3543960
34. Webber QMR Ferraro KM Hendrix JG Vander WE What do caribou eat? A review of the literature on caribou diet Can J Zool 2022 100 197 206 10.1139/cjz-2021-0162
Webber QMR, Ferraro KM, Hendrix JG, Vander WE. What do caribou eat? A review of the literature on caribou diet. Can J Zool. 2022;100:197–206.10.1139/cjz-2021-0162
35. Peignier M Webber QMR Koen EL Laforge MP Robitaille AL Vander WE Space use and social association in a gregarious ungulate: Testing the conspecific attraction and resource dispersion hypotheses Ecol Evol 2019 9 5133 5145 10.1002/ece3.5071 31110667
Peignier M, Webber QMR, Koen EL, Laforge MP, Robitaille AL, Vander WE. Space use and social association in a gregarious ungulate: Testing the conspecific attraction and resource dispersion hypotheses. Ecol Evol. 2019;9:5133–45.31110667 10.1002/ece3.5071
36. Boertje RD An Energy model for adult female Caribou of the Denali herd Alaska J Range Manag 1985 38 468 473 10.2307/3899725
Boertje RD. An Energy model for adult female Caribou of the Denali herd. Alaska J Range Manag. 1985;38:468–73.10.2307/3899725
37. Duquette LS Klein DR Activity budgets and group size of caribou during spring migration Can J Zool 1987 65 1 164 168 10.1139/z87-023
Duquette LS, Klein DR. Activity budgets and group size of caribou during spring migration. Can J Zool. 1987;65(1):164–8.10.1139/z87-023
38. Fortin D Courtois R Etcheverry P Dussault C Gingras A Winter selection of landscapes by woodland caribou: behavioural response to geographical gradients in habitat attributes J Appl Ecol 2008 45 1392 1400 10.1111/j.1365-2664.2008.01542.x
Fortin D, Courtois R, Etcheverry P, Dussault C, Gingras A. Winter selection of landscapes by woodland caribou: behavioural response to geographical gradients in habitat attributes. J Appl Ecol. 2008;45:1392–400.10.1111/j.1365-2664.2008.01542.x
39. Integrated-Informatics. Sustainable Development and Strategic Science Branch Land Cover Classification. St. John’s, NL; 2014. p. 3–19.
40. Bjørneraas K Van Moorter B Rolandsen CM Herfindal I Screening global positioning system location data for errors using animal movement characteristics J Wildl Manag 2010 74 1361 1366
Bjørneraas K, Van Moorter B, Rolandsen CM, Herfindal I. Screening global positioning system location data for errors using animal movement characteristics. J Wildl Manag. 2010;74:1361–6.
41. Börger L Franconi N De Michele G Gantz A Meschi F Manica A Effects of sampling regime on the mean and variance of home range size estimates J Anim Ecol 2006 75 1393 1405 10.1111/j.1365-2656.2006.01164.x 17032372
Börger L, Franconi N, De Michele G, Gantz A, Meschi F, Manica A, et al. Effects of sampling regime on the mean and variance of home range size estimates. J Anim Ecol. 2006;75:1393–405.17032372 10.1111/j.1365-2656.2006.01164.x
42. Frair JL Fieberg J Hebblewhite M Cagnacci F DeCesare NJ Pedrotti L Resolving issues of imprecise and habitat-biased locations in ecological analyses using GPS telemetry data Philos Trans Royal Soc B: Biol Sci 2010 365 2187 2200 10.1098/rstb.2010.0084
Frair JL, Fieberg J, Hebblewhite M, Cagnacci F, DeCesare NJ, Pedrotti L. Resolving issues of imprecise and habitat-biased locations in ecological analyses using GPS telemetry data. Philos Trans Royal Soc B: Biol Sci. 2010;365:2187–200.10.1098/rstb.2010.0084
43. Avgar T Potts JR Lewis MA Boyce MS Integrated step selection analysis: Bridging the gap between resource selection and animal movement Methods Ecol Evol 2016 7 619 630 10.1111/2041-210X.12528
Avgar T, Potts JR, Lewis MA, Boyce MS. Integrated step selection analysis: Bridging the gap between resource selection and animal movement. Methods Ecol Evol. 2016;7:619–30.10.1111/2041-210X.12528
44. Duchesne T Fortin D Rivest LP Equivalence between step selection functions and biased correlated random walks for statistical inference on animal movement PLoS ONE 2015 10 1 12 10.1371/journal.pone.0122947
Duchesne T, Fortin D, Rivest LP. Equivalence between step selection functions and biased correlated random walks for statistical inference on animal movement. PLoS ONE. 2015;10:1–12.10.1371/journal.pone.0122947
45. Basille M Fortin D Dussault C Bastille-Rousseau G Ouellet J-P Courtois R Plastic response of fearful prey to the spatio-temporal dynamics of predator distribution Ecology 2015 96 150511123704000 10.1890/14-1706.1
Basille M, Fortin D, Dussault C, Bastille-Rousseau G, Ouellet J-P, Courtois R. Plastic response of fearful prey to the spatio-temporal dynamics of predator distribution. Ecology. 2015;96:150511123704000.10.1890/14-1706.1
46. Fortin D Beyer HL Boyce MS Smith DW Duchesne T Mao JS Wolves influence elk movements: behavior shapes a trophic cascade in Yellowstone National Park Ecology 2005 86 1320 1330 10.1890/04-0953
Fortin D, Beyer HL, Boyce MS, Smith DW, Duchesne T, Mao JS. Wolves influence elk movements: behavior shapes a trophic cascade in Yellowstone National Park. Ecology. 2005;86:1320–30.10.1890/04-0953
47. Prokopenko CM Boyce MS Avgar T Characterizing wildlife behavioural responses to roads using integrated step selection analysis J Appl Ecol 2017 54 470 479 10.1111/1365-2664.12768
Prokopenko CM, Boyce MS, Avgar T. Characterizing wildlife behavioural responses to roads using integrated step selection analysis. J Appl Ecol. 2017;54:470–9.10.1111/1365-2664.12768
48. Signer J Fieberg J Avgar T Animal movement tools (amt): R package for managing tracking data and conducting habitat selection analyses Ecol Evol 2019 9 880 890 10.1002/ece3.4823 30766677
Signer J, Fieberg J, Avgar T. Animal movement tools (amt): R package for managing tracking data and conducting habitat selection analyses. Ecol Evol. 2019;9:880–90.30766677 10.1002/ece3.4823
49. Roberts DR Bahn V Ciuti S Boyce MS Elith J Guillera-Arroita G Cross-validation strategies for data with temporal, spatial, hierarchical, or phylogenetic structure Ecography 2017 40 913 929 10.1111/ecog.02881
Roberts DR, Bahn V, Ciuti S, Boyce MS, Elith J, Guillera-Arroita G, et al. Cross-validation strategies for data with temporal, spatial, hierarchical, or phylogenetic structure. Ecography. 2017;40:913–29.10.1111/ecog.02881
50. Fortin D Fortin ME Beyer HL Duchesne T Courant S Dancose K Group-size-mediated habitat selection and group fusion-fission dynamics of bison under predation risk Ecology. 2009 90 9 2480 90 10.1890/08-0345.1 19769126
Fortin D, Fortin ME, Beyer HL, Duchesne T, Courant S, Dancose K. Group-size-mediated habitat selection and group fusion-fission dynamics of bison under predation risk. Ecology. 2009;90(9):2480–90. 10.1890/08-0345.1.19769126 10.1890/08-0345.1
51. Kasozi H Montgomery RA Variability in the estimation of ungulate group sizes complicates ecological inference Ecol Evol 2020 10 14 6881 6889 10.1002/ece3.6463 32760498
Kasozi H, Montgomery RA. Variability in the estimation of ungulate group sizes complicates ecological inference. Ecol Evol. 2020;10(14):6881–9.32760498 10.1002/ece3.6463
52. Robitaille AL Webber QMR Turner JW Vander WE The problem and promise of scale in multilayer animal social networks Curr Zool 2021 67 113 123 10.1093/cz/zoaa052 33654495
Robitaille AL, Webber QMR, Turner JW, Vander WE. The problem and promise of scale in multilayer animal social networks. Curr Zool. 2021;67:113–23.33654495 10.1093/cz/zoaa052
53. Franks DW Ruxton GD James R Sampling animal association networks with the gambit of the group Behav Ecol Sociobiol 2010 64 493 503 10.1007/s00265-009-0865-8
Franks DW, Ruxton GD, James R. Sampling animal association networks with the gambit of the group. Behav Ecol Sociobiol. 2010;64:493–503.10.1007/s00265-009-0865-8
54. Cairns SJ Schwager SJ A comparison of association indices Anim Behav 1987 35 1454 1469 10.1016/S0003-3472(87)80018-0
Cairns SJ, Schwager SJ. A comparison of association indices. Anim Behav. 1987;35:1454–69.10.1016/S0003-3472(87)80018-0
55. Farine DR Whitehead H Constructing, conducting and interpreting animal social network analysis J Anim Ecol 2015 84 1144 1163 10.1111/1365-2656.12418 26172345
Farine DR, Whitehead H. Constructing, conducting and interpreting animal social network analysis. J Anim Ecol. 2015;84:1144–63.26172345 10.1111/1365-2656.12418
56. R Core Team. R: A language and environment for statistical computing. 2019. p. https://www.r-project.org/.
57. Robitaille AL Webber QMR Vander WE Conducting social network analysis with animal telemetry data: applications and methods using spatsoc Methods Ecol Evol 2019 10 1203 1211 10.1111/2041-210X.13215
Robitaille AL, Webber QMR, Vander WE. Conducting social network analysis with animal telemetry data: applications and methods using spatsoc. Methods Ecol Evol. 2019;10:1203–11.10.1111/2041-210X.13215
58. Csárdi G Nepusz T The igraph software package for complex network research InterJournal Complex Syst 2006 1695 1 9
Csárdi G, Nepusz T. The igraph software package for complex network research. InterJournal Complex Syst. 2006;1695:1–9.
59. Northrup JM Vander Wal E Bonar M Fieberg J Laforge MP Leclerc M Conceptual and methodological advances in habitat-selection modeling: guidelines for ecology and evolution Ecol Appl 2022 32 1 31 10.1002/eap.2470
Northrup JM, Vander Wal E, Bonar M, Fieberg J, Laforge MP, Leclerc M, et al. Conceptual and methodological advances in habitat-selection modeling: guidelines for ecology and evolution. Ecol Appl. 2022;32:1–31.10.1002/eap.2470
60. Fieberg J Rieger RH Zicus MC Schildcrout JS Regression modelling of correlated data in ecology: Subject-specific and population averaged response patterns J Appl Ecol 2009 46 1018 1025 10.1111/j.1365-2664.2009.01692.x
Fieberg J, Rieger RH, Zicus MC, Schildcrout JS. Regression modelling of correlated data in ecology: Subject-specific and population averaged response patterns. J Appl Ecol. 2009;46:1018–25.10.1111/j.1365-2664.2009.01692.x
61. Hebblewhite M Merrill E Modelling wildlife-human relationships for social species with mixed-effects resource selection models J Appl Ecol 2008 45 834 844 10.1111/j.1365-2664.2008.01466.x
Hebblewhite M, Merrill E. Modelling wildlife-human relationships for social species with mixed-effects resource selection models. J Appl Ecol. 2008;45:834–44.10.1111/j.1365-2664.2008.01466.x
62. Muff S Signer J Fieberg J Accounting for individual-specific variation in habitat-selection studies: efficient estimation of mixed-effects models using Bayesian or frequentist computation J Anim Ecol 2020 89 80 92 10.1111/1365-2656.13087 31454066
Muff S, Signer J, Fieberg J. Accounting for individual-specific variation in habitat-selection studies: efficient estimation of mixed-effects models using Bayesian or frequentist computation. J Anim Ecol. 2020;89:80–92.31454066 10.1111/1365-2656.13087
63. Burnham KP Anderson DR Model selection and multimodel inference [Internet] 2002 New York Springer
Burnham KP, Anderson DR. Model selection and multimodel inference [Internet]. New York: Springer; 2002.
64. Brooks ME Kristensen K van Benthem KJ Magnusson A Berg CW Nielsen A glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling R J 2017 9 378 400 10.32614/RJ-2017-066
Brooks ME, Kristensen K, van Benthem KJ, Magnusson A, Berg CW, Nielsen A, et al. glmmTMB balances speed and flexibility among packages for zero-inflated generalized linear mixed modeling. R J. 2017;9:378–400.10.32614/RJ-2017-066
65. Avgar T Lele SR Keim JL Boyce MS Relative selection strength: quantifying effect size in selection inference Ecol Evol 2017 7 5322 5330 10.1002/ece3.3122 28770070
Avgar T, Lele SR, Keim JL, Boyce MS. Relative selection strength: quantifying effect size in selection inference. Ecol Evol. 2017;7:5322–30.28770070 10.1002/ece3.3122
66. Jacobs A Group cohesiveness during collective movements: Travelling apart together Behav Proc 2010 84 678 680 10.1016/j.beproc.2010.03.004
Jacobs A. Group cohesiveness during collective movements: Travelling apart together. Behav Proc. 2010;84:678–80.10.1016/j.beproc.2010.03.004
67. Merkle JA Sigaud M Fortin D To follow or not? How animals in fusion-fission societies handle conflicting information during group decision-making Ecol Lett 2015 18 799 806 10.1111/ele.12457 26013202
Merkle JA, Sigaud M, Fortin D. To follow or not? How animals in fusion-fission societies handle conflicting information during group decision-making. Ecol Lett. 2015;18:799–806.26013202 10.1111/ele.12457
68. Sigaud M Merkle JA Cherry SG Fryxell JM Berdahl A Fortin D Collective decision-making promotes fitness loss in a fusion-fission society Ecol Lett 2017 20 33 40 10.1111/ele.12698 27873440
Sigaud M, Merkle JA, Cherry SG, Fryxell JM, Berdahl A, Fortin D. Collective decision-making promotes fitness loss in a fusion-fission society. Ecol Lett. 2017;20:33–40.27873440 10.1111/ele.12698
69. Courant S Fortin D Time allocation of bison in meadow patches driven by potential energy gains and group size dynamics Oikos 2012 121 1163 1173 10.1111/j.1600-0706.2011.19994.x
Courant S, Fortin D. Time allocation of bison in meadow patches driven by potential energy gains and group size dynamics. Oikos. 2012;121:1163–73.10.1111/j.1600-0706.2011.19994.x
70. Morrell LJ Ruxton GD James R Spatial positioning in the selfish herd Behav Ecol 2011 22 16 22 10.1093/beheco/arq157
Morrell LJ, Ruxton GD, James R. Spatial positioning in the selfish herd. Behav Ecol. 2011;22:16–22.10.1093/beheco/arq157
71. Heesen M Macdonald S Ostner J Schülke O Ecological and Social determinants of group cohesiveness and within-group spatial position in wild assamese macaques Ethology 2015 121 270 283 10.1111/eth.12336
Heesen M, Macdonald S, Ostner J, Schülke O. Ecological and Social determinants of group cohesiveness and within-group spatial position in wild assamese macaques. Ethology. 2015;121:270–83.10.1111/eth.12336
72. Barrette C Vandal D Social rank, dominance, antler size, and access to food in snow-bound wild woodland caribou Behaviour 1986 97 118 145 10.1163/156853986X00342
Barrette C, Vandal D. Social rank, dominance, antler size, and access to food in snow-bound wild woodland caribou. Behaviour. 1986;97:118–45.10.1163/156853986X00342
73. Hamilton WD Geometry for the selfish herd J Theor Biol 1971 31 295 311 10.1016/0022-5193(71)90189-5 5104951
Hamilton WD. Geometry for the selfish herd. J Theor Biol. 1971;31:295–311.5104951 10.1016/0022-5193(71)90189-5
74. Mayor SJ Schaefer JA Schneider DC Mahoney SP The spatial structure of habitat selection: a caribou’s-eye-view Acta Oecologica 2009 35 253 260 10.1016/j.actao.2008.11.004
Mayor SJ, Schaefer JA, Schneider DC, Mahoney SP. The spatial structure of habitat selection: a caribou’s-eye-view. Acta Oecologica. 2009;35:253–60.10.1016/j.actao.2008.11.004
75. Webber QMR Vander WE Context-dependent group size: effects of population density, habitat, and season Behav Ecol 2021 32 5 970 981 10.1093/beheco/arab070
Webber QMR, Vander WE. Context-dependent group size: effects of population density, habitat, and season. Behav Ecol. 2021;32(5):970–81.10.1093/beheco/arab070
76. Schaefer JA Mahoney SP Antlers on female caribou: Biogeography of the bones of contention Ecology 2001 82 3556 3560 10.1890/0012-9658(2001)082[3556:AOFCBO]2.0.CO;2
Schaefer JA, Mahoney SP. Antlers on female caribou: Biogeography of the bones of contention. Ecology. 2001;82:3556–60.10.1890/0012-9658(2001)082[3556:AOFCBO]2.0.CO;2
77. Aureli F Schaffner CM Boesch C Bearder SK Call J Chapman CA Fission-fusion dynamics: new research frameworks Curr Anthropol 2008 49 627 654 10.1086/586708
Aureli F, Schaffner CM, Boesch C, Bearder SK, Call J, Chapman CA, et al. Fission-fusion dynamics: new research frameworks. Curr Anthropol. 2008;49:627–54.10.1086/586708
78. Sueur C King AJ Conradt L Kerth G Lusseau D Mettke-Hofmann C Collective decision-making and fission-fusion dynamics: a conceptual framework Oikos 2011 120 1608 1617 10.1111/j.1600-0706.2011.19685.x
Sueur C, King AJ, Conradt L, Kerth G, Lusseau D, Mettke-Hofmann C, et al. Collective decision-making and fission-fusion dynamics: a conceptual framework. Oikos. 2011;120:1608–17.10.1111/j.1600-0706.2011.19685.x
79. Street GM Potts JR Börger L Beasley JC Demarais S Fryxell JM Solving the sample size problem for resource selection functions Methods Ecol Evol 2021 2021 1 11
Street GM, Potts JR, Börger L, Beasley JC, Demarais S, Fryxell JM, et al. Solving the sample size problem for resource selection functions. Methods Ecol Evol. 2021;2021:1–11.
80. Webber QMR Albery GF Farine DR Pinter-Wollman N Sharma N Spiegel O Behavioural ecology at the spatial–social interface Biol Rev 2023 98 868 886 10.1111/brv.12934 36691262
Webber QMR, Albery GF, Farine DR, Pinter-Wollman N, Sharma N, Spiegel O, et al. Behavioural ecology at the spatial–social interface. Biol Rev. 2023;98:868–86.36691262 10.1111/brv.12934
81. Webber QMR Vander WE An evolutionary framework outlining the integration of individual social and spatial ecology J Anim Ecol 2018 87 113 127 10.1111/1365-2656.12773 29055050
Webber QMR, Vander WE. An evolutionary framework outlining the integration of individual social and spatial ecology. J Anim Ecol. 2018;87:113–27.29055050 10.1111/1365-2656.12773
