Approximate reduction of linear population models governed by stochastic differential equations: application to multiregional models

In this work we develop approximate aggregation techniques in the context of slow-fast linear population models governed by stochastic differential equations and apply the results to the treatment of populations with spatial heterogeneity. Approximate aggregation techniques allow one to transform a complex system involving many coupled variables and in which there are processes with different time scales, by a simpler reduced model with a fewer number of ‘global’ variables, in such a way that the dynamics of the former can be approximated by that of the latter. In our model we contemplate a linear fast deterministic process together with a linear slow process in which the parameters are affected by additive noise, and give conditions for the solutions corresponding to positive initial conditions to remain positive for all times. By letting the fast process reach equilibrium we build a reduced system with a lesser number of variables, and provide results relating the asymptotic behaviour of the first- and second-order moments of the population vector for the original and the reduced system. The general technique is illustrated by analysing a multiregional stochastic system in which dispersal is deterministic and the rate growth of the populations in each patch is affected by additive noise.     

Extinction in relation to demographic and environmental stochasticity in age-structured models

The demographic variance of an age-structured population is defined. This parameter is further split into components generated by demographic stochasticity in each vital rate. The applicability of these parameters are investigated by checking how an age-structured population process can be approximated by a diffusion with only three parameters. These are the deterministic growth rate computed from the expected projection matrix and the environmental and demographic variances. We also consider age-structured populations where the fecundity at any stage is either zero or one, and there is neither environmental stochasticity nor dependence between individual fecundity and survival. In this case the demographic variance is uniquely determined by the vital rates defining the projection matrix. The demographic variance for a long-lived bird species, the wandering albatross in the southwestern part of the Indian Ocean, is estimated. We also compute estimates of the age-specific contributions to the total demographic variance from survival, fecundity and the covariance between survival and fecundity.     

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.     

种群统计随机性和环境随机性对种群绝灭的影响

影响种群绝灭的随机干扰可分为种群统计随机性、环境随机性和随机灾害三大类。在相对稳定的环境条件下和相对较短的时间内,以前两类随机干扰对种群绝灭的影响为生态学家关注的焦点。但是,由于自然种群动态及其影响因子的复杂特征,进一步深入研究随机干扰对种群绝灭的作用在理论上和实践上都必须发展新的技术手段。本文回顾了种群统计随机性与环境随机性的概念起源与发展,系统阐述了其分析方法。归纳了两类随机性在种群绝灭研究中的应用范围、作用方式和特点的异同和区别方法。各类随机作用与种群动态之间关系的理论研究与对种群绝灭机理的实践研究紧密相关。根据理论模型模拟和自然种群实际分析两方面的研究现状,作者提出了进一步深入研究随机作用与种群非线性动态方法的策略。指出了随机干扰影响种群绝灭过程的研究的方向：更多的研究将从单纯的定性分析随机干扰对种群动力学简单性质的作用,转向结合特定的种群非线性动态特征和各类随机力作用特点具体分析绝灭极端动态的成因,以期做出精确的预测。     

Incorporating environmental stochasticity within a biological population model

The birth and death transition rates for a population are modelled as functions of both the population size and the environmental condition. An assortment of important theoretical results and techniques that can be utilized to analyze such a population’s behaviour is covered. Consequently, these results and techniques are used to study two examples. Firstly, we study a population with a stable equilibrium state, whose per capita birth and death rates are linearly related to the environmental condition. (The environmental condition in turn is modelled as an Ornstein–Uhlenbeck process.) Secondly, we study a population affected by two interdependent environmental factors.     

Phenological spread in plants: A result of adaptations to environmental stochasticity?

Summary Plant species co-inhabiting a given geographical region often have distinetly different times of flowering. It is shown that such phenological spread, duc to short-term stochastic variation in weather variables, relaxes competition for empty sites to be colonized by diaspores. For sufficiently large spreads stable coexistence becomes possible. The applicability of the proposed hypothesis to the observed instances of phenological spread is discussed and shown to extend beyond that of other current theories.We thank L.-E. Liljelund, C. Solbreck and C. Wiklund for helpful comments. This work was carried out within the Swedish Coniferous Forest Project, supported by the Swedish Natural Science Research Council, the Swedish Environmental Protection Board, the Swedish Council of Forestry and Agricultural Research, and the Wallenberg foundation.     

Population models with environmental stochasticity

Two discrete population models, one with stochasticity in the carrying capacity and one with stochasticity in the per capita growth rate, are investigated. Conditions under which the corresponding Markov processes are null recurrent and positively recurrent are derived.     

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.     

A neutral theory with environmental stochasticity explains static and dynamic properties of ecological communities

