Estimating population dynamics without population data

We develop a biologically correct cost system for production systems facing invasive pests that allows the estimation of population dynamics without a priori knowledge of their true values. We apply that model to a data set for olive producers in Crete and derive from it predictions about the underlying population dynamics. Those dynamics are compared to information on population dynamics obtained from pest sampling with extremely favorable results.     

Estimating population dynamics without population data

We develop a biologically correct cost system for production systems facing invasive pests that allows the estimation of population dynamics without a priori knowledge of their true values. We apply that model to a data set for olive producers in Crete and derive from it predictions about the underlying population dynamics. Those dynamics are compared to information on population dynamics obtained from pest sampling with extremely favorable results.     

Ungulate population dynamics and optimization models

Optimal harvesting strategies for an ungulate population are estimated using stochastic dynamic programming. Data on the Llano Basin white-tailed deer (Odocoileus virginianus) population were used to construct a 2-variable population dynamics model. The model provided the basis for estimating optimal harvesting strategies as a feedback function of the current values of the state variables (prefawning older deer and juveniles). Optimal harvest strategies were insensitive to assumptions about the probability distributions of the stochastic variable (rainfall). The response of the population components to harvesting and the returns obtained from applying optimal strategies were explored through simulation. Mean annual harvest is about 15% of the population. Simplified harvesting strategies based on age-ratios as well as a simplified version based on optimal strategies—but assuming persisting equilibrium juvenile deer density—were compared to optimal strategies through examining values of information. Simplified harvesting strategies lead to a lower harvest over a 50-year simulation period.     

Hidden process models for animal population dynamics.

Hidden process models are a conceptually useful and practical way to simultaneously account for process variation in animal population dynamics and measurement errors in observations and estimates made on the population. Process variation, which can be both demographic and environmental, is modeled by linking a series of stochastic and deterministic subprocesses that characterize processes such as birth, survival, maturation, and movement. Observations of the population can be modeled as functions of true abundance with realistic probability distributions to describe observation or estimation error. Computer-intensive procedures, such as sequential Monte Carlo methods or Markov chain Monte Carlo, condition on the observed data to yield estimates of both the underlying true population abundances and the unknown population dynamics parameters. Formulation and fitting of a hidden process model are demonstrated for Sacramento River winter-run chinook salmon (Oncorhynchus tshawytsha).     

Switching state-space models for modeling penguin population dynamics

Environmental and Ecological Statistics - Tracking individual animals through time using mark-recapture methods is the gold standard for understanding how environmental conditions influence...     

Fitting complex population models by combining particle filters with Markov chain Monte Carlo

We show how a recent framework combining Markov chain Monte Carlo (MCMC) with particle filters (PFMCMC) may be used to estimate population state-space models. With the purpose of utilizing the strengths of each method, PFMCMC explores hidden states by particle filters, while process and observation parameters are estimated using an MCMC algorithm. PFMCMC is exemplified by analyzing time series data on a red kangaroo (Macropus rufus) population in New South Wales, Australia, using MCMC over model parameters based on an adaptive Metropolis-Hastings algorithm. We fit three population models to these data; a density-dependent logistic diffusion model with environmental variance, an unregulated stochastic exponential growth model, and a random-walk model. Bayes factors and posterior model probabilities show that there is little support for density dependence and that the random-walk model is the most parsimonious model. The particle filter Metropolis-Hastings algorithm is a brute-force method that may be used to fit a range of complex population models. Implementation is straightforward and less involved than standard MCMC for many models, and marginal densities for model selection can be obtained with little additional effort. The cost is mainly computational, resulting in long running times that may be improved by parallelizing the algorithm.     

Fitting the truncated negative binomial distribution to count data

Modeling empirical distributions of repeated counts with parametric probability distributions is a frequent problem when studying species abundance. One must choose a family of distributions which is flexible enough to take into account very diverse patterns and possess parameters with clear biological/ecological interpretations. The negative binomial distribution fulfills these criteria and was selected for modeling counts of marine fish and invertebrates. This distribution depends on a vector \(\left( K,\mathfrak {P}\right) \) of parameters, and ranges from the Poisson distribution (when \(K\rightarrow +\infty \)) to Fisher’s log-series, when \(K\rightarrow 0\). Moreover, these parameters have biological/ecological interpretations which are detailed in the literature and in this study. We compared three estimators of K, \(\mathfrak {P}\) and the parameter \(\alpha \) of Fisher’s log-series, following the work of Rao CR (Statistical ecology. Pennsylvania State University Press, University Park, 1971) on a three-parameter unstandardized variant of the negative binomial distribution. We further investigated the coherence underlying parameter values resulting from the different estimators, using both real count data collected in the Mauritanian Exclusive Economic Zone (MEEZ) during the period 1987–2010 and realistic simulations of these data. In the case of the MEEZ, we first built homogeneous lists of counts (replicates), by gathering observations of each species with respect to “typical environments” obtained by clustering the sampled stations. The best estimation of \(\left( K,\mathfrak {P}\right) \) was generally obtained by penalized minimum Hellinger distance estimation. Interestingly, the parameters of most of the correctly sampled species seem compatible with the classical birth-and-dead model of population growth with immigration by Kendall (Biometrika 35:6–15, 1948).     

