Skip to main content

Combined effects of landscape fragmentation and sampling frequency of movement data on the assessment of landscape connectivity

Abstract

Background

Network theory is largely applied in real-world systems to assess landscape connectivity using empirical or theoretical networks. Empirical networks are usually built from discontinuous individual movement trajectories without knowing the effect of relocation frequency on the assessment of landscape connectivity while theoretical networks generally rely on simple movement rules. We investigated the combined effects of relocation sampling frequency and landscape fragmentation on the assessment of landscape connectivity using simulated trajectories and empirical high-resolution (1 Hz) trajectories of Alpine ibex (Capra ibex). We also quantified the capacity of commonly used theoretical networks to accurately predict landscape connectivity from multiple movement processes.

Methods

We simulated forager trajectories from continuous correlated biased random walks in simulated landscapes with three levels of landscape fragmentation. High-resolution ibex trajectories were reconstructed using GPS-enabled multi-sensor biologging data and the dead-reckoning technique. For both simulated and empirical trajectories, we generated spatial networks from regularly resampled trajectories and assessed changes in their topology and information loss depending on the resampling frequency and landscape fragmentation. We finally built commonly used theoretical networks in the same landscapes and compared their predictions to actual connectivity.

Results

We demonstrated that an accurate assessment of landscape connectivity can be severely hampered (e.g., up to 66% of undetected visited patches and 29% of spurious links) when the relocation frequency is too coarse compared to the temporal dynamics of animal movement. However, the level of landscape fragmentation and underlying movement processes can both mitigate the effect of relocation sampling frequency. We also showed that network topologies emerging from different movement behaviours and a wide range of landscape fragmentation were complex, and that commonly used theoretical networks accurately predicted only 30–50% of landscape connectivity in such environments.

Conclusions

Very high-resolution trajectories were generally necessary to accurately identify complex network topologies and avoid the generation of spurious information on landscape connectivity. New technologies providing such high-resolution datasets over long periods should thus grow in the movement ecology sphere. In addition, commonly used theoretical models should be applied with caution to the study of landscape connectivity in real-world systems as they did not perform well as predictive tools.

Background

Studying landscape connectivity is crucial for better understanding many ecological processes including animal foraging behaviour [1], population distribution [2, 3], gene or pathogen flow [4], migration behaviour [5], species interactions [6, 7], or metapopulation persistence [8]. Analysis of landscape connectivity can also be useful to identify critical areas favouring movements [9,10,11] or to evaluate and compare the effect of land planning actions (e.g., site restoration) [12, 13]. Consequently, an accurate assessment and prediction of landscape connectivity is needed to avoid spurious understanding of multiple theoretical and applied landscape connectivity-related processes.

Network (or graph) theory has been largely used in multiple taxa to model and predict connectivity in fragmented landscapes [14,15,16,17,18]. In such landscapes, habitat patches are represented using nodes while observed or predicted movements among patches are represented using links [19]. Links can either mirror the Euclidean distance between patches (i.e., straight lines) [20] or account for matrix composition among patches (i.e., least cost path) [6]. The popularity of spatial networks relies on their easy implementation and the large amount of connectivity indices that can be derived from a network providing quantitative information on landscape connectivity at different levels (i.e., node, group of nodes or network level) [21,22,23]. Their adequatedness in real landscapes for correctly inferring ecological processes is, however, less immediately clear.

A landscape can be described using different network topologies including simple planar graphs (e.g., minimum spanning tree or minimum planar graph, MPG) [24] and more complex ones (e.g., scale-free or small-world networks) [25]. This result in different patterns of connectivity within the graph (Fig. 1A) [26]. For example, a MPG assumes that individuals move in a stepping stone fashion among resource patches, such that links in a MPG never cross (i.e., no shortcut). By contrast, in a scale-free network, few nodes (thereafter called hubs) are highly connected to the rest of the patches (thereafter called peripheral nodes) that have few connections [27]. Consequently, a scale-free network is highly sensitive to the removal of hubs but resistant to random removal of nodes, whereas a MPG shows similar effects of a disturbance independently of its spatial location [3, 28]. Network topology also has implications for processes occurring in spatial networks. Indeed, a disease would not spread much if it appears in a peripheral node while it would spread faster and to a larger extent if appearing in a hub or in a network having shortcuts among patches (such as in small-world networks) [20]. An accurate representation of actual landscape connectivity is thus of paramount importance to correctly identify network topology and, consequently, to understand and predict many landscape connectivity-related processes.

Fig. 1
figure 1

General description of network topologies and design of the study. Ibex photo ©: Franck Merlier

Application of spatial network theory to real-world landscapes generally relies on the modelling of either theoretical or empirical networks. Theoretical networks commonly include a MPG built from the distribution of resource patches [6, 15, 29], or a distance-based network built from both distribution of resource patches and empirically-estimated or expert-based dispersal distances [4, 9, 30]. Empirical graphs are generated by overlapping trajectories of remote-tracked individuals with resource patches to identify network links [17, 18, 31]. With the development of Global Positioning System (GPS) technology, many free-ranging species are now tracked over long periods of time [32, 33]. However, the weight constraints of collars fitted to animals limit battery capacities, generally resulting in a trade-off between relocation frequency and tracking duration depending on scientific questions [34]. For example, empirical networks have been built from GPS relocations collected every hour in plains bison (Bison bison bison) [3] or in African elephants (Loxodonta africana) [17], every two hours in woodland caribou (Rangifer tarandus caribou) [35], every four hours in brown bear (Ursus arctos) [36] or even every 24 h in several species of sea turtle [18]. However, the potential effect of relocation frequency on observed patterns of connectivity remains unclear. More generally, if the temporal dynamics of animal movement is too rapid compared to relocation frequency, detected patterns of landscape connectivity using networks might not reflect actual species movement behaviour. Indeed, discontinuous trajectories like the one generated from GPS monitoring could result in inaccurate estimations of inter-patch movements and, for example, the detection of artificial shortcuts due to missed stepping stone-like movements. In other words, an individual could stay 30 min in a feeding patch, then move to another feeding patch by making a quick stop along the way in a patch to drink but that series of inter-patch movements would not be detected if the individual relocation frequency is every hour or two hours. Consequently, relocation frequency could impact the assessment of connectivity in the sense that networks are perceived as more complex, while being actually simpler (e.g., artificial shortcuts could generate a small-world topology). On one hand, if networks are actually simpler than detected, it assures the relevant use of simple theoretical networks to predict landscape connectivity (i.e., such as the MPG). On the other hand, if networks are actually complex, there is an urgent need to develop theoretical network models that fit with complex functional patterns of connectivity.

In patchy landscapes, the level of habitat fragmentation could also impact observed patterns of landscape connectivity from remote-tracked individuals. Indeed, habitats show different levels of fragmentation starting from simple patch dissection to patch attrition with intermediate patterns such as patch shrinkage [37], resulting in habitat patches having different size, organisation and composition [38]. While patches can provide different resources to animals (e.g., food, shelter), their use (e.g., residency time, number of visits) can vary according to patch characteristics [3, 6, 35, 39]. For example, plains bison and woodland caribou were shown to stay longer in larger and highly profitable patches [35, 40]. Such variations in residency time within habitat patches could result in missed inter-patch movements (e.g., stepping-stones), notably when the sampling schedule is too coarse compared to the temporal dynamic of animal movement [41]. While the interplay between landscape fragmentation and relocation frequency could have noteworthy effects on the assessment of landscape connectivity using networks, it has yet to be appraised.

Finally, individuals could also adjust their movements as a response to space use of congeners, competitors or predators [42,43,44,45], such that landscape connectivity could result in different patterns according to the movement process at stake. For example, individuals can have limited access to the landscape when conspecifics are territorial and consequently have to restrain their movement to their own territory [46]. Similarly, space use and movements can be modified by preys in response to a threat (e.g., predator or human presence) [45, 47]. Consequently, we could expect network topology to be impacted by movement processes and, if so, these movement processes should be accounted for when building theoretical network models to predict connectivity in heterogeneous landscapes.

