A framework for studying transient dynamics of population projection matrix models

Empirical models are central to effective conservation and population management, and should be predictive of real-world dynamics. Available modelling methods are diverse, but analysis usually focuses on long-term dynamics that are unable to describe the complicated short-term time series that can arise even from simple models following ecological disturbances or perturbations. Recent interest in such transient dynamics has led to diverse methodologies for their quantification in density-independent, time-invariant population projection matrix (PPM) models, but the fragmented nature of this literature has stifled the widespread analysis of transients. We review the literature on transient analyses of linear PPM models and synthesise a coherent framework. We promote the use of standardised indices, and categorise indices according to their focus on either convergence times or transient population density, and on either transient bounds or case-specific transient dynamics. We use a large database of empirical PPM models to explore relationships between indices of transient dynamics. This analysis promotes the use of population inertia as a simple, versatile and informative predictor of transient population density, but criticises the utility of established indices of convergence times. Our findings should guide further development of analyses of transient population dynamics using PPMs or other empirical modelling techniques.     

Integral projection models perform better for small demographic data sets than matrix population models: a case study of two perennial herbs

1. Matrix population models are widely used to describe population dynamics, conduct population viability analyses and derive management recommendations for plant populations. For endangered or invasive species, management decisions are often based on small demographic data sets. Hence, there is a need for population models which accurately assess population performance from such small data sets.
2. We used demographic data on two perennial herbs with different life histories to compare the accuracy and precision of the traditional matrix population model and the recently developed integral projection model (IPM) in relation to the amount of data.
3. For large data sets both matrix models and IPMs produced identical estimates of population growth rate (λ). However, for small data sets containing fewer than 300 individuals, IPMs often produced smaller bias and variance for λ than matrix models despite different matrix structures and sampling techniques used to construct the matrix population models.
4. Synthesis and applications . Our results suggest that the smaller bias and variance of λ estimates make IPMs preferable to matrix population models for small demographic data sets with a few hundred individuals. These results are likely to be applicable to a wide range of herbaceous, perennial plant species where demographic fate can be modelled as a function of a continuous state variable such as size. We recommend the use of IPMs to assess population performance and management strategies particularly for endangered or invasive perennial herbs where little demographic data are available.     

Using statistics to design and estimate vital rates in matrix population models for a perennial herb

Matrix population models are widely used to assess population status and to inform management decisions. Despite existing theories for building such models, model construction is often partially based on expert opinion. So far, model structure has received relatively little attention, although it may affect estimates of population dynamics. Here, we assessed the consequences of two published matrix structures (a 4 × 4 matrix based on expert opinion and a 10 × 10 matrix based on statistical modeling) for estimates of vital rates and stochastic population dynamics of the long-lived herb Astragalus scaphoides. We explored the ways in which choice of model structure alters the accuracy (i.e., mean) and precision (i.e., variance) of predicted population dynamics. We found that model structure had a negligible effect on the accuracy and precision of vital rates and stochastic stage distribution. However, the 10 × 10 matrix produced lower estimates of stochastic population growth rates than the 4 × 4 matrix, and more accurately predicted the observed trends in population abundance for three out of four study populations. Moreover, estimates of realized variation in population growth rate due to fluctuations in population stage structure over time were occasionally sensitive to matrix structure, suggesting differential roles of transient dynamics. Our study indicates that statistical modeling for choosing categories in matrix models might be preferable over expert opinion to accurately predict population trends and can provide a more objective way for model construction when the biological knowledge of the species is limited.     

Demography_Lab,an educational application to evaluate population growth: Unstructured and matrix models

Training in Population Ecology asks for scalable applications capable of embarking students on a trip from basic concepts to the projection of populations under the various effects of density dependence and stochasticity. Demography_Lab is an educational tool for teaching Population Ecology aspiring to cover such a wide range of objectives. The application uses stochastic models to evaluate the future of populations. Demography_Lab may accommodate a wide range of life cycles and can construct models for populations with and without an age or stage structure. Difference equations are used for unstructured populations and matrix models for structured populations. Both types of models operate in discrete time. Models can be very simple, constructed with very limited demographic information or parameter‐rich, with a complex density‐dependence structure and detailed effects of the different sources of stochasticity. Demography_Lab allows for deterministic projections, asymptotic analysis, the extraction of confidence intervals for demographic parameters, and stochastic projections. Stochastic population growth is evaluated using up to three sources of stochasticity: environmental and demographic stochasticity and sampling error in obtaining the projection matrix. The user has full control on the effect of stochasticity on vital rates. The effect of the three sources of stochasticity may be evaluated independently for each vital rate. The user has also full control on density dependence. It may be included as a ceiling population size controlling the number of individuals in the population or it may be evaluated independently for each vital rate. Sensitivity analysis can be done for the asymptotic population growth rate or for the probability of extinction. Elasticity of the probability of extinction may be evaluated in response to changes in vital rates, and in response to changes in the intensity of density dependence and environmental stochasticity.     

