Abstract
Outdoor recreation is increasing and affects habitat use and selection by wildlife. These effects are challenging to study, especially for elusive species with large spatial requirements, as it is hard to obtain reliable proxies of recreational intensity over extensive areas. Commonly used proxies, such as the density of, or distance to, hiking paths, ignore outdoor recreation occurring on other linear feature types. Here we utilized crowdsourced data from the Strava training app to obtain a large-scale proxy for pedestrian outdoor recreation intensity in southeast Norway. We used the proxy and GPS-tracking data from collared Eurasian lynx (Lynx lynx) to investigate how recreation affects habitat selection at the home range scale and local scale by lynx during summer. We fitted resource selection functions at the two scales using conditional logistic regression. Our analysis revealed that lynx avoided areas of recreational activity at the local scale, but not at home range scale. Nonetheless, lynx frequently used areas associated with recreation, and to a greater degree at night than during the day. Our results suggest that local-scale avoidance of recreation and temporal adjustments of habitat use by lynx mitigate the need for a home range-scale response towards recreation. Scale-dependent responses and temporal adjustments in habitat use may facilitate coexistence between humans and large carnivores.
Similar content being viewed by others
Introduction
Large carnivores living in human-dominated landscapes face widespread land-use changes, infrastructure, and the presence of humans. The ability of large carnivores to share the landscape with humans will be increasingly challenged, as the human population is growing1. Land-use changes and human infrastructure impact the behaviour, habitat use, and habitat selection of large carnivores2. Although less conspicuous, outdoor recreation constitutes another disturbance pressure on wildlife in the human-dominated landscape. Despite being more benign than, for example, hunting and habitat conversion, recreation also affects behaviour, habitat use, and habitat selection3. For some wildlife species the negative impact from recreation is even considered an important conservation issue4 (e.g. Barbary Macaque Macaca sylvanus5 and cheetah Acinonyx jubatus5). Hence, understanding the effects of recreation on wildlife is crucial for mitigating its potential negative impacts at present and in the future.
The vast extent and intensity of outdoor recreation is illustrated by an estimate of 8 billion visits solely to protected areas per year globally and 4 billion visits per year in Europe alone6. Outdoor recreation is rising7 and, in combination with increasing urbanisation8,9, the patterns will likely change in the future, with growing levels of outdoor recreation especially around urban areas. These changes will have implications for wildlife. The effects of recreation on wildlife have received increased research attention throughout recent decades10 and a wide range of effects of recreation on wildlife has been documented3. Generally, animals tend to respond to humans in the same way they respond to a predator11. The detection of a human is usually followed by a behavioural and/or physiological response12,13. These responses might in turn influence energy budgets, habitat use and fitness. Hence, recreation can in the worst case influence population trends and the distribution of animals3.
Habitat selection, i.e. the disproportional use of a habitat type relative to its availability, can be inferred at multiple scales based on how availability is defined14. Animals may respond to the same habitat feature differently based on the scale of habitat selection15. For instance animals may avoid large patches of habitat that contain many linear features or receive high levels of recreation at the home range-scale16 (cf. 3rd order of habitat selection14). Animals may also avoid specific linear features that are associated with high levels of recreation at a local scale17,18 (cf. 4th order of habitat selection14), even if not avoided at a larger scale. To understand at which scale animals respond to recreation is therefore important for coexistence of humans and wildlife, as the ability to respond to recreation at the local scale19 may mitigate disturbance effects at the home range scale.
Large carnivores are now recovering in Europe from past persecution and returning to a human-dominated landscape20. Despite large carnivores being apex predators, humans influence their behaviour21 and ecosystem function22. Numerous studies have shown behavioural and physiological responses in large carnivores when humans directly approach on foot (e.g. brown bears Ursus arctos13,23, black bears Ursus americanus24, wolf Canis lupus25 and puma Puma concolor26), but less attention has been paid to how habitat selection and habitat use by large carnivores are affected by human recreational activity and not only the infrastructure associated with recreation. Recreation can reduce habitat quality27 and spatial avoidance of suitable habitat due to recreation can be functionally equivalent to habitat loss16,28. Large-scale segregation due to recreation has been documented for group-living herbivores in open landscapes29 and for wolverines (Gulo gulo) in mountain areas16 due to functional habitat loss. However, it is unclear how large carnivores, with large space requirements, respond to outdoor recreation over different scales in forested habitats.
The paucity of studies on the effect of recreation on habitat selection by large carnivores thus far is likely in part due to the lack of measures of human activity or recreation for the large areas in which these species live. Path density or distance to paths have been commonly used as proxies for human recreation10,30. These proxies ignore the fact that different path segments are associated with different levels of recreation and that other linear features than paths are also associated with recreation, such as forest roads and public roads. Estimating the level of recreation along linear features with, for example, human counters31 or camera traps requires vast resources if it is to be done at a scale relevant for large carnivores. Today, ubiquitous smartphones and smartwatches with built-in GPS-loggers have opened up new possibilities to obtain proxies for human activity at large spatial scales32. Users of certain software applications (apps), like the training app Strava (www.strava.com), agree to share their spatial locations with the company. For apps with large userbases, such data can provide relative proxies for the spatial distribution of recreational activity33,34. Here we utilize data from the Strava app as a proxy for pedestrian outdoor activity (walking, running, or hiking) during summer in southeast Norway and investigate how habitat selection and habitat use by Eurasian lynx (Lynx lynx) are influenced by recreation. Furthermore, we included density of hiking paths as a proxy for recreation and tested whether lynx responded differently to these two proxies of human recreation.
The Scandinavian lynx population has recovered after being hunted to the edge of extirpation in the mid-twentieth century35. The population continues to be exposed to extensive legal hunting, and human-caused sources of mortality are high (poaching and vehicle collisions in addition to hunting)36. Hence, there is a potential for strong avoidance of humans in this population. We studied the effect of summer recreation on habitat selection by lynx with resource selection functions (RSF37) at two spatial scales with different availability definitions; within the home range (home range-scale) and within a buffer of 1–2 km around locations used by the lynx (local-scale). If the lynx manage to coexist with humans by locally avoiding features with high recreational activity, they can still occupy human-dominated landscapes, which will facilitate coexistence compared to a total avoidance of humans. Furthermore, we investigated whether habitat selection was influenced by the time of day. Lastly, we explored how lynx habitat use of areas associated with recreation changed throughout the day.
Methods
Study area
The study area (approximately 43,000 km2) is located in southeast Norway (centroid coordinate: N 59.96982, E 9.693853), in Innlandet, Vestfold og Telemark, Oslo, and Viken counties (Fig. 1). The study area includes the most heavily populated areas in Norway, and has an overall human population of 2 million (www.ssb.no). In the southeast part of the study area, the landscape consists of forest fragmented with agricultural land and settlements, and a rolling topography. In the northwest, the topography is characterized by steep slopes, valley systems and some agricultural land along the valley floors. Most of the forests are heavily exploited by commercial forest industry and associated clearcut practices. The main tree species are Norway spruce (Picea abies), Scots pine (Pinus sylvestris) and birch (Betula sp.). Roe deer (Capreolus capreolus), free-ranging domestic sheep (Ovis aries), red deer (Cervus elaphus), and small prey species such as mountain hare (Lepus timidus), tetraonids and other birds comprise the diet of the lynx in the area38,39.
Animal capture and GPS-data
From 2008 to 2014, 25 lynx (11 females and 14 males) were captured in foot snares or wooden box traps and equipped with GPS-collars following a pre-established protocol40. All capture and handling procedures were approved by the Norwegian Experimental Animal Ethics Committee (permit numbers 2012/206992, 2010/161554, 2010/161563, 08/127430, 07/81885, 07/7883). This study was conducted in compliance with the ARRIVE guidelines, and all our methods were performed in accordance with relevant guidelines and regulations. In this study, we only analysed data from adult lynx in years when they were settled in a home range. GPS-collars were programmed to take between 1–19 GPS-locations per day. Fix schedules varied across individuals and alternated between intensive predation study periods and less intensive monitoring periods39. As the aim of this study was to explore how recreation affects habitat selection by resident lynx during the snow-free period, we included observations from 1st May to 31st October during all all years. To increase the robustness of the analyses, we only included combinations of lynx individuals and years (i.e. lynx-years) with ≥ 200 GPS-locations. This resulted in a final dataset of 13 611 GPS-locations for 22 lynx-years from 20 individual lynx (8 females and 12 males). Two females were included during two years each.
Recreation
We used crowdsourced human mobility data from Strava as a proxy for recreation. Strava is an app for smartphones and smartwatches, used primarily to record and upload georeferenced human training activities (hereafter referred to as “activity event”). Activity events can also be uploaded directly to the Strava webpage. Strava stores this data and processed version of it and can be accessed from Strava Metro. To maintain anonymity and conform with privacy regulations, access is limited to data processed by Strava through removal of personal identifiers and spatial and/temporal aggregation (Fig. 2). The processing involves linking individual activity events to nearby linear features (paths, roads etc.) in OpenStreetMap (OSM, www.openstreetmap.org). Hence, in the absence of OSM linear features close to an activity event (or parts of it), that event (or parts of it) is not included in the aggregated version of the dataset (Fig. 2). Further, to conform with privacy legislations, linear features with less than three unique users are removed and the number of activity events are rounded up to the nearest multiple of five. Based on our initial inspections of the Strava data, pedestrian activity (walking, running, hiking etc.) appeared to extend further into the forest compared to biking. Hence, we decided to focus our analysis only on pedestrian activities. We used Strava activities from pedestrians at a temporal resolution of one year, i.e., data where the timespan for at least three unique Strava users on a linear feature was a year, to maximise the spatial coverage. We included Strava data from 2016 to 2019 and summed this over all years to result in a single static covariate. A more detailed description of how we processed the Strava data can be found in Supplementary material S1 and Figure S1.
We also used the density of hiking paths as a proxy for recreation to facilitate comparison of a more traditional method with the Strava data. We obtained data on hiking paths from the OSM (key = highway and value = track and path) using the osmdata package41 in R (version 4.0.342). To derive our proxies for recreation, we divided our study area into a grid of 50 × 50 m (a compromise between computational time and highest possible resolution). For hiking paths, we calculated the line lengths of hiking paths inside each grid cell. For the Strava data, we summed up all the pedestrian activity inside each grid cell. We refer to this covariate as the Strava index. The value of this index in a given grid cell can be viewed as the total number of activity events taking place on all linear features within that grid cell during the entire period when the Strava data were collected. The Strava index accounts for different levels of human activity associated with all types of linear features (not only hiking paths) inside a grid cell, while the hiking path density is only a measure of the length of the hiking path inside a grid cell and does not account for the level of associated human activity.
Resource selection functions at two scales
We used resource selection functions (RSF37) to estimate habitat selection by lynx. RSF depends upon a use-available design, where locations used by GPS-collared lynx are compared to locations not used but considered as available (“null model”). For each used GPS-location, we randomly sampled 30 available locations (excluding open water). Based on our sensitivity analysis43 the choice of 30 available locations per used location appeared to be sufficient (see Figure S2–S3). We had two different availability definitions (Fig. 3). For the home range-scale habitat selection we sampled available locations inside the home range. Home ranges were estimated with the Brownian bridge kernel method44 with the kernelbb function in the adehabitathr package45 in R. We used the 95% isopleths. For the local-scale we sampled available locations inside a buffer with the mean step lengths (distance between two consecutive GPS-locations) of the individual as radius from the used GPS-location (Table S1). We then fitted conditional logistic regression to obtain the RSFs by using the coxph function in the survival package46 in R. The response was whether a location was used (1) or available (0). Conditional logistic regression relies on a matched design, where groups of observations are matched with given grouping IDs. For the home range-scale we matched used and available locations by including the lynx-year ID as strata (not to be confused with the Strava index) and for the local-scale we used an unique identifier for the buffer as the strata. This ensured that used locations in a home range were only compared to available locations sampled from the same home range, and also that the used location in a buffer was only compared to the available locations in the same buffer. Furthermore, we included the lynx-year ID as a cluster variable in the models to obtain robust standard error estimates for the coefficients47,48. We interpreted robust confidence intervals overlapping 0 as a lack of evidence for either avoidance or selection.
Covariates
Based on findings previously reported for lynx30,49,50,51, we considered the following covariates: slope, agricultural land (hereafter referred to as fields), forest, forest roads, public roads, and houses. In addition, we also considered the Strava index and hiking path density as proxies for recreation. Spatial data on fields and forest were obtained from an open access land cover map (AR5052, 1:50,000); and vector data on public roads, forest roads and houses from the Norwegian Mapping Authority (www.geonorge.no). We defined a house as a building approved as a residential building for the entire year. This excludes recreational cabins, which are associated with a variable degree of human activity. Slope was calculated with the terrain function (based on the 8 nearest neighbouring raster cells) in the raster package in R53, based on a digital elevation model with 50 m resolution from the Norwegian Mapping Authority (www.geonorge.no). All covariates were rasterized to a 50 × 50 m resolution (see Table S2 for simple summary of the covariates).
We prepared the covariates differently for the home range- and local-scale (see Table S2), to account for the different definitions of available habitat. For the home range-scale analysis we calculated the density of houses, forest roads and public roads, the proportion of forest, and the proportion of fields in a buffer of 1 km radius around each cell (50 × 50 m). For hiking path density we calculated the sum of all cells in a buffer of 1 km radius around the focal cell, while for the Strava index and slope we used the mean. We denote the hiking path density, Strava index, and slope with the subscript 1000 (e.g., Strava index1000) to separate them from the covariates with similar names used for the local-scale. For the local-scale analysis we used distance to house, distance to forest road, distance to public road and distance to field, instead of density or proportion metrics. We included forest as a binary variable (1 = forest and 0 = not forest). For the Strava index and hiking path density we calculated the mean (Strava index) or sum (hiking path) of the covariate in the four closest neighbouring cells and the focal cell. This ensured that the hiking path density and Strava index extended at least 50 m outside the linear feature they were associated with. For slope on the local-scale we used the original calculation as previously described. We denote the Strava index, hiking path density and slope used for the local-scale analysis with the subscript 50 (e.g., Strava index50) to separate them from the covariates used in the home range-scale analysis. In addition, we also considered day vs. night as a covariate in models at the local-scale. Night was defined as the time between sunset and sunrise, obtained by the sunriset function in R-package maptools54. We standardized all the continuous covariates by subtracting the mean and dividing by the standard deviation of all used and available locations prior to fitting the models. The distance to feature covariates were log-transformed, to make the effect of the covariate decrease with the distance from the feature (we added 1 m to all distances prior to log-transformation).
Candidate models and model selection
We considered four different candidate models for the home range-scale and eight different models for the local-scale analysis (see Table 1). The simplest model (“core model”) contained parameters that have previously been shown to be important for habitat selection by lynx. For both scales, we tested whether including recreation covariates improved our core model by adding the density of hiking paths (“path model”), the Strava index (“Strava model”) or both the density of hiking paths and the Strava index (“full model”). In the core model for the home range-scale we included the quadratic term of the proportion of forest as the results of Bouyer et al.51 suggested selection for an optimum less than 100% forest cover. In addition, we tested if local-scale habitat selection was influenced by time of day by including an interaction with night (true/false) for slope and forest as Filla, et al.30 showed that lynx select gentler slopes and spend more time in open habitats during the night. Given the diurnal activity of humans, we also included an interaction with night for the following human related covariates: distance to public and forest roads, hiking path density50, and the Strava index50. We did not include a main effect of night in the model as this variable was constant for any given stratum (single used location and associated sample of available locations) and could therefore not be meaningfully evaluated. We used Akaike’s Information Criterion (AIC55) to rank models with the strongest support, and based inferences on models ≤ 2 delta AIC of the top ranked model56. We interpreted that including recreational activity improved our models when either the full-, Strava-, or path model ranked highest or ≤ 2 delta AIC off the top-ranked model.
Time dependent habitat use of areas associated with Strava use
In addition to habitat selection, we explored how the habitat use of areas associated with the Strava index50 changed throughout the day. This was done by fitting a generalized additive mixed model (GAMM57) to only the GPS locations used by the lynx. As response, we used the GPS-locations of the lynx and coded them as 1 if the location had Strava activity (Strava index50 > 0) and 0 if not. We used a binomial distribution with a logit link to model the response. The predictors included hour of the day with a cyclic spline as a smoothing term and lynx ID as a random intercept. We used the gamm function in mgcv package58 in R to fit the model. The time of day was corrected for differences in day length using two anchors (one at sunrise and one at sunset) and the average method59 in the activity package60 in R.
Results
Mean human density within lynx home ranges was 26 km-2 (range: 4.8–166 km-2). On average, 12% (range: 0.8–31%) of the used GPS-locations per lynx, and 16% (range: 5.7–37%) of lynx home range areas, were located in grid cells with Strava index50 > 0. Corresponding values for Strava index1000 > 0 were 82% (range: 40–99.5%) of the used GPS-locations per lynx and 79% (46–98%) of their home range area.
Of the Strava activity events inside the lynx home ranges, 58% were located in forest, 24% occurred in built-up areas (urban, sub-urban, small towns etc.), 11% on fields, and 5% in alpine areas (open areas above the forest). Inside the lynx home ranges, roads were the linear features that had the highest levels of activity events, with hiking paths receiving fewer activity events (Figure S4).
Home range-scale habitat selection
The four candidate models differed only moderately from each other in terms of AIC (Table S4). The model including hiking paths ranked highest, but the full model was within 2 ΔAIC. We therefore present the results from the full model. In addition, the direction of effects (when included) was similar in both models.
We did not detect evidence that hiking path density1000 or Strava index1000 significantly influenced habitat selection at the home range-scale (Fig. 4, P value 0.502 and 0.483, respectively). Lynx selection increased with steeper slope (slope1000, P value < 2*10–16), higher proportion of fields (P value 2.56*10–4), higher forest road density (P value 3.74*10–4), lower house density (P value 3.96*10–4) and they selected forest cover with an optimum around 77% forest (see Fig. 6A, P value 0.458 and 4.86*10–6 for first and second order term, respectively). We did not find evidence that public road density significantly influenced habitat selection by lynx at the home range-scale (P value 0.153).
Local-scale habitat selection
The full model including interaction with time of day emerged as the top model for local-scale habitat selection (Table S4) and was clearly the most supported model. Models including “night” and/or the Strava index performed better than those without. Lynx avoidance increased with higher Strava index50 during both day (P value 9.51*10–4) and night (P value 0.017, Fig. 5). The effect of hiking path density50 was not significant (P value 0.079 and 0.920 for day and night, respectively), regardless of the time of day, although there was a trend towards avoidance during the day.
Lynx selection decreased with distance to fields (Fig. 6B P value 0.001). We did not detect an effect of distance to house on local-scale habitat selection (P value 0.063). Lynx selection increased with steeper slope (P value day < 2*10–16 and night < 2*10–16) and in forest (P value day: 1.63*10–6 and night: 3.61*10–4), but this effect was weaker during the night. During the day, lynx selection increased with distance from public roads (avoidance, P value 3.95*10–4) but showed no response to forest roads (P value 0.195), whereas at night, lynx selection decreased with distance to forest roads (selection for, P value 1.17*10–4) but showed no response to public roads (P value 0.318).
Time dependent habitat use of areas associated with the Strava index50
Despite lynx selection decreasing with the Strava index50 at the local-scale both during day and night, lynx still used areas with a Strava index50 > 0. The GAMM revealed a time-dependent use of these areas (Fig. 7). The proportion of lynx locations having a Strava index50 higher than 0 was lowest during the day, from around 08:00 to 16:00 when predicted proportion was in the range of 0.06 to 0.07. From 16:00 to 00:00 the proportion increased and reached a peak of 0.14 around 01:00, after the peak the proportion declined until 08:00.
Discussion
Our study revealed that lynx exhibit local-scale avoidance of areas with high levels of recreation in summer. Interestingly, while this effect was pronounced at the local-scale habitat selection, we did not detect any effect of recreation on home range-scale habitat selection. These results suggest that lynx are capable of adjusting their habitat selection and temporally adjusting their habitat use to recreation in a way that allows them to occupy human-dominated landscapes. This study also illustrates the added value presented by crowdsourced human mobility data as a reliable proxy for human activity in ecological studies, and the importance of accounting for the level of recreation associated with linear features of all types.
Accounting for the intensity of recreational use is important when studying the impact of recreation on wildlife, as the level of recreation can affect the animals’ responses17. The high level of recreational activity observed on other linear features than hiking paths suggests that using hiking paths as proxies for summer recreation fails to cover the entire spectrum of recreation. The inclusion of the crowdsourced human mobility data in our study revealed new details about the habitat selection and habitat use of lynx, similar to a recent study on brown bears (Ursus arctos) in Italy33. Lynx appeared to be relatively tolerant towards recreation, as we detected comparatively high use of the areas associated with the Strava index50 (11% of the used locations had values higher than 0). Due to the lack of avoidance at the home range-scale, lynx do not seem reluctant to occupy the same areas that humans use for recreation, but they do avoid the immediate surroundings of linear features associated with high levels of recreation (local-scale avoidance).
Previous studies have investigated the effect of recreation on home range-scale habitat selection and habitat use of large carnivores. These studies have reported selection for areas with nonmotorized winter recreation by Canada lynx (Lynx canadensis)61,62, avoidance of areas with higher intensity of winter recreation (both motorized and non-motorized) by wolverines16, avoidance of areas with higher recreational intensity by brown bears33 and daybed selection for areas assumed to receive less recreational activity by lynx63. In this study, we did not detect any spatial avoidance of hiking path density nor Strava index1000 (in a 1 km buffer) by lynx at the home range-scale habitat selection. Lynx have been reported to have relatively short flight initiation distances in forests, with a median distance at 50 m64. Hence, a substantial reduction in the need to initiate a flight response is likely achieved by local-scale avoidance of areas with high recreation levels at local-scale habitat selection. This local-scale avoidance might mitigate the need to exhibit larger scale avoidance; instead of avoiding a large forest area associated with high levels of recreation, lynx can still use it and can reduce the risk of encountering humans by local-scale avoidance of recreationists and temporal adjustments of habitat use.
During the day, local-scale avoidance of recreation likely reflects a selection of resting sites away from areas with high recreational use, which has also been reported for lynx in southern-Europe63. We found that lynx avoidance of areas with higher Strava index50 persisted throughout the night, at times when lynx are most active65 and humans are not. A lack of temporal adjustment in habitat selection towards areas associated with recreation (non-motorized and motorized recreation during winter) has also been reported for Canada lynx61. However, our results show that lynx, despite locally avoiding areas associated with higher levels of recreation, used areas associated with recreation quite often. Especially at night, when habitat use of Strava index50 was twice as high than during the day, showing some temporal adaptations towards recreation. Higher habitat use of areas associated with the Strava index50 could be explained by cost-effective transportation, as linear features have been shown to facilitate movement for other carnivores15,66.
Lynx in our study area occupy a human-dominated landscape and are thus capable of adjusting to human infrastructure67 and, as our results suggest, also recreation. The level of recreation in our study area might not be high enough to force lynx to adjust their habitat selection at larger scales. Nonetheless, our study area contains one of the most heavily used recreational areas in the immediate proximity of the capital of Norway (Oslo), and lynx still used this area. Additionally, the Scandinavian lynx population has been, and still is, subject to strong selection pressure to avoid humans due to hunting and poaching35,36. In this context, a lack of avoidance at the home range-scale habitat selection suggests that an area needs to receive substantially high levels of recreation before lynx start to avoid it at large scale, and that the other spatial and behavioural adaptations are sufficient. Our study area is forested, with abundant hiding cover and widespread access to rugged terrain and/or boulders. Dense horizontal cover has been shown to reduce the flight initiation distance for lynx64 and other large carnivores12. Hence, the effect of recreation might be less pronounced in forested landscapes with access to cover.
Strava data have proven useful in previous studies, and high correlations with ground truth data have been reported from cities in Norway68, the UK69, the USA70 and Australia71 as well as in rural areas in Austria34 and in Italy33. As our Strava index is an index of pedestrian recreation and the app is not used by everyone engaging in recreational activities in a defined area, true recreational activity is bound to be higher. For example, Venter, et al.68 found the ratio between the Strava data and human counters to range between 1:30 and 1:40 in Oslo, Norway, meaning that for each person using the Strava app there are an additional 30 to 40 people on the same track or road during the same time period. However, as long as Strava users are not using different areas than non-Strava users, this proxy for human activity should reliably represent relative recreational activity in our study area. We believe spatially crowdsourced data on human mobility or activity can open a range of new possibilities for wildlife research and inform management in the future. In cases like ours, with a study area of approximately 43,000 km2, crowdsourced data is currently the most feasible, and maybe the only, option for deriving a proxy for human activity across the whole area. The temporal mismatch (in years) between the collection of the lynx data and the Strava data have likely not impacted our results. For this mismatch to have an impact on our results, it has to be a drastic change in recreation patterns between the two periods. We argue that this is unlikely because there was a high correlation among the different years in the yearly Strava data from the period it was collected (Figure S1). However, there are other limitations of the data, e.g. the privacy legislations (minimum 3 unique users on a segment to be reported) impede the use of high resolution temporal Strava data (e.g. activity counts per hour or day) in areas receiving few activity events per temporal unit.
Conclusions
This study gives an example of the added-value from the application of crowdsourced human mobility data for ecological studies. Our results suggest that lynx reduce their direct interaction with pedestrian recreationists through local-scale avoidance and temporal adjustments in habitat use. The consequences of recreation for lynx in Norway are therefore likely minor as the impact of recreation appears to be spatially restricted to the immediate surroundings of linear features in which the recreation occurs. We believe the levels of recreation in our study area are not high enough to impede lynx from sharing the landscape with humans. Instead, spatial avoidance at local-scales and temporal adjustments in habitat use may facilitate coexistence between humans and large carnivores.
Data availability
Data used in this analysis are available from the corresponding author upon request.
References
United Nations. (Department of Economic and Social Affairs, Population Division, 2019).
Tucker, M. A. et al. Moving in the Anthropocene: Global reductions in terrestrial mammalian movements. Science 359, 466–469. https://doi.org/10.1126/science.aam9712 (2018).
Tablado, Z. & Jenni, L. Determinants of uncertainty in wildlife responses to human disturbance. Biol. Rev. 92, 216–233. https://doi.org/10.1111/brv.12224 (2017).
IUCN. IUCN Programme 2017–2020. (2016).
IUCN. The IUCN Red List of Threatened Species. Version 2021–3. (2021).
Balmford, A. et al. Walk on the wild side: estimating the global magnitude of visits to protected areas. PLoS Biol. 13, 6. https://doi.org/10.1371/journal.pbio.1002074 (2015).
Balmford, A. et al. A global perspective on trends in nature-based tourism. PLoS Biol. 7, 6. https://doi.org/10.1371/journal.pbio.1000144 (2009).
Seto, K. C., Guneralp, B. & Hutyra, L. R. Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools. Proc. Natl. Acad. Sci. U. S. A. 109, 16083–16088. https://doi.org/10.1073/pnas.1211658109 (2012).
Chen, G. Z. et al. Global projections of future urban land expansion under shared socioeconomic pathways. Nat. Commun. 11, 12. https://doi.org/10.1038/s41467-020-14386-x (2020).
Larson, C. L., Reed, S. E., Merenlender, A. M. & Crooks, K. R. Effects of recreation on animals revealed as widespread through a global systematic review. PLoS ONE 11, 21. https://doi.org/10.1371/journal.pone.0167259 (2016).
Frid, A. & Dill, L. Human-caused disturbance stimuli as a form of predation risk. Conserv. Ecol. 6, 16 (2002).
Moen, G. K., Stoen, O. G., Sahlen, V. & Swenson, J. E. Behaviour of solitary adult scandinavian brown bears (Ursus arctos) when approached by humans on foot. PLoS ONE 7, 7. https://doi.org/10.1371/journal.pone.0031699 (2012).
Le Grand, L. et al. Behavioral and physiological responses of scandinavian brown bears (ursus arctos) to dog hunts and human encounters. Front. Ecol. Evol. 7, 9. https://doi.org/10.3389/fevo.2019.00134 (2019).
Johnson, D. H. The comparison of usage and availability measurements for evaluating resource preference. Ecol. (Washington D C) 61, 65–71. https://doi.org/10.2307/1937156 (1980).
Zimmermann, B., Nelson, L., Wabakken, P., Sand, H. & Liberg, O. Behavioral responses of wolves to roads: scale-dependent ambivalence. Behav. Ecol. 25, 1353–1364. https://doi.org/10.1093/beheco/aru134 (2014).
Heinemeyer, K. et al. Wolverines in winter: indirect habitat loss and functional responses to backcountry recreation. Ecosphere https://doi.org/10.1002/ecs2.2611 (2019).
Ladle, A. et al. Grizzly bear response to spatio-temporal variability in human recreational activity. J. Appl. Ecol. 56, 375–386. https://doi.org/10.1111/1365-2664.13277 (2019).
Coppes, J., Burghardt, F., Hagen, R., Suchant, R. & Braunisch, V. Human recreation affects spatio-temporal habitat use patterns in red deer (Cervus elaphus). PLoS ONE 12, 19. https://doi.org/10.1371/journal.pone.0175134 (2017).
Kautz, T. M. et al. Large carnivore response to human road use suggests a landscape of coexistence. Global Ecol. Conserv. 30, e01772 (2021).
Chapron, G. et al. Recovery of large carnivores in Europe’s modern human-dominated landscapes. Science 346, 1517–1519. https://doi.org/10.1126/science.1257553 (2014).
Ordiz, A., Bischof, R. & Swenson, J. E. Saving large carnivores, but losing the apex predator?. Biol. Conserv. 168, 128–133. https://doi.org/10.1016/j.biocon.2013.09.024 (2013).
Estes, J. A. et al. Trophic downgrading of planet earth. Science 333, 301–306. https://doi.org/10.1126/science.1205106 (2011).
Ordiz, A. et al. Habituation, sensitization, or consistent behavioral responses? Brown bear responses after repeated approaches by humans on foot. Biol. Conserv. 232, 228–237. https://doi.org/10.1016/j.biocon.2019.01.016 (2019).
Smith, T. S., Oyster, J., Partridge, S. D., Martin, I. E. & Sisson, A. Assessing American black bear response to human activity at Kenai Fjords National Park, Alaska. Ursus 23, 179–191. https://doi.org/10.2192/ursus-d-11-00020.1 (2012).
Wam, H. K., Eldegard, K. & Hjeljord, O. Minor habituation to repeated experimental approaches in Scandinavian wolves. Eur. J. Wildl. Res. 60, 839–842. https://doi.org/10.1007/s10344-014-0841-0 (2014).
Sweanor, L. L., Logan, K. A. & Hornocker, M. G. Puma responses to close approaches by researchers. Wildl. Soc. Bull. 33, 905–913. https://doi.org/10.2193/0091-7648(2005)33[905:Prtcab]2.0.Co;2 (2005).
Coppes, J., Ehrlacher, J., Thiel, D., Suchant, R. & Braunisch, V. Outdoor recreation causes effective habitat reduction in capercaillie Tetrao urogallus: a major threat for geographically restricted populations. J. Avian Biol. 48, 1583–1594. https://doi.org/10.1111/jav.01239 (2017).
Coppes, J. et al. Habitat suitability modulates the response of wildlife to human recreation. Biol. Conserv. 227, 56–64. https://doi.org/10.1016/j.biocon.2018.08.018 (2018).
Gundersen, V., Vistad, O. I., Panzacchi, M., Strand, O. & van Moorter, B. Large-scale segregation of tourists and wild reindeer in three Norwegian national parks: Management implications. Tourism Manage. 75, 22–33. https://doi.org/10.1016/j.tourman.2019.04.017 (2019).
Filla, M. et al. Habitat selection by Eurasian lynx (Lynx lynx) is primarily driven by avoidance of human activity during day and prey availability during night. Ecol. Evol. 7, 6367–6381. https://doi.org/10.1002/ece3.3204 (2017).
Andersen, O., Gundersen, V., Wold, L. C. & Stange, E. Monitoring visitors to natural areas in wintertime: issues in counter accuracy. J. Sustain. Tour. 22, 550–560. https://doi.org/10.1080/09669582.2013.839693 (2014).
Marion, S. et al. A systematic review of methods for studying the impacts of outdoor recreation on terrestrial wildlife. Glob. Ecol. Conserv. 22, e00917 (2020).
Corradini, A. et al. Effects of cumulated outdoor activity on wildlife habitat use. Biol. Conserv. 253, 8. https://doi.org/10.1016/j.biocon.2020.108818 (2021).
Jager, H., Schirpke, U. & Tappeiner, U. Assessing conflicts between winter recreational activities and grouse species. J. Environ. Manage. 276, 9. https://doi.org/10.1016/j.jenvman.2020.111194 (2020).
Linnell, J. D. C., Broseth, H., Odden, J. & Nilsen, E. B. Sustainably harvesting a large Carnivore? Development of Eurasian Lynx populations in Norway during 160 years of shifting policy. Environ. Manage. 45, 1142–1154. https://doi.org/10.1007/s00267-010-9455-9 (2010).
Andren, H. et al. Survival rates and causes of mortality in Eurasian lynx (Lynx lynx) in multi-use landscapes. Biol. Conserv. 131, 23–32. https://doi.org/10.1016/j.biocon.2006.01.025 (2006).
Manly, B., McDonald, L., Thomas, D., McDonald, T. & Erickson, W. Resource Selection by Animals (Dordrecht: Kluwer Academic Publishers, 2002).
Odden, J., Linnell, J. D. C. & Andersen, R. Diet of Eurasian lynx, Lynx lynx, in the boreal forest of southeastern Norway: the relative importance of livestock and hares at low roe deer density. Eur. J. Wildl. Res. 52, 237–244. https://doi.org/10.1007/s10344-006-0052-4 (2006).
Gervasi, V., Nilsen, E. B., Odden, J., Bouyer, Y. & Linnell, J. D. C. The spatio-temporal distribution of wild and domestic ungulates modulates lynx kill rates in a multi-use landscape. J. Zool. 292, 175–183. https://doi.org/10.1111/jzo.12088 (2014).
Arnemo, J. M. & Evans, A. Biomedical protocols for free-ranging brown bears, wolves, wolverines and lynx (Hedmark University College Evenstad, 2017).
Padgham, M., Lovelace, R., Salmon, M. & Rudis, B. osmdata. J. Open Source Softw. 2, 305 (2017).
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. (2020).
Northrup, J. M., Hooten, M. B., Anderson, C. R. & Wittemyer, G. Practical guidance on characterizing availability in resource selection functions under a use-availability design. Ecology 94, 1456–1463. https://doi.org/10.1890/12-1688.1 (2013).
Horne, J. S., Garton, E. O., Krone, S. M. & Lewis, J. S. Analyzing animal movements using Brownian bridges. Ecology 88, 2354–2363. https://doi.org/10.1890/06-0957.1 (2007).
Calenge, C. The package “adehabitat” for the R software: a tool for the analysis of space and habitat use by animals. Ecol. Model. 197, 516–519. https://doi.org/10.1016/j.ecolmodel.2006.03.017 (2006).
Therneau, T. A Package for Survival Analysis in R. R package version 3.2–7. (2020).
Fay, M. P., Graubard, B. I., Freedman, L. S. & Midthune, D. N. Conditional logistic regression with sandwich estimators: Application to a meta-analysis. Biometrics 54, 195–208. https://doi.org/10.2307/2534007 (1998).
Prima, M.-C., Duchesne, T. & Fortin, D. Robust inference from conditional logistic regression applied to movement and habitat selection analysis. PLoS ONE 12, e0169779 (2017).
Basille, M. et al. Selecting habitat to survive: the impact of road density on survival in a large carnivore. PLoS ONE 8, 11. https://doi.org/10.1371/journal.pone.0065493 (2013).
Basille, M. et al. What shapes Eurasian lynx distribution in human dominated landscapes: selecting prey or avoiding people?. Ecography 32, 683–691. https://doi.org/10.1111/j.1600-0587.2009.05712.x (2009).
Bouyer, Y. et al. Eurasian lynx habitat selection in human-modified landscape in Norway: effects of different human habitat modifications and behavioral states. Biol. Conserv. 191, 291–299. https://doi.org/10.1016/j.biocon.2015.07.007 (2015).
Heggem, E. S. F., Mathisen, H. & Frydenlund, J. J. N. R. AR50–Arealressurskart i målestokk 1: 50 000. Et heldekkende arealressurskart for jord-og skogbruk. (2019).
Hijmans, R. J. raster: Geographic Data Analysis and Modeling. R package version 3.4–5. (2020).
Bivand, R. & Lewin-Koh, N. maptools: Tools for Handling Spatial Objects. R package version 1.1–1. https://CRAN.R-project.org/package=maptools. (2021).
Akaike, H. in IEEE Transactions on Automatic Control Vol. 19 716–723 (1974).
Burnham, K. P. & Anderson, D. R. Model selection and multimodel inference—a practical information-theoretic approach.2nd edn. Springer, New York (2002).
Hastie, T. J. & Tibshirani, R. J. Generalized additive models. Vol. 43 (CRC press, 1990).
Wood, S. N. Stable and efficient multiple smoothing parameter estimation for generalized additive models. J. Am. Stat. Assoc. 99, 673–686. https://doi.org/10.1198/016214504000000980 (2004).
Vazquez, C., Rowcliffe, J. M., Spoelstra, K. & Jansen, P. A. Comparing diel activity patterns of wildlife across latitudes and seasons: Time transformations using day length. Methods Ecol. Evol. 10, 2057–2066. https://doi.org/10.1111/2041-210x.13290 (2019).
Rowcliffe, M. activity: Animal Activity Statistics. R package version 1.3.1. (2021).
Olson, L. E., Squires, J. R., Roberts, E. K., Ivan, J. S. & Hebblewhite, M. Sharing the same slope: behavioral responses of a threatened mesocarnivore to motorized and nonmotorized winter recreation. Ecol. Evol. 8, 8555–8572. https://doi.org/10.1002/ece3.4382 (2018).
Squires, J. R., Olson, L. E., Roberts, E. K., Ivan, J. S. & Hebblewhite, M. Winter recreation and Canada lynx: reducing conflict through niche partitioning. Ecosphere 10, 22. https://doi.org/10.1002/ecs2.2876 (2019).
Belotti, E., Mayer, K., Kreisinger, J., Heurich, M. & Bufka, L. Recreational activities affect resting site selection and foraging time of Eurasian lynx (Lynx lynx). Hystrix 29, 181–189. https://doi.org/10.4404/hystrix-00053-2018 (2018).
Sunde, P., Stener, S. O. & Kvam, T. Tolerance to humans of resting lynxes Lynx lynx in a hunted population. Wildlife Biol. 4, 177–183 (1998).
Heurich, M. et al. Activity patterns of Eurasian Lynx are modulated by light regime and individual traits over a wide latitudinal range. PLoS ONE 9, 20. https://doi.org/10.1371/journal.pone.0114143 (2014).
Bischof, R., Gjevestad, J. G. O., Ordiz, A., Eldegard, K. & Milleret, C. High frequency GPS bursts and path-level analysis reveal linear feature tracking by red foxes. Sci. Rep. 9, 13. https://doi.org/10.1038/s41598-019-45150-x (2019).
Bouyer, Y. et al. Tolerance to anthropogenic disturbance by a large carnivore: the case of Eurasian lynx in south-eastern Norway. Anim. Conserv. 18, 271–278. https://doi.org/10.1111/acv.12168 (2015).
Venter, Z. S., Barton, D. N., Gundersen, V., Figari, H. & Nowell, M. Urban nature in a time of crisis: recreational use of green space increases during the COVID-19 outbreak in Oslo, Norway. Environ. Res. Lett. 15, 11. https://doi.org/10.1088/1748-9326/abb396 (2020).
Sun, Y. R., Du, Y. Y., Wang, Y. & Zhuang, L. Y. Examining associations of environmental characteristics with recreational cycling behaviour by street-level strava data. Int. J. Environ. Res. Public Health 14, 12. https://doi.org/10.3390/ijerph14060644 (2017).
Griffin, G. P. & Jiao, J. Where does bicycling for health happen? analysing volunteered geographic information through place and plexus. J. Transp. Health 2, 238–247. https://doi.org/10.1016/j.jth.2014.12.001 (2015).
Conrow, L., Wentz, E., Nelson, T. & Pettit, C. Comparing spatial patterns of crowdsourced and conventional bicycling datasets. Appl. Geogr. 92, 21–30. https://doi.org/10.1016/j.apgeog.2018.01.009 (2018).
Acknowledgements
The study was conducted within the frames of the Scandinavian Lynx Project, Scandlynx (http://scandlynx.nina.no/), and was supported by the Norwegian Environment Agency, the Research Council of Norway (grants 251112 and 281092), and the Nature Protection Division of the County Governor’s Office for Innlandet and Viken Counties. We would like to thank the hundreds of people who have helped us in the field during the years of this study. This report includes aggregated and de-identified data from Strava Metro.
Author information
Authors and Affiliations
Contributions
N.H.T., J.O. and J.L. conceived the idea. All authors contributed to the design of the study. J.O. collected the data. N.H.T. did the analysis and consulted with R.B., T.H. and J.M. J.O. and J.L. obtained the funding. N.H.T. wrote the manuscript. All authors discussed the results and contributed significantly to the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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/.
About this article
Cite this article
Thorsen, N.H., Bischof, R., Mattisson, J. et al. Smartphone app reveals that lynx avoid human recreationists on local scale, but not home range scale. Sci Rep 12, 4787 (2022). https://doi.org/10.1038/s41598-022-08468-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-08468-7
This article is cited by
-
Prerequisites for coexistence: human pressure and refuge habitat availability shape continental-scale habitat use patterns of a large carnivore
Landscape Ecology (2023)
-
Human recreation impacts seasonal activity and occupancy of American black bears (Ursus americanus) across the anthropogenic-wildland interface
Scientific Reports (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.