Bottom-up rather than top-down mechanisms determine mesocarnivore interactions in Norway.

Interactions among coexisting mesocarnivores can be influenced by different factors such as the presence of large carnivores, land-use, environmental productivity, or human disturbance. Disentangling the relative importance of bottom-up and top-down processes can be challenging, but it is important for biodiversity conservation and wildlife management. The aim of this study was to assess how the interactions among mesocarnivores (red fox Vulpes vulpes, badger Meles meles, and pine marten Martes martes) were affected by large carnivores (Eurasian lynx Lynx lynx and wolf Canis lupus), land cover variables (proportion of agricultural land and primary productivity), and human disturbance, as well as how these top-down and bottom-up mechanisms were influenced by season. We analyzed 3 years (2018-2020) of camera trapping observations from Norway and used structural equation models to assess hypothesized networks of causal relationships. Our results showed that land cover variables were more strongly associated with mesocarnivore detection rates than large carnivores in Norway. This might be caused by a combination of low density of large carnivores in an unproductive ecosystem with strong seasonality. Additionally, detection rates of all mesocarnivores showed positive associations among each other, which were stronger in winter. The prevalence of positive interactions among predators might indicate a tendency to use the same areas and resources combined with weak interference competition. Alternatively, it might indicate some kind of facilitative relationship among species. Human disturbance had contrasting effects for different species, benefiting the larger mesocarnivores (red fox and badger) probably through food subsidization, but negatively affecting apex predators (wolf and lynx) and smaller mesocarnivores (pine marten). In a human-dominated world, this highlights the importance of including anthropogenic influences in the study of species interactions.

Altitude, latitude and climate zone as determinants of mountain hare (Lepus timidus) coat colour change.

Local adaptation to annually changing environments has evolved in numerous species. Seasonal coat colour change is an adaptation that has evolved in multiple mammal and bird species occupying areas that experience seasonal snow cover. It has a critical impact on fitness as predation risk may increase when an individual is mismatched against its habitat's background colour. In this paper, we investigate the correlation between landscape covariates and moult timing in a native winter-adapted herbivore, the mountain hare (Lepus timidus), throughout Norway. Data was collected between 2011 and 2019 at 678 camera trap locations deployed across an environmental gradient. Based on this data, we created a Bayesian multinomial logistic regression model that quantified the correlations between landscape covariates and coat colour phenology and analysed among season and year moult timing variation. Our results demonstrate that mountain hare moult timing is strongly correlated with altitude and latitude with hares that live at higher latitudes and altitudes keeping their winter white coats for longer than their conspecifics that inhabit lower latitudes and altitudes. Moult timing was also weakly correlated with climate zone with hares that live in coastal climates keeping their winter white coats for longer than hares that live in continental climates. We found evidence of some among year moult timing variation in spring, but not in autumn. We conclude that mountain hare moult timing has adapted to local environmental conditions throughout Norway.

Predicting kill sites of an apex predator from GPS data in different multiprey systems.

Kill rates are a central parameter to assess the impact of predation on prey species. An accurate estimation of kill rates requires a correct identification of kill sites, often achieved by field-checking GPS location clusters (GLCs). However, there are potential sources of error included in kill-site identification, such as failing to detect GLCs that are kill sites, and misclassifying the generated GLCs (e.g., kill for nonkill) that were not field checked. Here, we address these two sources of error using a large GPS dataset of collared Eurasian lynx (Lynx lynx), an apex predator of conservation concern in Europe, in three multiprey systems, with different combinations of wild, semidomestic, and domestic prey. We first used a subsampling approach to investigate how different GPS-fix schedules affected the detection of GLC-indicated kill sites. Then, we evaluated the potential of the random forest algorithm to classify GLCs as nonkills, small prey kills, and ungulate kills. We show that the number of fixes can be reduced from seven to three fixes per night without missing more than 5% of the ungulate kills, in a system composed of wild prey. Reducing the number of fixes per 24 h decreased the probability of detecting GLCs connected with kill sites, particularly those of semidomestic or domestic prey, and small prey. Random forest successfully predicted between 73%-90% of ungulate kills, but failed to classify most small prey in all systems, with sensitivity (true positive rate) lower than 65%. Additionally, removing domestic prey improved the algorithm's overall accuracy. We provide a set of recommendations for studies focusing on kill-site detection that can be considered for other large carnivore species in addition to the Eurasian lynx. We recommend caution when working in systems including domestic prey, as the odds of underestimating kill rates are higher.

