Likelihood based population viability analysis in the presence of observation error

Population viability analysis (PVA) entails calculation of extinction risk, as defined by various extinction metrics, for a study population. These calculations strongly depend on the form of the population growth model and inclusion of demographic and/or environmental stochasticity. Form of the model and its parameters are determined based on observed population time series data. A typical population time series, consisting of estimated population sizes, inevitably has some observation error and likely has missing observations. In this paper, we present a likelihood based PVA in the presence of observation error and missing data. We illustrate the importance of incorporation of observation error in PVA by reanalyzing the population time series of song sparrow Melospiza melodia on Mandarte Island, British Columbia, Canada from 1975-1998. Using Akaike information criterion we show that model with observation error fits the data better than the one without observation error. The extinction risks predicted by with and without observation error models are quite different. Further analysis of possible causes for observation error revealed that some component of the observation error might be due to unreported dispersal. A complete analysis of such data, thus, would require explicit spatial models and data on dispersal along with observation error. Our conclusions are, therefore, two-fold: 1) observation errors in PVA matter and 2) integrating these errors in PVA is not always enough and can still lead to important biases in parameter estimates if other processes such as dispersal are ignored.

[1]  Brian Dennis,et al.  Data cloning: easy maximum likelihood estimation for complex ecological models using Bayesian Markov chain Monte Carlo methods. , 2007, Ecology letters.

[2]  Ray Hilborn,et al.  The Influence of Model Structure on Conclusions about the Viability and Harvesting of Serengeti Wildebeest , 1997 .

[3]  J N Smith,et al.  Estimating the time to extinction in an island population of song sparrows , 2000, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[4]  Michael A. Larson,et al.  Incorporating parametric uncertainty into population viability analysis models , 2011 .

[5]  Henrik Madsen,et al.  Estimation methods for nonlinear state-space models in ecology , 2011 .

[6]  M. Shaffer.  Minimum Population Sizes for Species Conservation , 1981 .

[7]  G. Grigg,et al.  Of sheep and rain: large-scale population dynamics of the red kangaroo , 2005 .

[8]  W. Hochachka,et al.  A metapopulation approach to the population biology of the Song Sparrow Melospiza melodia , 2008 .

[9]  M. Gilpin,et al.  Global models of growth and competition. , 1973, Proceedings of the National Academy of Sciences of the United States of America.

[10]  M. Taper,et al.  Environmental Variation and the Persistence of Small Populations. , 1992, Ecological applications : a publication of the Ecological Society of America.

[11]  J. Maynard Smith,et al.  The Stability of Predator‐Prey Systems , 1973 .

[12]  B. Taylor,et al.  The Reliability of Using Population Viability Analysis for Risk Classification of Species , 1995 .

[13]  S. Lele,et al.  ESTIMATING DENSITY DEPENDENCE, PROCESS NOISE, AND OBSERVATION ERROR , 2006 .

[14]  Brian Dennis,et al.  Estimation of Growth and Extinction Parameters for Endangered Species , 1991 .

[15]  H. Akçakaya,et al.  Assessing human impact despite uncertainty:viability of the northern spotted owl metapopulation in the northwestern USA , 1998, Biodiversity and Conservation.

[16]  Brian Dennis,et al.  ESTIMATING POPULATION TREND AND PROCESS VARIATION FOR PVA IN THE PRESENCE OF SAMPLING ERROR , 2004 .

[17]  Ian R. Harris.  Predictive fit for natural exponential families , 1989 .

[18]  Brian Dennis,et al.  Hierarchical models in ecology: confidence intervals, hypothesis testing, and model selection using data cloning. , 2009, Ecology.

[19]  G. Kitagawa.  Non-Gaussian State—Space Modeling of Nonstationary Time Series , 1987 .

[20]  P. Arcese.  How Fit are Floaters? Consequences of Alternative Territorial Behaviors in a Nonmigratory Sparrow , 1989, The American Naturalist.

[21]  Robert P Freckleton,et al.  Census error and the detection of density dependence. , 2006, The Journal of animal ecology.