Dynamics of age- and size-structured populations in fluctuating environments: Applications of stochastic matrix models to natural populations

Recent developments of the theory of stochastic matrix modeling have made it possible to estimate general properties of age- and size-structured populations in fluctuating environments. However, applications of the theory to natural populations are still few. The empirical studies which have used stochastic matrix models are reviewed here to examine whether predictions made by the theory can be generally found in wild populations. The organisms studied include terrestrial grasses and herbs, a seaweed, a fish, a reptile, a deer and some marine invertebrates. In all the studies, the stochastic population growth rate (ln λ _s ) was no greater than the deterministic population growth rate determined using average vital rates, suggesting that the model based only on average vital rates may overestimate growth rates of populations in fluctuating environments. Factors affecting ln λ _s include the magnitude of variation in vital rates, probability distribution of random environments, fluctuation in different types of vital rates, covariances between vital rates, and autocorrelation between successive environments. However, comprehensive rules were hardly found through the comparisons of the empirical studies. Based on shortcomings of previous studies, I address some important subjects which should be examined in future studies.     

Evaluation of the population dynamics of the forage legume (Lotus corniculatus), using matrix population models

The population dynamics of perennial crop plants are influenced by numerous factors, including management practices. Conditions in the field vary from year to year, and matrix population models are useful for evaluating population behaviour in relation to environmental variability. In Missouri, the stand persistence of birdsfoot trefoil ( Lotus corniculatus ), a perennial legume, is often limited by disease and poor seed production. A stage-based, matrix population model was developed to evaluate the population dynamics of birdsfoot trefoil in relation to clipping treatment. The plant growth stages represented in the model were seeds, seedlings, mature vegetative and reproductive plants. Two phases of population growth were evaluated in clipped and unclipped stands. Establishment-phase populations were characterized by relatively high mortality and low reproduction. Elasticity analysis indicated that growth of these populations was most sensitive to the survival of vegetative plants. Mature vegetative plants and seeds comprised the majority of surviving individuals in clipped and unclipped populations, respectively; however, establishment-phase populations under both management treatments tended toward extinction. Populations in the post-establishment phase of growth were characterized by relatively low mortality and high reproduction. Population growth in this phase of growth was most sensitive to seed production, and most individuals in these populations were at the seed stage.     

On the reproductive value and the spectrum of a population projection matrix with implications for dynamic population models

The eigenvalues of a population projection matrix-except for the Lotka coefficient-are uniquely determined by the reproductive values and the survival. This relation (proposed earlier, but not really well known in western literature) follows from another useful relation between fertility, reproductive values, survival, and Lotka’s coefficient. These results are applied to provide demographic interpretations to the intrinsically dynamic and metastable population models by Schoen and co-workers.     

一种模拟昆虫种群动态的改进的变维矩阵模型

提出了一种模拟昆虫种群动态的改进的变维矩阵模型,该模型以发有历期为维数,采用分解与合成的方法变维,并考虑了个体间的发育差异。经模拟检验,模型模拟结果略优于徐汝梅等（1981）变维矩阵模型的结果。     

Open removal models with temporary emigration and population dynamics to inform invasive animal management

