Expected relative fitness and the adaptive topography of fluctuating selection

Wright's adaptive topography describes gene frequency evolution as a maximization of mean fitness in a constant environment. I extended this to a fluctuating environment by unifying theories of stochastic demography and fluctuating selection, assuming small or moderate fluctuations in demographic rates with a stationary distribution, and weak selection among the types. The demography of a large population, composed of haploid genotypes at a single locus or normally distributed phenotypes, can then be approximated as a diffusion process and transformed to produce the dynamics of population size, N, and gene frequency, p, or mean phenotype, . The expected evolution of p or is a product of genetic variability and the gradient of the long-run growth rate of the population, , with respect to p or . This shows that the expected evolution maximizes , the mean Malthusian fitness in the average environment minus half the environmental variance in population growth rate. Thus, as a function of p or represents an adaptive topography that, despite environmental fluctuations, does not change with time. The haploid model is dominated by environmental stochasticity, so the expected maximization is not realized. Different constraints on quantitative genetic variability, and stabilizing selection in the average environment, allow evolution of the mean phenotype to undergo a stochastic maximization of . Although the expected evolution maximizes the long-run growth rate of the population, for a genotype or phenotype the long-run growth rate is not a valid measure of fitness in a fluctuating environment. The haploid and quantitative character models both reveal that the expected relative fitness of a type is its Malthusian fitness in the average environment minus the environmental covariance between its growth rate and that of the population.     

Estimating population growth rates from stochastic Leslie matrices

Stochastic versions of exponential growth models predict that even when r or λ values calculated from mean vital statistics indicate exponential growth, most of the individual populations may become extinct. Several recent papers have considered this problem and some misunderstanding has arisen due to the difference between mathematical expectation of population size and most likely course of population growth. We replicated Boyce's (1977, 1979) simulations of population growth with age structure and a single randomly varying vital statistic, and reconciled some of these differences. Mean number can be projected using the dominant eigenvalue of the mean Leslie matrix, but the modal number may be considerably lower. We compared several measures of the rate of growth of the geometric mean or median of numbers and conclude that Tuljapurkar's α is an acceptable measure.     

Population dynamics in variable environments. II. Correlated environments,sensitivity analysis and dynamics

Mean and mean square number are studied for age-structured populations with serially correlated temporally fluctuating vital rates. Results are that (1) Moments of population number can be used effectively to analyse growth rates of the coefficient of variation and an approximate median population number. (2) Analytical approximations to the growth rates of moments reveal dynamic consequences of covarying phenotypic traits and of temporal correlation along environmental sequences. (3) Dynamic properties can be explicitly related to the static sensitivity of an average vital rate matrix. (4) The use of (1), (2) and (3) allows an extension of many applications of static vital rate theory to dynamics with fluctuating rates.     

Adaptive topography of fluctuating selection in a Mendelian population

An adaptive topography is derived for a large randomly mating diploid population under weak density-independent selection in a fluctuating environment. Assuming a stationary distribution of environmental states with no temporal autocorrelation, a diffusion approximation for population size and allele frequency, p, reveals that the expected change in p involves the gradient with respect to p of the stochastic intrinsic rate of increase (the density-independent long-run growth rate), r = r - sigma 2 e/2, where r is the mean Malthusian fitness in the average environment and is the environmental variance in population growth rate. The expected relative fitness of a genotype is its Malthusian fitness in the average environment minus the covariance of its fitness with population growth rate. The influence of fitness correlation between genotypes is illustrated by an analysis of the Haldane-Jayakar model of fluctuating selection on a single diallelic locus, and on two loci with additive effects on a quantitative character.     

Estimating stochastic elasticities directly from longitudinal data

