Evolutionary rescue under demographic and environmental stochasticity

Populations suffer two types of stochasticity: demographic stochasticity, from sampling error in offspring number, and environmental stochasticity, from temporal variation in the growth rate. By modelling evolution through phenotypic selection following an abrupt environmental change, we investigate how genetic and demographic dynamics, as well as effects on population survival of the genetic variance and of the strength of stabilizing selection, differ under the two types of stochasticity. We show that population survival probability declines sharply with stronger stabilizing selection under demographic stochasticity, but declines more continuously when environmental stochasticity is strengthened. However, the genetic variance that confers the highest population survival probability differs little under demographic and environmental stochasticity. Since the influence of demographic stochasticity is stronger when population size is smaller, a slow initial decline of genetic variance, which allows quicker evolution, is important for population persistence. In contrast, the influence of environmental stochasticity is population-size-independent, so higher initial fitness becomes important for survival under strong environmental stochasticity. The two types of stochasticity interact in a more than multiplicative way in reducing the population survival probability. Our work suggests the importance of explicitly distinguishing and measuring the forms of stochasticity during evolutionary rescue.     

Demographic stochasticity and temporal autocorrelation in the dynamics of structured populations

Populations can show temporal autocorrelation in the dynamics arising from different mechanisms, including fluctuations in the demographic structure. This autocorrelation is often treated as a complicating factor in the analyses of stochastic population growth and extinction risk. However, it also reflects important information about the demographic structure. Here, we consider how temporal autocorrelation is related to demographic stochasticity in structured populations. Demographic stochasticity arises from inherent randomness in the demographic processes of individuals, like survival and reproduction, and the resulting impact on population growth is measured by the demographic variance. Earlier studies have shown that population structure have positive or negative effects on the demographic variance compared to a model where the structure is ignored. Here, we derive a new expression for the demographic variance of a structured population, using the temporal autocorrelation function of the population growth rate. We show that the relative difference in demographic variance when the structure is included or ignored (the effect of structure on demographic variance) is approximately twice the sum of the autocorrelations. We demonstrate the result for a simple hypothetical example, as well as a set of empirical examples using age‐structured models of 24 mammals from the demographic database COMADRE. In the empirical examples, the sum of the autocorrelation function was negative in all cases, indicating that age structure generally has a negative effect on the demographic variance (i.e. the demographic variance is lower compared to that of a model where the structure is ignored). Other kinds of structure, such as spatial heterogeneity affecting fecundity, can have positive effects on the demographic variance, and the sum of the autocorrelations will then be positive. These results yield new insights into the complex interplay between population structure, demographic variance, and temporal autocorrelation, that shapes the population dynamics and extinction risk of populations.     

EVOLUTION IN FLUCTUATING ENVIRONMENTS: DECOMPOSING SELECTION INTO ADDITIVE COMPONENTS OF THE ROBERTSON–PRICE EQUATION

We analyze the stochastic components of the Robertson–Price equation for the evolution of quantitative characters that enables decomposition of the selection differential into components due to demographic and environmental stochasticity. We show how these two types of stochasticity affect the evolution of multivariate quantitative characters by defining demographic and environmental variances as components of individual fitness. The exact covariance formula for selection is decomposed into three components, the deterministic mean value, as well as stochastic demographic and environmental components. We show that demographic and environmental stochasticity generate random genetic drift and fluctuating selection, respectively. This provides a common theoretical framework for linking ecological and evolutionary processes. Demographic stochasticity can cause random variation in selection differentials independent of fluctuating selection caused by environmental variation. We use this model of selection to illustrate that the effect on the expected selection differential of random variation in individual fitness is dependent on population size, and that the strength of fluctuating selection is affected by how environmental variation affects the covariance in Malthusian fitness between individuals with different phenotypes. Thus, our approach enables us to partition out the effects of fluctuating selection from the effects of selection due to random variation in individual fitness caused by demographic stochasticity.     

Life history evolution and demographic stochasticity