Removal sampling data are the primary source of monitoring information for many populations (e.g., invasive species, fisheries). Population dynamics, temporary emigration, and imperfect detection are common sources of variation in monitoring data and are key parameters for informing management. We developed two open robust‐design removal models for simultaneously modeling population dynamics, temporary emigration, and imperfect detection: a random walk linear trend model (estimable without ancillary information), and a 2‐age class informed population model (InfoPM, closely related to integrated population models) that incorporated prior information for age‐structured vital rates and relative juvenile availability. We applied both models to multiyear, removal trapping time‐series of a large invasive lizard (Argentine black and white tegu, Salvator merianae) in three management areas of South Florida to evaluate the effectiveness of management programs. Although estimates of the two models were similar, the InfoPMs generally returned more precise estimates, partitioned dynamics into births, deaths, net migration, and provided a decision support tool to predict population dynamics under different effort scenarios while accounting for uncertainty. Trends in tegu superpopulation abundance estimates were increasing in two management areas despite generally high removal rates. However, tegu abundance appeared to decline in the Core management area, where trapping density was the highest and immigration the lowest. Finally, comparing abundance predictions of no‐removal scenarios to those estimated in each management area suggested significant population reductions due to management. These results suggest that local tegu population control via systematic trapping may be feasible with high enough trap density and limited immigration; and highlights the value of these trapping programs. We provided the first estimates of tegu abundance, capture probabilities, and population dynamics, which is critical for effective management. Furthermore, our models are applicable to a wide range of monitoring programs (e.g., carcass recovery or removal point‐counts).     

Relation of transition probabilities in the Leslie model to those from experimental cumulative distributions

We explore the relationship between transition probabilities in the Leslie model and those derived from experimental cumulative distributions. The nature of the two kinds of probabilities are discussed, and a formula derived for converting from one to the other. A numerical example is given to illustrate the differences.     

Forecasting species range dynamics with process‐explicit models: matching methods to applications

Knowing where species occur is fundamental to many ecological and environmental applications. Species distribution models (SDMs) are typically based on correlations between species occurrence data and environmental predictors, with ecological processes captured only implicitly. However, there is a growing interest in approaches that explicitly model processes such as physiology, dispersal, demography and biotic interactions. These models are believed to offer more robust predictions, particularly when extrapolating to novel conditions. Many process–explicit approaches are now available, but it is not clear how we can best draw on this expanded modelling toolbox to address ecological problems and inform management decisions. Here, we review a range of process–explicit models to determine their strengths and limitations, as well as their current use. Focusing on four common applications of SDMs – regulatory planning, extinction risk, climate refugia and invasive species – we then explore which models best meet management needs. We identify barriers to more widespread and effective use of process‐explicit models and outline how these might be overcome. As well as technical and data challenges, there is a pressing need for more thorough evaluation of model predictions to guide investment in method development and ensure the promise of these new approaches is fully realised.     

The reliability of approximate reduction techniques in population models with two time scales

As a result of the complexity inherent in some natural systems, mathematical models employed in ecology are often governed by a large number of variables. For instance, in the study of population dynamics we often find multiregional models for structured populations in which individuals are classified regarding their age and their spatial location. Dealing with such structured populations leads to high dimensional models. Moreover, in many instances the dynamics of the system is controlled by processes whose time scales are very different from each other. For example, in multiregional models migration is often a fast process in comparison to the growth of the population.Approximate reduction techniques take advantage of the presence of different time scales in a system to introduce approximations that allow one to transform the original system into a simpler low dimensional system. In this way, the dynamics of the original system can be approximated in terms of that of the reduced system. This work deals with the study of that approximation. In particular, we work with a non-autonomous discrete time model previously presented in the literature and obtain different bounds for the error we incur when we describe the dynamics of the original system in terms of the reduced one.The results are illustrated by some numerical simulations corresponding to the reduction of a Leslie type model for a population structured in two age classes and living in a two patch system.     

Improved viability of populations with diverse life-history portfolios

A principle shared by both economists and ecologists is that a diversified portfolio spreads risk, but this idea has little empirical support in the field of population biology. We found that population growth rates (recruits per spawner) and life-history diversity as measured by variation in freshwater and ocean residency were negatively correlated across short time periods (one to two generations), but positively correlated at longer time periods, in nine Bristol Bay sockeye salmon populations. Further, the relationship between variation in growth rate and life-history diversity was consistently negative. These findings strongly suggest that life-history diversity can both increase production and buffer population fluctuations, particularly over long time periods. Our findings provide new insights into the importance of biocomplexity beyond spatio-temporal aspects of populations, and suggest that maintaining diverse life-history portfolios of populations may be crucial for their resilience to unfavourable conditions like habitat loss and climate change.     

Linear discrete models with different time scales