The elasticities of long-run population growth rate with respect to vital rates are useful in studying selection on vital rates, and in evaluating management policy that aims to control vital rates. In temporally varying environments, elasticity is often calculated from simulations that assume a probability distribution for the environmental states. Here we develop a method to estimate elasticities directly from demographic data. Using a time-series of demographic matrices and age-structure we construct a consistent statistical estimator of elasticity that converges to the correct limiting value as the sample length increases. We also construct confidence intervals for elasticities from temporal data and suggest tools for testing hypotheses about the strength of selection. We use data on a natural population to show that our method can indeed accurately estimate elasticities using relatively short time series.     

A diffusion model for population growth in random environment

The growth of a population in a randomly varying environment is modeled by replacing the Malthusian growth rate with a delta-correlated normal process. The population size is then shown to be a random process, lognormally distributed, obeying a diffusion equation of the Fokker-Planck type. The first passage time p.d.f. through any arbitrarily assigned value and the probability of absorption are derived. The asymptotic behavior of the population size is investigated.     

Is the impact of environmental noise visible in the dynamics of age-structured populations?

Climate change has ignited lively research into its impact on various population-level processes. The research agenda in ecology says that some of the fluctuations in population size are accountable for by the external noise (e.g. weather) modulating the dynamics of populations. We obeyed the agenda by assuming population growth after a resource-limited Leslie matrix model in an age-structured population. The renewal process was disturbed by superimposing noise on the development of numbers in one or several age groups. We constructed models for iteroparous and semelparous breeders so that, for both categories, the population growth rate was matching. We analysed how the modulated population dynamics correlates with the noise signal with different time-lags. No significant correlations were observed for semelparous breeders, whereas for iteroparous breeders high correlations were frequently observed with time-lags of 71 year or longer. However, the latter occurs under red-coloured noise and for low growth rates when the disturbance is on the youngest age group only. It is laborious to find any clear signs of the (red) noise- and age group-specific fluctuations if the disturbance influences older age groups only. These results cast doubts on the possibility of detecting the signature of external disturbance after it has modulated temporal fluctuations in age-structured populations.     

Linking vital rates to invasiveness of a perennial herb

Invaders generally show better individual performance than non-invaders and, therefore, vital rates (survival, growth, fecundity) could potentially be used to predict species invasiveness outside their native range. Comparative studies have usually correlated vital rates with the invasiveness status of species, while few studies have investigated them in relation to population growth rate. Here, I examined the influence of five vital rates (plant establishment, survival, growth, flowering probability, seed production) and their variability (across geographic regions, habitat types, population sizes and population densities) on population growth rate (λ) using data from 37 populations of an invasive, iteroparous herb (Lupinus polyphyllus) in a part of its invaded range in Finland. Variation in vital rates was often related to habitat type and population density. The performance of the populations varied from declining to rapidly increasing independently of habitat type, population size or population density, but differed between regions. The population growth rate increased linearly with plant establishment, and with the survival and growth of vegetative individuals, while the survival of flowering individuals and annual seed production were not related to λ. The vital rates responsible for rapid population growth varied among populations. These findings highlight the importance of both regional and local conditions to plant population dynamics, demonstrating that individual vital rates do not necessarily correlate with λ. Therefore, to understand the role of individual vital rates in a species ability to invade, it is necessary to quantify their effect on population growth rate.     

Bayesian Estimates of Transition Probabilities in Seven Small Lithophytic Orchid Populations: Maximizing Data Availability from Many Small Samples

Predicting population dynamics for rare species is of paramount importance in order to evaluate the likelihood of extinction and planning conservation strategies. However, evaluating and predicting population viability can be hindered from a lack of data. Rare species frequently have small populations, so estimates of vital rates are often very uncertain due to lack of data. We evaluated the vital rates of seven small populations from two watersheds with varying light environment of a common epiphytic orchid using Bayesian methods of parameter estimation. From the Lefkovitch matrices we predicted the deterministic population growth rates, elasticities, stable stage distributions and the credible intervals of the statistics. Populations were surveyed on a monthly basis between 18–34 months. In some of the populations few or no transitions in some of the vital rates were observed throughout the sampling period, however, we were able to predict the most likely vital rates using a Bayesian model that incorporated the transitions rates from the other populations. Asymptotic population growth rate varied among the seven orchid populations. There was little difference in population growth rate among watersheds even though it was expected because of physical differences as a result of differing canopy cover and watershed width. Elasticity analyses of Lepanthes rupestris suggest that growth rate is more sensitive to survival followed by growth, shrinking and the reproductive rates. The Bayesian approach helped to estimate transition probabilities that were uncommon or variable in some populations. Moreover, it increased the precision of the parameter estimates as compared to traditional approaches.     