Summary Can demographic stochasticity bias the evolution of life history traits? Under a neutral version of the Cole-Charnov-Schaffer model, variance in offspring number for both annuals and perennials depends on the precise values of fitness components. Either annuals or perennials may have the larger variance (for equal ), depending on the importance of random survivalversus fixed reproduction. By extension, the variance in offspring number should generally depend on whether is mainly composed of highly variable elements or elements with limited variation. Thus, data about the variability of demographic parameters may be as important as data about their mean values.This result concerns only one source of demographic stochasticity, the probabilistic nature of demographic processes like survival. The other source of demographic stochasticity is the fact that populations are composed of whole numbers of individuals (integer arithmetic). Integer arithmetic without probabilistic demography (or environmental variation) can make it difficult for rare invaders to persist in populations even when selection would favour the invaders in a deterministic model. Integer arithmetic can also cause population coexistence when the equivalent deterministic model leads to exclusion. This effect disappears when demography is probabilistic, and probably also when there is environmental variation. Thus probabilistic demography and environmental variation may make some population patterns more, rather than less, understandable.     

The extended Moran effect and large-scale synchronous fluctuations in the size of great tit and blue tit populations

1. Synchronous fluctuations of geographically separated populations are in general explained by the Moran effect, i.e. a common influence on the local population dynamics of environmental variables that are correlated in space. Empirical support for such a Moran effect has been difficult to provide, mainly due to problems separating out effects of local population dynamics, demographic stochasticity and dispersal that also influence the spatial scaling of population processes. Here we generalize the Moran effect by decomposing the spatial autocorrelation function for fluctuations in the size of great tit Parus major and blue tit Cyanistes caeruleus populations into components due to spatial correlations in the environmental noise, local differences in the strength of density regulation and the effects of demographic stochasticity. 2. Differences between localities in the strength of density dependence and nonlinearity in the density regulation had a small effect on population synchrony, whereas demographic stochasticity reduced the effects of the spatial correlation in environmental noise on the spatial correlations in population size by 21.7% and 23.3% in the great tit and blue tit, respectively. 3. Different environmental variables, such as beech mast and climate, induce a common environmental forcing on the dynamics of central European great and blue tit populations. This generates synchronous fluctuations in the size of populations located several hundred kilometres apart. 4. Although these environmental variables were autocorrelated over large areas, their contribution to the spatial synchrony in the population fluctuations differed, dependent on the spatial scaling of their effects on the local population dynamics. We also demonstrate that this effect can lead to the paradoxical result that a common environmental variable can induce spatial desynchronization of the population fluctuations. 5. This demonstrates that a proper understanding of the ecological consequences of environmental changes, especially those that occur simultaneously over large areas, will require information about the spatial scaling of their effects on local population dynamics.     

Estimating the growth of a newly established moose population using reproductive value

Estimating the population growth rate and environmental stochasticity of long-lived species is difficult because annual variation in population size is influenced by temporal autocorrelations caused by fluctuations in the age-structure. Here we use the dynamics of the reproductive value to estimate the long-term growth rate s and the environmental variance of a moose population that recently colonized the island of Vega in northern Norway. We show that the population growth rate was high (ŝ=0.26). The major stochastic influences on the population dynamics were due to demographic stochasticity, whereas the environmental variance was not significantly different from 0. This supports the suggestion that population growth rates of polytocous ungulates are high, and that demographic stochasticity must be assessed when estimating the growth of small ungulate populations.     

Experimental support of the scaling rule for demographic stochasticity

A scaling rule of ecological theory, accepted but lacking experimental confirmation, is that the magnitude of fluctuations in population densities due to demographic stochasticity scales inversely with the square root of population numbers. This supposition is based on analyses of models exhibiting exponential growth or stable equilibria. Using two quantitative measures, we extend the scaling rule to situations in which population densities fluctuate due to nonlinear deterministic dynamics. These measures are applied to populations of the flour beetle Tribolium castaneum that display chaotic dynamics in both 20-g and 60-g habitats. Populations cultured in the larger habitat exhibit a clarification of the deterministic dynamics, which follows the inverse square root rule. Lattice effects, a deterministic phenomenon caused by the discrete nature of individuals, can cause deviations from the scaling rule when population numbers are small. The scaling rule is robust to the probability distribution used to model demographic variation among individuals.     

Fluctuating Environments and Phytoplankton Community Structure: A Stochastic Model

Spatial heterogeneity in organism and resource distributions can generate temporal heterogeneity in resource access for simple organisms like phytoplankton. The role of temporal heterogeneity as a structuring force for simple communities is investigated via models of phytoplankton with contrasting life histories competing for a single fluctuating resource. A stochastic model in which environmental and demographic stochasticity are treated separately is compared with a model with deterministic resource variation to assess the importance of stochasticity. When compared with the deterministic model, the stochastic model allows for coexistence over a wider range of parameter values (or life-history types). The model suggests that demographic stochasticity alone is far more important in increasing the possibility of coexistence than environmental stochasticity alone. However, the combined effects of both types of stochasticity produce the largest likelihood of coexistence. Finally, the influence of relative nutrient levels and nutrient pulse frequency on these results is addressed. We relate our findings to variable environment theory with evidence for both relative nonlinearity and the storage effect acting in this model. We show for the first time that temporal dynamics generated by demographic stochasticity may operate like the storage effect at particular spatial scales.     

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

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