Understanding the forces shaping ecological communities is crucial to basic science and conservation. Neutral theory has made considerable progress in explaining static properties of communities, like species abundance distributions (SADs), with a simple and generic model, but was criticised for making unrealistic predictions of fundamental dynamic patterns and for being sensitive to interspecific differences in fitness. Here, we show that a generalised neutral theory incorporating environmental stochasticity may resolve these limitations. We apply the theory to real data (the tropical forest of Barro Colorado Island) and demonstrate that it much better explains the properties of short‐term population fluctuations and the decay of compositional similarity with time, while retaining the ability to explain SADs. Furthermore, the predictions are considerably more robust to interspecific fitness differences. Our results suggest that this integration of niches and stochasticity may serve as a minimalistic framework explaining fundamental static and dynamic characteristics of ecological communities.     

Linear discrete population models with two time scales in fast changing environments II: non-autonomous case

As the result of the complexity inherent in nature, mathematical models employed in ecology are often governed by a large number of variables. For instance, in the study of population dynamics we often deal with models for structured populations in which individuals are classified regarding their age, size, activity or location, and this structuring of the population leads to high dimensional systems. In many instances, the dynamics of the system is controlled by processes whose time scales are very different from each other. Aggregation techniques take advantage of this situation to build a low dimensional reduced system from which behavior we can approximate the dynamics of the complex original system.In this work we extend aggregation techniques to the case of time dependent discrete population models with two time scales where both the fast and the slow processes are allowed to change at their own characteristic time scale, generalizing the results of previous studies. We propose a non-autonomous model with two time scales, construct an aggregated model and give relationship between the variables governing the original and the reduced systems. We also explore how the properties of strong and weak ergodicity, regarding the capacity of the system to forget initial conditions, of the original system can be studied in terms of the reduced system.     

Sources of stochasticity in models of sex allocation in spatially structured populations

There are many ways to include stochastic effects in models of sex allocation evolution. These include variability in the number of mating partners and fecundity in a rich literature that goes back 20 years. The effects of variance in the fecundity and number of mating partners have typically been considered separately from the stochastic effects of mortality. However, I show that these processes produce mathematically equivalent models with subtly different biological details. These scenarios differ in the way that information becomes available to individuals because the parents often have information on mating partners while they are making sex allocation decisions, but must make these decisions before brood mortality takes place. This makes it possible to test which mechanism, stochastic mortality or variation in mating partners, is responsible for observed sex ratios. Alternatively, asymmetric variance between sexual functions can cause skewed sex allocation, even in the absence of local mate competition. This allows the evolution of either female- or male-biased sex ratios depending on which sexual function is more variable.     

Univariate Community Assembly Analysis (UniCAA): Combining hierarchical models with null models to test the influence of spatially restricted dispersal,environmental filtering,and stochasticity on community assembly

Identifying the influence of stochastic processes and of deterministic processes, such as dispersal of individuals of different species and trait‐based environmental filtering, has long been a challenge in studies of community assembly. Here, we present the Univariate Community Assembly Analysis (UniCAA) and test its ability to address three hypotheses: species occurrences within communities are (a) limited by spatially restricted dispersal; (b) environmentally filtered; or (c) the outcome of stochasticity—so that as community size decreases—species that are common outside a local community have a disproportionately higher probability of occurrence than rare species. The comparison with a null model allows assessing if the influence of each of the three processes differs from what one would expect under a purely stochastic distribution of species. We tested the framework by simulating “empirical” metacommunities under 15 scenarios that differed with respect to the strengths of spatially restricted dispersal (restricted vs. not restricted); habitat isolation (low, intermediate, and high immigration rates); and environmental filtering (strong, intermediate, and no filtering). Through these tests, we found that UniCAA rarely produced false positives for the influence of the three processes, yielding a type‐I error rate ≤5%. The type‐II error rate, that is, production of false negatives, was also acceptable and within the typical cutoff (20%). We demonstrate that the UniCAA provides a flexible framework for retrieving the processes behind community assembly and propose avenues for future developments of the framework.     

Predicting fluctuations of reintroduced ibex populations: the importance of density dependence, environmental stochasticity and uncertain population estimates

1. Development of population projections requires estimates of observation error, parameters characterizing expected dynamics such as the specific population growth rate and the form of density regulation, the influence of stochastic factors on population dynamics, and quantification of the uncertainty in the parameter estimates. 2. Here we construct a Population Prediction Interval (PPI) based on Bayesian state space modelling of future population growth of 28 reintroduced ibex populations in Switzerland that have been censused for up to 68 years. Our aim is to examine whether the interpopulation variation in the precision of the population projections is related to differences in the parameters characterizing the expected dynamics, in the effects of environmental stochasticity, in the magnitude of uncertainty in the population parameters, or in the observation error. 3. The error in the population censuses was small. The median coefficient of variation in the estimates across populations was 5.1%. 4. Significant density regulation was present in 53.6% of the populations, but was in general weak. 5. The width of the PPI calculated for a period of 5 years showed large variation among populations, and was explained by differences in the impact of environmental stochasticity on population dynamics. 6. In spite of the high accuracy in population estimates, the uncertainty in the parameter estimates was still large. This uncertainty affected the precision in the population predictions, but it decreased with increasing length of study period, mainly due to higher precision in the estimates of the environmental variance in the longer time-series. 7. These analyses reveal that predictions of future population fluctuations of weakly density-regulated populations such as the ibex often become uncertain. Credible population predictions require that this uncertainty is properly quantified.     

