We present an individual-based model of the Norwegian wolf Canis lupus population, which is used to evaluate the effectiveness of current and potential management policies in fulfilling the Norwegian Government's stated aim of maintaining three breeding packs within a designated wolf zone. The model estimates the ‘functional extinction rate’ of the population, defined as the proportion of years in which breeding wolf packs are absent. Under the current conditions according to estimates from Scandinavia, with observed values of natural survival rates (0.903) and unauthorised mortality (0.203) and allowing for immigration from Sweden, the model predicts that the probability of functional extinction is as low as 0.07. This output variable is highly sensitive to the demographic parameters, and if alternative estimates of natural survival rates (0.73 for cubs and 0.83 for adults) reported for wolves elsewhere and higher rates of unauthorised mortality (0.4) are utilised, the functional extinction rate is 0.67 ± 0.15, and the population is dependent on maintenance by immigration from Sweden. The main determinants of the functional extinction rate are the unauthorised mortality rate and the immigration rate from Sweden. The Scandinavian population as a whole shows a rapid non-linear increase in probability of extinction at unauthorised mortality rates >0.10. Varying levels of the current management interventions (increasing the size of the wolf zone and target number of packs) are ineffective; only when the unauthorised mortality rate falls below 0.30 is a self-sustaining population in Norway able to establish. An adaptive harvest policy with culls targeted only at dispersing animals, or taking place only when the population exceeds a threshold level, could be sustainable if the unauthorised mortality rate is reduced. The fact that the Norwegian Government has been explicit about its management strategy and objectives has allowed us to test the ability of this strategy to meet the objectives, and we have shown that it is dependant upon maintaining the current circumstances alongside a high adult survival rate to be able to do so. Given the potentially critical role of the Swedish population in sustaining Norway's wolves, there is a strong case for joint management of the Scandinavian population. These insights are likely to be relevant for the management of other species living across geopolitical boundaries.
Large carnivore management regimes are often controversial. The extremely mobile and wide-ranging nature of these species means that they can rarely be conserved just within protected areas (Linnell et al. 2005b). Successful conservation of large carnivore populations involves management of entire landscapes, often across administrative boundaries. The conservation of trans-boundary populations is particularly challenging, as management regimes may vary hugely as animals move over borders. For example, in Scandinavia, Sweden has a large population of brown bears Ursus arctos, whereas the Norwegian population is dominated by immigrant young males from Sweden (Swensonetal. 1998). For the wolverine Gulo gulo, the entire northern population in Sweden and Norway is very sensitive to off-take on the Norwegian side of the border (Sæther et al. 2005). Nevertheless, the populations are managed independently.
Throughout much of the northern hemisphere, the most controversial large carnivore species, with a particularly extensive range, is the grey wolf Canis lupus. Wolves have for thousands of years lived in close proximity to mankind, and the relationship has often been volatile (Mech 1970). In recent centuries, the perception of the wolf has been as an agricultural pest, a competitor with hunters, and a direct threat to human safety. Consequently, the wolf has experienced a sharp decline globally, and extirpation from much of its former range (Mech & Boitani 2003). However, in common with other large charismatic carnivores, many people now recognise the wolf as a vulnerable component of global biodiversity, as a tool for ecological management, and as a potent ‘flagship’ symbol of the wild (Linnell et al. 2000). The wolf thus epitomises many of the issues faced by policy-makers aiming to enable coexistence of people and wildlife (Woodroffe et al. 2005).
Norway's wolves
On the Scandinavian peninsula, having faced deliberate persecution for centuries, wolves were considered functionally extinct by the 1960s (Wabakken et al. 2001b). However, in 1983, a single breeding pack was discovered in southeastern Norway, close to the Swedish border, which had probably immigrated from Finland or Russia (Wabakken et al. 2001b, Linnell et al. 2005a). Their numbers slowly increased, but faltered under the effects of genetic depression and demographic stochasticity (Liberg et al. 2005) until apparently being sustained by the arrival of a third immigrant in 1991 (Vila et al. 2003). Since then, the population has expanded rapidly; by 2005, the total Scandinavian population was generally held to have grown to 119–143, with the Norwegian subpopulation numbering ∼ 20 wolves (Wabakken et al. 2005).
Since its re-establishment, the Norwegian wolf population has attracted substantial controversy, as well as intensive research and population management (see e.g. Nilsson 2003, Skogen 2001, Wabakken et al. 2001a). Politically, the wolf became a focus for the more general conflict between rural and urban communities, the former leaning toward eradication of the wolf as a threat to humans, livestock and big game, and the latter toward its conservation (Skogen 2001). In the face of strong demands from the anti-carnivore lobby, in 2000, the Norwegian government sanctioned a culling programme. One of the more recent culls, in 2005, killed five of the 22–24 individuals in the country, comprising an alpha female and two scent-marking pairs, leaving 16 wolves and only one functional breeding pack. This led to major international concern about the future of Norway's population (information from reports by animal-interest groups such as Proact Mammal Campaigns 2006, and the WWF and BBC news websites; confirmed through expert consultation).
The stated aim of the Norwegian government is to maintain three breeding pairs within a designated wolf zone, as a compromise between wolf conservation and livestock rearing in the region (Miljøverndepartementet 2004). This three-pack limit does not include packs that straddle the border. From a biological point of view, Norway's wolves form one contiguous population with Sweden's wolves. In Sweden, government policy at the time of writing is to have 20 breeding packs within the country, but the whole wolf population is likely to be affected by the Norwegian culling given that the current population is still only around 140 animals. This has led to correspondence between Sweden and Norway, regarding their joint responsibility for wolf conservation under the Convention on Biological Diversity and the ‘Bern’ Convention; however, the Minister of the Environment of Norway recently stated that there will be no joint management of large carnivores with Sweden (Erik Solheim, official statement on 25 October 2007).
Modelling wolf behaviour
There have been several Individual Based Models (IBMs) developed for wolves (e.g. Vucetich et al. 1997, Haight et al. 1998, Chapron et al. 2003, Nilsen et al. 2007), while Nilsson (2003) used a structured matrix approach to investigate the effects of catastrophes, hunting and inbreeding on the Scandinavian wolf population. Population modelling is a useful tool to guide the management of endangered species, allowing the exploration of the impact of different management strategies on population viability for a given set of limiting factors (Coulson et al. 2001). Socio-political considerations determine the acceptability of particular strategies, but modelling can evaluate the likely impact of different strategies on the viability of Norway's wolves, and hence on the government's ability to fulfill its stated management objectives. Models can also be used to explore the biological consequences of the objectives themselves. In this paper, we develop a stochastic individual-based model for wolf population dynamics, and use it to explore the effects of key demographic parameters on the viability of the Norwegian wolf population. The demographic parameters include variation in mortality, dispersal, pack number and pack establishment rates linked to different management strategies, as well as the extent of unauthorised killing. We then ask what effect the target of three breeding packs and the current size of the wolf zone have on the viability of both the Norwegian and the wider Scandinavian population, and hence whether Norway's management aim is compatible with Sweden's. This modelling framework is applicable to wolf population dynamics more broadly, while the exploration of management options is both an interesting case study of the interactions between social and biological factors in carnivore management, and of the quantitative evaluation of the implications of a government's stated conservation objective.
Methods
Due to the importance of the wolf's social structure for its population dynamics and the small size of the Scandinavian population, we used an individualbased model (IBM). We modified the model developed by Nilsen et al. (2007) to reflect the biology of the Scandinavian population. Other than Nilsson's (2003) model, which is not an IBM, previous models of wolf population dynamics have mainly focussed on the interaction between wolves and their prey. In our study, we changed the focus towards assessing options for wolf population management. Given this focus, and the fact that the Scandinavian population is not likely to be prey-limited due to its small size and the abundance of suitable prey in the area (Nilsen et al. 2005), we treat wolf-prey interactions in a simplistic manner and focus on the dynamics of the wolf population itself. Also, we consider only the demography of the population, and we do not include management for genetic variation; Nilsson (2003) warns that inbreeding is a potentially serious threat to the long-term viability of the Scandinavian population, but our focus is on short-term management for population viability, and hence, genetic considerations are outside our remit. Nevertheless, it is clear that genetics have an important role to play in such a small and isolated population of mammals, and hence, the results of our study should be considered in this light (see also Liberg et al. 2005).
The Norwegian subpopulation was treated as an isolated population with a trickle of immigrant wolves from an unstructured outside source (i.e. Sweden); there was no spatial demographic substructuring of the population as a whole. Each individual within the IBM is uniquely identified and has an age, sex, stage and pack membership associated with it. Transitions between stages are related to age and social status, and all transitions take place probabilistically for each individual in random order. Each pack consists of a dominant pair and their offspring, and only the dominant pair reproduce. The initial population was assumed to be 21 wolves, comprising three packs of six individuals with three dispersers, to simulate the current Norwegian population: each pack contained a 4-yearold alpha pair with four descendants (two males and two females) of ages 1–3, and the dispersers (two males and one female) were aged between two and three years.
Our model has four compartments, reflecting different wolf life-stages (Fig. 1A), and six demographic processes are modelled, in the following order: survival, dispersal, establishment of a territory, anthropogenic mortality, reproduction and immigration from Sweden (see Fig. 1B and Table 1). Each year, the previous year's cubs can either disperse, remain in the natal pack as a subadult, or die. The existing subadults can disperse, die or remain in the pack as subdominant adults. Subdominant adults each year have a probability of dispersal and mortality, and otherwise, they remain in their natal pack. Dominant animals have a probability of survival and a reproductive rate. If a dominant animal dies, it is replaced by the oldest subdominant animal of the right sex in the pack. If there is no suitable animal available within the pack, the position remains vacant until it is filled by a dispersing animal or until a young individual moves up into the subdominant adult stage. The pack cannot reproduce until a new dominant pair is established. If both dominant animals die, the rest of the pack disperses. We made this assumption primarily to simplify the model, but it is considered reasonable with regards to the literature: Fuller et al. (2003) support our assumption, and in Brainerd et al. (2008), the loss of both breeders was noted to cause pack dissolution in 85% of cases, and in only 9% of cases did a pack reproduce the following season under such circumstances.
Figure 1.
Model presentation of the wolf life cycle with symbols for the vital rates referring to the parameters listed in Table 1. Note that the transition from subdominant to dominant adults is not automatic, but depends on the death of an alpha animal. The processes represented are movements between stage classes, and survival rates are age specific, hence there can be several different survival rates operating in one stage. In B), the order of implementation of the demographic processes in the model is given on an annual basis.
![f01_412.eps](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/f01_412.jpg)
Table 1.
Parameter values used in the model. ‘Symbol’ refers to symbols used in Figure 1A. Baseline values are probabilities unless otherwise stated; values in brackets give the values recently published by Liberg et al. 2008. The standard deviations (SD) for the distributions around the mean values used in the sensitivity analyses are shown in the SD column. Parameter values were taken from the literature, except those marked ‘Expert’ which were conservative estimates reached during consultation with wolf experts in Norway.
![t01_412.gif](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/t01_412.gif)
The first step in the annual cycle is that each individual other than the dominant pair has a probability of joining the pool of dispersing animals. Each disperser > 2 years old then has a probability of establishing a territory. If an animal is establishing a territory, it first looks for an existing pack with a lone alpha animal of the opposite sex, and next, it establishes a new territory and waits for a partner. If neither is available, it remains in the disperser pool. It is assumed that the probability of any disperser being successful in establishing a new territory (and hence, the total number of territories available) is linearly related to the size of the designated wolf zone.
Once the territory establishment process is complete, dominant pairs reproduce, and produce a large or small litter with a given probability, the size of the litter being based on data in the literature. For instance, we took the number of cubs in a small litter from Mech (1970), which relates to data from a study of thinly distributed wolves in Alaska. Next, immigration from Sweden occurs, with the immigrants joining the disperser pool. These immigrants are assumed to be subadults and are assigned a sex with a probability based on empirical observation (Mech 1970); note that these observations relate to studies of a small population in Finland, which was repopulating at the time, and hence, it was far from carrying capacity, so it was appropriate for our model.
Next, we modelled survival of each individual as a function of overall wolf population density and prey availability. As we did not consider prey as limiting, but as fluctuating from year to year due to climatic conditions, we modelled the prey availability by adding an error term to each survival rate, drawn from a uniform distribution with a mean of 1.0 and a standard deviation of ± 10%, and varying on an annual basis. In the absence of explicit feedback between the wolf and prey populations, we modelled density dependence in a ratio-dependent manner; that is, beyond the threshold population size at which density dependence begins to act, all survival rates were multiplied by the ratio of the threshold population size and the current total wolf population (Eberhardt et al. 2003). We set this threshold at 100 wolves based on expert consultation; this is the number estimated to fit comfortably within the current wolf zone (Pedersen et al. 2005). Basic survival rates were age dependent, and not related to other factors such as the life stage of the wolf (e.g. survival rates were not modified based on whether the wolf was subdominant, dominant or dispersing). This is ecologically unrealistic; dispersing wolves in particular are likely to be at greater risk of mortality than wolves resident in a territory (Pletscher et al. 1997). However, in our model dispersing wolves do not materially affect the functional extinction rate until they find a territory, and therefore disperser mortality is subsumed into the other parameters.
Finally, we applied anthropogenic mortality of two kinds: legal culls and unauthorised mortality. The latter includes both deliberate and accidental killing, such as road deaths, shooting or poisoning. We assume that unauthorised mortality targets individuals from all stage classes with equal probability, whereas legal culls could potentially target individuals from particular stage classes. The parameters in Table 1 are the result of a full literature review and are referred to as the ‘baseline’ parameters throughout. Liberg et al. 's (2008) estimates of survival rates (including unauthorised mortality) for wolves in Sweden are substantially different to the values obtained from the literature. We explicitly considered the effect of using these values instead of the baseline values, due to the sensitivity of the model to these two parameters.
Our model was coded in R (R Core Development Team 2007), and each simulation was run 50 times over 100 years in order to generate the main results; we considered 50 the minimum number of repetitions necessary statistically to keep the confidence intervals (CIs) acceptably narrow, and this minimum was used due to the time taken to run the model. The output parameter of interest is the number of years in the 100 year simulation in which there are no functional breeding packs. This may occur when the total population size is above zero, and does not represent population extinction, partly because of internal recruitment of animals into breeding from the subdominant stages, and partly because of immigration from the contiguous Swedish component of the population. A true probability of extinction is not possible to calculate for a population receiving continuous immigration, but because the Government of Norway's stated aim is to maintain a breeding population of wolves, this metric is a policy-relevant measure of management failure. We call this metric the ‘functional extinction rate’; hence, the functional extinction rate is equivalent to the proportion of the 100 year time span of the model, in which no alpha pairs were left in the population.
Table 2.
Results of a sensitivity analysis in which a generalised linear model was fitted with a binomial link, and the functional extinction rate as the dependent variable. The table shows the coefficients of all the parameters in the full model, and the significant interactions. Significance is coded as: *** = P<0.001, ** = P<0.01, * = P<0.05, NS = P>0.05.
![t02_412.gif](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/t02_412.gif)
We ran sensitivity analyses by varying all parameters simultaneously both over a feasible range, given in Table 1, and by picking values from a distribution with a standard deviation of 10%, according to the method described by McCarthy et al. (1995). We ran the model for 500 random parameter combinations, and each was repeated 50 times. We used the wolf population's functional extinction rate as the dependent variable in a general linear model (GLM) with a binomial link, including all parameters and their interactions. The linearity of the relationship between the parameters and the functional extinction rate was examined by plotting the parameter values against the logit of the extinction risk (McCarthy et al. 1995), and the significance of each parameter as a determinant of extinction risk was inferred from the coefficients of the GLM.
Sensitivity analyses revealed that the parameters to which the model was most sensitive were the immigration rate from Sweden and the unauthorised mortality rate (Table 2). These are the two parameters with the highest uncertainty attached, but are also parameters that can, to some extent, be controlled by management interventions. Other important parameters were the probability that an immigrant wolf from Sweden was a male, the probability of establishing a territory and survival rates (see Table 2). When the sensitivity analysis was run with a ± 10% change in parameter values, rather than varying them over a feasible range, the results were qualitatively the same, except that the importance of the survival rates increased. As expected, there was no discernable relationship between most of the variables and the functional extinction rate. The relationship between the functional extinction rate and the immigration and unauthorised mortality rates was linear with an r2 of 0.61 and 0.27, respectively. The total uncertainty in extinction risk that stems from the uncertainty in the initial parameter values can be calculated based on the inherent uncertainty from the three most influential variables (unauthorised mortality levels ± 10%, immigration rate ± 10%, probability of success in establishing a territory ± 5%). From these three dominant sources of model sensitivity, the total uncertainty in extinction risk (due to uncertainty in initial parameter values) can be estimated at ± 15%, using the basic error formula:
where σy = uncertainty in the dependent variable, σxi = uncertainty in the ith independent variable, y = the dependent variable, xi = the ith independent variable, i = 1–3.This is a conservative estimate of the experimental uncertainty in the extinction risk, and assuming that the relationship between extinction risk and any independent variable is of constant form as parameter values vary, it was applied to the rest of the results where appropriate.
We carried out validation of overall model behaviour by applying the model to data from the Yellowstone wolf population (it was not possible to validate the model against data from the Norwegian population, as the relevant time series is too short). The Yellowstone population grew from a small initial population introduced in 1995, in an area of abundant prey, according to available data (taken from the official web-based data stores maintained by the US Centre for Biological Diversity, and Yellowstone National Park, accessed in 2006). Hence, the qualitative behaviour of the population in Yellowstone is likely to be similar to that of the Norwegian population in the absence of human interference, particularly given the constancy of wolf biology over a wide geographic range (Mech 1970). Therefore, starting with the known wolf population introduced to Yellowstone in 1995, the model was used to predict the subsequent growth of that population. The official data sources mentioned were used to establish upper and lower bounds of the estimated wolf population in Yellowstone from 1995 onwards, and the model was found to predict a wolf population that was comfortably within these bounds, thus validating general model behaviour. We carried out validation of the estimate for immigration rate by applying the model to the Scandinavian population as a whole and calculating the net immigration rate across the Swedish-Norwegian border (Scenario 1 b in Table 3). The model was then used to explore a number of management scenarios (see Table 3).
Table 3.
Different management scenarios explored using our model. Scenarios numbers 2–5 model the Norwegian population.
![t03_412.gif](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/t03_412.gif)
Results
The baseline model
The baseline parameter values given in Table 1 produced a population that fluctuated dramatically, with a mean number of wolves of 4.2 ± 3.4 (the standard deviation is used throughout), and an average of 67 ± 15% of years having no functional breeding packs (Scenario 2a, Fig. 2A). With these rates, the Norwegian wolf population is sustained by immigration from Sweden, and would otherwise quickly go extinct. If we assume no legal culling, the functional extinction rate is virtually unchanged, at 0.64 ± 0.15 (Scenario 2b). The results are, however, sensitive to parameter values, and the functional extinction rate (Scenario 2a, 2b) is significantly reduced to 0.07 ± 0.15 if the far higher natural/ unauthorised survival rates estimated by Liberg et al. (2008) are utilised. In the absence of either unauthorised mortality or legal culling, and with constant immigration from Sweden, the model predicts a zero functional extinction rate, and instead rapid growth to an equilibrium population fluctuating around 250 animals (Scenario 2c, see Fig. 2B). Given that our density dependence and prey dependence assumptions are crude, this result is not intended predicatively, but merely as an exploration of the model behaviour. The model assumptions are best suited to exploring dynamics at low population sizes when there is negligible resource limitation.
Figure 2.
A single run of the baseline model (see parameter values in Table 1; Scenario 2a in Table 3), showing the number of individual wolves and of packs in Norway over a 100-year period (A). In this run, the functional extinction rate (number of years with no breeding packs) was 0.66, and the mean size of the wolf population was 4.16. In B), a run of the baseline model with no culling or unauthorised mortality (Scenario 2c) shows that the population stabilises at a mean size of 263 wolves in 47 packs.
![f02_412.eps](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/f02_412.jpg)
The Norwegian population is only one part of the Scandinavian population, and so the model was run for the Scandinavian population as a whole, for comparative purposes. The following specifications were used (Scenario lb): there was no immigration (the Scandinavian population as a whole is effectively isolated), no culling (the Norwegian culls act on a small proportion of the population, so were ignored for simplicity), unauthorised mortality levels were set to 0.20 (as estimated by Liberg et al. 2008), and the initial population was 15 packs (approximately the 2005 size). The metric of extinction was the proportion of runs in which the population reached zero (because the population is isolated and assumed not recolonisable). With the baseline set of parameter values, the extinction rate was 0.92. This value was highly dependent on the unauthorised mortality rate (Fig. 3), with a probability of extinction of zero when unauthorised mortality was zero, and a rapid, non-linear increase in probability of extinction between unauthorised mortality rates of 0.10 and 0.20. The extinction probability was also far higher than the 0.05–0.29 calculated by Nilsson (2003) for various short-term scenarios using an agestructured model rather than an IBM. The discrepancy is likely to be due to a number of differences in the model, including the lower unauthorised mortality rate and higher initial population assumed by Nilsson (2003); and indeed, when rerun using the survival rates estimated by Liberg et al. (2008), our model does agree with Nilsson (with a mean functional extinction rate of 0.06).
Figure 3.
Proportion of runs in which the Scandinavian population as a whole is extirpated (Probability of Extinction), shown as a function of the unauthorised mortality rate (Scenario 3d; baseline value - 0.20).
![f03_412.eps](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/f03_412.jpg)
We used the model of the Scandinavian population as a whole to calculate the net number of dispersers moving from the Swedish to the Norwegian populations, on the assumption that Norway has three of the 15 packs. The value was calculated as 4.6 ± 2.7 wolves, not significantly different from the 4 ± 1 animals used in the model, which was based on expert opinion. This gave us further confidence in the validity of the model assumptions.
Management interventions
The current management regime in Norway revolves around a target maximum number of packs (currently three), with the presumption that any packs formed above this number will be broken up by culling the dominant pair. The other major component of the management regime is the designation of a ‘wolf zone’, within which wolves are allowed to exist, and outside which animals are culled. We modelled changes in the size of the wolf zone through the proxy of the probability of establishing a territory, on the assumption that the probability that a dispersing animal is able to establish a territory is linearly related to the area available. This assumption is questionable, but is taken in the absence of empirical evidence about the relationship between the area of the wolf zone and the wolf demography. On varying the target maximum number of packs (Scenario 3a) and the probability of territory establishment (Scenario 3b), we see that above a target number of packs of about five, the probability of functional extinction is unchanged as the target number of packs increases. The total number of wolves in the area (irrespective of whether they are in packs) is also insensitive to the target number of packs after the initial increase to five packs, and is very low, between five and seven animals (Fig. 4A). The effect of an increase in the probability of establishment is also minor. It is clear that the limiting factor for Norway's wolves is neither the size of the wolf zone nor the chosen target number of packs, and hence that varying either aspect of the current management regime will not improve the population's viability. Increasing the actual number of packs in the initial starting population rather than the maximum number of packs allowed, on the other hand, has a direct, linear relationship with functional extinction rate (Scenario 4a, Fig. 4B).
Figure 4.
Average number of wolves in the population after 100 years as a function of the target number of packs (A ; Scenario 3a; baseline value - 3). In B), the total size of the Norwegian population, as a proxy for re-establishment success, is plotted against the number of packs in the starting population (Scenario 4a; baseline value - 3).
![f04_412.eps](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/f04_412.jpg)
Figure 5.
Functional extinction rate of the Norwegian population plotted against the mean number of immigrants entering the population annually from Sweden (Scenario 3c; baseline value - 4).
![f05_412.eps](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/f05_412.jpg)
The functional extinction rate is highly sensitive to the immigration rate from Sweden, and as the immigration rate increases, the sensitivity of functional extinction rate to the target number of packs also increases (Scenario 3c, Fig. 5 and Table 4). This again demonstrates that under the current levels of unauthorised mortality, the main driver of the viability of Norway's wolves is the success of Sweden's management interventions. Similarly, looking at the Scandinavian population as a whole, the existence of immigration from Finland or Russia at the rate of just two animals per generation (the level required to ensure that inbreeding will not be an issue) reduced the overall extinction risk significantly, i.e. from 0.92 ± 0.10 to 0.60 ± 0.15 when using the baseline parameter values (Scenario 1c). Starting from a zero population (i.e. after a second extirpation of Norway's wolves), the probability of a viable population re-establishing through immigration from Sweden depends strongly on the unauthorised mortality rate (Scenario 4b, Fig. 6). The number of wolves in the population after 100 years increased to a high level when unauthorised mortality rates fell below 0.30 with other parameters at the baseline values, showing a self-sufficient population not reliant on continued immigration from Sweden. The same quantitative relationship is observed using the higher survival rates reported by Liberg et al. (2008), but with a higher unauthorised mortality threshold below which the population is relatively secure (0.40–0.50).
Table 4.
Threshold target number of packs at which the extinction probability levels out, and the value at which it levels out, for different immigration rates (number of wolves).
![t04_412.gif](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/t04_412.gif)
Figure 6.
Number of wolves left after 100 years, as a proxy for extinction risk, plotted against the unauthorised mortality rate. The initial population size is assumed to be zero and no culling takes place (Scenario 4b).
![f06_412.eps](ContentImages/Journals/wbio/15/4/08-010/graphic/WebImages/f06_412.jpg)
Public opinion may dictate that a legal wolf cull is required in order to ensure that wolves are accepted within the wolf zone. Hence, we explored the effects of a legal cull on wolf population viability. If the cull is a straight proportion of the population size, it leads to increases in functional extinction rate equivalent to the effects of raising the unauthorised mortality rate by the same amount (Scenario 5a, see Fig. 3). A decision rule was then tested by which a set number of dispersing wolves could be killed once the number of dispersers had reached a threshold level (Scenario 5b). This policy aims to maintain the stability of the breeding packs, and it has negligible effect on extinction probability. However, the cull is only possible in certain years, with the number of years dependent on the level at which the threshold is set. For example, if one wolf can be killed every time there are one or more disperser, 40 ± 5 wolves are killed in 100 years, whereas if two wolves can be killed when there are two or more dispersers, 10 ± 5 wolves are killed in 100 years. Again, this is dependant upon the assumed parameter values, and more wolves are culled if higher natural survival rates are assumed.
Given the low number of wolves actually killed, the feasibility of this strategy depends on whether the mere fact that there is a legal cull in principle is enough to satisfy the public, and upon whether it would be possible to implement such a policy in practice. We tested another variant of the legal hunting quota; a proportional threshold cull, which has been shown previously to be a robust hunting strategy under conditions of uncertainty (Scenario 5c; Engen et al. 1997). Under this strategy, a fixed proportion of the population, with individuals selected at random, is culled once the population size has exceeded a certain threshold. This was explored both under the baseline unauthorised mortality rate and on the assumption that unauthorised mortality levels are at half this level, following Liberg et al. (2008). Under baseline unauthorised mortality levels and for culling thresholds > 15 wolves, the population size was so low that hunting virtually never took place, and so the extinction risk was unaffected. However, with an unauthorised mortality rate of 0.2, the overall extinction risk was lower, and there were corresponding changes in the actual hunting rate and the extinction risk. For instance, with baseline parameter values, hunting beginning at a threshold of nine wolves (i.e. approximately two packs) and a proportional hunting rate of 0.18, on average 1-3 wolves could be legally culled a year, with a mean functional extinction rate of 0.29. By way of comparison, using the higher survival rates from Liberg et al. (2008), with a proportional hunting rate of 0.18, but a higher threshold of 15 wolves (approximately three packs), as many as 6–7 wolves could be culled on average per year, with a mean functional extinction rate of 0.09.
Discussion
Our model of Norway's wolf population estimates a functional extinction rate of 0.67 ± 0.15 under the baseline parameter values based on a literature survey, whereas it was 0.07 with the current high survival reported by Liberg et al. 2008. The former result is comparable to Vucetich et al.'s (1997) finding of a 0.70 extinction probability for an isolated population of 50 wolves over 100 years, also using an IBM. The functional extinction rate is completely dominated by the relative values of the immigration rate from Sweden and unauthorised mortality, and the outcome is insensitive to current management policies, such as the size of the wolf zone and the target number of wolf packs. This is a worrying finding, firstly because unauthorised mortality is difficult both to monitor and to control, and secondly because of the precarious state of the Scandinavian population as a whole. The extinction rate is sensitive to both natural and anthropogenic mortality rates. Currently, Norway appears likely to be a sink for the Swedish part of the population, completely reliant on Sweden to meet its management objectives. Norway also has a policy of higher management off-take for other carnivores, i.e. Norway has the dominant part of the trans-boundary population of wolverines, and Saether et al. (2005) modelled that current off-takes on the Norwegian side of the border make the northern population likely to go extinct. Wolves in Norway are not necessarily currently self-sustaining, and if the current situation remains unchanged, we would expect to see very low and highly variable population sizes (in line with empirical observations from recent wolf censuses in Norway; Wabakken et al. 2001a, 2002, 2004a, b, 2005), with regular functional extinction of the breeding packs followed by recolonisation from Sweden. The model does not account for any effects of low genetic variation (see the section Methods), which is likely to be a characteristic of such a small population. Although a discussion on genetics goes beyond the scope of our paper, it is likely that this would be a further factor contributing to a high probability of functional extinction of the Norwegian wolf population (see Nilsson 2003, Liberg et al. 2005).
However, our analysis has also revealed some potential ways forward. We show that if unauthorised mortality were brought under control, Norway's population could grow to a self-sustaining level within the current wolf zone, assuming they can utilise this area to the full. With only a few additional immigrants per year from Finland or Russia, and with reduced unauthorised mortality, the Scandinavian population is potentially viable both ecologically and genetically. If Liberg et al.'s (2008) estimates of unauthorised and natural mortality in Scandinavia as a whole are also applicable to the Norwegian subpopulation, this would be enough to reduce the functional extinction rate to a low level. Research to reduce uncertainty surrounding estimates of mortality rates in Norway, and actions to reduce these mortality rates, are key priorities for scientists and managers aiming to ensure a viable, self-sustaining wolf population in Norway, as envisaged in government policy.
A permit-based legal cull, of the proportional threshold type recommended by Engen et al. (1997), could provide a means of control by providing a mechanism for legal wolf hunts. By its precautionary nature, it would allow Norway to fulfil its conservation objectives and support a more robust wolf population; potentially also reducing reliance upon the ‘wolf zone’ approach to conservation, which can lead to social conflict. Given that many Norwegians are in favour of a continued wolf presence in Norway, and that the Swedish wolf population is currently indirectly supporting Norway's wolf population through cross-border movements, there is a case for increased cooperation between the two governments. This might take the form of a joint wolf management agreement, in which Sweden was recognised as the partner maintaining the majority of the wolf population and acting as a source for Norway's population, while Norway took the role of ensuring immigration from Finland or Russia to safeguard the genetic viability of the joint population. The economics of management of joint stocks have been explored for other species (e.g. Skonkoft 2005), and this body of theory could form the basis for an equitable agreement.
Our study is a useful contribution to the wolf management literature, partly because we have been explicit about testing the implications of the Norwegian government's stated management objectives. We have shown which parameters impact most strongly on the government's ability to fulfil these objectives, and which are most in need of improved monitoring in order to reduce uncertainty about outcomes. We have also been able to demonstrate the relative effectiveness of current and potential management interventions. It is important that governments state their objectives and management strategies in the transparent way that the Norwegian government has done, as this enables evaluations of this sort to be carried out.
Acknowledgements
we thank NERC and the Research Council of Norway (YFF to Atle Mysterud) for financial support of this study, and all experts and other consultants for their help and advice, including particularly Leif Egil Loe, Jos Milner, Torstein Storaas, Thomas Stromseth and Barbara Zimmerman.