Bias in estimating the Malthusian parameter for Leslie matrices

For Leslie matrices of order 3 × 3 or larger, conditions for concavity or convexity of the Malthusian parameter in each of the entries in the matrix are given. Both cases are possible so it follows that the expected population growth rate computed from a Leslie matrix whose entries are random variables can be either smaller or larger than the growth rate computed from the expected value of the matrix. Boyce [(1977) Theor. Pop. Biol.12] showed that in the 2 × 2 case this bias is always positive; we give a numerical example illustrating the magnitude of the bias in this case, and compare it with the sampling error of the parameter for the same example.     

How does temporal variability affect predictions of weed population numbers?

1. Population models that are used to predict weed population dynamics or the impact of control measures on weed abundance typically ignore temporal variability in life-history parameters and control measures, and utilize mean arithmetic population growth rates to predict population abundance.
2. We demonstrate that the persistence of weeds in a stochastically varying environment depends on the geometric mean population growth rate being greater than zero, rather than the arithmetic mean population growth rate being greater than zero.
3. In a stochastically varying environment we show that temporal variability in fecundity, germination and survivorship will tend to decrease population size, relative to predictions based on arithmetic means. Conversely, variability in competitive effects and weed control will tend to increase population size, relative to predictions based on arithmetic mean values. The distinction between these two sets of parameters is that increases in the former will increase population growth rate, whereas increases in the latter will decrease it.
4. We argue that population models based on arithmetic mean population growth rates will tend to over-estimate population size. Numerical simulations indicate that this bias may be considerable.
5. Since short-term studies cannot, in general, estimate the geometric mean growth rate of a population we suggest several approaches for estimating the degree of bias in the predictions of models owing to the effects of variability. Accounting for such variability is necessary since current models for the dynamics of weed populations are based on arithmetic mean measures of population growth and hence likely to be biased.     

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.     

Darwinian fitness

The term Darwinian fitness refers to the capacity of a variant type to invade and displace the resident population in competition for available resources. Classical models of this dynamical process claim that competitive outcome is a deterministic event which is regulated by the population growth rate, called the Malthusian parameter. Recent analytic studies of the dynamics of competition in terms of diffusion processes show that growth rate predicts invasion success only in populations of infinite size. In populations of finite size, competitive outcome is a stochastic process--contingent on resource constraints--which is determined by the rate at which a population returns to its steady state condition after a random perturbation in the individual birth and death rates. This return rate, a measure of robustness or population stability, is analytically characterized by the demographic parameter, evolutionary entropy, a measure of the uncertainty in the age of the mother of a randomly chosen newborn. This article appeals to computational and numerical methods to contrast the predictive power of the Malthusian and the entropic principles. The computational analysis rejects the Malthusian model and is consistent with of the entropic principle. These studies thus provide support for the general claim that entropy is the appropriate measure of Darwinian fitness and constitutes an evolutionary parameter with broad predictive and explanatory powers.     

The Invisible Cliff: Abrupt Imposition of Malthusian Equilibrium in a Natural-Fertility,Agrarian Society