Aggregation of variables allows to approximate a large scale dynamical system (the micro-system) involving many variables into a reduced system (the macro-system) described by a few number of global variables. Approximate aggregation can be performed when different time scales are involved in the dynamics of the micro-system. Perturbation methods enable to approximate the large micro-system by a macro-system going on at a slow time scale. Aggregation has been performed for systems of ordinary differential equations in which time is a continuous variable. In this contribution, we extend aggregation methods to time-discrete models of population dynamics. Time discrete micro-models with two time scales are presented. We use perturbation methods to obtain a slow macro-model. The asymptotic behaviours of the micro and macro-systems are characterized by the main eigenvalues and the associated eigenvectors. We compare the asymptotic behaviours of both systems which are shown to be similar to a certain order.     

Bayesian experimental design for nonlinear mixed-effects models with application to HIV dynamics

Bayesian experimental design is investigated for Bayesian analysis of nonlinear mixed-effects models. Existence of the posterior risk for parameter estimation is shown. When the same prior distribution is used for both design and inference, existence of the preposterior risk for design is also proven. If the prior distribution used in design is different from that used for inference, sufficient conditions are established for existence of the preposterior risk for design. A case study of design for an experiment in population HIV dynamics is provided.     

Predicting Aggregation Kinetics of DU 145 Prostate Cancer Cells in Liquid-Overlay Culture

The predictive capacity of a novel population-balance model to simulate aggregation kinetics of attachment-dependent cells at the resolution of one-cell increments has been evaluated. Using spheroid assembly of DU 145 human prostate cancer cells as a representative system, the mathematical model proved to be robust in simulating aggregation over a 5-fold range of surface densities from 5×10³ to 2.5×10⁴ cells/cm² with a single matrix of rate constants. For cultures at 1×10⁵ cells/cm², more than 75% of simulated aggregate concentrations are within the standard deviation of measured concentrations. For the two extreme densities, at least two-thirds of model predictions are within 35% of the mean for experimental data. Error in model predictions is attributed to uncertainty in measurements and intrinsic changes in aggregation. The model has application to the rational design of spheroids in tissue engineering and bioseparation processes in pharmaceutical manufacturing.     

Frequency‐dependent two‐sex models: a new approach to sex ratio evolution with multiple maternal conditions

Mothers that experience different individual or environmental conditions may produce different proportions of male to female offspring. The Trivers‐Willard hypothesis, for instance, suggests that mothers with different qualities (size, health, etc.) will use different sex ratios if maternal quality differentially affects sex‐specific reproductive success. Condition‐dependent, or facultative, sex ratio strategies like these allow multiple sex ratios to coexist within a population. They also create complex population structure due to the presence of multiple maternal conditions. As a result, modeling facultative sex ratio evolution requires not only sex ratio strategies with multiple components, but also two‐sex population models with explicit stage structure. To this end, we combine nonlinear, frequency‐dependent matrix models and multidimensional adaptive dynamics to create a new framework for studying sex ratio evolution. We illustrate the applications of this framework with two case studies where the sex ratios depend one of two possible maternal conditions (age or quality). In these cases, we identify evolutionarily singular sex ratio strategies, find instances where one maternal condition produces exclusively male or female offspring, and show that sex ratio biases depend on the relative reproductive value ratios for each sex.     

Loss of animal seed dispersal increases extinction risk in a tropical tree species due to pervasive negative density dependence across life stages

Overhunting in tropical forests reduces populations of vertebrate seed dispersers. If reduced seed dispersal has a negative impact on tree population viability, overhunting could lead to altered forest structure and dynamics, including decreased biodiversity. However, empirical data showing decreased animal-dispersed tree abundance in overhunted forests contradict demographic models which predict minimal sensitivity of tree population growth rate to early life stages. One resolution to this discrepancy is that seed dispersal determines spatial aggregation, which could have demographic consequences for all life stages. We tested the impact of dispersal loss on population viability of a tropical tree species, Miliusa horsfieldii, currently dispersed by an intact community of large mammals in a Thai forest. We evaluated the effect of spatial aggregation for all tree life stages, from seeds to adult trees, and constructed simulation models to compare population viability with and without animal-mediated seed dispersal. In simulated populations, disperser loss increased spatial aggregation by fourfold, leading to increased negative density dependence across the life cycle and a 10-fold increase in the probability of extinction. Given that the majority of tree species in tropical forests are animal-dispersed, overhunting will potentially result in forests that are fundamentally different from those existing now.