Modelling crayfish population dynamics using catch data: A size-structured model

In this study we explore a rather unique time series (1979–2002) of catch data of the crayfish Astacus astacus in Lake Steinsfjorden (SE Norway) in combination with temperature data and data on Canadian pondweed Elodea canadensis coverage. In 1977, E. canadensis was for the first time observed in the lake. Over the following years, the plant established dense covers over large parts of the shallow areas, excluding the crayfish from these areas and causing a sudden drop in population size. A size-structured model with bi-stability including a range of observed stage-specific life-cycle attributes (e.g. growth, fecundity, fertility, sex-ratio), population specific parameters and density-dependant (shelters, cannibalism, unspecified predators, competition between individuals, catch, number of traps) as well as density independent factors (temperature and Elodea coverage) were constructed to evaluate the various drivers for the population dynamics, and as a predictive tool for assessing the effects of future changes. Our model revealed that the decline primarily was due to density-dependant effect of the Elodea expansion with reduced number of hides and thus increased risk for predation and cannibalism, but also that temperature played an important role related to recruitment. The model should be relevant for crayfish stock management in general, and by demonstrating the major role of temperature, it is also relevant for predicting population responses under a changing climate. The model should also be applicable to other crustaceans and species with discrete growth and late maturation.     

The scale transition: scaling up population dynamics with field data

Applying the recent developments of scale transition theory, we demonstrate a systematic approach to the problem of scaling up local scale interactions to regional scale dynamics with field data. Dynamics on larger spatial scales differ from the predictions of local dynamics alone because of an interaction between nonlinearity in population dynamics at the local scale and spatial variation in density and environmental factors over the regional population. Our systematic approach to scaling up involves the following five steps. First, define a model for dynamics on the local spatial scale. Second, apply scale transition theory to identify key interactions between nonlinearity and spatial variation that translate local dynamics to the regional scale. Third, measure local-scale model parameters to determine nonlinearities at local scales. Fourth, measure spatial variation. Finally, combine nonlinearity and variation measures to obtain the scale transition. Using field data for the dynamics of grazers and periphyton in a freshwater stream, we show that scale transition terms greatly reduce the growth and equilibrium density of the periphyton population at the stream scale compared to rock scale populations, confirming the importance of spatial mechanisms to stream-scale dynamics.     

Ecotoxicology and population dynamics: Using DEBtox models in a Leslie modeling approach

Although the ecological risks of toxic chemicals are usually assessed on the basis of individual responses, such as survival, reproduction or growth, ecotoxicologists are now attempting to assess the impact of environmental pollution on the dynamics of naturally exposed populations. The main issue is how to infer the likely impact on the population of the toxic effects observed at the individual level. Dynamic energy budget in toxicology (DEBtox) is the most user-friendly software currently available to analyze the experimental data obtained in toxicity tests performed on individuals. Because toxic effects are diverse and because the sensitivity of individuals varies considerably depending on life-cycle stage, Leslie models offer a convenient way of predicting toxicant effects on population dynamics.In the present study, we first show how parameter inputs, estimated from individual data using DEBtox, can be coupled using a Leslie matrix population model. Then, using experimental data obtained with Chironomus riparius, we show how the effects of a pesticide (methiocarb) on the population growth rate of a laboratory population can be estimated. Lastly, we perform a complex sensitivity analysis to pinpoint critical age classes within the population for the purposes of the field management of populations.     

A new probability model with application to heavy-tailed hydrological data

Environmental and Ecological Statistics - Because of the dramatic changes that are being observed in the climatic conditions of the world, such as excess of rains, drought and huge floods, we...     

Fitting a parabola to growth data of fishes and some applications to fisheries

