Continuous-discrete models of the dynamics of an isolated population and of 2 competing species

The paper presents the analysis of various mathematical models for dynamics of isolated population and for competition between two species. It is assumed that mortality is continuous and birth of individuals of new generations takes place in certain fixed moments. Influence of winter upon the population dynamics and conditions of classic discrete model "deduction" of population dynamics (in particular, Moran-Ricker and Hassel's models) are investigated. Dynamic regimes of models under various assumptions about the birth and death rates upon the population states are also examined. Analysis of models of isolated population dynamics with nonoverlapping generations showed the density changes regularly if the birth rate is constant. Moreover, there exists a unique global stable level and population size stabilizes asymptotically at this equilibrium, i.e. cycle and chaotic regimes in various discrete models depend on correlation between individual productivity and population state in previous time. When the correlation is exponential upon mean population size the discrete Hassel model is realized. Modification of basis model, based on the assumption that during winter survival/death changes are constant, showed that population size at global level is stable. Generally, the dependence of population rate upon "winter parameters" has nonlinear character. Nonparametric models of competition between two species does not vary if the individual productivity is constant. In a phase space there are several stable stationary states and population stabilizes at one or other level asymptotically. So, in discrete models of competition between two species oscillation can be explained by dependence of population growth rate on the population size at previous times.     

A model for the dynamics of an annual plant population

A model for the dynamics of a single species population of plants is proposed and its use demonstrated by the analysis of a simple example. The model incorporates the effects of microsite variation by allowing for individual differences in growth and death rates within each season. We demonstrate that an increase in the variance in individual growth rates may increase both the chances that a plant population will persist and the equilibrium size of that population. We also show that even if size-dependent death is occurring, it may not have a significant effect on the shape of the size frequency distribution. An extension of the model to multispecies communities of plants suggests an experimental procedure to determine whether competition is responsible for excluding a particular plant species from a community that appears otherwise to be suitable. A more detailed analysis of the model for a two-species community produces conditions for competitive coexistence reminiscent of those from the Lotka-Volterra competition equations. Another extension suggests that selection will favor those genotypes that maximize the product of germination probability and mass of seeds produced, if survivorship and growth are not substantially altered. Finally, an analog to r- and K-selection theory for animal populations is developed. Selection in low-density populations favors increasing growth rate, and in high-density populations favors minimizing the effect of neighbors on one's own growth rate.     

Density-dependence as a size-independent regulatory mechanism

The growth function of populations is central in biomathematics. The main dogma is the existence of density-dependence mechanisms, which can be modelled with distinct functional forms that depend on the size of the population. One important class of regulatory functions is the theta-logistic, which generalizes the logistic equation. Using this model as a motivation, this paper introduces a simple dynamical reformulation that generalizes many growth functions. The reformulation consists of two equations, one for population size, and one for the growth rate. Furthermore, the model shows that although population is density-dependent, the dynamics of the growth rate does not depend either on population size, nor on the carrying capacity. Actually, the growth equation is uncoupled from the population size equation, and the model has only two parameters, a Malthusian parameter rho and a competition coefficient theta. Distinct sign combinations of these parameters reproduce not only the family of theta-logistics, but also the van Bertalanffy, Gompertz and Potential Growth equations, among other possibilities. It is also shown that, except for two critical points, there is a general size-scaling relation that includes those appearing in the most important allometric theories, including the recently proposed Metabolic Theory of Ecology. With this model, several issues of general interest are discussed such as the growth of animal population, extinctions, cell growth and allometry, and the effect of environment over a population.     

Inter- and intraspecific relationships between performance and temperature in a cryptic species complex of the rotifer Brachionus plicatilis

The strategy of decreasing size with increasing temperature operates at regional and phenotypic scale and presents a puzzle to researchers. In this work, we studied two aspects of the temperature–performance relationship along a temperature gradient, (i) comparing the population growth rates of three cryptic Brachionus species differing in adult size, and (ii) assessing the phenotypic plasticity of adult size, in one clone per species. The working hypotheses were that (i) the bigger the species the lower its optimal temperature for population growth, and (ii) the higher the temperature the smaller the individual within each focal species. The results showed that (i) the optimal temperature for population growth is related to species size in a manner foreseen by Bergmanns’ rule for two of the three species (the third, biggest species, performed evenly well at all temperatures examined, what could be explained by its generally eurioic character), and that (ii) the strategy of body size adjustment to environmental temperature differs between species and may depend on the level of temperature specialization. This work demonstrated the usefulness of inter- and intraspecific comparisons for studying the role of growth strategies in adaptation to temperature.     

Effects of body size and temperature on population growth

For at least 200 years, since the time of Malthus, population growth has been recognized as providing a critical link between the performance of individual organisms and the ecology and evolution of species. We present a theory that shows how the intrinsic rate of exponential population growth, rmax, and the carrying capacity, K, depend on individual metabolic rate and resource supply rate. To do this, we construct equations for the metabolic rates of entire populations by summing over individuals, and then we combine these population-level equations with Malthusian growth. Thus, the theory makes explicit the relationship between rates of resource supply in the environment and rates of production of new biomass and individuals. These individual-level and population-level processes are inextricably linked because metabolism sets both the demand for environmental resources and the resource allocation to survival, growth, and reproduction. We use the theory to make explicit how and why rmax exhibits its characteristic dependence on body size and temperature. Data for aerobic eukaryotes, including algae, protists, insects, zooplankton, fishes, and mammals, support these predicted scalings for rmax. The metabolic flux of energy and materials also dictates that the carrying capacity or equilibrium density of populations should decrease with increasing body size and increasing temperature. Finally, we argue that body mass and body temperature, through their effects on metabolic rate, can explain most of the variation in fecundity and mortality rates. Data for marine fishes in the field support these predictions for instantaneous rates of mortality. This theory links the rates of metabolism and resource use of individuals to life-history attributes and population dynamics for a broad assortment of organisms, from unicellular organisms to mammals.     

Analysis of interspecific competition in perennial plants using Life Table Response Experiments

The impact of interspecific competition is usually measured by its effect upon plant growth, neglecting impacts upon other stages of the life cycle such as fecundity which have a direct influence upon individual fitness and the asymptotic population growth rate of a population (λ). We used parameterized matrix models for three perennial plant species grown with and without interspecific competition to illustrate how the methodology of Life Table Response Experiments (LTRE) can be used to link any change in population dynamics to changes in any part of the life cycle. Plants were herbaceous grassland species grown for two years in a field experiment at Rothamsted Experimental Station, England. Interspecific competition reduced λ by over 90% in all species. Survival and growth were slightly affected by competition whereas plant fecundity was greatly reduced. Nearly all of the observed difference in λ between the competition treatments was explained by the fecundity terms, and more precisely by a large difference in the number of seeds, and a high sensitivity of λ to the germination rate. Whereas most competition studies focus on the measurement of change in individual fitness, our study illustrates how informative it is to take account not only of the effect of competition upon vital rates but also of how different vital rates affect population growth rate.     

Multiple environmental changes interact to modify species dynamics and invasion rates

Multiple aspects of the environment often change at the same time, influencing populations directly by modifying their physiology, but also indirectly by influencing other interacting species. The impacts of each environmental change upon population dynamics are usually assumed be independent of the state of other aspects of the environment, despite evidence at the individual level indicating that the combined impacts are often non‐additive. The importance of indirect effects mediated through community interactions also has high uncertainty. We used experimental microcosms to determine whether environmental factors interact to affect species dynamics and the relative importance of direct and indirect effects on species dynamics. We factorially manipulated three aspects of the environment (temperature, food availability and salinity) and examined reciprocal invasions of competing protist species under each environment. Experimental observations were used to parameterize a dynamic model of the system. Using this model and a novel variance decomposition method, we examined the mechanisms by which environmental changes altered species invasion rates. The three environmental factors interacted when modifying species growth rates, intra‐ and interspecific competition, causing the impact of each environmental change on species dynamics to depend crucially on the state of other aspects of the environment. Indirect changes in the abundance of the resident competitor and its interspecific competitive ability were the main cause of environmental driven variation in invasion rates, whilst direct effects on species intrinsic growth rates were relatively unimportant. This indicates that, to understand and ultimately predict species and community responses to multiple environmental changes, we should consider their joint impacts and the mechanisms by which they interact to modify key ecological processes such as competition.     

Density dependent selection incorporating intraspecific competition 1. A haploid model

A haploid model is introduced and analyzed in which intraspecific competition is incorporated within a density dependent framework. It is assumed that each genotype has a unique carrying capacity corresponding to the equilibrium population size when fixed for that type. Each genotypic fitness at a single multi-allelic locus is a function of a distinctive effective population size formed by adding the numbers of each genotype present, weighted by an intraspecific competition coefficient. As a result, the fitnesses depend upon the relative frequencies of the various genotypes as well as the total population size. Intergenotypic interactions can have a profound effect upon the outcome of the population. In particular, when the density effect of one individual upon another depends upon their respective genotypes, a unique stable interior equilibrium is possible in which all alleles are present. This stands in contrast to the purely density dependent haploid system in which the only possible stable state corresponds to fixation for the type with the highest carrying capacity. In the present model selective advantage is determined by a balance between carrying capacity and sensitivity to density pressures from other genotypes. Fixation for the genotype with the highest carrying capacity, for instance, will not be stable if it exerts a sufficiently weak competitive effect upon the other genotypes. In the diallelic case, maintenance of both alleles at a stable equilibrium requires that the net intragenotypic competition between individuals of like genotype be stronger than that between unlike types. As for purely density regulated systems, there may be no stable equilibria and/or regular and chaotic cycling may occur. The results may also be interpreted in terms of a discrete time model of interspecific competition with each haplotype representing a different species.     

Complex population dynamics and the coalescent under neutrality

Estimates of the coalescent effective population size N(e) can be poorly correlated with the true population size. The relationship between N(e) and the population size is sensitive to the way in which birth and death rates vary over time. The problem of inference is exacerbated when the mechanisms underlying population dynamics are complex and depend on many parameters. In instances where nonparametric estimators of N(e) such as the skyline struggle to reproduce the correct demographic history, model-based estimators that can draw on prior information about population size and growth rates may be more efficient. A coalescent model is developed for a large class of populations such that the demographic history is described by a deterministic nonlinear dynamical system of arbitrary dimension. This class of demographic model differs from those typically used in population genetics. Birth and death rates are not fixed, and no assumptions are made regarding the fraction of the population sampled. Furthermore, the population may be structured in such a way that gene copies reproduce both within and across demes. For this large class of models, it is shown how to derive the rate of coalescence, as well as the likelihood of a gene genealogy with heterochronous sampling and labeled taxa, and how to simulate a coalescent tree conditional on a complex demographic history. This theoretical framework encapsulates many of the models used by ecologists and epidemiologists and should facilitate the integration of population genetics with the study of mathematical population dynamics.     

Modeling the growth of individuals in crowded plant populations

Aims We present an improved model for the growth of individuals in plant populations experiencing competition.Methods Individuals grow sigmoidally according to the Birch model, which is similar to the more commonly used Richards model, but has the advantage that initial plant growth is always exponential. The individual plant growth models are coupled so that there is a maximum total biomass for the population. The effects of size-asymmetric competition are modeled with a parameter that reflects the size advantage that larger individual have over smaller individuals. We fit the model to data on individual growth in crowded populations of Chenopodium album .Important findings When individual plant growth curves were not coupled, there was a negative or no correlation between initial growth rate and final size, suggesting that competitive interactions were more important in determining final plant size than were plants' initial growth rates. The coupled growth equations fit the data better than individual, uncoupled growth models, even though the number of estimated parameters in the coupled competitive growth model was far fewer, indicating the importance of modeling competition and the degree of size-asymmetric growth explicitly. A quantitative understanding of stand development in terms of the growth of individuals, as altered by competition, is within reach.     

Population modelling of Gelidium sesquipedale (Rhodophyta, Gelidiales)

A matrix population model of Gelidium sesquipedale, a commercial agarophyte from the Northeast Atlantic, was developed based on demographic data obtained during two years in a commercial stand of Cape Espichel, Portugal. G. sesquipedale individuals were classified into categories such as life cycle phase, spores, juveniles and adult frond size, because the species vital rates, fecundity, fertility, survival, growth and breakage depend on them. We also exemplify the use of a user-friendly modelling software, Stella, to develop a structured-population model. This is the first time this software has been used to model the demography of seaweed populations. The Stella model developed here behaved very similarly to the matrix model, because of its particular construction, which causes the forcing functions to be discrete rather than continuous. The relative importance of spore recruitment and vegetative growth of new fronds in both population growth and population structure was investigated. Elasticity analysis suggests that vegetative recruitment is the most important demographic parameter controlling population growth together with survival and transitions between juveniles (1–6 cm fronds) and class 1 fronds (6–9 cm fronds). On the other hand, sexual reproduction may, by itself, efficiently control the relative proportion of gametophytes and tetrasporophytes in the population, even though its contribution to recruitment is extremely small. A 40% difference in the growth rates of gametophyte and tetrasporophyte submatrices resulted from natural differences in spore recruitment rates. This revised version was published online in June 2006 with corrections to the Cover Date.     

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.     

A population biological approach to the collective dynamics of countries undergoing demographic transition

Birth rates have been declining in higher-income countries since the middle of the 19th century. A growing number of other countries have entered this demographic transition to lower fertility, as socioeconomic development continues. Analyses of this demographic transition vary widely, but most analyze individual populations in isolation from others, and most come from fields outside the biological sciences. Here, we develop a population biological model of population dynamics in higher-income countries. Individual countries evolve through density-regulated growth, where gradual evolution toward higher population densities boosts productivity (and hence socioeconomic growth) through economics of agglomeration and scale, in turn reducing birth rates. The exchange of technology and capital between countries can further boost productivity gains in any given country, thus contributing to its demographic transition. As a result, countries can down-regulate one another's population growth through mutual improvements in productivity. The model is fitted to time series data on population size, GDP per capita, and birth rates for the United States, the United Kingdom, and France. The metapopulation dynamics are also characterized across a range of parameter values close to the fitted values. This work may help advance population biological approaches to understanding the implications of the fertility demographic transition for modern human populations. This is relevant to developing long-term predictions of the earth's total population size, which must be based upon a model that incorporates underlying mechanisms.     

Mixture models for estimating the size of a closed population when capture rates vary among individuals

We develop a parameterization of the beta-binomial mixture that provides sensible inferences about the size of a closed population when probabilities of capture or detection vary among individuals. Three classes of mixture models (beta-binomial, logistic-normal, and latent-class) are fitted to recaptures of snowshoe hares for estimating abundance and to counts of bird species for estimating species richness. In both sets of data, rates of detection appear to vary more among individuals (animals or species) than among sampling occasions or locations. The estimates of population size and species richness are sensitive to model-specific assumptions about the latent distribution of individual rates of detection. We demonstrate using simulation experiments that conventional diagnostics for assessing model adequacy, such as deviance, cannot be relied on for selecting classes of mixture models that produce valid inferences about population size. Prior knowledge about sources of individual heterogeneity in detection rates, if available, should be used to help select among classes of mixture models that are to be used for inference.     

Species pool size and invasibility of island communities: a null model of sampling effects

The success of alien species on oceanic islands is considered to be one of the classic observed patterns in ecology. Explanations for this pattern are based on lower species richness on islands and the lower resistance of species‐poor communities to invaders, but this argument needs re‐examination. The important difference between islands and mainland is in the size of species pools, not in local species richness; invasibility of islands should therefore be addressed in terms of differences in species pools. Here I examine whether differences in species pools can affect invasibility in a lottery model with pools of identical native and exotic species. While in a neutral model with all species identical, invasibility does not depend on the species pool, a model with non‐zero variation in population growth rates predicts higher invasibility of communities of smaller pools. This is because of species sampling; drawing species from larger pools increases the probability that an assemblage will include fast growing species. Such assemblages are more likely to exclude random invaders. This constitutes a mechanism through which smaller species pools (such as those of isolated islands) can directly underlie differences in invasibility.     

Physiologically structured models – from versatile technique to ecological theory

A ubiquitous feature of natural communities is the variation in size that can be observed between organisms, a variation that to a substantial degree is intraspecific. Size variation within species by necessity implies that ecological interactions vary both in intensity and type over the life cycle of an individual. Physiologically structured population models (PSPMs) constitute a modelling approach especially designed to analyse these size‐dependent interactions as they explicitly link individual level processes such as consumption and growth to population dynamics. We discuss two cases where PSPMs have been used to analyse the dynamics of size‐structured populations. In the first case, a model of a size‐structured consumer population feeding on a non‐structured prey was successful in predicting both qualitative (mechanisms) and quantitative (individual growth, survival, cycle amplitude) aspects of the population dynamics of a planktivorous fish population. We conclude that single generation cycles as a result of intercohort competition is a general outcome of size‐structured consumer–resource interactions. In the second case, involving both cohort competition and cannibalism, we show that PSPMs may predict double asymptotic growth trajectories with individuals ending up as giants. These growth trajectories, which have also been observed in field data, could not be predicted from individual level information, but are emergent properties of the population feedback on individual processes. In contrast to the size‐structured consumer–resource model, the dynamics in this case cannot be reduced to simpler lumped stage‐based models, but can only be analysed within the domain of PSPMs. Parameter values used in PSPMs adhere to the individual level and are derived independently from the system at focus, whereas model predictions involve both population level processes and individual level processes under conditions of population feedback. This leads to an increased ability to test model predictions but also to a larger set of variables that is predicted at both the individual and population level. The results turn out to be relatively robust to specific model assumptions and thus render a higher degree of generality than purely individual‐based models. At the same time, PSPMs offer a much higher degree of realism, precision and testing ability than lumped stage‐based or non‐structured models. The results of our analyses so far suggest that also in more complex species configurations only a limited set of mechanisms determines the dynamics of PSPMs. We therefore conclude that there is a high potential for developing an individual‐based, size‐dependent community theory using PSPMs.     

Resource-mediated demographic variation during the midsummer succession of a cladoceran community

SUMMARY. 1. The midsummer succession of cladoceran zooplankton in Lake Texoma was examined by demographic characterization of lake populations and a parallel set of life table experiments in which the dominant cladoceran species were raised on naturally occurring lake phytoplankton at lake temperatures. Diaphanosoma leuchtenbergianum occurred throughout the study period (mid June-mid August 1981) but Daphnia parvula, D. galeata mendotae and Ceriodaphnia lacustris were rarely detected after 15 July.
2. Analyses of size structure and size-specific fecundity of the cladoceran populations showed that (1) size structure shifted to smaller individuals in mid July, (2) size-specific fecundity was relatively constant in the smaller size classes and was constant or declined in larger females of each species, and (3) size at maturation did not vary significantly during the study of any population.
3. Individual growth rates of the two Daphnia species and Ceriodaphnia in the life table experiments were significantly lower in mid July when populations were declining and when size structure of each population was more dominated by smaller Individuals. The reduction in individual growth rates was short-lived in the life table experiments; growth rates returned to pre-decline levels by the end of July.
4. Two effects of temporary reductions in individual growth rates when size-specific fecundity and size at maturation are constant are (1) that size structure is temporarily shifted to smaller individuals, and (2) the potential rate of population growth is reduced because of delays in reproduction. Because data on size structure and size-specific fecundity are commonly used to evaluate the relative effects of predators and resource changes on populations it is noteworthy that erroneous conclusions about the mechanism of population decline may be reached when individual growth rates vary.     

A mathematical model of population dynamics for Batesian mimicry system

We analyse a mathematical model of the population dynamics among a mimic, a corresponding model, and their common predator populations. Predator changes its search-and-attack probability by forming and losing its search image. It cannot distinguish the mimic from the model. Once a predator eats a model individual, it comes to omit both the model and the mimic species from its diet menu. If a predator eats a mimic individual, it comes to increase the search-and-attack probability for both model and mimic. The predator may lose the repulsive/attractive search image with a probability per day. By analysing our model, we can derive the mathematical condition for the persistence of model and mimic populations, and then get the result that the condition for the persistence of model population does not depend on the mimic population size, while the condition for the persistence of mimic population does depend the predator's memory of search image.     

Simplifying a physiologically structured population model to a stage-structured biomass model

We formulate and analyze an archetypal consumer-resource model in terms of ordinary differential equations that consistently translates individual life history processes, in particular food-dependent growth in body size and stage-specific differences between juveniles and adults in resource use and mortality, to the population level. This stage-structured model is derived as an approximation to a physiologically structured population model, which accounts for a complete size-distribution of the consumer population and which is based on assumptions about the energy budget and size-dependent life history of individual consumers. The approximation ensures that under equilibrium conditions predictions of both models are completely identical. In addition we find that under non-equilibrium conditions the stage-structured model gives rise to dynamics that closely approximate the dynamics exhibited by the size-structured model, as long as adult consumers are superior foragers than juveniles with a higher mass-specific ingestion rate. When the mass-specific intake rate of juvenile consumers is higher, the size-structured model exhibits single-generation cycles, in which a single cohort of consumers dominates population dynamics throughout its life time and the population composition varies over time between a dominance by juveniles and adults, respectively. The stage-structured model does not capture these dynamics because it incorporates a distributed time delay between the birth and maturation of an individual organism in contrast to the size-structured model, in which maturation is a discrete event in individual life history. We investigate model dynamics with both semi-chemostat and logistic resource growth.     

A mathematical model of population dynamics for Batesian mimicry system

We analyse a mathematical model of the population dynamics among a mimic, a corresponding model, and their common predator populations. Predator changes its search-and-attack probability by forming and losing its search image. It cannot distinguish the mimic from the model. Once a predator eats a model individual, it comes to omit both the model and the mimic species from its diet menu. If a predator eats a mimic individual, it comes to increase the search-and-attack probability for both model and mimic. The predator may lose the repulsive/attractive search image with a probability per day. By analysing our model, we can derive the mathematical condition for the persistence of model and mimic populations, and then get the result that the condition for the persistence of model population does not depend on the mimic population size, while the condition for the persistence of mimic population does depend the predator's memory of search image.