Establishment probability in newly founded populations
© Gusset et al.; licensee BioMed Central Ltd. 2012
Received: 16 February 2012
Accepted: 13 June 2012
Published: 20 June 2012
Skip to main content
© Gusset et al.; licensee BioMed Central Ltd. 2012
Received: 16 February 2012
Accepted: 13 June 2012
Published: 20 June 2012
Establishment success in newly founded populations relies on reaching the established phase, which is defined by characteristic fluctuations of the population’s state variables. Stochastic population models can be used to quantify the establishment probability of newly founded populations; however, so far no simple but robust method for doing so existed. To determine a critical initial number of individuals that need to be released to reach the established phase, we used a novel application of the “Wissel plot”, where –ln(1 – P0(t)) is plotted against time t. This plot is based on the equation , which relates the probability of extinction by time t, P0(t), to two constants: c1 describes the probability of a newly founded population to reach the established phase, whereas ω1 describes the population’s probability of extinction per short time interval once established.
For illustration, we applied the method to a previously developed stochastic population model of the endangered African wild dog (Lycaon pictus). A newly founded population reaches the established phase if the intercept of the (extrapolated) linear parts of the “Wissel plot” with the y-axis, which is –ln(c1), is negative. For wild dogs in our model, this is the case if a critical initial number of four packs, consisting of eight individuals each, are released.
The method we present to quantify the establishment probability of newly founded populations is generic and inferences thus are transferable to other systems across the field of conservation biology. In contrast to other methods, our approach disaggregates the components of a population’s viability by distinguishing establishment from persistence.
Trying to (re)establish populations by releasing individuals into suitable habitat is an important element of modern conservation practice. The success of such release attempts depends largely on two factors, namely the newly founded population reaching the established phase and, once this stage is reached, maintaining itself in the release area (on the importance of this distinction, see ). The established phase is defined by characteristic fluctuations of the population’s state variables (e.g. number of individuals or age structure), in which case population dynamics is no longer affected by initial conditions . Establishment success often depends on how many individuals or groups are released, but high post-release dispersal can create a disparity between release population size and the effective initial size of a newly founded population . Establishment does not imply persistence, as the ecological capacity of the release area may be small and environmental stochasticity large (for a visualization of the established phase, see ).
In a previous study , we used a stochastic population model of the endangered African wild dog (Lycaon pictus) to quantify the critical initial number of packs (two) and individuals per pack (six) necessary for a reintroduced population of this species to maintain itself in the release area. Persistence was virtually impossible, unless the population was frequently supplemented. Given sufficient prey, the intervals between artificially adding a pack seem to be the most important factor governing the persistence of a small, reintroduced wild dog population. However,  focused on persistence and post-release management, but not on the probability that the newly founded population reaches the established phase in the first place.
Two important questions arising in any release attempt  thus are: (1) How many individuals or groups should be released so that population establishment is more or less guaranteed? (2) In the case of no supplementation, how is persistence of an established population affected by post-release dispersal (i.e. loss of dispersers due to emigration from the release area)? Here, we present a method to quantify the establishment probability of newly founded populations, using attempts to reintroduce wild dogs for illustration. In contrast to other methods, our approach disaggregates the components of a population’s viability by distinguishing establishment from persistence.
To tackle these questions, we used our previously developed individual-based model for wild dogs . Individual-based models enable us to explore how population characteristics emerge from the ways in which individuals behave and interact with each other . The model includes social structure and behaviour, but nevertheless is conceptually simple. It is parameterized with data from a 25-year field study in Hluhluwe-iMfolozi Park, South Africa .
In short, the model was designed to predict the probability of small, reintroduced populations of wild dogs establishing themselves and persisting in the release area under various levels of reintroduction effort (for details, see ). In contrast to , no post-release management interventions occurred in the present application of the model and disperser groups were simulated to leave the release area with various probabilities.
Our validation procedure ensured that the model correctly captures internal relationships between variables and to some degree the internal organization of the real system (see ). This suggests that the model is appropriate for its intended purpose, as it could reproduce multiple output patterns observed at different hierarchical levels of the system , which were not imposed onto the model but emerged from interactions between the simulated individuals, packs and disperser groups.
Model parameters, reference values and results of the local sensitivity analysis for a reintroduced wild dog population (initial condition: number of packs = 4; probability for a disperser group to leave release area = 0.4)
+10% of parameter value
–10%of parameter value
Reproduction in newly formed packs (p)
Reproduction in established packs (p)
Pack size (v)
8.1 ± 1.1
Litter size (v)
7.9 ± 0.8
Primary sex ratio (p)
0.55 ± 0.06
Ecological capacity (v)
Density dependence threshold (v)
Dispersal in males (p)
Dispersal in females (p)
Disperser group size threshold (v)
Pack formation (p)
Dominant displacement (p)
Mortality in male pups (p)
0.07 ± 0.06
Mortality in female pups (p)
0.16 ± 0.14
Mortality in yearling males (p)
0.29 ± 0.14
Mortality in yearling females (p)
0.20 ± 0.20
Mortality in young adult males (p)
0.17 ± 0.08
Mortality in young adult females (p)
0.01 ± 0.01
Mortality in old adult males (p)
0.30 ± 0.16
Mortality in old adult females (p)
0.22 ± 0.16
Dispersal mortality in males (p)
Dispersal mortality in females (p)
Catastrophe occurrence (p)
Catastrophe severity (p)
Our model thus appears to capture the essential characteristics of a real wild dog population and to be relatively robust to parameter uncertainty (see ). Collectively, this suggests that the model is structurally realistic  enough to place confidence in inferences about real wild dog populations based on modelling results.
To determine a critical initial number of packs (consisting of eight individuals each; ) to ensure establishment, we used the “Wissel plot” (formerly referred to as “ln(1 – P0) plot”; ). This plot is based on the equation , which relates the probability of extinction by time t P0(t), to two constants, c1 and ω1. The former, c1, reflects the initial state of a population at time t = 0. If this state is in the range of states that can be observed in the established phase, c1 is equal to one; if a population initially is so small that it does not necessarily reach the established phase but might go extinct beforehand, c1 is smaller than one . The other constant, ω1, is independent of the initial state of a population and describes the probability of extinction per short time interval in the established phase, which is constant. The inverse of this risk, Tm = 1/ω1, can be defined as the “intrinsic mean time to extinction” . It describes the intrinsic persistence of a population in a given area and environment. For c1 = 1, Tm is equal to the arithmetic mean time to extinction that can be determined from repeated simulations starting from the same initial population .
If the initial number of packs is larger than one but still too small, c1 is smaller than one and thus the intercept of the (extrapolated) linear parts of the “Wissel plot” with the y-axis, which is –ln(c1), is positive (Figure 1). Figure 2 shows that the y-intercept is positive if the initial number of packs is two or three. Thus, to ensure establishment, four packs should be released. As expected, this is more than the number of packs (two) necessary for a reintroduced population to maintain itself in the release area with frequent supplementation . However, the extinction times presented in Figure 2 are very short. An intrinsic mean time to extinction Tm of 10,000 years corresponds to an extinction risk of 1% in 100 years , thus even the largest Tm obtained in our case results in an extinction risk exceeding 10%.
The asymptotic nature of our results for establishment (Figure 2) is suggestive of a critical threshold size typical for an Allee effect (cf. [3, 13]). From Figures 1 to 3 it becomes clear that it does not make sense to release more than four packs, because this neither increases population establishment nor promotes persistence. Strikingly, we previously established empirically that a critical minimum number of four packs, which simultaneously produce enough unrelated dispersers, are necessary for successful pack formation events to occur , and thus to maintain population viability. With a given number of individuals available for reintroduction, consecutive releases of several smaller packs (Figure 2 in ) result in higher population viability than a single release of a few larger packs (Figure 3), as this may buffer a newly founded population from environmental stochasticity, but such frequent supplementations may not be feasible.
Sensitivity analysis (Table 1) showed that reproduction (both the probability of producing a litter and litter size), and thus the production of future dispersers, most strongly affect population viability (cf. [14, 15]). Focusing conservation management on enhancing reproduction thus seems particularly advisable in this case, and behavioural traits were indeed found to most strongly affect the survival of reintroduced wild dogs . Variation in release pack size, on the other hand, had a small impact on population viability (Table 1), as we previously established empirically .
Release strategies are often based on intuition and trial-and-error rather than a critical appraisal of the available evidence . Reintroduction biology is a typical field where initial conditions are important: a population may not realize its intrinsic ability to persist because it is too small, and consequently goes extinct before establishment. The dichotomy between establishment and persistence is useful because newly founded populations can fail to reach the established phase in conditions that would enable persistence once the population is established .
Therefore, reintroduction biology and related disciplines (e.g. invasion and restoration ecology as well as the emerging field of assisted colonization) would likely benefit from using structurally realistic models as well as adopting the plot and concepts proposed by  to assess establishment separately from persistence (i.e. to differentiate between initial and intrinsic aspects). The method we present for doing this is generic and inferences thus are transferable to other systems. The plot of –ln(c1) in Figure 2 (intercept with y-axis) represents a novel application of the “Wissel plot”  to quantify the establishment probability of newly founded populations across the field of conservation biology. It should be emphasized, though, that the reliability of recommendations for real (re)introductions, which are based on the approach presented here, depends on the reliability of the underlying population model, which needs to be well documented, tested and validated.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.