Denoting a fish length or weight at age t by X _t , a reference age by t _m , and the corresponding fish length or weight by X _m , the relation between age and length or weight may be described by a parabola as follows: $$\left| {X_t } \right. - X_m \left| = \right.a + b(\left| {t - t_m } \right.\left| ) \right. + c(\left| t \right. - t_m \left| ) \right.^2$$ or $$X_t = A + b(\left| {t - t_m } \right.\left| ) \right. + c(\left| t \right. - t_m \left| ) \right.^2$$ where a, b and c are constants. Each of the above Eqs. describes one curve at ages older than t _m and another one at younger ages, which is made possible by means of the transformation of t to (|t-t _m |). In 2 cases out of 10, the parabola takes the form of a cubic equation. Evidence is given that, as the growth data become fewer, the better fit of the parabola or cubic equation will probably be less in comparison to the von Bertalanffy equation (1938, 1949) as developed by Beverton and Holt (1957) and the power-growth equation (Rafail, 1971), and vice versa. This growth equation is used to derive models for estimating the optimum age and yield for fish populations.     

Assessing the effects of disease and bleaching on Florida Keys corals by fitting population models to data

Coral diseases have increased in frequency over the past few decades and have important influences on the structure and composition of coral reef communities. However, there is limited information on the etiologies of many coral diseases, and pathways through which coral diseases are acquired and transmitted are still in question. Furthermore, it is difficult to assess the impacts of disease on coral populations because outbreaks often co-occur with temperature-induced bleaching and anthropogenic stressors. We developed spatially explicit population models of coral disease and bleaching dynamics to quantify the impact of six common diseases on Florida Keys corals, including aspergillosis, dark spots, white band, white plague, white patch, and Caribbean yellow band. Models were fit to an 8-year data set of coral abundance, disease prevalence, and bleaching prevalence. Model selection was used to assess alternative pathways for disease transmission, and the influence of environmental stressors, including sea temperature and human population density, on disease prevalence and coral mortality. Classic disease transmission from contagious to susceptible colonies provided the best-fit model only for aspergillosis. For other diseases, external disease forcing, such as through a vector or directly from pathogens in the environment, provided the best fit to observed data. Estimates of disease reproductive ratio values (R₀) were less than one for each disease, indicating coral colonies were below densities required for diseases to become established through contagious spread alone. Incidences of white band and white patch disease were associated with greater susceptibility or slower recovery of bleached colonies, and no disease outbreaks were associated with periods of elevated sea temperatures alone. Projections of best-fit models indicated that, atleast during the period of this study, disease and bleaching did not have substantial impacts on populations and impaired rates of population growth appeared to be attributable to other stressors. By applying epidemiological models to field data, our study gives qualitative insights into the dynamics of coral diseases, relative stressor impacts, and directions for future research.     

Fitting animal survival models with temporal random effects

Estimating temporal variance in animal demographic parameters is of particular importance in population biology. We implement the Schall’s algorithm for incorporating temporal random effects in survival models using recovery data. Our frequentist approach is based on a formulation of band-recovery models with random effects as generalized linear mixed models and a linearization of the link function conditional on the random effects. A simulation study shows that our procedure provides unbiased and precise estimates. The method is then implemented on two case studies using recovery data on fish and birds.     

Fitting fertilisation kinetics models for free-spawning marine invertebrates

 To determine how fertilisation varied with sperm concentration for two species of scallop, Chlamys (Equichlamys) bifrons (Lamarck) and C. asperrima (Lamarck), we performed a simple series of sperm dilution experiments, and measured egg size and sperm swimming speeds. C. bifrons eggs were much larger (average diam=116.5 μm), and sperm swimming speeds faster (209.8 μm s⁻¹), than C. asperrima (71.2 μm, 166.0 μm s⁻¹). In both species, maximum fertilisation occurred at an ambient sperm concentration of around 100 sperm μl⁻¹; the maximum proportion of eggs fertilised was less than 0.70 in the C. bifrons experiments, but nearer 1.0 with C. asperrima. At high sperm concentrations (>100 sperm μl⁻¹), fertilisation decreased (presumably due to polyspermy) with increasing sperm concentration, but decreased more rapidly in C. bifrons than C. asperrima. A polyspermy-adjusted fertilisation kinetics model could be fitted to the experimental data, but unique parameter estimates could not be determined. Received: 7 October 1999 / Accepted: 8 July 2000     

Transient population dynamics in periodic matrix models: methodology and effects of cyclic permutations