[22]  Benjamin Gompertz,et al.  XXIV. On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies. In a letter to Francis Baily, Esq. F. R. S. &c , 1825, Philosophical Transactions of the Royal Society of London.

[23]  Jon T. Schnute,et al.  A General Framework for Developing Sequential Fisheries Models , 1994 .

[24]  B. Danielson,et al.  Spatially Explicit Population Models: Current Forms and Future Uses , 1995 .

[25]  Perry de Valpine,et al.  Review of methods for fitting time-series models with process and observation error and likelihood calculations for nonlinear, non-Gaussian state-space models , 2002 .

[26]  Øyvind Bakke,et al.  Demographic and Environmental Stochasticity Concepts and Definitions , 1998 .

[27]  Christian Wissel,et al.  The intrinsic mean time to extinction: a unifying approach to analysing persistence and viability of populations , 2004 .

[28]  James S. Clark,et al.  POPULATION TIME SERIES: PROCESS VARIABILITY, OBSERVATION ERRORS, MISSING VALUES, LAGS, AND HIDDEN STATES , 2004 .

[29]  James D. Hamilton.  A standard error for the estimated state vector of a state-space model , 1986 .

[30]  Subhash R. Lele,et al.  Estimability and Likelihood Inference for Generalized Linear Mixed Models Using Data Cloning , 2010 .

[31]  Subhash R. Lele,et al.  STATISTICAL ANALYSIS OF POPULATION DYNAMICS INSPACE AND TIME USING ESTIMATING FUNCTIONS , 1998 .

[32]  Leo Polansky,et al.  Likelihood ridges and multimodality in population growth rate models. , 2009, Ecology.

[33]  Subhash R Lele,et al.  Sampling variability and estimates of density dependence: a composite-likelihood approach. , 2006, Ecology.

[34]  Brian Dennis,et al.  DENSITY DEPENDENCE IN TIME SERIES OBSERVATIONS OF NATURAL POPULATIONS: ESTIMATION AND TESTING' , 1994 .

[35]  Robert M. May,et al.  Patterns of Dynamical Behaviour in Single-Species Populations , 1976 .

[36]  M. Taper,et al.  JOINT DENSITY DEPENDENCE , 1998 .

[37]  K. Wiegand,et al.  The Role of Density Regulation in Extinction Processes and Population Viability Analysis , 2004, Biodiversity & Conservation.

[38]  Donald Ludwig,et al.  Stability, Regulation, and the Determination of Abundance in an Insular Song Sparrow Population , 1992 .

[39]  Péter Sólymos,et al.  dclone: Data Cloning in R , 2010, R J..

[40]  S. T. Buckland,et al.  Hidden process models for animal population dynamics. , 2006, Ecological applications : a publication of the Ecological Society of America.

[41]  Nicolas Schtickzelle,et al.  Metapopulation viability analysis of the bog fritillary butterfly using RAMAS/GIS , 2004 .

[42]  W. Ricker.  Stock and Recruitment , 1954 .

[43]  David R. Anderson,et al.  Multimodel Inference , 2004 .

[44]  R. Sibly,et al.  The effects of environmental perturbation and measurement error on estimates of the shape parameter in the theta-logistic model of population regulation , 2008 .

[45]  Steinar Engen,et al.  Predicting fluctuations of reintroduced ibex populations: the importance of density dependence, environmental stochasticity and uncertain population estimates. , 2007, The Journal of animal ecology.

[46]  D. Ludwig.  Is it meaningful to estimate a probability of extinction , 1999 .

[47]  Alan Hastings,et al.  FITTING POPULATION MODELS INCORPORATING PROCESS NOISE AND OBSERVATION ERROR , 2002 .

[48]  Brian Dennis,et al.  Joint effects of density dependence and rainfall on abundance of San Joaquin kit fox , 2000 .

[49]  Robin McCleery,et al.  Environmental Stochasticity and Extinction Risk in a Population of a Small Songbird, the Great Tit , 1998, The American Naturalist.