On the brink between extinction and persistence
 Cino Pertoldi^{1, 2}Email author,
 Lars A Bach^{3} and
 Volker Loeschcke^{1}
DOI: 10.1186/17456150347
© Pertoldi et al; licensee BioMed Central Ltd. 2008
Received: 18 November 2008
Accepted: 19 November 2008
Published: 19 November 2008
Abstract
The nature of size fluctuations is crucial in forecasting future population persistence, independently of whether the variability stems from external forces or from the dynamics of the population renewal process. The risk of intercepting zero is highly dependent on the way the variance of the population size relates to its mean. The minimum population size required for a population not to go extinct can be determined by a scaling equation relating the variance to the arithmetic mean. By the use of a derived expression for the harmonic mean defined by the parameters of the scaling equation we show how it is possible to separate the domains of persistence from those of extinction and to facilitate the identification of populations on the brink of extinction.
Reviewers
This article was reviewed by Mark W. Schwartz (nominated by Peter Olofsson), Josef Bryja (nominated by Aniko Szabo) and WaiYuanTan. For the full reviews, please go to the Reviewers' Comments section.
Background
Natural populations are affected by σ^{2}_{e} and σ^{2}_{d} which in turn affect the expected time to extinction [1]. Although the time to extinction is expected to increase with population size [2], other factors influence the dynamics of populations as e.g. density dependent mechanisms and population growth rate [3, 4]. The specific population growth rate is affecting the longterm persistence of populations, whereas the most immediate effects on the risk of extinction are mainly due to σ^{2}_{e}. Especially in small populations, population persistence is affected by σ^{2}_{d} [5]. Reproduction and survival of offspring can be dependent or independent of population size and in combination cause populations to fluctuate. Moreover, recent work points to the importance of altered σ^{2}_{e} on the variation of vital rates, which obviously feedback on the demographic dynamics [6]. Studies developing methods for estimating population viability are growing steadily [7, 8].
Presentation of hypothesis
Here we propose a simple model to estimate the risk of extinction and population persistence based upon a twoparameter description of the harmonic mean (HM), defined by the parameters of the scaling equation [9]. This method allows a separation of the domains of population persistence versus those of extinction and enables the identification of populations on the brink of extinction as it allows the estimation of the minimum population size required for population persistence. The HM has the special property that a single occurrence/generation of zero suffices to cause ultimate population extinction – much like the behaviour of real populations in the absence of migration.
Testing the hypothesis
where K is a measure of individual level variability, and β is the scaling exponent. Taylor's power law is well documented for animal populations and suggests that increasing σ^{2}_{d} increases individual reproductive variance while increasing σ^{2}_{e} increases reproductive covariance. For populations experiencing constant per capita σ^{2}_{e}, the regression of log σ^{2} versus log $\overline{\mu}$ gives a line with a slope of 2 for trivial mathematical reasons. Data from time series of natural populations suggest that β may lie anywhere in the range of 0.6 to 2.8 [10, 11]. The range of K has been estimated empirically in many populations to lie in the range 0.10 <K < 8.32 [10, 11].
The degree of reproductive covariance among individuals affects the scaling exponent β. Completely correlated reproduction results in β = 2, while independent reproduction results in β = 1 [12]. Organisms that have highly correlated responses to environmental fluctuations will exhibit less variable reproductive patterns than organisms that experience a high degree of σ^{2}_{e} on the individual scale.
HM can be considered a proxy for N_{E}, especially for populations with discrete generations, whereas when there are overlapping generations, the time scale becomes crucial and must be defined such that it relates to the discrete models. Applying the power law to HM allows us to study the separation of domains of attraction to zero from domains of nonzero HM. Subsequently, these domains can be interpreted in terms of persistence and extinction of populations and hence help identifying populations destined for extinction. To explore the domains we investigated various combinations of the parameters β, K and $\overline{\mu}$. The consequences of β > 2 have been discussed earlier [9], and it has been found that HM, for certain values of β and K, first increases with increasing $\overline{\mu}$ but then starts to decrease for higher $\overline{\mu}$.
which is the necessary condition for population persistence for population fluctuations between two fixed sizes [16]. Moreover, the inequality is valid for any distribution of population abundances following normal or lognormal distributions [16].
Implications of the hypothesis
The methodology used will improve estimates of population viability and the fact that the parameters used in the model can be obtained from ordinary data sets gives this method a potentially wide applicability for comparing the chances of survival for fluctuating populations when facing environmental stochasticity. However, we note some of the potential difficulties associated with the estimation of the parameters of the model. The sampling variance or error in the measurement of population abundance, which generally scales directly with $\overline{\mu}$, would by itself produce a powerlaw slope of less than 2. This problem however is not important in data sets where the sample is large and thus sampling variance small.
The model is not taking into account the consequences of σ^{2}_{d}, which certainly affects the dynamics of populations. Such effects are, however, most prevalent for small populations. For larger population sizes, σ^{2}_{e} in per capita population growth rates leads to a variance that tends to scale with the square of the mean and hence overshadows the variance from σ^{2}_{d} which scales directly with $\overline{\mu}$. Hence, larger population size gives a better capacity of the model to predict the minimum viable population size. The model will therefore add to the ongoing debate about the relative importance of σ^{2}_{d} versus σ^{2}_{e} in determining the extinction risk. In fact we can deduce from the model that for certain values of K and β a population will become extinct even if its population size is sufficiently large to restrict the impact of σ^{2}_{d}. For the same reason we can see that a population can go extinct even if N_{E} is high (and therefore harbour relatively high levels of genetic variability).
The fact that β is not only often different from 2 but also changes over time, has some implications for the association between population variability and extinction risk. Intuitively, for a given average abundance, one expects the risk of extinction to increase with temporal variability, however, many studies conducted on longterm data from natural populations have found the opposite result [2, 17, 18]. These studies use temporal variability as a direct proxy for population vulnerability, where population variability measures are calculated from time series data (std logN, CV). However, both the logarithmic transformation and the CV can be only properly applied if the variances scale proportionally to the square of the mean (β = 2), and given the fact that β is often different from 2 and is also changing with time it is quite evident that these correlation studies are potentially afflicted by this problem. The reasons for the discordant results obtained in these correlational studies have also been subject of a long debate and the relative importance of density dependent processes on population dynamics has been compared to the relative importance of environmental variability [19].
β is determined by the degree of correlation of the reproductive efforts among individuals within a population. In populations where the density is low we will expect a lack of correlation between reproductive efforts among individuals as the intraspecific interactions between individuals are very low, whereas if the population reaches higher density (for example when approaching the carrying capacity) we will expect a higher amount of interactions between individuals and an increased value of β which will increase the minimum population size necessary for avoiding extinction. Hence, a population near the carrying capacity with β near the value of 2 should be more prone to extinction, as when an environmental stochastic event is added β will become larger than 2, which means an increased risk of extinction. This finding could explain some of the metapopulation dynamics observed where a colonizing population initially increases in size only to suddenly go extinct. With environmental variability we do not only mean the amplitude of fluctuations of the environment but we also consider the spectrum of the fluctuations; one pole of this spectrum is 'white noise', where values of the variable are independent and variability entirely stationary, and at the other end of the spectrum there is the random walk or densityindependent drift, for which the expected population variability grows at a rate proportional to the length of the series [20]. This effect has been linked to a reddened spectrum and there is increasing evidence that longterm ecological data sets showa reddened spectrum [21]. Inchausti & Halley (2002) concluded therefore that traditional measures of population variability need to be supplemented by a spectrum analysis.
The reddening of population dynamics has been suggested to be due to shortterm fluctuations (of small amplitude) which are superimposed on ever larger longperiod variations. However, the fact that β depends on the density of the population makes it quite evident that β and K (which are both influencing the variance of the population size) should be considered when interpreting the fluctuations of the population, and their changes in time should be taken into account when estimating the spectrum of the population dynamics. The need for diagnostic methods in population biology and conservation management undoubtedly becomes accentuated in the years to come [22]. In particular there will be a need for detecting regime shifts in the dynamic behaviour of populations as changes of the global environment begin to accelerate. This model allows an estimation of the importance of σ^{2}_{e} on the two parameters β and K and on how much alteration of the parameters will push the population towards the extinction threshold [23]. Environmental stochasticity will in fact increase or reduce the amplitude of the population fluctuations depending on the sign of the correlation between population size and environmental fluctuations as:
σ^{2}_{tot} = σ^{2} + σ^{2}_{e} + 2r(σσ_{e}),
where σ^{2}_{tot} is the variance of the population size in presence of environmental noise, σ^{2} is the variance of the population size in absence of environmental noise, σ^{2}_{e} is the environmental noise and r(σσ_{e}) is the covariance between the environmental noise and the population fluctuation. The covariance is given by two times the product of r and the std of the population size and the environmental fluctuations (σ and σ_{e} respectively). Hence, a negative correlation (r < 0) between environmental stochasticity and population fluctuations will decrease the fluctuations of the population size, with σ^{2}_{tot} < σ^{2}, whereas in case of a positive correlation (r > 0) we will observe an increase in population fluctuations in the presence of environmental stochasticity (σ^{2}_{tot} > σ^{2}).
Reviewers' comments
Reviewer's report 1
Mark W. Schwartz, Department of Environmental Science & Policy, University of California, Davis Davis, CA 95616. Nominated by Peter Olofsson.
Reviewer comments
This paper proposes using the harmonic mean of population to determine a minimum population size necessary to avoid extinction. This is equivalent to what conservation biologists have referred to as a 'minimum viable population'. In 1987 Michael Soule edited a book on Minimum Viable Populations. Subsequent to Soule's book there has been much literature on why the philosophical approach of estimating an MVP for conservation is inappropriate. Basically, MVP's encourage managers to manage to the lowest possible N. Mistakes in this realm being particularly costly. This, in fact, is exactly analogous to the problems with fisheries management. Promoting a measure that implies a bulletproof estimate of a population size that is not threatened with extinction just begs for mistakes and failures if managers target/permit exploitation to reduce populations to those numbers. Thus, although Pertoldi and colleagues propose an advancement in terms of thinking parameters that may describe the minimum viable population, the very concept remains problematic with respect to management application.
Author's response
We agree with the Reviewer that the estimation of MVP for conservation purposes can be problematic and therefore we added some sentences to the Discussion section to point out these shortcomings. Furthermore, we explained how our model can make an important contribution to the ongoing debate about the changes of the extinction risk of populations in presence of different spectral noises.
I am concerned that the HM can empirically incorrect owing to nonequilibrium dynamics in ecological systems. That harmonic mean estimates an extinction threshold assuming two things: that the dynamics driving the mean are at equilibrium, which they are often not, and that we actually know the variance, which we often do not. Both problems are particularly acute in species of conservation concern. To illustrate the latter, see Pimm (book, ~1995), who shows variance increasing through time as populations are further sampled. Sampling to estimate a mean routinely underestimates the frequency of extreme events, and thus underestimates extinction risk. Thus, our harmonic mean will increase through time and additional sampling. The fact that the drivers that affect the mean are assumed to be constant is a dangerous assumption illustrated by climate change. We simply can't assume, for management purposes, that we live in an equilibrium world. Things change, the mean changes with it.
Author's response
We agree that in natural populations equilibrium situations are rather the exception than the rule. However in populations with a long generation time it is possible to assume that the populations have a relatively stable dynamics, at least in the time frame that is interesting from a conservation perspective (which is typically not longer than a century). In the Discussion section we also mentioned how changes of β and K with increasing or decreasing population size could affect the spectrum of noise, suggesting new ideas for future studies. Lastly we also discussed how deviations of β from 2 can invalidate the comparisons of population variability where the coefficient of variation or the logarithm of the standard deviation have been used.
Finally, this model appears to assume that populations are stable and near an extinction threshold. Imagine a population that is clearly on the brink of extinction. Surveys find, 4, 8, 6, 4, 12, 10 breeding females over a series of generations. The arithmetic mean of these numbers is 7.3, the harmonic mean is 6.15. I don't think that these authors, or any responsible conservation biologist, would want to be placed in a position of arguing that a harmonic mean of a population that was already in dire danger of extinction was secure with 7 breeding females. Rather, we would favor an interpretation that if inbreeding depression or demographic stochasticity doesn't doom this population, then catastrophic events (e.g., from hurricanes to habitat loss) surely will.
Author's response
We agree with the Reviewer that with small population size demographic stochasticity is playing a big role for the probability of extinction. Therefore both in the Introduction and the Discussion sections we discussed the concept of demographic stochasticity and we made clear that our model will only be applicable to populations which are large enough, i.e. populations where demographic stochasticity is not a major concern for their persistence. A reference dealing with the debate on the importance of density dependent processes on population dynamics compared to environmental variability has been added.
Reviewer's report 2
Josef Bryja, Department of Population Biology, Institute of Vertebrate Biology AS CR, 675 02 Studenec 122, Czech Republic. Nominated by Aniko Szabo
It seems to me that the paper may contain results that is worth publishing as it presents a novel approach which can be important for many scientific fields; Macroecology, Conservation Genetics and Population Ecology. The paper is overall well written and the maths are simply and clearly expressed. The authors do not explain all mechanisms that impact N_{e}, but I guess that it is because of limited space. If some factors impact both β and N_{e}, we expect things to be more complicated. It does not mean that the authors work is less relevant but, to the opposite, it open new perspectives that should be discussed. I think the introduction is sometimes difficult to follow due to the fact that the authors refer to their results in the introduction. Thank you for the opportunity to read about your inspiring work.
Author's response
We thank the Reviewer for comments and have modified the Introduction according to the suggestions given.
Reviewer's report 3
WaiYuanTan, The University of Memphis Department of Mathematical Sciences. Memphis, TN 38152 United States
The authors assume that the mean and the variance will determine the pattern. This is equivalent to assume normal probability distributions for the random variables. This may be true in some cases but may not be so in other cases. The authors should provide strong evidence to just this. (Simply to quote one paper may not be enough). Aside from the above, I do not have any objection to the paper.
Author's response
We decided to add the following sentence and references to the MS following the Reviewer's comments: "Population size is expected to follow a lognormal distribution (which becomes normally distributed if logtransformed), given that it is the product of temporally multiplicative renewal processes [13–15]. The assumption of the model (a normal or lognormal distribution) is valid both for populations under exponential growth or decline and for populations at equilibrium [13]. It has been demonstrated in a survey of 544 longterm time series of terrestrial and aquatic organisms that about onehalf of them are lognormally distributed, which implicates that our model can be applied also to longterm time series [15]."
Abbreviations
 $\overline{\mu}$ :