In this study, we address these outstanding methodological questions by investigating the combined effects of relocation sampling frequency and landscape fragmentation on the assessment of landscape connectivity from networks using (i) simulated trajectories accounting for four movement processes (foraging, foraging + avoidance of an elusive predator, foraging + avoidance of a stalking predator, foraging and territoriality) in simulated habitats with variable levels of fragmentation and (ii) empirical high-resolution (1 Hz) trajectories of Alpine ibex (Capra ibex) tracked in the northern French Alps (Belledonne massif) (Fig. 1C- panels 1-3). For both simulated and empirical trajectories, we also compared predictions from commonly used theoretical networks (i.e., MPG and distance-based graphs) to actual connectivity obtained from high-resolution data (Fig. 1C- panel 4).

Methods

Simulation study

Modelling of movement processes

We used a recently developed modelling framework to simulate individual trajectories in fragmented landscapes from continuous correlated biased random walks [48,49,50]. This framework simulates an animal grazing across stationary resources that deplete and regenerate, based on three processes: consumption and regeneration of resources, resource memory of the grazer, and state-specific biased correlated movement process [48]. The framework can also accommodate predator avoidance during movements by simulating predator appearance in the landscape and a resulting fleeing-behaviour by the grazer that memorizes an encounter location and avoids it for a parameterized period of time [49]. The model can finally include multiple territorial foragers that regularly scent-mark their territory and avoid scent-marks of conspecifics [50]. Detailed descriptions of foraging, predator avoidance and territoriality models are available in [48, 49] and [50], respectively, such that we only provide a brief description of each modelled process in Appendix S1.

Modelling of landscape fragmentation

We generated three types of landscapes of 50 × 50 cells with different levels of fragmentation using an exponential variogram model. This model was characterized by four parameters including the nugget (i.e., variogram’s intercept), the sill (i.e., variogram’s asymptote), the range (r, i.e., distance beyond which variables are no longer autocorrelated) and the trend (t, i.e., average predicted value over the landscape). Constant parameters of the model were: nugget = 0 and sill = 1. We set range and trend coefficients to obtain three levels of fragmentation: low (r = 5, t = 0.3), intermediate (r = 2, t = − 0.5) and high (r = 2, t = − 1.5) [51, 52]. We produced 100 landscapes for each set of parameters to account for stochasticity. Following [48], negative values were truncated to 0 and landscapes were normalized to sum to one. Thus, cells with positive values reflected resources of gradual quality while null values corresponded to non-resource cells. We used the R package gstat (v2.1-1) to produce the landscapes [54]. We then calculated four indices of landscape fragmentation for each simulated landscape to make sure fragmentation characteristics changed between the different landscapes. To do so, we first truncated cells with positive values to one to delimit resource patches. We then calculated the aggregation index (i.e., the ratio between the amounts of realized shared edges and maximum possible shared edges, range: [0–100]), the patch cohesion index (i.e., equivalent of ratio of patch perimeters over patch areas, range: [0–100]), the division index (i.e., the probability that two randomly selected cells are not located in the same patch, range: [0–1]) and the proportion of cells occupied by resource patches (range: [0–100]) [37]. We selected those indices because they range between 0 and 1 such that they are comparable between different landscapes. We used the R package landscapemetrics to calculate the fragmentation indices [55].

Model runs and construction of spatial networks

