Recolonization following past persecution questions the importance of persistent snow cover as a range limiting factor for wolverines.

Globally, climate is changing rapidly, which causes shifts in many species' distributions, stressing the need to understand their response to changing environmental conditions to inform conservation and management. Northern latitudes are expected to experience strongest changes in climate, with milder winters and decreasing snow cover. The wolverine (Gulo gulo) is a circumpolar, threatened carnivore distributed in northern tundra, boreal, and subboreal habitats. Previous studies have suggested that wolverine distribution and reproduction are constrained by a strong association with persistent spring snow cover. We assess this hypothesis by relating spatial distribution of 1589 reproductive events, a fitness-related proxy for female reproduction and survival, to snow cover over two decades. Wolverine distribution has increased and number of reproductive events increased 20 times in areas lacking spring snow cover during our study period, despite low monitoring effort where snow is sparse. Thus, the relationship between reproductive events and persistent spring snow cover weakened during this period. These findings show that wolverine reproductive success and hence distribution are less dependent on spring snow cover than expected. This has important implications for projections of future habitat availability, and thus distribution, of this threatened species. Our study also illustrates how past persecution, or other factors, that have restricted species distribution to remote areas can mask actual effects of environmental parameters, whose importance reveals when populations expand beyond previously restricted ranges. Overwhelming evidence shows that climate change is affecting many species and ecological processes, but forecasting potential consequences on a given species requires longitudinal data to revisit hypotheses and reassess the direction and magnitude of climate effects with new data. This is especially important for conservation-oriented management of species inhabiting dynamic systems where environmental factors and human activities interact, a common scenario for many species in different ecosystems around the globe.

Evaluating expert-based habitat suitability information of terrestrial mammals with GPS-tracking data.

Aim: Macroecological studies that require habitat suitability data for many species often derive this information from expert opinion. However, expert-based information is inherently subjective and thus prone to errors. The increasing availability of GPS tracking data offers opportunities to evaluate and supplement expert-based information with detailed empirical evidence. Here, we compared expert-based habitat suitability information from the International Union for Conservation of Nature (IUCN) with habitat suitability information derived from GPS-tracking data of 1,498 individuals from 49 mammal species. Location: Worldwide. Time period: 1998-2021. Major taxa studied: Forty-nine terrestrial mammal species. Methods: Using GPS data, we estimated two measures of habitat suitability for each individual animal: proportional habitat use (proportion of GPS locations within a habitat type), and selection ratio (habitat use relative to its availability). For each individual we then evaluated whether the GPS-based habitat suitability measures were in agreement with the IUCN data. To that end, we calculated the probability that the ranking of empirical habitat suitability measures was in agreement with IUCN's classification into suitable, marginal and unsuitable habitat types. Results: IUCN habitat suitability data were in accordance with the GPS data (> 95% probability of agreement) for 33 out of 49 species based on proportional habitat use estimates and for 25 out of 49 species based on selection ratios. In addition, 37 and 34 species had a > 50% probability of agreement based on proportional habitat use and selection ratios, respectively. Main conclusions: We show how GPS-tracking data can be used to evaluate IUCN habitat suitability data. Our findings indicate that for the majority of species included in this study, it is appropriate to use IUCN habitat suitability data in macroecological studies. Furthermore, we show that GPS-tracking data can be used to identify and prioritize species and habitat types for re-evaluation of IUCN habitat suitability data.

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.

Northward range expansion of rooting ungulates decreases detritivore and predatory mite abundances in boreal forests.