Analysis of a natural fertility agrarian society with a multi-variate model of population ecology isolates three distinct phases of population growth following settlement of a new habitat: (1) a sometimes lengthy copial phase of surplus food production and constant vital rates; (2) a brief transition phase in which food shortages rapidly cause increased mortality and lessened fertility; and (3) a Malthusian phase of indefinite length in which vital rates and quality of life are depressed, sometimes strikingly so. Copial phase duration declines with increases in the size of the founding group, maximum life expectancy and fertility; it increases with habitat area and yield per hectare; and, it is unaffected by the sensitivity of vital rates to hunger. Transition phase duration is unaffected by size of founding population and area of settlement; it declines with yield, life expectancy, fertility and the sensitivity of vital rates to hunger. We characterize the transition phase as the Malthusian transition interval (MTI), in order to highlight how little time populations generally have to adjust. Under food-limited density dependence, the copial phase passes quickly to an equilibrium of grim Malthusian constraints, in the manner of a runner dashing over an invisible cliff. The three-phase pattern diverges from widely held intuitions based on standard Lotka-Verhulst approaches to population regulation, with implications for the analysis of socio-cultural evolution, agricultural intensification, bioarchaeological interpretation of food stress in prehistoric societies, and state-level collapse.     

濒危植物长柄双花木自然种群数量动态

 运用种群生命表、生殖力表、Leslie矩阵模型和时间序列预测分析方法,研究了濒危植物长柄双花木(Disanthus cercidifolius var. longipes)种群的动态变化过程,揭示了长柄双花木各龄级植株的数量动态规律。结果表明：长柄双花木为缓慢负增长型种群;种群的净增殖率、内禀增长率以及周限增长率都较低,世代平均周期较长;Leslie矩阵模型和时间序列预测分析均表明在未来50年内长柄双花木各年龄级种群数量会出现波动性的消长,但其种群总数将逐步下降。导致种群下降的主要原因可能是人为砍伐及由此造成的生境破碎化等。     

The structure and dynamics of woodreed Calamagrostis canescens population: a modelling approach

A scale of ontogenetic states has been developed for woodreed Calamagrostis canescens, a perennial species dominating the grass layer of fell forest areas. The population structure is considered as a set of age-stage groups of individuals differing both in the ontogenetic stage and the chronological age measured in years. to describe the dynamics through years a special kind of matrix formalism has been proposed which is reducible neither to the classic Leslie matrix for an age-structured population, nor to the well-known Lefkovitch matrix for a stage-structured one, and which does not suffer from excessiveness of the "two-dimensional" representation for the structure implying the projection matrix of a block pattern. It has been shown however that the protection matrix corresponding to C. canescens life-history graph embodies the canonical features of matrix formalism for structured population dynamics, such as the exponential population growth or decline, the convergence to a stable equilibrium structure, the calculable indicator of growth/decline/equilibrium (i.e., a measure of the population reproductive potential) as well as possibility to determine the relative reproductive value of each group. On the other hand, "left-sidedness of the age spectrum", a property that is often observed in real populations and is inherent in Leslie models of growing populations, may fail in the age-stage-structured model. The aggregation of age-stage groups into the age classes is possible only under special strict relationship among the age-stage-specific vital rates of the population. The both circumstances serve a methodical indication that an additional dimension such as the stages, for example, ought to be introduced into the age structure of the model population.     

Dynamics of African ungulate populations with fluctuating, density-independent calf survival

Leslie matrix simulations were performed, employing idealized ungulate life histories and density-independent, either random or periodic, variation in calf survival, p₀. Long-term average rate of increase was reduced whenever variation was added to p₀. The iteroparous life histories strongly buffered high-frequency vital rate variations, so that large fluctuations in p₀ resulted in relatively small fluctuations in population size. The buffering effect was greatest for species with the longest lifespans. The results suggest that wild ungulate populations will not fluctuate greatly in response to density-independent, environmental-driven vital rates. Such populations will reach high densities, thus experiencing strong density-dependent regulation, only a few times in many generations.     

Demographic consequences of population subdivision on the long-furred woolly mouse opossum (Micoureus paraguayanus) from the Atlantic Forest