Many biological populations are subject to periodically changing environments such as years with or without fire, or rotation of crop types. The dynamics and management options for such populations are frequently investigated using periodic matrix models. However the analysis is usually limited to long-term results (asymptotic population growth rate and its sensitivity to perturbations of vital rates). In non-periodic matrix models it has been shown that long-term results may be misleading as populations are rarely in their stable structure. We therefore develop methods to analyze transient dynamics of periodic matrix models. In particular, we show how to calculate the effects of perturbations on population size within and at the end of environmental cycles. Using a model of a weed population subject to a crop rotation, we show that different cyclic permutations produce different patterns of sensitivity of population size and different population sizes. By examining how the starting environment interacts with the initial conditions, we explain how different patterns arise. Such understanding is critical to developing effective management and monitoring strategies for populations subject to periodically recurring environments.     

On the latent state estimation of nonlinear population dynamics using Bayesian and non-Bayesian state-space models

Nonlinear state-space models have been increasingly applied to study population dynamics and data assimilation in environmental sciences. State-space models can account for process error and measurement error simultaneously to correct for the bias in the estimates of system state and model parameters. However, few studies have compared the performance of different nonlinear state-space models for reconstructing the state of population dynamics from noisy time series. This study compared the performance of the extended Kalman filter (EKF), unscented Kalman filter (UKF) and Bayesian nonlinear state-space models (BNSSM) through simulations. Synthetic population time series were generated using the theta logistic model with known parameters, and normally distributed process and measurement errors were introduced using the Monte Carlo simulations. At higher levels of nonlinearity, the UKF and BNSSM had lower root mean square error (RMSE) than the EKF. The BNSSM performed reliably across all levels of nonlinearity, whereas increased levels of nonlinearity resulted in higher RMSE of the EKF. The Metropolis–Hastings algorithm within the Gibbs algorithm was used to fit the theta logistic model to synthetic time series to estimate model parameters. The estimated posterior distribution of the parameter θ indicated that the 95% credible intervals included the true values of θ (=0.5 and 1.5), but did not include 1.0 and 0.0. Future studies need to incorporate the adaptive Metropolis algorithm to estimate unknown model parameters for broad applications of Bayesian nonlinear state-space models in ecological studies.     

Zero-inflated models with application to spatial count data

Count data arises in many contexts. Here our concern is with spatial count data which exhibit an excessive number of zeros. Using the class of zero-inflated count models provides a flexible way to address this problem. Available covariate information suggests formulation of such modeling within a regression framework. We employ zero-inflated Poisson regression models. Spatial association is introduced through suitable random effects yielding a hierarchical model. We propose fitting this model within a Bayesian framework considering issues of posterior propriety, informative prior specification and well-behaved simulation based model fitting. Finally, we illustrate the model fitting with a data set involving counts of isopod nest burrows for 1649 pixels over a portion of the Negev desert in Israel.     

Non-linear biological responses to disturbance: consequences on population dynamics

We assessed how non-linear biological responses to environmental noise, or “noise filtering”, impact the spectra of density-dependent population dynamics, and the correlation between noise and population dynamics. The noise was assumed to affect population growth rate in a discrete-time population model by Hassell [J. Anim. Ecol. 44 (1975) 283–295] where the population growth rate was linked to the environment with an optimum type filter. When compared to unfiltered noise, the filtered noise can distort the stationary distribution of population values. The optimum type filter can make cyclic population dynamics more regular and low population values can become more frequent or rare depending on the strength of density dependence. Filtering can cause blue shifted and red shifted population dynamics and determine the strength of correlation between environmental noise and population size. In most cases, optimum type filtering makes linear correlation between population dynamics and noise weaker. The filter effect on population spectra and noise versus population correlation is sensitive to changes in population model parameters, the location where noise hits the filter, and noise colour.     

The pattern of species turnover resulting from stochastic population dynamics: The model and field data

The model of random population dynamics in a sampling site returns geometric distribution of longevities of continuous presence (=persistence) and Poisson distribution of the presence–absence transitions. This discrete-time stochastic process describes the presence–absence pattern observed in the beetles surveyed 6 years on Mount Carmel, Israel. Homogeneous pools of species mostly on the Families rank, exhibit the predicted by the model patterns. Conformity to an ergodic hypothesis is the criterion of ecological homogeneity. This criterion assumes the equivalence of short-term behavior of entire pool and long-term behavior of any species from this pool. The pool of all 801 species of Order Coleoptera does not match the model. Thus a taxon of an arbitrary rank may not be considered a priory as a unit of ecological study. Determined from field data parameters of the model are biased and magnitude of the bias depends on longevity of the survey. Parameter of distribution depends also on species tolerance, which is the level adaptation of given species to given environment in given time interval. Random process of species turnover may be considered as a game of species to gain their presence against deteriorative fluctuations of environmental conditions.