Quantifying the checks and balances of collaborative governance systems for adaptive carnivore management.

Recovering or threatened carnivore populations are often harvested to minimise their impact on human activities, such as livestock farming or game hunting. Increasingly, harvest quota decisions involve a set of scientific, administrative and political institutions operating at national and sub-national levels whose interactions and collective decision-making aim to increase the legitimacy of management and ensure population targets are met. In practice, however, assessments of how quota decisions change between these different actors and what consequences these changes have on population trends are rare.We combine a state-space population modelling approach with an analysis of quota decisions taken at both regional and national levels between 2007 and 2018 to build a set of decision-making models that together predict annual harvest quota values for Eurasian lynx (Lynx lynx) in Norway.We reveal a tendency for administrative decision-makers to compensate for consistent quota increases by political actors, particularly when the lynx population size estimate is above the regional target. Using population forecasts based on the ensemble of decision-making models, we show that such buffering of political biases ensures lynx population size remains close to regional and national targets in the long term.Our results go beyond the usual qualitative assessment of collaborative governance systems for carnivore management, revealing a system of checks and balances that, in the case of lynx in Norway, ensures both multi-stakeholder participation and sustainable harvest quotas. Nevertheless, we highlight important inter-regional differences in decision-making and population forecasts, the socio-ecological drivers of which need to be better understood to prevent future population declines. Synthesis and applications. Our work analyses the sequence of decisions leading to yearly quotas for lynx harvest in Norway, highlighting the collaborative and structural processes that together shape harvest sustainability. In doing so, we provide a predictive framework to evaluate participatory decision-making processes in wildlife management, paving the way for scientists and decision-makers to collaborate more widely in identifying where decision biases might lie and how institutional arrangements can be optimised to minimise them. We emphasise, however, that this is only possible if wildlife management decisions are documented and transparent.

Timing and synchrony of birth in Eurasian lynx across Europe.

The ecology and evolution of reproductive timing and synchrony have been a topic of great interest in evolutionary ecology for decades. Originally motivated by questions related to behavioral and reproductive adaptation to environmental conditions, the topic has acquired new relevance in the face of climate change. However, there has been relatively little research on reproductive phenology in mammalian carnivores. The Eurasian lynx (Lynx lynx) occurs across the Eurasian continent, covering three of the four main climate regions of the world. Thus, their distribution includes a large variation in climatic conditions, making it an ideal species to explore reproductive phenology. Here, we used data on multiple reproductive events from 169 lynx females across Europe. Mean birth date was May 28 (April 23 to July 1), but was ~10 days later in northern Europe than in central and southern Europe. Birth dates were relatively synchronized across Europe, but more so in the north than in the south. Timing of birth was delayed by colder May temperatures. Severe and cold weather may affect neonatal survival via hypothermia and avoiding inclement weather early in the season may select against early births, especially at northern latitudes. Overall, only about half of the kittens born survived until onset of winter but whether kittens were born relatively late or early did not affect kitten survival. Lynx are strict seasonal breeders but still show a degree of flexibility to adapt the timing of birth to surrounding environmental conditions. We argue that lynx give birth later when exposed to colder spring temperatures and have more synchronized births when the window of favorable conditions for raising kittens is shorter. This suggests that lynx are well adapted to different environmental conditions, from dry and warm climates to alpine, boreal, and arctic climates. This variation in reproductive timing will be favorable in times of climate change, as organisms with high plasticity are more likely to adjust to new environmental conditions.

Smartphone app reveals that lynx avoid human recreationists on local scale, but not home range scale.

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.

Camera trap records of leucistic Eurasian badgers (Meles meles) in central Norway.

Coat coloration plays an important role in communication, camouflage, and sexual selection in animals. Genetic mutations can lead to anomalous colorations such as melanism and leucism, where animals appear, respectively, darker or lighter than normal. Reporting abnormal coloration in wild animals is an important first step to understand the distribution, prevalence, and potential fitness consequences of these rare events. Here, we report several records of suspected leucism in the Eurasian badger (Meles meles) in a population in central Norway. Several camera traps recorded at least two leucistic individuals between 2017 and 2020. It took considerable effort, almost 400,000 camera trap nights over a period of 10 years all over Norway, to obtain a total of eleven records of leucistic badgers, indicating the rarity of this phenotype. It is unclear what has caused the presence of multiple leucistic badgers in a single population, but recent colonization and lack of predators might have played a role. Due to our observations, future studies can now be developed to study the underlying mechanisms and potential consequences of leucism in this badger population. The increasing use of networks of camera traps in wildlife research will provide new opportunities to record rare coloration in wild animals.