In the last few decades wild boar populations have expanded northwards, colonizing boreal forests. The soil disturbances caused by wild boar rooting may have an impact on soil organisms that play a key role in organic matter turnover. However, the impact of wild boar colonization on boreal forest ecosystems and soil organisms remains largely unknown. We investigated the effect of natural and simulated rooting on decomposer and predatory soil mites (total, adult and juvenile abundances; and adult-juvenile proportion). Our simulated rooting experiment aimed to disentangle the effects of (i) bioturbation due to soil mixing and (ii) removing organic material (wild boar food resources) on soil mites. Our results showed a decline in the abundance of adult soil mites in response to both natural and artificial rooting, while juvenile abundance and the relative proportion of adults and juveniles were not affected. The expansion of wild boar northwards and into new habitats has negative effects on soil decomposer abundances in boreal forests which may cascade through the soil food web ultimately affecting ecosystem processes. Our study also suggests that a combined use of natural and controlled experimental approaches is the way forward to reveal any subtle interaction between aboveground and belowground organisms and the ecosystem functions they drive.

Environmental and seasonal correlates of capercaillie movement traits in a Swedish wind farm.

Animals continuously interact with their environment through behavioral decisions, rendering the appropriate choice of movement speed and directionality an important phenotypic trait. Anthropogenic activities may alter animal behavior, including movement. A detailed understanding of movement decisions is therefore of great relevance for science and conservation alike. The study of movement decisions in relation to environmental and seasonal cues requires continuous observation of movement behavior, recently made possible by high-resolution telemetry. We studied movement traits of 13 capercaillie (Tetrao urogallus), a mainly ground-moving forest bird species of conservation interest, over two summer seasons in a Swedish windfarm using high-resolution GPS tracking data (5-min sampling interval). We filtered and removed unreliable movement steps using accelerometer data and step characteristics. We explored variation in movement speed and directionality in relation to environmental and seasonal covariates using generalized additive mixed models (GAMMs). We found evidence for clear daily and seasonal variation in speed and directionality of movement that reflected behavioral adjustments to biological and environmental seasonality. Capercaillie moved slower when more turbines were visible and faster close to turbine access roads. Movement speed and directionality were highest on open bogs, lowest on recent clear-cuts (<5 y.o.), and intermediate in all types of forest. Our results provide novel insights into the seasonal and environmental correlates of capercaillie movement patterns and supplement previous behavioral observations on lekking behavior and wind turbine avoidance with a more mechanistic understanding.

Behavioral effects of wolf presence on moose habitat selection: testing the landscape of fear hypothesis in an anthropogenic landscape.

Landscape of fear refers to the spatial variation in prey perception of predation risk, that under certain conditions, may lead to changes in their behavior. Behavioral responses of prey in relation to large carnivore predation risk have mainly been conducted in areas with low anthropogenic impact. We used long-term data on the distribution of moose in different habitat types in a system characterized by intensive management of all three trophic levels (silviculture, harvest of wolves and moose) to study effects on moose habitat selection resulting from the return of an apex predator, the wolf. We assumed that coursing predators such as wolves will cause an increased risk for moose in some habitat types and tested the hypotheses that moose will avoid open or young forest habitats following wolf establishment. After wolf recolonization, moose reduced their use of one type of open habitat (bog) but there was neither change in the use of the other open habitat type (clear-cut), nor in their use of young forest. Wolf establishment did not influence the use of habitat close to dense habitat when being in open habitats. Thus, the effect of wolves varied among habitat types and there was no unidirectional support for a behavioral effect of wolves' establishment on moose habitat use. Human-driven habitat heterogeneity, concentration of moose forage to certain habitat types, and the effects of a multiple predator guild on moose may all contribute to the results found. We conclude that the landscape of fear is likely to have weak ecological effects on moose in this system.

Ungulate management in European national parks: Why a more integrated European policy is needed.