Seed dormancy and delayed flowering in monocarpic plants: selective interactions in a stochastic environment

We explore the effects of temporal variation in multiple demographic rates on the joint evolution of delayed reproduction and seed dormancy using integral projection models (IPMs). To do this, we extend the standard IPM to include a discrete state variable representing the number of seeds in the seed bank, density-dependent recruitment, and temporal variation in demography. Parameter estimates for Carlina vulgaris and Carduus nutans are obtained from long-term studies. Carlina is relatively long lived and has a short-lived seed bank, whereas most Carduus plants flower in their first year and the seed bank is long lived. Using the evolutionarily stable strategy (ESS) approach, we predict the observed flowering and germination strategies. There is excellent agreement between the predictions and the field observations. The effects of temporal variation on the joint ESS are partitioned into components arising from nonlinear averaging (systematic changes in the mean resulting from the interaction between variability and nonlinearity) and nonequilibrium dynamics (fluctuations in fitness caused by temporal variation). This shows that temporal variation can have substantial effects on the observed flowering and germination strategies and that covariance between demographic processes is important. We extend the models to include spatial population structure and assess the robustness of the results from the nonspatial models.     

Durable Resistance to Crop Pathogens: An Epidemiological Framework to Predict Risk under Uncertainty

Increasing the durability of crop resistance to plant pathogens is one of the key goals of virulence management. Despite the recognition of the importance of demographic and environmental stochasticity on the dynamics of an epidemic, their effects on the evolution of the pathogen and durability of resistance has not received attention. We formulated a stochastic epidemiological model, based on the Kramer-Moyal expansion of the Master Equation, to investigate how random fluctuations affect the dynamics of an epidemic and how these effects feed through to the evolution of the pathogen and durability of resistance. We focused on two hypotheses: firstly, a previous deterministic model has suggested that the effect of cropping ratio (the proportion of land area occupied by the resistant crop) on the durability of crop resistance is negligible. Increasing the cropping ratio increases the area of uninfected host, but the resistance is more rapidly broken; these two effects counteract each other. We tested the hypothesis that similar counteracting effects would occur when we take account of demographic stochasticity, but found that the durability does depend on the cropping ratio. Secondly, we tested whether a superimposed external source of stochasticity (for example due to environmental variation or to intermittent fungicide application) interacts with the intrinsic demographic fluctuations and how such interaction affects the durability of resistance. We show that in the pathosystem considered here, in general large stochastic fluctuations in epidemics enhance extinction of the pathogen. This is more likely to occur at large cropping ratios and for particular frequencies of the periodic external perturbation (stochastic resonance). The results suggest possible disease control practises by exploiting the natural sources of stochasticity.     

Geographical gradients in the population dynamics of North American prairie ducks