GPS collars have an apparent positive effect on the survival of a large carnivore.

Are instrumented animals representative of the population, given the potential bias caused by selective sampling and the influence of capture, handling and wearing bio-loggers? The answer is elusive owing to the challenges of obtaining comparable data from individuals with and without bio-loggers. Using non-invasive genetic data of a large carnivore, the wolverine (Gulo gulo) in Scandinavia, and an open-population spatial capture-recapture model, we found a 16 (credible interval: 4-30) percentage points lower mortality probability for GPS-collared individuals compared with individuals without GPS collars. While the risk of dying from legal culling was comparable for collared and non-collared wolverines, the former experienced lower probability of mortality due to causes other than legal culling. The aforementioned effect was pronounced despite a potentially lower age-and therefore likely higher natural mortality-of collared individuals. Reports of positive effects of bio-loggers on the survival of individuals are uncommon and we argue that GPS collars could shield animals from poaching. Our results highlight the challenges of drawing population-level inferences for populations subjected to poaching when using data from instrumented individuals.

Extreme home range sizes among Eurasian lynx at the northern edge of their biogeographic range.

Eurasian lynx (Lynx lynx) have a wide distribution across Eurasia. The northern edge of this distribution is in Norway, where they reach up to 72 degrees north. We conducted a study of lynx space use in this region from 2007 to 2013 using GPS telemetry. The home range sizes averaged 2,606 (± 438 SE) km2 for males (n = 9 ranges) and 1,456 (± 179 SE) km2 for females (n = 24 ranges). These are the largest home ranges reported for any large felid, and indeed are only matched by polar bears, arctic living wolves, and grizzly bears among all the Carnivora. The habitat occupied was almost entirely treeless alpine tundra, with home ranges only containing from 20% to 25% of forest. These data have clear implications for the spatial planning of lynx management in the far north as the current management zones are located in unsuitable habitats and are not large enough to encompass individual lynx movements.

Harvest models of small populations of a large carnivore using Bayesian forecasting.

Harvesting large carnivores can be a management tool for meeting politically set goals for their desired abundance. However, harvesting carnivores creates its own set of conflicts in both society and among conservation professionals, where one consequence is a need to demonstrate that management is sustainable, evidence-based, and guided by science. Furthermore, because large carnivores often also have high degrees of legal protection, harvest quotas have to be carefully justified and constantly adjusted to avoid damaging their conservation status. We developed a Bayesian state-space model to support adaptive management of Eurasian lynx harvesting in Scandinavia. The model uses data from the annual monitoring of lynx abundance and results from long-term field research on lynx biology, which has provided detailed estimates of key demographic parameters. We used the model to predict the probability that the forecasted population size will be below or above the management objectives when subjected to different harvest quotas. The model presented here informs decision makers about the policy risks of alternative harvest levels. Earlier versions of the model have been available for wildlife managers in both Sweden and Norway to guide lynx harvest quotas and the model predictions showed good agreement with observations. We combined monitoring data with data on vital rates and were able to estimate unobserved additional mortality rates, which are most probably due to poaching. In both countries, the past quota setting strategy suggests that there has been a de facto threshold strategy with increasing proportion, which means that there is no harvest below a certain population size, but above this threshold there is an increasing proportion of the population harvested as the population size increases. The annual assessment of the monitoring results, the use of forecasting models, and a threshold harvest approach to quota setting will all reduce the risk of lynx population sizes moving outside the desired goals. The approach we illustrate could be adapted to other populations of mammals worldwide.

Eurasian lynx fitness shows little variation across Scandinavian human-dominated landscapes.