1. Primary objectives of national parks usually include both, the protection of natural processes and species conservation. When these objectives conflict, as occurs because of the cascading effects of large mammals (i.e., ungulates and large carnivores) on lower trophic levels, park managers have to decide upon the appropriate management while considering various local circumstances. 2. To analyse if ungulate management strategies are in accordance with the objectives defined for protected areas, we assessed the current status of ungulate management across European national parks using the naturalness concept and identified the variables that influence the management. 3. We collected data on ungulate management from 209 European national parks in 29 countries by means of a large-scale questionnaire survey. Ungulate management in the parks was compared by creating two naturalness scores. The first score reflects ungulate and large carnivore species compositions, and the second evaluates human intervention on ungulate populations. We then tested whether the two naturalness score categories are influenced by the management objectives, park size, years since establishment, percentage of government-owned land, and human impact on the environment (human influence index) using two generalized additive mixed models. 4. In 67.9% of the national parks, wildlife is regulated by culling (40.2%) or hunting (10.5%) or both (17.2%). Artificial feeding occurred in 81.3% of the national parks and only 28.5% of the national parks had a non-intervention zone covering at least 75% of the area. Furthermore, ungulate management differed greatly among the different countries, likely because of differences in hunting traditions and cultural and political backgrounds. Ungulate management was also influenced by park size, human impact on the landscape, and national park objectives, but after removing these variables from the full model the reduced models only showed a small change in the deviance explained. In areas with higher anthropogenic pressure, wildlife diversity tended to be lower and a higher number of domesticated species tended to be present. Human intervention (culling and artificial feeding) was lower in smaller national parks and when park objectives followed those set by the International Union for the Conservation of Nature (IUCN). 5. Our study shows that many European national parks do not fulfil the aims of protected area management as set by IUCN guidelines. In contrast to the USA and Canada, Europe currently has no common ungulate management policy within national parks. This lack of a common policy together with differences in species composition, hunting traditions, and cultural or political context has led to differences in ungulate management among European countries. To fulfil the aims and objectives of national parks and to develop ungulate management strategies further, we highlight the importance of creating a more integrated European ungulate management policy to meet the aims of national parks.

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.

Games as Tools to Address Conservation Conflicts.

Conservation conflicts represent complex multilayered problems that are challenging to study. We explore the utility of theoretical, experimental, and constructivist approaches to games to help to understand and manage these challenges. We show how these approaches can help to develop theory, understand patterns in conflict, and highlight potentially effective management solutions. The choice of approach should be guided by the research question and by whether the focus is on testing hypotheses, predicting behaviour, or engaging stakeholders. Games provide an exciting opportunity to help to unravel the complexity in conflicts, while researchers need an awareness of the limitations and ethical constraints involved. Given the opportunities, this field will benefit from greater investment and development.

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.

Tracking neighbours promotes the coexistence of large carnivores.

The study of competition and coexistence among similar interacting species has long been considered a cornerstone in evolutionary and community ecology. However, understanding coexistence remains a challenge. Using two similar and sympatric competing large carnivores, Eurasian lynx and wolverines, we tested the hypotheses that tracking among heterospecifics and reactive responses to potential risk decreases the probability of an agonistic encounter when predators access shared food resources, thus facilitating coexistence. Lynx and wolverines actively avoided each other, with the degree of avoidance being greater for simultaneous than time-delayed predator locations. Wolverines reacted to the presence of lynx at relatively short distances (mean: 383 m). In general, lynx stayed longer, and were more stationary, around reindeer carcasses than wolverines. However, when both predators were present at the same time around a carcass, lynx shortened their visits, while wolverine behavior did not change. Our results support the idea that risk avoidance is a reactive, rather than a predictive, process. Since wolverines have adapted to coexist with lynx, exploiting lynx-killed reindeer carcasses while avoiding potential encounters, the combined presence of both predators may reduce wolverine kill rate and thus the total impact of these predators on semi-domestic reindeer in Scandinavia. Consequently, population management directed at lynx may affect wolverine populations and human-wolverine conflicts.

Large impact of Eurasian lynx predation on roe deer population dynamics.