Stochastic models of some endemic infections

Stochastic models are established and studied for several endemic infections with demography. Approximations of quasi-stationary distributions and of times to extinction are derived for stochastic versions of SI, SIS, SIR, and SIRS models. The approximations are valid for sufficiently large population sizes. Conditions for validity of the approximations are given for each of the models. These are also conditions for validity of the corresponding deterministic model. It is noted that some deterministic models are unacceptable approximations of the stochastic models for a large range of realistic parameter values.     

Evolution of dispersal in metapopulations with local density dependence and demographic stochasticity

In this paper, we predict the outcome of dispersal evolution in metapopulations based on the following assumptions: (i) population dynamics within patches are density-regulated by realistic growth functions; (ii) demographic stochasticity resulting from finite population sizes within patches is accounted for; and (iii) the transition of individuals between patches is explicitly modelled by a disperser pool. We show, first, that evolutionarily stable dispersal rates do not necessarily increase with rates for the local extinction of populations due to external disturbances in habitable patches. Second, we describe how demographic stochasticity affects the evolution of dispersal rates: evolutionarily stable dispersal rates remain high even when disturbance-related rates of local extinction are low, and a variety of qualitatively different responses of adapted dispersal rates to varied levels of disturbance become possible. This paper shows, for the first time, that evolution of dispersal rates may give rise to monotonically increasing or decreasing responses, as well as to intermediate maxima or minima.     

Variance in animal longevity: contributions of heterogeneity and stochasticity

Variance in longevity among individuals may arise as an effect of heterogeneity (differences in mortality rates experienced at the same age or stage) or as an effect of individual stochasticity (the outcome of random demographic events during the life cycle). Decomposing the variance into components due to heterogeneity and stochasticity is crucial for evolutionary analyses.In this study, we analyze longevity from ten studies of invertebrates in the laboratory, and use the results to partition the variance in longevity into its components. To do so, we fit finite mixtures of Weibull survival functions to each data set by maximum likelihood, using the EM algorithm. We used the Bayesian Information Criterion to select the most well supported model. The results of the mixture analysis were used to construct an age × stage-classified matrix model, with heterogeneity groups as stages, from which we calculated the variance in longevity and its components. Almost all data sets revealed evidence of some degree of heterogeneity. The median contribution of unobserved heterogeneity to the total variance was 35%, with the remaining 65% due to stochasticity. The differences among groups in mean longevity were typically on the order of 30% of the overall life expectancy. There was considerable variation among data sets in both the magnitude of heterogeneity and the proportion of variance due to heterogeneity, but no clear patterns were apparent in relation to sex, taxon, or environmental conditions.     

Tuberculosis models with fast and slow dynamics: the role of close and casual contacts

Models that incorporate local and individual interactions are introduced in the context of the transmission dynamics of tuberculosis (TB). The multi-level contact structure implicitly assumes that individuals are at risk of infection from close contacts in generalized household (clusters) as well as from casual (random) contacts in the general population. Epidemiological time scales are used to reduce the dimensionality of the model and singular perturbation methods are used to corroborate the results of time-scale approximations. The concept and impact of optimal average cluster or generalized household size on TB dynamics is discussed. We also discuss the potential impact of our results on the spread of TB.     

Modeling baseline conditions of ecological indicators: Marine renewable energy environmental monitoring

Ecological indicators are often collected to detect and monitor environmental change. Statistical models are used to estimate natural variability, pre-existing trends, and environmental predictors of baseline indicator conditions. Establishing standard models for baseline characterization is critical to the effective design and implementation of environmental monitoring programs. An anthropogenic activity that requires monitoring is the development of Marine Renewable Energy sites. Currently, there are no standards for the analysis of environmental monitoring data for these development sites. Marine Renewable Energy monitoring data are used as a case study to develop and apply a model evaluation to establish best practices for characterizing baseline ecological indicator data. We examined a range of models, including six generalized regression models, four time series models, and three nonparametric models. Because monitoring data are not always normally distributed, we evaluated model ability to characterize normal and non-normal data using hydroacoustic metrics that serve as proxies for ecological indicator data. The nonparametric support vector regression and random forest models, and parametric state-space time series models generally were the most accurate in interpolating the normal metric data. Support vector regression and state-space models best interpolated the non-normally distributed data. If parametric results are preferred, then state-space models are the most robust for baseline characterization. Evaluation of a wide range of models provides a comprehensive characterization of the case study data, and highlights advantages of models rarely used in Marine Renewable Energy environmental monitoring. Our model findings are relevant for any ecological indicator data with similar properties, and the evaluation approach is applicable to any monitoring program.