Despite extensive research on the ecology and behavioural adaptations of large carnivores in human-dominated landscapes, information about the fitness consequences of sharing landscapes is still limited. We assessed the variation in three consecutive components of female fitness: the probability of reproduction, litter size and juvenile survival in relation to environmental and human factors in a solitary carnivore, the Eurasian lynx (Lynx lynx), occurring in human-dominated landscapes in Scandinavia. We used demographic data from 57 radio-collared adult females between 1995-2011 (126 radio-years). Overall, the yearly probability of female reproduction was 0.80, mean litter size was 2.34 (range 1-4) and the probability to find a female that reproduced in the spring being accompanied by at least one offspring during the subsequent winter was 0.70. We did not find evidence that food availability was a key factor influencing female fitness. Female lynx may adapt to food availability when establishing their home ranges by adopting an obstinate strategy, ensuring a minimum amount of prey necessary for survival and reproduction even during periods of prey scarcity. In human-dominated landscapes, where sufficient prey are available for lynx, mortality risk may have a larger influence on lynx population dynamics compared to food availability. Our results suggest that lynx population dynamics in human-dominated landscapes may be mainly driven by human impacts on survival.

Framing pictures: A conceptual framework to identify and correct for biases in detection probability of camera traps enabling multi-species comparison.

Obtaining reliable species observations is of great importance in animal ecology and wildlife conservation. An increasing number of studies use camera traps (CTs) to study wildlife communities, and an increasing effort is made to make better use and reuse of the large amounts of data that are produced. It is in these circumstances that it becomes paramount to correct for the species- and study-specific variation in imperfect detection within CTs. We reviewed the literature and used our own experience to compile a list of factors that affect CT detection of animals. We did this within a conceptual framework of six distinct scales separating out the influences of (a) animal characteristics, (b) CT specifications, (c) CT set-up protocols, and (d) environmental variables. We identified 40 factors that can potentially influence the detection of animals by CTs at these six scales. Many of these factors were related to only a few overarching parameters. Most of the animal characteristics scale with body mass and diet type, and most environmental characteristics differ with season or latitude such that remote sensing products like NDVI could be used as a proxy index to capture this variation. Factors that influence detection at the microsite and camera scales are probably the most important in determining CT detection of animals. The type of study and specific research question will determine which factors should be corrected. Corrections can be done by directly adjusting the CT metric of interest or by using covariates in a statistical framework. Our conceptual framework can be used to design better CT studies and help when analyzing CT data. Furthermore, it provides an overview of which factors should be reported in CT studies to make them repeatable, comparable, and their data reusable. This should greatly improve the possibilities for global scale analyses of (reused) CT data.

The range of the mange: Spatiotemporal patterns of sarcoptic mange in red foxes (Vulpes vulpes) as revealed by camera trapping.

Sarcoptic mange is a widely distributed disease that affects numerous mammalian species. We used camera traps to investigate the apparent prevalence and spatiotemporal dynamics of sarcoptic mange in a red fox population in southeastern Norway. We monitored red foxes for five years using 305 camera traps distributed across an 18000 km2 area. A total of 6581 fox events were examined to visually identify mange compatible lesions. We investigated factors associated with the occurrence of mange by using logistic models within a Bayesian framework, whereas the spatiotemporal dynamics of the disease were analysed with space-time scan statistics. The apparent prevalence of the disease fluctuated over the study period with a mean of 3.15% and credible interval [1.25, 6.37], and our best logistic model explaining the presence of red foxes with mange-compatible lesions included time since the beginning of the study and the interaction between distance to settlement and season as explanatory variables. The scan analyses detected several potential clusters of the disease that varied in persistence and size, and the locations in the cluster with the highest probability were closer to human settlements than the other survey locations. Our results indicate that red foxes in an advanced stage of the disease are most likely found closer to human settlements during periods of low wild prey availability (winter). We discuss different potential causes. Furthermore, the disease appears to follow a pattern of small localized outbreaks rather than sporadic isolated events.

Intensity of space use reveals conditional sex-specific effects of prey and conspecific density on home range size.

