
MBW:Density Dependence in SingleSpecies PopulationsFrom MathBioThis article is a review of the paper Density Dependence in SingleSpecies Populations by M.P Hassell published in Journal of Animal Ecology, Vol. 44 No. 1, 283295 (1975). by Eric Greenwald and Scott Tatum
ContentsDensity Dependence in SingleSpecies PopulationsAbstractThere have been many attempts to describe density dependent populations throughout history. Most models can only describe portions of the density dependence, such as the high density portion of the population, but have failed to be able to fully describe any data. The objective of this new model it to combine features of previous models to get a "complete" model of density dependence in populations. This article discusses the old models and the drawbacks of each of these, then discusses the new model and its stability, and finally discusses the application and results of the new model. ContextThere have been several models that describe how singlespecies populations grow based on their density. In general populations exhibit very different dynamics for high density and low density populations. This relationship can be seen in the Density Dependent Relationships to the right. Each figure is a plot of vs where is the initial population density and is the surviving population density. Figure 1 below shows two cases from the Density Dependent Relationships figure. Figure 1: Two examples of density dependent populations.^{[1]} As one can see the populations growth is greatly effected by the starting densities. There is a critical density where the kvalues change slope drastically from low to high. HistoryDensity Dependent ModelThe above figures show several density dependent relationships that arise from intraspecific competition for a fixed amount of food. Varley & Gradwell ^{[2]} and Haldane^{[3]} express the mortality as (also known as the kvalue) which is a ratio of the starting density to the surviving density. A simple model to describe a density dependent relationship is given by Varley & Gradwell^{[4]} and Morris^{[5]} is where and b are constants that relate the mortality with the population density. Putting this into logarithmic form and rearranging we can get This model is only good for the linear portions of the data and cannot model the transition between low and high density relationships seen in Figure 1. May^{[6]} examined Equation 1 using a net rate of increase to give the equation where and are population sizes at successive generations. May^{[6]} also examined the stability and showed that it is governed solely by the parameter b. The population is stable when . When the population decays exponentially to equilibrium while will "overcompensate" and exhibit damped oscillations back to equilibrium. If the model predicts increasing oscillations. The model in Equation 3 has a limitation of being linear on on log scale. This creates a problem when the line intercepts the x axis. This would give a threshold population density, (seen in the Density Dependent Relationships figure f). This is only possible when is negative meaning that there is a larger surviving population than the initial population, which would require either reproduction or immigration prior to the sampling of the survivors. A corrected model was developed by Varley, Gradwell & Hassell^{[7]} with a simple correction This model solves the problem of but if the population will oscillate irregularly around the equilibrium with a pattern dependent on the initial condition which is not something that would occur in nature. The expected outcome is one of stable limit cycles around the equilibrium point which can be achieved when "smoothing" is introduced. May^{[6]} created the model where replaces the density dependence term b. The function is defined with the limiting forms Scramble and ContestIntraspecific competition can be defined in two ways: scramble and contest^{[8]}. Scramble is when there is a equal partitioning of resources and there is an abrupt change from from complete survival to 100% mortality. This can be defined as In contest each successful animal get all the resources it requires and unsuccessful animals do not get sufficient resources for survival^{[7]}. This is usually demonstrated by a fixed number of refugees where the surviving number remains constant while the population increases. This leads to a density dependence given by The figure below shows an example of each case just described. Figure 2: Density Dependent Relationships for Extremes: a) Scramble, b) Contest. This model aproaches the logistic equation except for the discontinuity at . Since these two examples are the extremes the normal situation will lay somewhere in the middle with when would represent different amounts of competition and when <N_t > N_c </math> represents a combination of scramble and contest. Figure 1 shows this range between the multiple figures. A model that describes the extremes of scramble has been presented by May^{[6]} as a limiting case of a population of periodical cicadas has the form The mortality rate increases exponentially as increases and has the limiting forms given by Equation 5 of The stability properties were shown by May^{[6]} to be governed by alone. As the value of increases the population exhibits exponential damping, damped oscillations, and stable limit cycles respectively. The drawback of this model is that it does not fit models where density dependence is linear at high population densities. The goal is to combine these above models to get one that will fit linearly at high densities but model the transition around properly. Mathematical ModelHassell proposes a new density dependent model that smooths the transition between across the critical population boundary. This new model, of the following form , has to have the following constraints. To be able to satisfy these conditions as well as constraint that when a logistic function was used by Hassell with the following form.
Where and , the same same as above, are variables representing the net rate of increase and the high density slope, respectively. The new parameter, defines the threshold density so . The Importance of this model change can be seen in the "Example of New Model" figure (right) that the extremes follow the defined limit requirements but the intermediate region is a smooth transitional mix between the two. Stability AnalysisFirst, the equilibrium point is defined where . The stability of this system is determined by looking at the slope of with respect to at the equilibrium point . May et al. showed that when the system experiences exponential damping, the system experiences damped oscillations and when the system has stable limit cycles , where ^{[6]}
The above equation can be plugged into the stability limit to be able to derive a stability relationship between and . These relationships are shown in the figure below. Figure 3: Stability dependence on model parameters ^{[1]} The figure below shows for a fixed a and b when you vary you can see the different types of stability. Figure 4: Plots of N against generation where a=0.01 and b=4 and varying ^{[1]} To see an example of the implications of stability analysis, see this article on a reactiondiffusion model of cancer diffusion. Model ApplicationHassell then applied his model to some data for different insect populations and showed that the model worked for multiple cases even though some did not fit very well. Hassell showcases some of these models in the figure to the right (Application of Model to Data"). The main issue with much of the experimental data is that the range of population densities that are often studied are not broad enough so the data may not cross this critical density boundary. None the less the model can still be used to fit these data fairly well. DiscussionThe primary bonus of this new model is that it has a smooth transition across the critical density. One major weakness of this model is that the abrupt transition that has been observed in some insect types as can be seen in Figure 1 a. This weakness is considered minor because "under field conditions, such sharp discontinuities as seen in Fig. 1 are much less likely" ^{[1]} (see "Density Dependent Relationships" for reference to "Fig 1") The stability analysis also leads to the conclusion that the model result significantly depends on how you approximate the variable because as shown in Figure 3, as long as the stability regime will depend on . Therefore Hassell says that you cant just assume is equal to the fecundity, instead one needs to take into account all of the mortalities during the life cycle. One example proposed by Hassell is when looking at for pupae.
Another useful area that this model can be used is looking at how perturbations. Normally, when a system is perturbed past an equilibrium solution, there tends to be overcompensation but this model can see both the oscillation that occurs with over compensation or see smooth growth to the equilibrium. This is due to the fact that this model can have both low and high density solutions, and that the stability varies for both and . Thus for low values of even with a pertubation we will still be in the exponential damping region thus no oscillations will occur Project categorizationMathematicsThis paper develops a density dependent population equation which is a function of successive generations, various constants, and a parameter referred to as the net rate of increase. For analysis purpose, the author performed stability analysis and formed logistic equations and plots. Model typeThe article starts with a brief history of prior density dependent population models and then combined and corrected them to help accomplish the goal of establishing a "complete" model of density dependence in populations. Thus, the model is composed of important features of previous density dependent models with the ultimate goal of a smooth transition across the critical density boundary. Biological systemThe biological system studied is density dependent populations. Examples include certain insect populations and radioactive decay. Due to the generality of this article, it can be applied to many biological systems. Additional paper discussionThe paper "On the mechanistic underpinning of discretetime population models with complex dynamics" by Stefan A.H. Geritz and Eva Kisdi (Journal of Theoretical Biology 228 (2004) 261–269) cites Hassell's paper. This paper presents the mechanistic basis for various discretetime population models (such as Hassell's) that can produce limit cycles and chaotic dynamics. The authors specifically reference Hassell's model as an example of a logistic model that was used as a descriptive model that can be interpreted mechanistically. The paper begins with the development of a simple resource consumption model with a certain population that is living off a limited food source and producing eggs at a rate proportional to food intake. After looking at the dynamics within a year, the authors derive a model for the betweenyear dynamics, but they transition from continuoustime to discretetime accordingly. They then notice that the betweenyears dynamics are always overcompensated due to a lack of resources coming from outside of the system, thus possibly leading to complex dynamics or even extinction. In order to emphasize the importance of these findings, popular models such as the Beverton and Holt (1975) model, the Ricker model, and the Hassell model are all analyzed and shown to exhibit this exact problem. Finally, their model generalizes to account for multiple resources and consumer species in addition to multiple trophic levels. Density can also effect the populations of multiple species in an environment where species populations are dependent on each other. For an example of how density will effect population in a predatorprey model view MBW:Wolf Territoriality, WolfDeer Interaction and Survival. References