Arithmetic mean
 CV:

Coefficient of variation
 r:

Correlation coefficient
 σ^{2}_{d}:

Demographic stochasticity
 N_{E}:

: Effective population size
 σ^{2}_{e}:

Environmental stochasticity
 HM:

Harmonic mean
 log:

Logarithm
 N :

Population size
 std:

Standard deviation
 σ^{2}:

Variance.
Declarations
Acknowledgements
This study has been partly supported by a Marie Curie Transfer of Knowledge Fellowship BIORESC of European Community's Sixth Framework Program (contract number MTKDCT2005029957). Furthermore we wish to thank the ConGen program (funded by the European Science Foundation), the Danish Natural Science Research Council for financial support to CP (grant number: #21010526 and #21030125) and the Nordic council of ministers for supporting Ecoclim, a Nordic center of excellence (LAB).
Authors’ Affiliations
References
 Drake JM, Lodge DM: Effects of environmental variability on extinction and establishment. Ecol Lett. 2004, 7: 2630. 10.1046/j.14610248.2003.00546.x.View ArticleGoogle Scholar
 Lande R: Risks of population extinction from demographic and environmental stochasticity and random catastrophes. Am Nat. 1993, 142: 911927. 10.1086/285580.View ArticleGoogle Scholar
 Sæther BE, Engen S: Routes to extinction. Macroecology. Edited by: Blackburn T, Gaston K. 2003, Blackwell Publishing, Oxford, UK, 218236.Google Scholar
 Henle K, Sarre S, Wiegand K: The role of density regulation in extinction processes and population viability analysis. Biodivers Conserv. 2004, 13: 952. 10.1023/B:BIOC.0000004312.41575.83.View ArticleGoogle Scholar
 Engen S, Bakke Ø, Islam A: Demographic and environmental stochasticity: concepts and definitions. Biometrics. 1998, 54: 840846. 10.2307/2533838.View ArticleGoogle Scholar
 Boyce MS, Haridas CV, Lee CT, NCEAS: Stochastic Demography Working Group: Demography in an increasing variable world. Trends Ecol Evol. 2006, 21: 141148. 10.1016/j.tree.2005.11.018.PubMedView ArticleGoogle Scholar
 Akçakaya HR: Estimating the variance of survival rates and fecundities. Anim Conserv. 2002, 5: 333336.View ArticleGoogle Scholar
 Kohlmann SG, Schmidt GA, Garcelon DK: A population viability analysis for the island fox on Santa Catalina Island, California. Ecol Model. 2005, 183: 7794. 10.1016/j.ecolmodel.2004.07.022.View ArticleGoogle Scholar
 Pertoldi C, Bach LA, Barker JSF, Lundberg P, Loeschcke V: The consequences of the variancemean rescaling effect on effective population size. Oikos. 2007, 116: 769774. 10.1111/j.00301299.2007.15672.x.View ArticleGoogle Scholar
 Hanski I: On patterns of temporal and spatial variation in animal populations. Ann Zool Fenn. 1982, 19: 2138.Google Scholar
 Taylor LR, Woiwod IP: Comparative synoptic dynamics. 1. Relationships between interspecific and interspecific spatial and temporal variance mean population parameters. J Anim Ecol. 1982, 51: 879906. 10.2307/4012.View ArticleGoogle Scholar
 Ballantyne F, Kerkhoff AJ: The observed range for temporal meanvariance scaling exponents can be explained by reproductive correlation. Oikos. 2007, 116: 174180. 10.1111/j.2006.00301299.15383.x.View ArticleGoogle Scholar
 May RM: Patterns of species abundances and diversity. A practical guide to heavy tails: statistical techniques for analysing heavy tailed distributions. Edited by: Cody M, Diamond J. 1975, Birkhauser, 283258.Google Scholar
 Dennis B, Patil GP: Applications in Ecology. Lognormal distributions: theory and applications. Edited by: EL, Shimizu K. 1988, Marcel Dekker, 303330.Google Scholar
 Halley J, Inchausti P: Lognormality in ecological time series. Oikos. 2002, 99: 518530. 10.1034/j.16000706.2002.11962.x.View ArticleGoogle Scholar
 Groom M, Pascual MA: The analysis of population persistence: an outlook on the practice of viability analysis. Edited by: Fiedler PL, Kareiva P. 1998, Chapman and Hall, New York, USA, 427.Google Scholar
 Pimm S, Jones H, Diamond J: On the risk of extinction. Am Nat. 1988, 132: 757785. 10.1086/284889.View ArticleGoogle Scholar
 Pimm S: Life on an intermittent edge. Trends Ecol Evol. 1993, 8: 4546. 10.1016/01695347(93)90156J.PubMedView ArticleGoogle Scholar
 Turchin P: Population regulation: a synthetic view. Oikos. 1999, 84: 153159. 10.2307/3546876.View ArticleGoogle Scholar
 Halley J, Kunin W: Extinction risk and the 1/f family of noise models. Theor Popul Biol. 1999, 56: 215239. 10.1006/tpbi.1999.1424.PubMedView ArticleGoogle Scholar
 Inchausti P, Halley J: On the relation between temporal variability and persistence time in animal populations. J Anim Ecol. 2003, 72: 899908. 10.1046/j.13652656.2003.00767.x.View ArticleGoogle Scholar
 Pertoldi C, Bijlsma R, Loeschcke V: Conservation genetics in a globally changing environment: present problems, paradoxes and future challenges. Biodivers Conserv. 2007, 16: 41474163. 10.1007/s1053100792124.View ArticleGoogle Scholar
 Pertoldi C, Bach LA: Evolutionary aspects of climate induced changes and the need for multidisciplinarity. J Therm Biol. 2007, 32: 118124. 10.1016/j.jtherbio.2007.01.011.View ArticleGoogle Scholar
Copyright
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.