Home range (HR) size variation is often linked to resource abundance, with sex differences expected to relate to sex-specific fitness consequences. However, studies generally fail to disentangle the effects of the two main drivers of HR size variation, food and conspecific density, and rarely consider how their relative influence change over spatiotemporal scales. We used location data from 77 Eurasian lynx (Lynx lynx) from a 16-year Scandinavian study to examine HR sizes variation relative to prey and conspecific density at different spatiotemporal scales. By varying the isopleth parameter (intensity of use) defining the HR, we show that sex-specific effects were conditional on the spatial scale considered. Males had larger HRs than females in all seasons. Females' total HR size declined as prey and conspecific density increased, whereas males' total HR was only affected by conspecific density. However, as the intensity of use within the HR increased (from 90% to 50% isopleth), the relationship between prey density and area showed opposing patterns for females and males; for females, the prey density effect was reduced, while for males, prey became increasingly important. Thus, prey influenced the size of key regions within male HRs, despite total HR size being independent of prey density. Males reduced their HR size during the mating season, likely to remain close to individual females in estrous. Females reduced their HR size postreproduction probably because of movement constrains imposed by dependent young. Our findings highlight the importance of simultaneously considering resources and intraspecific interactions as HR size determinants. We show that sex-specific demands influence the importance of prey and conspecific density on space use at different spatiotemporal scales. Thus, unless a gradient of space use intensity is examined, factors not related to total HR size might be disregarded despite their importance in determining size of key regions within the HR.

Activity patterns of Eurasian lynx are modulated by light regime and individual traits over a wide latitudinal range.

The activity patterns of most terrestrial animals are regarded as being primarily influenced by light, although other factors, such as sexual cycle and climatic conditions, can modify the underlying patterns. However, most activity studies have been limited to a single study area, which in turn limit the variability of light conditions and other factors. Here we considered a range of variables that might potentially influence the activity of a large carnivore, the Eurasian lynx, in a network of studies conducted with identical methodology in different areas spanning latitudes from 49°7'N in central Europe to 70°00'N in northern Scandinavia. The variables considered both light conditions, ranging from a day with a complete day-night cycle to polar night and polar day, as well as individual traits of the animals. We analysed activity data of 38 individual free-ranging lynx equipped with GPS-collars with acceleration sensors, covering more than 11,000 lynx days. Mixed linear additive models revealed that the lynx activity level was not influenced by the daily daylight duration and the activity pattern was bimodal, even during polar night and polar day. The duration of the active phase of the activity cycle varied with the widening and narrowing of the photoperiod. Activity varied significantly with moonlight. Among adults, males were more active than females, and subadult lynx were more active than adults. In polar regions, the amplitude of the lynx daily activity pattern was low, likely as a result of the polycyclic activity pattern of their main prey, reindeer. At lower latitudes, the basic lynx activity pattern peaked during twilight, corresponding to the crepuscular activity pattern of the main prey, roe deer. Our results indicated that the basic activity of lynx is independent of light conditions, but is modified by both individual traits and the activity pattern of the locally most important prey.

Physiologically persistent Corpora lutea in Eurasian lynx (Lynx lynx) - longitudinal ultrasound and endocrine examinations intra-vitam.

Felids generally follow a poly-estrous reproductive strategy. Eurasian lynx (Lynx lynx) display a different pattern of reproductive cyclicity where physiologically persistent corpora lutea (CLs) induce a mono-estrous condition which results in highly seasonal reproduction. The present study was based around a sono-morphological and endocrine study of captive Eurasian lynx, and a control-study on free-ranging lynx. We verified that CLs persist after pregnancy and pseudo-pregnancy for at least a two-year period. We could show that lynx are able to enter estrus in the following year, while CLs from the previous years persisted in structure and only temporarily reduced their function for the period of estrus onset or birth, which is unique among felids. The almost constant luteal progesterone secretion (average of 5 ng/ml serum) seems to prevent folliculogenesis outside the breeding season and has converted a poly-estrous general felid cycle into a mono-estrous cycle specific for lynx. The hormonal regulation mechanism which causes lynx to have the longest CL lifespan amongst mammals remains unclear. The described non-felid like ovarian physiology appears to be a remarkably non-plastic system. The lynx's reproductive ability to adapt to environmental and anthropogenic changes needs further investigation.

One size fits all: Eurasian lynx females share a common optimal litter size.

Lack proposed that the average clutch size of altricial species should be determined by the average maximum number of young the parents can raise such that all females in a given population should share a common optimal clutch size. Support for this model remains equivocal and recent studies have suggested that intra-population variation in clutch size is adaptive because each female has its own optimal clutch size associated with its intrinsic ability to raise offspring. Although Lack litter size and condition-dependent litter size are presented as two competing models, both are based on the concept of individual optimization. We propose a unified optimal litter size model (called 'adaptive litter size') and identify a set of conditions under which a common vs. a state-dependent optimal litter size should be observed. We test whether females of Eurasian lynx (Lynx lynx) have a common optimal litter size, or whether they adjust their litter size according to their state. We used a detailed individual-based data set collected from contrasting populations of Eurasian lynx in Scandinavia. Observed reproductive patterns in female lynx provide strong support for the existence of a common optimal litter size. Litter size did not vary according to female body mass or reproductive category, or among contrasted populations and years. A litter size of 2 was associated with a higher fitness than both smaller and larger litters, and thus corresponded to the 'adaptive litter size' for female lynx. We suggest that the reproductive pattern of female lynx might correspond to a risk avoidance tactic common to all individuals, which has evolved in response to strong environmental constraints generated by a highly unpredictable food supply during lactation.