1. Geographic gradients in population dynamics may occur because of spatial variation in resources that affect the deterministic components of the dynamics (i.e. carrying capacity, the specific growth rate at small densities or the strength of density regulation) or because of spatial variation in the effects of environmental stochasticity. To evaluate these, we used a hierarchical Bayesian approach to estimate parameters characterizing deterministic components and stochastic influences on population dynamics of eight species of ducks (mallard, northern pintail, blue-winged teal, gadwall, northern shoveler, American wigeon, canvasback and redhead (Anas platyrhynchos, A. acuta, A. discors, A. strepera, A. clypeata, A. americana, Aythya valisineria and Ay. americana, respectively) breeding in the North American prairies, and then tested whether these parameters varied latitudinally. 2. We also examined the influence of temporal variation in the availability of wetlands, spring temperature and winter precipitation on population dynamics to determine whether geographical gradients in population dynamics were related to large-scale variation in environmental effects. Population variability, as measured by the variance of the population fluctuations around the carrying capacity K, decreased with latitude for all species except canvasback. This decrease in population variability was caused by a combination of latitudinal gradients in the strength of density dependence, carrying capacity and process variance, for which details varied by species. 3. The effects of environmental covariates on population dynamics also varied latitudinally, particularly for mallard, northern pintail and northern shoveler. However, the proportion of the process variance explained by environmental covariates, with the exception of mallard, tended to be small. 4. Thus, geographical gradients in population dynamics of prairie ducks resulted from latitudinal gradients in both deterministic and stochastic components, and were likely influenced by spatial differences in the distribution of wetland types and shapes, agricultural practices and dispersal processes. 5. These results suggest that future management of these species could be improved by implementing harvest models that account explicitly for spatial variation in density effects and environmental stochasticity on population abundance.     

Scaling of the mean and variance of population dynamics under fluctuating regimes

Theoretical ecologists have long sought to understand how the persistence of populations depends on the interactions between exogenous (biotic and abiotic) and endogenous (e.g., demographic and genetic) drivers of population dynamics. Recent work focuses on the autocorrelation structure of environmental perturbations and its effects on the persistence of populations. Accurate estimation of extinction times and especially determination of the mechanisms affecting extinction times is important for biodiversity conservation. Here we examine the interaction between environmental fluctuations and the scaling effect of the mean population size with its variance. We investigate how interactions between environmental and demographic stochasticity can affect the mean time to extinction, change optimal patch size dynamics, and how it can alter the often-assumed linear relationship between the census size and the effective population size. The importance of the correlation between environmental and demographic variation depends on the relative importance of the two types of variation. We found the correlation to be important when the two types of variation were approximately equal; however, the importance of the correlation diminishes as one source of variation dominates. The implications of these findings are discussed from a conservation and eco-evolutionary point of view.     

Demographic heterogeneity impacts density-dependent population dynamics

Among-individual variation in vital parameters such as birth and death rates that is unrelated to age, stage, sex, or environmental fluctuations is referred to as demographic heterogeneity. This kind of heterogeneity is prevalent in ecological populations, but is almost always left out of models. Demographic heterogeneity has been shown to affect demographic stochasticity in small populations and to increase growth rates for density-independent populations. The latter is due to ??cohort selection,?? where the most frail individuals die out first, lowering the cohort??s average mortality as it ages. The importance of cohort selection to population dynamics has only recently been recognized. We use a continuous-time model with density dependence, based on the logistic equation, to study the effects of demographic heterogeneity in mortality and reproduction. Reproductive heterogeneity is introduced in three ways: parent fertility, offspring viability, and parent?Coffspring correlation. We find that both the low-density growth rate and the equilibrium population size increase as the magnitude of mortality heterogeneity increases or as parent?Coffspring phenotypic correlation increases. Population dynamics are affected by complex interactions among the different types of heterogeneity, and trade-off scenarios are examined which can sometimes reverse the effect of increased heterogeneity. We show that there are a number of different homogeneous approximations to heterogeneous models, but all fail to capture important parts of the dynamics of the full model.     

Density dependence and stochastic variation in a newly established population of a small songbird

Models describing fluctuations in population size should include both density dependence and stochastic effects. We examine the relative contribution of variation in parameters of the expected dynamics as well as demographic and environmental stochasticity to fluctuations in a population of a small passerine bird, the pied flycatcher, that was newly established in a Dutch study area. Using the theta-logistic model of density regulation, we demonstrate that the estimated quasi-stationary distribution including demographic stochasticity is close to the stationary distribution ignoring demographic stochasticity, indicating a long expected time to extinction. We also show that the variance in the estimated quasi-stationary distribution is especially sensitive to variation in the density regulation function. Reliable population projections must therefore account for uncertainties in parameter estimates which we do by using the population prediction interval (PPI). After 2 years the width of the 90% PPI was already larger than the corresponding estimated range of variation in the quasi-stationary distribution. More precise prediction of future population size than can be derived from the quasi-stationary distribution could only be made for a time span less than about five years.     

Paternal but not maternal age influences early-life performance of offspring in a long-lived seabird

Variability in demographic traits between individuals within populations has profound implications for both evolutionary processes and population dynamics. Parental effects as a source of non-genetic inheritance are important processes to consider to understand the causes of individual variation. In iteroparous species, parental age is known to influence strongly reproductive success and offspring quality, but consequences on an offspring fitness component after independence are much less studied. Based on 37 years longitudinal monitoring of a long-lived seabird, the wandering albatross, we investigate delayed effects of parental age on offspring fitness components. We provide evidence that parental age influences offspring performance beyond the age of independence. By distinguishing maternal and paternal age effects, we demonstrate that paternal age, but not maternal age, impacts negatively post-fledging offspring performance.     

Impacts of environmental variability in open populations and communities: "inflation" in sink environments

Ecological communities are typically open to the immigration and emigration of individuals, and also variable through time. In this paper we argue that interesting and potentially important effects arise when one splices together spatial fluxes and temporal variability. The particular system we examine is a sink habitat, where a species faces deterministic extinction but is rescued by recurrent immigration. We have shown, using a simple extension of the canonical exponential growth model in a time-varying environment, that variation "inflates" the average abundance of sink populations. We can analytically quantify the magnitude of this effect in several special cases (square-wave temporal variation and Gaussian stochastic variation). The inflationary effect can be large in "intermittent" sinks (where there are periods with positive growth), and when temporal variation is strongly autocorrelated. The effect appears to be robust to incorporation of demographic stochasticity (due to discrete birth-death-immigration processes), and to direct density dependence. With discrete generations, however, one can observe a wide range of effects of temporal variation, including depression as well as inflation. We argue that the inflationary effect of temporal variation in sink habitats can have important implications for community structure, because it can increase the average abundance (and hence local impacts) of species that on average are being excluded from a local community. We illustrate the latter effect using a familiar model of exploitative competition for a single limiting resource. We demonstrate that temporal variation can reverse local competitive dominance, even to the extent of allowing an inferior competitor maintained by immigration to exclude a competing species that would be locally superior in a constant environment.     

Distance decay of similarity,effects of environmental noise and ecological heterogeneity among species in the spatio‐temporal dynamics of a dispersal‐limited community

The joint spatial and temporal fluctuations in community structure may be due to dispersal, variation in environmental conditions, ecological heterogeneity among species and demographic stochasticity. These factors are not mutually exclusive, and their relative contribution towards shaping species abundance distributions and in causing species fluctuations have been hard to disentangle. To better understand community dynamics when the exchange of individuals between localities is very low, we studied the dynamics of the freshwater zooplankton communities in 17 lakes located in independent catchment areas, sampled at end of summer from 2002 to 2008 in Norway. We analysed the joint spatial and temporal fluctuations in the community structure by fitting the two‐dimensional Poisson lognormal model under a two‐stage sampling scheme. We partitioned the variance of the distribution of log abundance for a random species at a random time and location into components of demographic stochasticity, ecological heterogeneity among species, and independent environmental noise components for the different species. Non‐neutral mechanisms such as ecological heterogeneity among species (20%) and spatiotemporal variation in the environment (75%) explained the majority of the variance in log abundances. Overdispersion relative to Poisson sampling and demographic stochasticity had a small contribution to the variance (5%). Among a set of environmental variables, lake acidity was the environmental variable that was most strongly related to decay of community similarity in space and time.     

Coevolution at multiple spatial scales: Linum marginale–Melampsora lini – from the individual to the species

Coevolutionary processes are intrinsically spatial as well as temporal, and occur at many different scales. These range from single populations dominated by demographic and genetic stochasticity, to metapopulations in which colonisation/extinction dynamics have a large influence, and larger geographic regions where phylogenetic patterns and historical events become important. We present data for the genetically and demographically well-characterised plant host–pathogen interaction, the Linum marginale–Melampsora lini system, and use this to demonstrate the varying nature of resistance and virulence structure across these spatial scales. At the within population level, our results indicate considerable variability in resistance and virulence, but little evidence of coordinated changes in host and pathogen. Studies involving comparisons among multiple demes within a single metapopulation show that adjacent populations often have asynchronous disease dynamics and large differences in diversity and frequency of resistance and virulence phenotypes. Nevertheless, at this scale, there is also evidence of spatial structure in that more closely adjacent host populations are significantly more likely to have similar resistance phenotypes and mean levels of resistance. At larger scales, comparisons among adjacent metapopulations indicate that quantitative differences in host mating system and other life history features can have further major consequences for how host and pathogen variation is packaged. Finally, comparisons at continental and among host-species levels show variation consistent with specialisation and speciation in the pathogen. This revised version was published online in July 2006 with corrections to the Cover Date.     

Consequences of Environmental Fluctuations on Taylor’s Power Law and Implications for the Dynamics and Persistence of Populations

Conservation Biologists have found that demographic stochasticity causes the mean time to extinction to increase exponentially with population size. This has proved helpful in analyses determining extinction times and characterizing the pathway to extinction. The aim of this investigation is to explore the possible interactions between environmental/demographic noises and the scaling effect of the mean population size with its variance, which is expected to follow Taylor’s power law relationship. We showed that the combined effects of environmental/demographic noises and the scaling of population size variability interact with the population dynamics and affect the mean time to extinction.