We ran each model (F: Foraging-only, F + Pe: Foraging and avoidance of an elusive predator, F + Ps: Foraging and avoidance of a stalking predator, F + T: Foraging and territoriality) in each simulated landscape resulting in 1200 simulations (100 landscapes of each type × 4 movement processes × 3 levels of landscape fragmentation) (Fig. 1C- panel 1). Continuous time models were implemented in Java, with time discretized with small regular intervals Δt approximating dt (available code at https://zenodo.org/records/6104214). For each simulation, we ran the model during 10 000 timesteps. As for the F + T model, territories have to emerge first, we ran 20 000 time steps and discarded the first 10 000 to only keep movement behaviour when territories were in place [50]. Model parameters for each simulated process were initialized based on [48,49,50], and are reported in Table S1.

All of the subsequent analyses were performed using the R software (v4.3.0, [53]). We then built a network for each complete simulated trajectory (i.e., 1200 networks). To do so, we identified all movements that occurred between two resource patches and used unique inter-patch movements to build the network (i.e., unweighted and undirectional network) (Fig. 1C- panels 2-3). Then, we regularly resampled each trajectory at various frequencies to mimic discontinuous trajectories like the one generated from GPS monitoring. Notably, we drew one location out of two, one out of three, etc.…, until one location out of fifty, resulting in 4999, 3333, etc.…, down to 199 movement steps, respectively, for each trajectory (Fig. 1C- panel 2). We then built the spatial networks using these resampled trajectories (Fig. 1C- panel 3).

Fig. 2
figure 2

Degree distributions of simulated networks in three types of landscape fragmentation from four movement processes. F: Foraging, F + Pe: Foraging and avoidance of an elusive predator, F + Ps: Foraging and avoidance of a stalking predator, F + T: Foraging and territoriality. All distributions decayed as a power function with scale parameter being higher than 1 (Table S3)

Evaluation of simulated spatial network

Effects of landscape fragmentation and sampling frequency on network topology

Multiple network metrics are available to describe landscape connectivity and to distinguish between network types [22]. Here, we focused on the degree distribution (i.e., the distribution of the number of links that each node has), the global clustering coefficient C (i.e., the number of connected triplets of nodes over the total number of triplets in the graph) and the characteristic path length L (i.e., the number of links in the shortest path between two nodes, averaged over all pairs of nodes) (Fig. 1B). These metrics allow evaluating whether a network topology is complex and similar to classic complex networks such as scale-free or small-world networks (Fig. 1A) [25]. In a scale-free network, the degree distribution follows a power-law distribution with scale parameter \(\alpha\) > 1 (i.e., the probability of a node to be connected to \(k\) other nodes, \(P\left(X=k\right)\propto {k}^{-\alpha }, \alpha >1, k > \text{ 0}\)), reflecting the presence of hubs and peripheral nodes within the network [27, 56]. In a small-world network, many triplets of nodes are connected resulting in a high clustering coefficient and some links provide long range connection between distant nodes (e.g., shortcuts) leading to a short characteristic path length (i.e., fast traversability of the graph) [56, 57]. The ad hoc method to evaluate whether a network is small-world is to compare C and L to the clustering coefficient and characteristic path length of an associated random graph built from the same number of links and nodes than the one of interest [58, 59]. Notably, in a small-world network, the clustering coefficient is much larger than the clustering coefficient of its associated random graph and the characteristic path length is similar to the one of its associated random graph.

To evaluate the effects of landscape fragmentation, sampling frequency and movement process driving forager behaviour on the network topology, we first extracted the node degree of all visited patches of each network. We then fitted a linear mixed-effects model to the log-transformed node degree as a function of landscape fragmentation, the number of relocations used to build the network (correlated to sampling frequency, Fig. 1C-2) and the interaction of both (Fig. 1C-3). We included a random effect of landscape ID to account for autocorrelation in the degrees of patches coming from the same landscape. We ran a model for each simulated movement process (i.e., F, F + Pe, F + Ps and F + T). We also computed network degree distribution for each level of landscape fragmentation and each sampling frequency and fitted the log-transformed frequency of node degrees to the log-transformed degrees using a linear model to estimate whether networks were scale-free (i.e., in such case estimated scale parameter should be higher than 1). We then extracted C and L of each network and built 100 random networks for each simulated network based on the same number of links and visited patches and compared their clustering coefficient (CR) and characteristic path length (LR) to C and L to assess whether the empirical network had the small-world property. More precisely, we followed [60] and calculated the ratio σ = (C/CR) / (L/LR) for each random graph meaning that we had one hundred calculation of σ for each simulated network. As these ratios should be higher than 1 if the simulated network is small-world [60], we defined a simulated network to be small-world if min(σ) > 1. Finally, we calculated the proportion of simulated networks being small-world for each simulated movement process as a function of landscape fragmentation and the number of relocations used to build the network. We used the R packages igraph (v1.4.3, [61]) to build the networks and calculate their connectivity indices, and lme4 (v1.1–33, [62]) and lmerTest (v3.1–3, [63]) to fit the statistical models.

Adequacy with commonly used theoretical networks

We first built the minimum planar graph in each simulated landscape using the R package grainscape [53]. We also simulated two distance-based networks and a random graph in each landscape for each simulated movement process (Fig. 1C- panel 4). Notably, we built distance-based networks using two thresholds for link length: one corresponding to the median length (Lmed) of realized links obtained from unique inter-patch movements of the complete trajectory (DIST50), and one corresponding to the 95% quantile (L95%) of the distribution of realized link lengths (DIST95). Consequently, only patches that were closer than Lmed or L95% were thus connected in DIST50 and DIST95, respectively. For the random graph, we randomly assigned the same number of realized links among patches of the landscape. We then used the Cohen’s Kappa (κ) to assess the performance of theoretical spatial networks (MPG, DIST50, DIST95 and Random) to correctly predict connectivity emerging from different movement models in landscapes having various levels of fragmentation [64, 65]. We specifically calculated κ by comparing predicted links (0/1) of theoretical networks to realized links (0/1) from simulated movement processes in each landscape.

Empirical study—movements of Alpine ibex in the Belledonne massif

Study area and Alpine ibex population monitoring

The empirical study took place in the Belledonne massif located in the northern French Alps (45° 13′ N, 6° 4′ E). The tagged Alpine ibex ranged over a 35km2 area at an altitude around 2100 m within the core of the massif. Ibex are notably adapted to movements in steep terrains, which provide them refuge from perceived predation risk, in particular the steepest rocky cliffs. However, cliffs generally provide few foraging resources, which are rather concentrated in Alpine grasslands [66, 67], such that ibex move frequently among various habitat types to get both food and cover. In total, 10 male ibex (between 7 and 12 years old) were monitored in 2017 using Lotek (3300S and Litetrack models) or Vectronic (Vertex Plus model) GPS collars. Relocation frequency varied between 1 location per hour to 1 location every 2 h and ibex were tracked most of the time nearly one year. Among these 10 ibex, 6 were additionally equipped with Daily Diary multi-sensor biologgers (Wildbyte Technologies 2020), which include tri-axial accelerometer and tri-axial magnetometer sensors recording at 25 Hz/8 Hz, respectively. Data from these biologgers were combined with GPS data collected every one or two hours to reconstruct the movement trajectory every second using the dead-reckoning technique (DR; hereafter GPS-enabled dead-reckoning) [68]. Given a starting point, DR uses the accelerometry and magnetometry data to calculate the speed and heading of an individual animal every second and thereby reconstruct the full high-frequency movement path [69]. To successfully use the DR method, in addition to correctly calibrating the sensors and correct for hard and soft iron bias (see [69]), it is important to recalibrate the reconstructed DR track as errors can accumulate and increase if not accounted for—this is typically done using GPS locations collected at a lower frequency, as detailed in [70]. We focused our analysis on the period for which we had very high-resolution trajectories for all 6 ibex (i.e., June) resulting in a total of 3412 GPS relocations (341 locations on average per individual, range: 204–395) from the 10 GPS-collared ibex and 18 070 914 DR relocations (3 011 819 locations on average per individual, range: 3 010 683–3 012 356) from the 6 ibex that were also equipped with biologgers. Variation in the number of relocations between ibex from GPS-enabled dead-reckoning came from a few missing relocation estimates but since this phenomenon was marginal (i.e., less than 0.5% of the dataset) and the time difference between two successive relocations never exceeded 4 s, we considered missing location to be the same as immediately preceding location.

We used slope and aspect layers provided by the R package terra to describe terrain characteristics in the ibex range at a resolution of 13 × 13 m. We also used description of habitats from the vegetation map provided by the Conservatoire Botanique National Alpin [71]. Notably, 19 habitats were identified within the ibex range from the vegetation map, which we reclassified into 8 classes including: closed forest, grassland, herbaceous—mineral, low ligneous, mineral, open forest, snow and other.

Identification of resource patches for Alpine ibex

To identify ibex resource patches, we first quantified their habitat selection using a resource selection function (RSF) and geolocations of the 10 GPS-collared ibex. In an RSF, characteristics of observed GPS points of an individual are compared to random locations [72]. For each observed location of Alpine ibex, we thus generated 5 random locations within the 99% contour of the utilization distribution (UD) of each individual. UD was estimated using a brownian bridge kernel calculated from the kernelbb function of the R package adehabitatHR with sig2 = 25 (i.e., GPS measurement error) and sig1 being estimated from the liker function of adehabitatHR [73,74,75]. Observed and random locations were then characterized by land cover types using the reclassified vegetation map (8 classes), the log-transformed distance to steep slope (i.e., slope > 40°, refuge effect) and the northness (i.e., sinus of aspect, potential thermal cover) [75]. We then used a generalized linear mixed-effects model with a binomial distribution to estimate RSF parameters. We included a random intercept for individual’s ID to control for the non-independence and the unbalanced number of observations per individual. Model robustness was assessed using k-fold cross validation [76]. More specifically, 80% of the individuals were randomly selected to fit the model. We then used the model to calculate predicted values for observed and random locations of the remaining 20% of individuals. Predicted values were divided into 10 equal bins such that observed and random locations obtained each a bin rank between 1 and 10. Bin rank (1–10) was then compared to the frequency of observed locations in each bin using Spearman’s rank correlation (\(\overline{{r}_{s}}\)). The procedure was repeated 100 times to obtain mean and range of (\(\overline{{r}_{s}}\)).

We then calculated from the fitted RSF a selection probability for each 13 × 13 m pixel of the study area and retained pixels with a probability above the 75% quantile of probability distribution. We then aggregated adjacent selected pixels to generate resource patches for ibex. We filtered patches on polygons larger than 2000 m2 to take into account the radius of GPS measurement error (i.e., π*252 m2). Finally, once ibex resource patches were identified, we calculated the four indices of landscape fragmentation (i.e., the aggregation, patch cohesion and division indices, and the proportion of cells occupied by habitat patch) to identify to which fragmentation scenario from the simulated landscapes the ibex actual landscape was more similar. We used the R packages adehabitatHR (v0.4.21), sp (v1.6-1), stars (v0.6-1) and raster (v3.6-20) to perform all geographical information system work and lme4 (v1.1-33) and lmerTest (v3.1-3) to run the statistical analysis [52, 62, 63, 73, 77, 78].

Construction of Alpine ibex spatial network at various sampling designs

We used GPS-enabled dead-reckoning relocations obtained every second (1 Hz) for 6 ibex to first build a very high-resolution spatial network. As for the simulation study, we identified all movement steps that occurred between two resource patches and used unique inter-patch movements to build the network. We then regularly resampled each trajectory to mimic various sampling designs and built the resulting spatial networks. Notably, we resampled each trajectory by joining one location out of 30 (30 s), 60 (1 min), 300 (5 min), 600 (10 min), 1200 (20 min), 1800 (30 min), 3600 (1 h), 7200 (2 h) and 21,600 (6 h). We randomly selected the first relocation to start the resampling for each ibex and each sampling design to generate stochasticity. For example, when relocations were sampled every 30 s, the relocation to start with was randomly selected between the first 30 available relocations while for the sampling at 1 h, the relocation to start with was selected between the first 3600. In addition, for the 2 and 6 h sampling design, we also shifted the sampling time by one hour every day to have a complete sample of all possible hours. For example, for the 6 h sampling design, we sampled relocations day 1 at 00:00, 06:00, 12:00, 18:00 and 24:00, then day 2 at 01:00, 07:00, 13:00, 19:00, 01:00, then day 3 at 02:00, 08:00 and so on.

Evaluation of Alpine ibex spatial network

We performed the same analysis on ibex spatial networks as in the simulation study. We first evaluated the effect of sampling frequency on node degree by fitting a linear model to all log-transformed degrees as a function of the log-transformed number of relocations used to build the networks. We also assessed whether the networks had the scale-free and small-world properties by estimating the scale parameter of each degree distribution and by comparing their clustering coefficient and characteristic path length to the ones of random graphs, respectively (see section Effects of landscape fragmentation and sampling frequency on network topology). Secondly, we built the four commonly used theoretical spatial networks (i.e., MPG, DIST50, DIST95 and random graph) to predict connectivity among ibex resource patches. We restricted the predictions to patches that were encompassed within the 100% minimum convex polygon of the very high-resolution relocations (i.e., every second) of the 6 ibex. We then compared predictions of each theoretical network to the empirical observation of ibex spatial network derived from the very high-resolution trajectory (1 Hz), using Cohen’s Kappa.

In addition, we also calculated, for each network obtained from the resampled trajectories, the proportion of visited patches detected, the proportion of realized and spurious links detected. To do so, we used the spatial network built from the 1 Hz trajectory as the reference and extracted from it the lists of visited patches and realized links. We then evaluated whether these visited patches and links were found again in resampled networks. We defined spurious links as links that were observed in the resampled networks but missing from the reference network.

Results

Simulation study

Effects of landscape fragmentation and sampling frequency on simulated networks

Simulated landscapes showed various fragmentation characteristics according to the four indices of landscape fragmentation (described in Table S2). Each degree distribution calculated from 100 simulated networks decayed as a power function with scale parameter being greater than 1 (Table S3), indicating a strong heterogeneous pattern of connectivity among patches in all simulated networks. Few patches had many links (i.e., the hubs) and most of the remaining patches had few links, typical of scale-free networks (Figs. 2, S1 and S2). The scale-free property persisted even when sampling frequency decreased (i.e., fewer relocations) for all types of landscape fragmentation and all simulated movement processes (Fig. 2, Table S3). However, node degree depended on both sampling frequency and landscape fragmentation for all four simulated movement processes (Fig. 3). Node degree globally increased with the number of relocations but at a stronger rate in highly fragmented landscapes for the movement processes reflecting only foraging (i.e., F) or foraging and predator avoidance (i.e., F + Pe and F + Ps, Fig. 3), implying higher information loss on landscape connectivity in fragmented landscapes when sampling frequency gets coarser. This effect, however, dissipated when foragers had to constrain their movements within their own territory (Fig. 3, F + T).

Fig. 3
figure 3

Estimates and their 95% confidence interval of average node degree of simulated networks. One model was run for each simulated movement process: F: Foraging, F + Pe: Foraging and avoidance of an elusive predator, F + Ps: Foraging and avoidance of a stalking predator, F + T: Foraging and territoriality. Each model included as covariates the level of landscape fragmentation, the number of relocations used to build the spatial network (correlated to the sampling frequency) and the interaction of both. Pseudo-R2 of the models were: F: 0.14; F + Pe: 0.18; F + Ps: 0.18; F + T: 0.26

The emergence of the small-world property was linked to both the simulated movement process and the level of landscape fragmentation (Fig. 4). Indeed, from the full trajectory (i.e., 10,000 relocations), all simulated networks were small-world in highly and medium fragmented landscapes for the movement processes reflecting only foraging (i.e., F) or foraging and predator avoidance (i.e., F + Pe and F + Ps, Fig. 4). This proportion decreased though when landscapes had low levels of fragmentation or when territoriality was at stake in the movement process (F + T, Fig. 4). In addition, a decrease in the number of relocations generally led to a decrease in the proportion of small-world networks at a rate that depended on both the movement process and the fragmentation level (Fig. 4), meaning that a coarser sampling frequency impeded the detection of complex network topologies.

Fig. 4
figure 4

Pourcentage of simulated networks having the small-world property. The small-world property is defined as a relatively high clustering coefficient and a similar characteristic path length compared to a random graph. The number of relocations is correlated to the sampling frequency with higher number of relocations indicating fine sampling frequency. The combination of three levels of landscape fragmentation and four movement processes were simulated: F: Foraging, F + Pe: Foraging and avoidance of an elusive predator, F + Ps: Foraging and avoidance of a stalking predator, F + T: Foraging and territoriality

Adequacy with commonly used theoretical spatial networks

Commonly used theoretical networks showed various levels of accuracy to predict connectivity emerging from different movement processes in landscapes having various levels of fragmentation (Table 1) and never predicted correctly more than half of the landscape connectivity. Indeed, Cohen’s Kappa varied between 0 and 0.49 with lowest values for the random graph and highest values for the DIST50 in all four simulated movement processes and all levels of landscape fragmentation. The accuracy provided by the MPG was quite similar to the one of DIST50 although the DIST50 outperformed the MPG when landscape fragmentation was low and the foraging or foraging and territoriality processes were modelled, or when landscape fragmentation was high and the foraging process was simulated (Table 1). Note that the accuracy of both DIST50 and MPG dropped around 0.25 in highly fragmented landscape when foraging and predator avoidance or foraging and territoriality were simulated (Table 1).

Table 1 Adequacy (Cohen’s Kappa) between theoretical spatial networks and simulated networks from four movement processes

Case study

Identification of resource patches for Alpine ibex

Habitat selection by male ibex in June depended on both the distance to steep slopes and land cover class (Table S4). Males selected notably low levels of ligneous cover and avoided mineral and snow covered areas when compared to grasslands (i.e., reference habitat). Besides, they selected cover that was close to steep slopes (Table S4). Ibex resource patches were highly aggregated as the aggregation index was 78.4 and the patch cohesion index reached 96.2. However, the division index was equal to 0.90 indicating that despite being aggregated, patches were also highly divided in the landscape. The proportion of cells occupied by habitat patches totalled 23.8%. Consequently, the level of fragmentation in ibex landscape was relatively similar to the one of landscapes defined as medium or low fragmention in our simulations, according to these four indices of landscape fragmentation.

Evaluation of Alpine ibex spatial network

Similarly to the simulation results, topology of ibex spatial network was impacted by sampling frequency as node degree increased significantly with the number of relocations used to build individual trajectory (\(\text{log}\left(D\right)=0.043 \left(\pm 0.010\right)*\text{log}\left(N\right), p<0.001, {R}^{2}=0.03\), Fig. 5). The amount of variance explained (\({R}^{2}\)) was however quite low indicating that other factors should impact patch connectivity more strongly. In addition, most ibex networks were scale-free according to the estimation of scale parameter (i.e., \(\widehat{\alpha }>1\), range for scale-free networks [min–max]: 1.03–1.30, \({R}^{2}\) of fitted models [min–max]: 0.47–0.75). However, two networks (i.e., sampling relocations every 30 s and 1 min) did not show scale-free properties (i.e., \({\widehat{\alpha }}_{every 30 sec.}\) = 0.94, R2 = 0.49 and \({\widehat{\alpha }}_{every 1 min.}\)= 0.93, R2 = 0.47) even though their scale parameter were both very close to one. Furthermore, most ibex networks had the small-world property according to the comparison of their clustering coefficient and characteristic path length to the one of random graphs (Table 2). This property was however not detected when sampling frequency was too coarse (i.e., relocations every 6 h, Table 2).

Fig. 5
figure 5

Degree distributions of Alpine ibex spatial networks in Belledonne massif (French Alps) in June 2017. Colors indicate the relocation frequency used to build individual trajectories that generated the spatial network. Lines indicate estimated degree distribution from a power law ajusted to the observed degree distribution (points)

Table 2 Properties of Alpine ibex spatial networks in the Belledonne massif (French Alps) in June 2017

A strong loss of information on animal movement appeared quickly when sampling frequency became coarser (Figs. 6, 7). For example, when relocations were collected every 30 min, we detected 76% of visited patches and 52% of realized links and when relocations were collected every 2 h, these values dropped to 51% and 26%, respectively. This loss of information did not evolve linearly with sampling frequency as a significant amount of information was lost within the range of one relocation every 5 min to one relocation every 30 min (Fig. 7). Sampling frequency also impacted the detected pattern of connectivity as many identified links were actually spurious (Fig. 7). For example, 20% of detected links were spurious when relocations were sampled every 20 min and it raised to 27% when relocations were collected every hour (Fig. 7).

Fig. 6
figure 6

Alpine ibex spatial networks for three sampling schedules in Belledonne massif (French Alps) in June 2017. Resources patches are shown in light grey and black lines indicate inter-patch movements identified from ibex trajectories. Red patches and red dotted links are patches and links, respectively, that are detected from the full trajectory but not with the resampled trajectory. Yellow dotted links represent spurious links that were only detected with the resampled trajectory

Fig. 7
figure 7

Effect of sampling frequency on information loss in Alpine ibex spatial networks. Proportion of visited patches detected (green), proportion of realized links detected (yellow) and proportion of spurious links detected (purple) in Alpine ibex spatial networks in Belledonne massif (French Alps) in June 2017 as a function of relocation sampling frequency

Commonly used theoretical networks showed relatively poor capacity to accurately predict landscape connectivity derived from our empirical dataset on Alpine ibex. Indeed, Cohen’s Kappa equalled 0.29, 0.10, 0.26 and 0.00 for the DIST50, DIST95, MPG and random graphs, respectively. Consequently, the DIST50 predicted more accurately inter-patch movements for ibex than the MPG, even though the difference is weak, but it only predicted successfully about one-third of landscape connectivity.

Discussion

Using both simulated and empirical movement trajectories at very high-resolution of individuals showing different movement behaviours in landscapes with varying levels of fragmentation, we first demonstrated that an accurate assessment of landscape connectivity patterns can be hampered when the relocation sampling frequency is too coarse compared to the temporal dynamics of animal movement. However, we also revealed that the level of landscape fragmentation and the processes driving animal movement can both mitigate the effect of relocation sampling frequency on the detection of accurate landscape connectivity patterns. In addition, we showed that network topologies emerging from different movement behaviours and a wide range of landscape fragmentation were generally complex, and that commonly used theoretical networks, usually of simple topologies, did not perform well to accurately predict landscape connectivity in such environments. Finally, we also pinpointed that a significant amount of information on landscape connectivity is rapidly lost when sampling frequency of individual relocations becomes coarser.

Empirical spatial networks have been built from individual trajectories obtained at varying sampling frequencies without knowing so far the potential effect of trajectory resolution on the assessment of landscape connectivity using networks [6, 17, 18, 35, 36, 79, 80]. Here, we demonstrated that the sampling frequency of individual relocations can notably impede the detection of network small-world property but the strength of this effect was strongly linked to the level of landscape fragmentation. Indeed, a higher number of relocations (i.e., finer sampling frequency) was necessary to detect the small-world property of simulated networks in low fragmented than in highly or medium fragmented landscapes. The characteristics of small-world network result in many connected triplets of nodes (i.e., potential local functional use) with efficient movement within the network (i.e., through shortcuts) [57, 59]. Consequently, the missed detection of such patterns could result in inaccurate understanding of how landscapes are connected and for example, could lead to missed identification of patches or connections used for specific functional attributes (e.g., stopovers or stepping-stones) [31, 41, 81]. Similarly, we did not detect the small-world property of Alpine ibex empirical network when sampling frequency was too coarse (i.e., one location every 6 h). On the other hand, the scale-free property of simulated networks was always detected independently of the grain of relocation sampling frequency. However, the skewness of network degree distribution was impacted by both the level of landscape fragmentation and sampling frequency. Indeed, the local connectivity (i.e., node degree) was generally less heterogeneous (i.e., less skewed) in patchy landscapes when full trajectories were used to infer landscape connectivity but when sampling frequency became coarser, the skewness of degree distribution became stronger. These results indicate that commonly used discontinuous trajectories could lack robustness to infer landscape connectivity in real-world systems.

While habitat patch characteristics can potentially reduce the effect of sampling frequency on detected patterns of connectivity, the mitigation could yet be effective in opposite landscapes (i.e., high vs low fragmentation) depending on network topology (i.e., scale-free or small-world). Notably, in low fragmented landscapes, fine sampling frequency should be set to correctly detect small-world property while in highly fragmented landscapes, fine sampling frequency should be set to accurately assess scale-free property, according to our simulations. Many studies demonstrated that animals adjust their behaviour regarding patch selection and intensity of use as a response to patch environmental attributes [40, 82,83,84,85]. We indeed demonstrated that movements in landscape having different attributes can generate different connectivity patterns and that such patterns actually require fine sampling frequency to be correctly appraised. For example, we showed that about 50% of information on ibex movements in a low to medium fragmented environment was lost when ibex were relocated every 30 min and about 25% of their inferred movements were spurious with this sampling rate. While vertebrate movements are generally drawn from discontinuous individual trajectories collected at a coarser rate than one location per hour [18, 36, 80, 86, 87], our empirical study raised a warning on the use of such movement inference to the assessment of landscape connectivity from spatial networks.

The effect of relocation sampling frequency on detected patterns of landscape connectivity was also largely mitigated when simulated individuals were territorial in all types of fragmented landscapes. In the simulation, territorial individuals constrained their movements to their own territory due to the presence of conspecific scent-markings that were avoided [50]. Consequently, territorial individuals had access to a smaller part of the landscape than non-territorial individuals (i.e., F, F + Pe and F + Ps scenarios) such that movements within their territory were probably more redundant resulting in less information loss from coarser sampling frequencies. While our simulations provide new insights on the interplay between animal movement behaviour, landscape fragmentation and relocation sampling frequency on detected patterns of landscape connectivity, it has yet to be appraised using more empirical studies to provide additional support to our results. It could notably be helpful to further improve the design of studies aiming at evaluating landscape connectivity from empirical networks.

Most of our results, either from empirical or simulated cases, conclude that very high-resolution trajectories should be used to correctly infer patterns of landscape connectivity from empirical networks. New technological developments in animal tracking should help meeting this need. Indeed, a wide range of species, from small to large sizes, can now be tracked using GPS tags that collect more frequent data and over longer periods that in the past (e.g., by relying on solar energy) [88]. GPS relocations can be collected as frequently as every second, but generally at a cost of limited duration of tracking due to the high battery needs [89]. The application of dead-reckoning techniques to biologging and GPS data could also help collecting very high-resolution trajectories over longer periods [90]. DR relies on accelerometry and magnetometry data to calculate speed and heading along individual trajectories that can further be translated into very high-resolution coordinates (e.g., location every second, [68]). As battery constraints are less limited for magnetometers and accelerometers, sensor data can be collected much more frequently than solely GPS data and, when combined together to correct for magnetic drift, proved to be quite performant [68, 90]. In this study, Alpine ibex were equipped with tri-axial accelerometer and tri-axial magnetometer sensors recording at 25 Hz/ 8 Hz, respectively, and we combined those data with GPS data collected every two hours to produce locations every second [68]. DR has largely been used on domestic and aquatic wild species but its application remains rare for terrestrial wild species while it could provide detailed information about the movement paths of animals between GPS fixes and open new opportunities in movement ecology [91,92,93,94,95].

Spatial networks reflecting movements among resource patches have been identified as having complex topologies in a large body of litterature for various species [3, 31, 35, 79, 96]. We suggested that network topologies could be simpler than actually detected due to missed stepping-stone movements from coarse sampling frequency and consequently be more similar to simple topologies like the MPG or distance-based networks. Here, we refuted this hypothesis as all networks obtained from very high-resolution trajectories were actually complex (i.e., always scale-free and most of the time small-world). Observations of complex patterns of landscape connectivity could actually be expected as animal movement results from the complex interplay of animal state, navigation and motion capacity and external factors [43]. Consequently, commonly used theoretical networks did not succeed to accurately predict landscape connectivity in our simulated or empirical landscapes. These network models make general and relatively simple hypotheses on the drivers of animal movement among resource patches and are actually quite commonly used as decision-making tools [29, 97,98,99]. However, we demonstrated that they should be used with caution as they only succeeded to correctly predict between 30 and 50% of landscape connectivity when tested for various movement behaviours in different landscapes. Other tools such as individual based models (IBMs), have been combined with network theory to assess landscape connectivity [100, 101]. An IBM models individual movement in a landscape from parameterized rules of movements and consequently allows to generate mechanistic-based stochastic movement trajectories in heterogeneous environments [96, 102, 103]. Many IBMs have been developped to simulate dispersal behaviour in heterogeneous landscapes (i.e., individual based dispersal models [101]) notably because IBMs can address questions even when empirical knowledge is insufficient for linking individual-level processes to landscape-level patterns [104]. However, the application of IBM to the study of landscape connectivity using networks in the context of home-ranging behaviour (e.g., including recursive movements [105, 106]) remains underexploited and could cope with the urgent need to develop theoretical network models that fit with complex but more realistic functional patterns of connectivity [104].

Conclusions

Animal behaviour, the level of landscape fragmentation and the frequency of individual relocations were shown to have an effect on the assessment of landscape connectivity using spatial networks. In general, very high-resolution trajectories (e.g., relocations collected every 1 min in Alpine ibex) were necessary to accurately identify complex network topologies and avoid the generation of spurious information on landscape connectivity. New technologies providing such high-resolution datasets over long periods, as for example dead-reckoning techniques, should thus grow in the movement ecology sphere to help unravelling many ecological questions and conservation challenges. In addition, commonly used theoretical models should be applied with caution to the study of landscape connectivity in real-world systems as they did not perform well as predictive tools. Instead, other modelling tools should be challenged to assess whether they can correctly generate more realistic functional patterns of landscape connectivity in actual environments.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

DIST50:

Distance-based network where patches closer than the median realized link length are connected

DIST95:

Distance-based network where patches closer than the 95% quantile of realized link length are connected

DR:

Dead-reckoning

GPS:

Global positioning system

IBM:

Individual-based model

MPG:

Minimum planar graph

References

  1. Henry M, Pons J-M, Cosson J-F. Foraging behaviour of a frugivorous bat helps bridge landscape connectivity and ecological processes in a fragmented rainforest. J Anim Ecol. 2007;76:801–13.

    Article  PubMed  Google Scholar 

  2. Jacoby DMP, Brooks EJ, Croft DP, Sims DW. Developing a deeper understanding of animal movements and spatial dynamics through novel application of network analyses. Methods Ecol Evol. 2012;3:574–83.

    Article  Google Scholar 

  3. Prima M-C, Duchesne T, Fortin A, Rivest L-P, Fortin D. Combining network theory and reaction–advection–diffusion modelling for predicting animal distribution in dynamic environments. Methods Ecol Evol. 2018;9:1221–31.

    Article  Google Scholar 

  4. Ferrari JR, Preisser EL, Fitzpatrick MC. Modeling the spread of invasive species using dynamic network models. Biol Invasions. 2014;16:949–60.

    Article  Google Scholar 

  5. Kauffman MJ, Aikens EO, Esmaeili S, Kaczensky P, Middleton A, Monteith KL, et al. Causes, consequences, and conservation of ungulate migration. Ann Rev Ecol Evol Syst. 2021;52:453–78.

    Article  Google Scholar 

  6. Courbin N, Fortin D, Dussault C, Courtois R. Logging-induced changes in habitat network connectivity shape behavioral interactions in the wolf–caribou–moose system. Ecol Monogr. 2014;84:265–85.

    Article  Google Scholar 

  7. Santos M, Cagnolo L, Roslin T, Marrero HJ, Vázquez DP. Landscape connectivity explains interaction network patterns at multiple scales. Ecology. 2019;100:e02883.

    Article  PubMed  Google Scholar 

  8. Hanski I, Schulz T, Wong SC, Ahola V, Ruokolainen A, Ojanen SP. Ecological and genetic basis of metapopulation persistence of the Glanville fritillary butterfly in fragmented landscapes. Nat Commun. 2017;8:14504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Ziółkowska E, Ostapowicz K, Kuemmerle T, Perzanowski K, Radeloff VC, Kozak J. Potential habitat connectivity of European bison (Bison bonasus) in the Carpathians. Biol Conserv. 2012;146:188–96.

    Article  Google Scholar 

  10. Belote RT, Barnett K, Zeller K, Brennan A, Gage J. Examining local and regional ecological connectivity throughout North America. Landsc Ecol. 2022;37:2977–90.

    Article  Google Scholar 

  11. Chauveau V, Garel M, Toïgo C, Anderwald P, Bunz Y, Bouche M, Cagnacci F, Canut M, Cavailhes J, Champly I, Filli F, Frey-Roos A, Gressmann G, Herfindal I, Jurgeit F, Martinelli L, Papet R, Petit E, Ramanzin M, Semenzato P, Vannard E, Loison A, Coulon A, Marchand P. Identifying the environmental drivers of corridors and predicting connectivity between seasonal ranges in multiple populations of Alpine ibex (Capra ibex) as tools for conserving migration. Divers Distrib. 2024;e13894.

  12. Lookingbill T, Carter S, Gorsira B, Kingdon C. Using landscape analysis to evaluate ecological impacts of battlefield restoration. Park Sci. 2008;25:60–5.

    Google Scholar 

  13. Bergsten A, Zetterberg A. To model the landscape as a network: a practitioner’s perspective. Landsc Urban Plan. 2013;119:35–43.

    Article  Google Scholar 

  14. Urban D, Keitt T. Landscape connectivity: a graph-theoretic perspective. Ecology. 2001;82:1205–18.

    Article  Google Scholar 

  15. O’Brien D, Manseau M, Fall A, Fortin M-J. Testing the importance of spatial configuration of winter habitat for woodland caribou: an application of graph theory. Biol Conserv. 2006;130:70–83.

    Article  Google Scholar 

  16. Cook Z, Franks DW, Robinson EJH. Efficiency and robustness of ant colony transportation networks. Behav Ecol Sociobiol. 2013;68:509–17.

    Article  Google Scholar 

  17. Wittemyer G, Keating LM, Vollrath F, Douglas-Hamilton I. Graph theory illustrates spatial and temporal features that structure elephant rest locations and reflect risk perception. Ecography. 2017;40:598–605.

    Article  Google Scholar 

  18. Kot CY, Åkesson S, Alfaro-Shigueto J, Amorocho Llanos DF, Antonopoulou M, Balazs GH, et al. Network analysis of sea turtle movements and connectivity: a tool for conservation prioritization. Divers Distrib. 2022;28:810–29.

    Article  Google Scholar 

  19. Fall A, Fortin M-J, Manseau M, O’Brien D. Spatial graphs: principles and applications for habitat connectivity. Ecosystems. 2007;10:448–61.

    Article  Google Scholar 

  20. Minor ES, Urban DL. A graph-theory framework for evaluating landscape connectivity and conservation planning. Conserv Biol. 2008;22:297–307.

    Article  PubMed  Google Scholar 

  21. Saura S, Pascual-Hortal L. A new habitat availability index to integrate connectivity in landscape conservation planning: comparison with existing indices and application to a case study. Landsc Urban Plan. 2007;83:91–103.

    Article  Google Scholar 

  22. Rayfield B, Fortin M-J, Fall A. Connectivity for conservation: a framework to classify network measures. Ecology. 2011;92:847–58.

    Article  PubMed  Google Scholar 

  23. Martensen AC, Saura S, Fortin M-J. Spatio-temporal connectivity: assessing the amount of reachable habitat in dynamic landscapes. Methods Ecol Evol. 2017;8:1253–64.

    Article  Google Scholar 

  24. Reunanen P, Fall A, Nikula A. Spatial graphs as templates for habitat networks in boreal landscapes. Biodivers Conserv. 2012;21:3569–84.

    Article  Google Scholar 

  25. Boccaletti S, Latora V, Moreno Y, Chavez M, Hwang D-U. Complex networks: Structure and dynamics. Phys Rep. 2006;424:175–308.

    Article  Google Scholar 

  26. Barthélemy M. Spatial networks. Phys Rep. 2011;499:1–101.

    Article  Google Scholar 

  27. Barabasi AL, Albert R. Emergence of scaling in random networks. Science. 1999;286:509–12.

    Article  CAS  PubMed  Google Scholar 

  28. Albert R, Jeong H, Barabasi A-L. Error and attack tolerance of complex networks. Nature. 2000;406:378–82.

    Article  CAS  PubMed  Google Scholar 

  29. Leblond M, Rudolph T, Boisjoly D, Dussault C, St-Laurent M-H. Science-informed policy decisions lead to the creation of a protected area for a wide-ranging species at risk. Conserv Sci Pract. 2022;4:e12833.

    Article  Google Scholar 

  30. Minor ES, Lookingbill TR. A multiscale network analysis of protected-area connectivity for mammals in the United States. Conserv Biol. 2010;24:1549–58.

    Article  PubMed  Google Scholar 

  31. Fox RJ, Bellwood DR. Herbivores in a small world: network theory highlights vulnerability in the function of herbivory on coral reefs. Funct Ecol. 2014;28:642–51.

    Article  Google Scholar 

  32. Cagnacci F, Boitani L, Powell RA, Boyce MS. Animal ecology meets GPS-based radiotelemetry: a perfect storm of opportunities and challenges. Philos Trans R Soc Lond B Biol Sci. 2010;365:2157–62.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Wilson RR, Horne JS, Rode KD, Regehr EV, Durner GM. Identifying polar bear resource selection patterns to inform offshore development in a dynamic and changing Arctic. Ecosphere. 2014;5:art136.

    Article  Google Scholar 

  34. Hofman MPG, Hayward MW, Heim M, Marchand P, Rolandsen CM, Mattisson J, et al. Right on track? Performance of satellite telemetry in terrestrial wildlife research. PLoS ONE. 2019;14:e0216223.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Prima M-C, Duchesne T, Fortin A, Rivest L-P, Drapeau P, St-Laurent M-H, et al. A landscape experiment of spatial network robustness and space-use reorganization following habitat fragmentation. Funct Ecol. 2019;33:1663–73.

    Article  Google Scholar 

  36. Almpanidou V, Mazaris AD, Mertzanis Y, Avraam I, Antoniou I, Pantis JD, et al. Providing insights on habitat connectivity for male brown bears: a combination of habitat suitability and landscape graph-based models. Ecol Model. 2014;286:37–44.

    Article  Google Scholar 

  37. Jaeger JAG. Landscape division, splitting index, and effective mesh size: new measures of landscape fragmentation. Landsc Ecol. 2000;15:115–30.

    Article  Google Scholar 

  38. Sánchez-de-Jesús HA, Arroyo-Rodríguez V, Andresen E, Escobar F. Forest loss and matrix composition are the major drivers shaping dung beetle assemblages in a fragmented rainforest. Landsc Ecol. 2016;31:843–54.

    Article  Google Scholar 

  39. Cresswell JE, Osborne JL. The effect of patch size and separation on bumblebee foraging in oilseed rape: implications for gene flow. J Appl Ecol. 2004;41:539–46.

    Article  Google Scholar 

  40. Courant S, Fortin D. Search efficiency of free-ranging plains bison for optimal food items. Anim Behav. 2012;84:1039–49.

    Article  Google Scholar 

  41. Saura S, Bodin Ö, Fortin M-J. Stepping stones are crucial for species’ long-distance dispersal and range expansion through habitat networks. J Appl Ecol. 2014;51:171–82.

    Article  Google Scholar 

  42. Moorcroft PR, Lewis MA, Crabtree RL. Mechanistic home range models capture spatial patterns and dynamics of coyote territories in Yellowstone. Proc R Soc Lond B Biol Sci. 2006;273:1651–9.

    Google Scholar 

  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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Laundré JW. Behavioral response races, predator-prey shell games, ecology of fear, and patch use of pumas and their ungulate prey. Ecology. 2010;91:2995–3007.

    Article  PubMed  Google Scholar 

  45. Courbin N, Loveridge AJ, Macdonald DW, Fritz H, Valeix M, Makuwe ET, et al. Reactive responses of zebras to lion encounters shape their predator–prey space game at large scale. Oikos. 2016;125:829–38.

    Article  Google Scholar 

  46. Bateman AW, Lewis MA, Gall G, Manser MB, Clutton-Brock TH. Territoriality and home-range dynamics in meerkats, Suricata suricatta: a mechanistic modelling approach. J Anim Ecol. 2015;84:260–71.

    Article  PubMed  Google Scholar 

  47. Courbin N, Garel M, Marchand P, Duparc A, Debeffe L, Börger L, et al. Interacting lethal and nonlethal human activities shape complex risk tolerance behaviors in a mountain herbivore. Ecol Appl. 2022;32:e2640.

    Article  PubMed  Google Scholar 

  48. Bracis C, Gurarie E, Van Moorter B, Goodwin RA. Memory effects on movement behavior in animal foraging. PLoS ONE. 2015;10:e0136057.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Bracis C, Gurarie E, Rutter JD, Goodwin RA. Remembering the good and the bad: memory-based mediation of the food–safety trade-off in dynamic landscapes. Theor Ecol. 2018;11:305–19.

    Article  Google Scholar 

  50. Theng M, Milleret C, Bracis C, Cassey P, Delean S. Confronting spatial capture–recapture models with realistic animal movement simulations. Ecology. 2022;103:e3676.

    Article  PubMed  Google Scholar 

  51. Deutsch CV, Journel AG. GSLIB: Geostatistical software library and user’s guide. Oxford: Oxford University Press; 1992.

    Google Scholar 

  52. Bivand RS, Pebesma E, Gómez-Rubio V. Applied spatial data analysis with R. New York: Springer; 2013.

    Book  Google Scholar 

  53. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing; 2023.

  54. Pebesma E. Multivariable geostatistics in S: the gstat package. Comput Geosci. 2004;30:683–91.

    Article  Google Scholar 

  55. Hesselbarth MHK, Sciaini M, With KA, Wiegand K, Nowosad J. landscapemetrics: an open-source R tool to calculate landscape metrics. Ecography. 2019;42:1648–57.

    Article  Google Scholar 

  56. Lee U, Mashour GA. The role of network science in the study of anesthetic state transitions. Anesthesiology. 2018;129:1029–44.

    Article  PubMed  Google Scholar 

  57. Watts DJ, Strogatz SH. Collective dynamics of ‘small-world’ networks. Nature. 1998;393:440–2.

    Article  CAS  PubMed  Google Scholar 

  58. Newman MEJ. The structure and function of complex networks. SIAM Rev. 2003;45:167–256.

    Article  Google Scholar 

  59. Telesford QK, Joyce KE, Hayasaka S, Burdette JH, Laurienti PJ. The ubiquity of small-world networks. Brain Connect. 2011;1:367–75.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Humphries MD, Gurney K. Network ‘small-world-ness’: A quantitative method for determining canonical network equivalence. PLoS ONE. 2008;3:e0002051.

    Article  PubMed  Google Scholar 

  61. Csardi G, Nepusz T. The igraph software package for complex network research. InterJournal. 2005;Complex Systems:1695.

  62. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.

    Article  Google Scholar 

  63. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest Package: tests in linear mixed effects models. J Stat Softw. 2017;82:1–26.

    Article  Google Scholar 

  64. Cohen J. A coefficient of agreement for nominal scales. Educ Psychol Meas. 1960;20:37–46.

    Article  Google Scholar 

  65. Fletcher RJ, Acevedo MA, Reichert BE, Pias KE, Kitchens WM. Social network models predict movement and connectivity in ecological landscapes. PNAS. 2011;108:19282–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Villaret JC, Bon R, Rivet A. Sexual segregation of habitat by the Alpine ibex in the French Alps. J Mammal. 1997;78:1273–81.

    Article  Google Scholar 

  67. Grignoli S, Rossi I, Bertolotto E, Bassano B, Apollonio M. Influence of the kid on space use and habitat selection of female Alpine ibex. J Wildl Manag. 2007;71:713–9.

    Article  Google Scholar 

  68. Redcliffe J. Alpine ungulate movement: Quantification of spatiotemporal environmental energetics and social interaction [Degree of Doctor of Philosophy in Biological Sciences]. Swansea University; 2021.

  69. Gunner RM, Holton MD, Scantlebury MD, van Schalkwyk OL, English HM, Williams HJ, et al. Dead-reckoning animal movements in R: a reappraisal using Gundog. Tracks. Anim Biotelemetry. 2021;9:23.

    Article  Google Scholar 

  70. Gunner RM, Holton MD, Scantlebury DM, Hopkins P, Shepard ELC, Fell AJ, et al. How often should dead-reckoned animal movement paths be corrected for drift? Anim Biotelemetry. 2021;9:43.

    Article  PubMed  PubMed Central  Google Scholar 

  71. Sanz T, Villaret JC. Catalogue des végétations de l’Isère. Classification physionomique et phytosociologique avec clés de détermination. Conservatoire botanique national alpin, Ministère de la Transition écologique et solidaire; 2018.

  72. Manly BFJ, McDonald LL, Thomas D, McDonald TL, Erickson WP. Resource selection by animals: statistical design and analysis for field studies. Kluwer Academic Dordrecht; 2002.

  73. Calenge C. The package “adehabitat” for the R software: A tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197:516–9.

    Article  Google Scholar 

  74. Horne JS, Garton EO, Krone SM, Lewis JS. Analyzing animal movements using Brownian bridges. Ecology. 2007;88:2354–63.

    Article  PubMed  Google Scholar 

  75. Margaillan L. Simultaneous GPS monitoring during summer reveals habitat selection in male Alpine ibex is shaped by resource and interference competition with sheep herds. OFB, EPHE, PSL; 2021.

  76. Boyce MS, Vernier PR, Nielsen SE, Schmiegelow FKA. Evaluating resource selection functions. Ecol Model. 2002;157:281–300.

    Article  Google Scholar 

  77. Hijmans R, van Etten J. Raster: Raster: Geographic data analysis and modeling. R Package Version. 2014;517:2–12.

    Google Scholar 

  78. Bivand EP Roger. Spatial data science: With applications in R. Chapman and Hall/CRC; 2023.

  79. Rhodes M, Wardell-Johnson GW, Rhodes MP, Raymond B. Applying network analysis to the conservation of habitat trees in urban environments: a case study from Brisbane, Australia. Conser Biol. 2006;20:861–70.

    Article  Google Scholar 

  80. Pasquaretta C, Dubois T, Gomez-Moracho T, Delepoulle VP, Le Loc’h G, Heeb P, et al. Analysis of temporal patterns in animal movement networks. Methods Ecol Evol. 2021;12:101–13.

    Article  Google Scholar 

  81. Monteith KL, Hayes MM, Kauffman MJ, Copeland HE, Sawyer H. Functional attributes of ungulate migration: landscape features facilitate movement and access to forage. Ecol Appl. 2018;28:2153–64.

    Article  PubMed  Google Scholar 

  82. McMahon LA, Rachlow JL, Shipley LA, Forbey JS, Johnson TR. Habitat selection differs across hierarchical behaviors: selection of patches and intensity of patch use. Ecosphere. 2017;8:e01993.

    Article  Google Scholar 

  83. Dzialak MR, Olson CV, Webb SL, Harju SM, Winstead JB. Incorporating within- and between-patch resource selection in identification of critical habitat for brood-rearing greater sage-grouse. Ecol Process. 2015;4:5.

    Article  Google Scholar 

  84. Lesmerises R, Ouellet J-P, Dussault C, St-Laurent M-H. The influence of landscape matrix on isolated patch use by wide-ranging animals: conservation lessons for woodland caribou. Ecol Evol. 2013;3:2880–91.

    Article  PubMed  PubMed Central  Google Scholar 

  85. Searle KR, Thompson Hobbs N, Shipley LA. Should I stay or should I go? Patch departure decisions by herbivores at multiple scales. Oikos. 2005;111:417–24.

    Article  Google Scholar 

  86. Prima M-C, Duchesne T, Merkle JA, Chamaillé-Jammes S, Fortin D. Multi-mode movement decisions across widely ranging behavioral processes. PLoS ONE. 2022;17:e0272538.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Courbin N, Loveridge AJ, Fritz H, Macdonald DW, Patin R, Valeix M, et al. Zebra diel migrations reduce encounter risk with lions at night. J Anim Ecol. 2019;88:92–101.

    Article  PubMed  Google Scholar 

  88. Kays R, Crofoot MC, Jetz W, Wikelski M. Terrestrial animal tracking as an eye on life and planet. Science. 2015;348:aaa2478.

    Article  PubMed  Google Scholar 

  89. Li J, Brugere I, Ziebart B, Berger-Wolf T, Crofoot M, Farine D. Social information improves location prediction in the wild. Workshops at the Twenty-Ninth AAAI Conference on Artificial Intelligence. 2015.

  90. Bidder OR, Walker JS, Jones MW, Holton MD, Urge P, Scantlebury DM, et al. Step by step: reconstruction of terrestrial animal movement paths by dead-reckoning. Mov Ecol. 2015;3:23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  91. Mitani Y, Sato K, Ito S, Cameron MF, Siniff DB, Naito Y. A method for reconstructing three-dimensional dive profiles of marine mammals using geomagnetic intensity data: results from two lactating Weddell seals. Polar Biol. 2003;26:311–7.

    Article  Google Scholar 

  92. Wensveen PJ, Thomas L, Miller PJO. A path reconstruction method integrating dead-reckoning and position fixes applied to humpback whales. Mov Ecol. 2015;3:31.

    Article  PubMed  PubMed Central  Google Scholar 

  93. Dewhirst OP, Evans HK, Roskilly K, Harvey RJ, Hubel TY, Wilson AM. Improving the accuracy of estimates of animal path and travel distance using GPS drift-corrected dead reckoning. Ecol Evol. 2016;6:6210–22.

    Article  PubMed  PubMed Central  Google Scholar 

  94. Bidder OR, Soresina M, Shepard ELC, Halsey LG, Quintana F, Gómez-Laich A, et al. The need for speed: testing acceleration for estimating animal travel rates in terrestrial dead-reckoning systems. Zoology. 2012;115:58–64.

    Article  PubMed  Google Scholar 

  95. Munden R, Börger L, Wilson RP, Redcliffe J, Loison A, Garel M, et al. Making sense of ultrahigh-resolution movement data: A new algorithm for inferring sites of interest. Ecol Evol. 2019;9:265–74.

    Article  PubMed  Google Scholar 

  96. Lookingbill TR, Gardner RH, Ferrari JR, Keller CE. Combining a dispersal model with network theory to assess habitat connectivity. Ecol Appl. 2010;20:427–41.

    Article  PubMed  Google Scholar 

  97. Bourdouxhe A, Duflot R, Radoux J, Dufrêne M. Comparison of methods to model species habitat networks for decision-making in nature conservation: The case of the wildcat in southern Belgium. J Nat Conserv. 2020;58:125901.

    Article  Google Scholar 

  98. Martinez L. L’aménagement urbain : (En) quête sur la préservation des espaces ouverts : L’exemple de la métropole Grenoble-Alpes [These de doctorat]. Université Grenoble Alpes; 2022.

  99. Albert CH, Rayfield B, Dumitru M, Gonzalez A. Applying network theory to prioritize multispecies habitat networks that are robust to climate and land-use change. Conserv Biol. 2017;31:1383–96.

    Article  PubMed  Google Scholar 

  100. Grimm V, Revilla E, Berger U, Jeltsch F, Mooij WM, Railsback SF, et al. Pattern-oriented modeling of agent-based complex systems: Lessons from ecology. Science. 2005;310:987–91.

    Article  PubMed  Google Scholar 

  101. Diniz MF, Cushman SA, Machado RB, De Marco JP. Landscape connectivity modeling from the perspective of animal dispersal. Landsc Ecol. 2020;35:41–58.

    Article  Google Scholar 

  102. da Rocha ÉG, Brigatti E, Niebuhr BB, Ribeiro MC, Vieira MV. Dispersal movement through fragmented landscapes: the role of stepping stones and perceptual range. Landsc Ecol. 2021;36:3249–67.

    Article  Google Scholar 

  103. Vanlandeghem V, Drapeau P, Prima M-C, St-Laurent M-H, Fortin D. Management-mediated predation rate in the caribou–moose–wolf system: spatial configuration of logging activities matters. Ecosphere. 2021;12:e03550.

    Article  Google Scholar 

  104. Pe’er G, Henle K, Dislich C, Frank K. Breaking functional connectivity into components: A novel approach using an individual-based model, and first outcomes. PLoS ONE. 2011;6:e22355.

    Article  PubMed  PubMed Central  Google Scholar 

  105. Berger-Tal O, Bar-David S. Recursive movement patterns: review and synthesis across species. Ecosphere. 2015;6:art149.

    Article  Google Scholar 

  106. Merkle JA, Potts JR, Fortin D. Energy benefits and emergent space use patterns of an empirically parameterized model of memory-based patch selection. Oikos. 2016;126:185–95.

    Google Scholar 

Download references

Acknowledgements

The authors thank all participants for helping with the capture and marking of ibex and for the recovery of GPS collars: volunteers and professionals from the Office Français de la Biodiversité (Carole Toïgo, Mathieu Beurier and SD38).

Funding

M.-C.P. was funded by the French ministry of environment.

Author information

Authors and Affiliations

Authors

Contributions

M.-C.P., M.G., P.M. and F.B conceived the ideas and designed the methodology; M.-C.P. generated simulated datasets; P.M. and M.G. provided Alpine ibex dataset; L.B. and J.R. generated very high-resolution trajectories of Alpine ibex from dead-reckoning techniques; M.-C.P. performed GIS and statistical analysis of both simulated and empirical datasets; M.G., P.M., L.B. and F.B. secured fundings for the project; M.-C.P wrote first draft of the manuscript and all authors contributed substantially to revisions and gave final approval for publications.

Corresponding author

Correspondence to Marie-Caroline Prima.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Prima, MC., Garel, M., Marchand, P. et al. Combined effects of landscape fragmentation and sampling frequency of movement data on the assessment of landscape connectivity. Mov Ecol 12, 63 (2024). https://doi.org/10.1186/s40462-024-00492-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40462-024-00492-8

Keywords