Density of wild prey modulates lynx kill rates on free-ranging domestic sheep.

Understanding the factors shaping the dynamics of carnivore-livestock conflicts is vital to facilitate large carnivore conservation in multi-use landscapes. We investigated how the density of their main wild prey, roe deer Capreolus capreolus, modulates individual Eurasian lynx Lynx lynx kill rates on free-ranging domestic sheep Ovis aries across a range of sheep and roe deer densities. Lynx kill rates on free-ranging domestic sheep were collected in south-eastern Norway from 1995 to 2011 along a gradient of different livestock and wild prey densities using VHF and GPS telemetry. We used zero-inflated negative binomial (ZINB) models including lynx sex, sheep density and an index of roe deer density as explanatory variables to model observed kill rates on sheep, and ranked the models based on their AICc values. The model including the effects of lynx sex and sheep density in the zero-inflation model and the effect of lynx sex and roe deer density in the negative binomial part received most support. Irrespective of sheep density and sex, we found the lowest sheep kill rates in areas with high densities of roe deer. As roe deer density decreased, males killed sheep at higher rates, and this pattern held for both high and low sheep densities. Similarly, females killed sheep at higher rates in areas with high densities of sheep and low densities of roe deer. However, when sheep densities were low females rarely killed sheep irrespective of roe deer density. Our quantification of depredation rates can be the first step towards establishing fairer compensation systems based on more accurate and area specific estimation of losses. This study demonstrates how we can use ecological theory to predict where losses of sheep will be greatest, and can be used to identify areas where mitigation measures are most likely to be needed.

Selecting habitat to survive: the impact of road density on survival in a large carnivore.

Habitat selection studies generally assume that animals select habitat and food resources at multiple scales to maximise their fitness. However, animals sometimes prefer habitats of apparently low quality, especially when considering the costs associated with spatially heterogeneous human disturbance. We used spatial variation in human disturbance, and its consequences on lynx survival, a direct fitness component, to test the Hierarchical Habitat Selection hypothesis from a population of Eurasian lynx Lynx lynx in southern Norway. Data from 46 lynx monitored with telemetry indicated that a high proportion of forest strongly reduced the risk of mortality from legal hunting at the home range scale, while increasing road density strongly increased such risk at the finer scale within the home range. We found hierarchical effects of the impact of human disturbance, with a higher road density at a large scale reinforcing its negative impact at a fine scale. Conversely, we demonstrated that lynx shifted their habitat selection to avoid areas with the highest road densities within their home ranges, thus supporting a compensatory mechanism at fine scale enabling lynx to mitigate the impact of large-scale disturbance. Human impact, positively associated with high road accessibility, was thus a stronger driver of lynx space use at a finer scale, with home range characteristics nevertheless constraining habitat selection. Our study demonstrates the truly hierarchical nature of habitat selection, which aims at maximising fitness by selecting against limiting factors at multiple spatial scales, and indicates that scale-specific heterogeneity of the environment is driving individual spatial behaviour, by means of trade-offs across spatial scales.

Patterns of variation in reproductive parameters in Eurasian lynx (Lynx lynx).

Detailed knowledge of the variation in demographic rates is central for our ability to understand the evolution of life history strategies and population dynamics, and to plan for the conservation of endangered species. We studied variation in reproductive output of 61 radio-collared Eurasian lynx females in four Scandinavian study sites spanning a total of 223 lynx-years. Specifically, we examined how the breeding proportion and litter size varied among study areas and age classes (2-year-old vs. >2-year-old females). In general, the breeding proportion varied between age classes and study sites, whereas we did not detect such variation in litter size. The lack of differences in litter sizes among age classes is at odds with most findings in large mammals, and we argue that this is because the level of prenatal investment is relatively low in felids compared to their substantial levels of postnatal care.