Habitat destruction and fragmentation severely affected the Atlantic Forest. Formerly contiguous populations may become subdivided into a larger number of smaller populations, threatening their long-term persistence. The computer package VORTEX was used to simulate the consequences of habitat fragmentation and population subdivision on Micoureus paraguayanus, an endemic arboreal marsupial of the Atlantic Forest. Scenarios simulated hypothetical populations of 100 and 2000 animals being partitioned into 1–10 populations, linked by varying rates of inter-patch dispersal, and also evaluated male-biased dispersal. Results demonstrated that a single population was more stable than an ensemble of populations of equal size, irrespective of dispersal rate. Small populations (10–20 individuals) exhibited high instability due to demographic stochasticity, and were characterized by high rates of extinction, smaller values for metapopulation growth and larger fluctuations in population size and growth rate. Dispersal effects on metapopulation persistence were related to the size of the populations and to the sexes that were capable of dispersing. Male-biased dispersal had no noticeable effects on metapopulation extinction dynamics, whereas scenarios modelling dispersal by both sexes positively affected metapopulation dynamics through higher growth rates, smaller fluctuations in growth rate, larger final metapopulation sizes and lower probabilities of extinction. The present study highlights the complex relationships between metapopulation size, population subdivision, habitat fragmentation, rate of inter-patch dispersal and sex-biased dispersal and indicates the importance of gaining a better understanding of dispersal and its interactions with correlations between disturbance events.     

The units of selection and measures of fitness

This paper presents a unified account of the properties of the measures, Malthusian parameter and entropy in predicting evolutionary change in populations of macromolecules, cells and individuals. The Malthusian parameter describes the intrinsic rate of increase of the population. The entropy describes the intrinsic variability in populations: it characterizes the variability in mutation and replication rates in populations of macromolecules; the rate of decay of synchrony in populations of cells; and the degree of iteroparity in populations of individuals. The Malthusian parameter determines ultimate population numbers: under constant environmental conditions, it is the rate of increase when equilibrium conditions are attained. Entropy determines population stability: the gain in the Malthusian parameter due to small fluctuations in the life-cycle variables is determined by entropy. These properties, which are valid for populations of macromolecules, cells and individuals, show that the Malthusian parameter and entropy act as complimentary fitness indices in understanding evolutionary change in populations.     

Stochastic LTRE analysis of the effects of herbivory on the population dynamics of a perennial grassland herb

Herbivores can have strong deleterious effects on vital rates (growth, reproduction, and survival) and thus negatively impact the population dynamics of plant species. In practice, however, these effects might be strongly correlated, for example as a result of tradeoffs between vital rates. To get better insights into the effects of herbivory on the population dynamics of the long‐lived grassland plant Primula veris population projection matrices were constructed from demographic data collected between 1999 and 2008 (nine annual transitions). Data were collected in two large grassland populations, each of which was subjected to two treatments (grazing by cattle versus a mowing treatment), yielding a total of 36 matrices. We applied a lower‐level vital rate life table response experiment (LTRE) using the small noise approximation (SNA) of the stochastic population growth rate to disentangle the contributions of changes in mean vital rates, variability in vital rates, correlations between vital rates and vital rate elasticities to the difference in the stochastic growth rate. Stochastic growth rates (a= log λ_S) were significantly lower in grazed than in mown plots (a= 0.0185 and 0.1019, respectively). SNA LTRE analysis showed that contributions of mean vital rates by far made the largest contribution to the observed difference in a between grazed and control plots. In particular, changes in sexual reproduction rates made the largest contributions to lower the stochastic growth rate in grazed plots: both adult flowering probabilities and flower and seed production were importantly lower in grazed populations, but these negative effects were largely buffered by increased establishment and seedling survival rates. Among the stochastic terms of the SNA decomposition, contributions of covariance and correlations between vital rates had the largest impact, whereas contributions of elasticities were smaller. The strongest correlation driver was the association between adult survival and seedling establishment, suggesting that environmental conditions favouring adult survival also are beneficial for seedling establishment. Overall, our results show that herbivory had a strong negative effect on the long‐term population growth rate of P. veris that was primarily mediated by differences in fecundity (flower and seed production) and germination.