The effects of predation on ungulate populations depend on several factors. One of the most important factors is the proportion of predation that is additive or compensatory respectively to other mortality in the prey, i.e., the relative effect of top-down and bottom-up processes. We estimated Eurasian lynx (Lynx lynx) kill rate on roe deer (Capreolus capreolus) using radio-collared lynx. Kill rate was strongly affected by lynx social status. For males it was 4.85 ± 1.30 S.E. roe deer per 30 days, for females with kittens 6.23 ± 0.83 S.E. and for solitary females 2.71 ± 0.47 S.E. We found very weak support for effects of prey density (both for Type I (linear) and Type II (non-linear) functional responses) and of season (winter, summer) on lynx kill rate. Additionally, we analysed the growth rate in a roe deer population from 1985 to 2005 in an area, which lynx naturally re-colonized in 1996. The annual roe deer growth rate was lower after lynx re-colonized the study area, but it was also negatively influenced by roe deer density. Before lynx colonized the area roe deer growth rate was λ = 1.079 (± 0.061 S.E.), while after lynx re-colonization it was λ = 0.94 (± 0.051 S.E.). Thus, the growth rate in the roe deer population decreased by Δλ = 0.14 (± 0.080 S.E.) after lynx re-colonized the study area, which corresponded to the estimated lynx predation rate on roe deer (0.11 ± 0.042 S.E.), suggesting that lynx predation was mainly additive to other mortality in roe deer. To conclude, this study suggests that lynx predation together with density dependent factors both influence the roe deer population dynamics. Thus, both top-down and bottom-up processes operated at the same time in this predator-prey system.

Recovery of large carnivores in Europe's modern human-dominated landscapes.

The conservation of large carnivores is a formidable challenge for biodiversity conservation. Using a data set on the past and current status of brown bears (Ursus arctos), Eurasian lynx (Lynx lynx), gray wolves (Canis lupus), and wolverines (Gulo gulo) in European countries, we show that roughly one-third of mainland Europe hosts at least one large carnivore species, with stable or increasing abundance in most cases in 21st-century records. The reasons for this overall conservation success include protective legislation, supportive public opinion, and a variety of practices making coexistence between large carnivores and people possible. The European situation reveals that large carnivores and people can share the same landscape.

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.

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.

Habitat selection and risk of predation: re-colonization by lynx had limited impact on habitat selection by roe deer.

Risk of predation is an evolutionary force that affects behaviors of virtually all animals. In this study, we examined how habitat selection by roe deer was affected by risk of predation by Eurasian lynx - the main predator of roe deer in Scandinavia. Specifically, we compared how habitat selection by roe deer varied (1) before and after lynx re-established in the study area and (2) in relation to habitat-specific risk of predation by lynx. All analyses were conducted at the spatial and temporal scales of home ranges and seasons. We did not find any evidence that roe deer avoided habitats in which the risk of predation by lynx was greatest and information-theoretic model selection showed that re-colonization by lynx had limited impact on habitat selection by roe deer despite lynx predation causing 65% of known mortalities after lynx re-colonized the area. Instead we found that habitat selection decreased when habitat availability increased for 2 of 5 habitat types (a pattern referred to as functional response in habitat selection). Limited impact of re-colonization by lynx on habitat selection by roe deer in this study differs from elk in North America altering both daily and seasonal patterns in habitat selection at the spatial scales of habitat patches and home ranges when wolves were reintroduced to Yellowstone National Park. Our study thus provides further evidence of the complexity by which animals respond to risk of predation and suggest that it may vary between ecosystems and predator-prey constellations.

Higher levels of multiple ecosystem services are found in forests with more tree species.

Forests are of major importance to human society, contributing several crucial ecosystem services. Biodiversity is suggested to positively influence multiple services but evidence from natural systems at scales relevant to management is scarce. Here, across a scale of 400,000 km(2), we report that tree species richness in production forests shows positive to positively hump-shaped relationships with multiple ecosystem services. These include production of tree biomass, soil carbon storage, berry production and game production potential. For example, biomass production was approximately 50% greater with five than with one tree species. In addition, we show positive relationships between tree species richness and proxies for other biodiversity components. Importantly, no single tree species was able to promote all services, and some services were negatively correlated to each other. Management of production forests will therefore benefit from considering multiple tree species to sustain the full range of benefits that the society obtains from forests.