Environmental carrying capacity - Math Insight

Math Insight.

Environmental carrying capacity. Suggested background. Video introduction. The discrete logistic equation. Environmental limits decrease growth rate. One of the simplest types of discrete dynamical systems describes the exponential growth of a population, where reproduction in each time step is proportional to the number of individuals. In such exponential growth, the population grows faster and faster, continuing to double in size at regular intervals. Exponential growth may be a good model for early stages of population growth, such as we can observe with a model for the initial growth of a bacteria population. However, as one can readily see by exploring some of the problems involving exponential growth, exponential growth cannot continue for long or a population will reach astronomical values. One reason that exponential growth cannot continue is because the environment in which a population lives cannot support an indefinitely large population of a given species. Based on the amount of available food, space, water, and other necessities, an environment will have a finite carrying capacity , which is the maximum population size that the environment can sustain indefinitely. As populations near their environmental limits, the growth rates decrease. The logistic model. One way to develop a model that incorporates a carrying capacity is to make the reproduction be proportional not only to the population size but also to the fraction of unused carrying capacity. The resulting model is the logistic model for population growth. Let $P_t$ denote the population size in time period $t$. Let $M$ be the carrying capacity, which is the population size that the environment can support. Then, $\displaystyle \frac $ is the fraction of the carrying capacity that is used, and $\displaystyle 1 - \frac $ is the unused fraction of the carrying capacity. Because of the double proportionality, the annual increase is proportional to the product the population size and the unused carrying capacity: \begin P_ - P_t = r \, P_t \, \times \left( 1 - \frac \right) \label \end Equation \eqref is the discrete logistic equation. Logistic growth is similar to exponential growth when the population is small. The quantity $r$ is called the low density growth rate and plays the same role as the parameter $r$ defined in exponential growth page. In fact, for small population sizes, we can show that we get the equation for exponential growth back. When $P_t$ is “small” compared to $M$, then the factor $\frac $ is close to zero and the factor $\left( 1 - \frac \right)$ is close to 1. If we replace the factor $\left( 1 - \frac \right)$ with 1, the logistic equation \eqref becomes \begin P_ - P_t \approx r \, P_t. \end which is the same as equation (2) of the exponential growth page. Once the population size $P_t$ becomes appreciable compared carrying capacity $M$, then the growth rate slows down becomes the factor $\left( 1 - \frac \right)$ becomes substantially smaller than 1. As $P_t$ approaches the value of $M$, the factor $\left(1 - \frac \right)$ approaches 0. In this case, the logistic equation \eqref for the population change becomes \begin P_ - P_t \approx r \times P_t \times 0 = 0. \end As the population change $P_ -P_t$ is zero when $P_t=M$, the carrying capacity is an equilibrium of the logistic equation. To be concrete, consider a low density growth rate $r = 0.4$ and a carrying capacity $M = 1000$. Then equation \ref becomes \begin P_ - P_t = 0.4 \times P_t \times \left(1 - \frac \right) \end When $P_t = 10$ (small compared to 1000), then $1 - \frac = 0.99 \approx 1$ and \begin P_ - P_t \approx 0.4 \times P_t \times 1 = 0.4 \times P_t \end The population increases approximately 40% each time period, and there is exponential growth. This growth continues until the population size compared to $1000$ is significant. When $P_t = 500$, $1 - \frac = 0.5$ and \begin P_ - P_t \approx 0.4 \times P_t \times 0.5 = 0.2 \times P_t \end The population is still growing, but only at 20% per time period, one-half of that at low density. When the populations size, $P_t$ reaches 990 (high density, almost to carrying capacity), $1 - \frac = 0.01$ and \begin P_ - P_t \approx 0.4 \times P_t \times 0.01 = 0.004 \times P_t \end and the growth rate has fallen to 0.4% per time period. The following applet shows the solutions to the dynamical system \begin P_0 &= 10 \\ P_ & = P_t + 0.4 P_t \left( 1 - \frac \right) \end in the blue circles and \begin W_0 &= 10\\ W_ & = W_t + 0.4 W_t \end in the red X's. The exponential growth $W_t$ follows the logistic growth $P_t$ until $P_t \approx 200$ and $1 - \frac \approx 0.8$. Logistic and exponential growth. Illustration of how logistic and exponential growth agree for small population sizes and diverge as the population gets larger.

The blue circles show the solution to the dynamical system \begin P_ & = P_t + r P_t \left( 1 - \frac \right), \end and the red X's show the solution to the dynamical system \begin W_ & = W_t + r W_t. \end Both dynamical systems use the same initial condition $W_0=P_0$, which you can change by typing in a value in the box or dragging the purple diamond with your mouse. The carrying capacity $M$ is illustrated by the horizontal line. You can change $M$ and the low density growth rate $r$ by typing in values in the corresponding boxes. You can use the buttons at the top to zoom in and out as well as pan the view. Alternate forms. The form of the logistic equation \eqref that we have studied, \begin P_ - P_t = r \, P_t \, \times \left( 1 - \frac \right), \end gives the change $P_ -P_t$ of the population as a function of the low density growth rate $r$ and the carrying capacity $M$. Such equation is the form of a difference equation, since the change is the difference between the population sizes. Sometimes, it is convenient to write the equation in alternate forms. We could normalize the difference equation to remove the carrying capacity parameter $M$. We could also rewrite the equation by solving it for the new population size $P_t$, putting it in function iteration form. Normalized difference equation. The carrying capacity $M$ is an important quantity as it determines the population size where the population begins to decline. However, the actual value of the carrying capacity doesn't fundamentally change the behavior of the system. No matter what (positive) value of $M$ we start with, we can normalize the equation to remove the dependence on $M$, giving us the same equation but with the number 1 instead of $M$. To perform this normalization, we divide each term of the above equation by $M$ to obtain \begin \frac > - \frac = r \times \frac \times \left(1 - \frac \right). \end Then, we can defined a new variable \begin Q_t = \frac , \end which gives the population size relative to the carrying capacity. $Q_t=1$ means the population has reached carrying capacity; $Q_t=1/2$ means the population is at half the carrying capacity, etc. In terms of $Q_ $, the dynamical system is simply \begin Q_ - Q_t = r \times Q_t \times \left(1 - Q_t \right), \end which is the original equation with the number 1 for the carrying capacity. In this way, we can justify that the carrying capacity doesn't fundamentally change the behavior of the system. For any value of $M$, the rescaled population size $Q_t$ follows the same equation. This normalized equation is simpler as it has only one parameter $r$. When going from the original system to the normalized system, sometimes one doesn't both to change the notation from $P_t$ to $Q_t$. The same variable name, such as $P_t$, can used ambiguously for both the original population size (what we wrote as $P_t$) and the normalized population size (what we wrote as $Q_t$). Function iteration form. The original model \eqref is not the traditional way of representing the discrete time version logistic equation. A more traditional way is to solve the equation for $P_ $ by adding $P_t$ to both sides, yielding \begin P_ = r P_t \left( 1 - \frac \right) + P_t. \end This right hand side of the new equation is still a quadratic equation in $P_t$. We can multiply it out and then refactor the equation to put it into a similar form as the original one. \begin P_ &= r P_t - \frac + P_t\\ &= (r+1) P_t - \frac \\ &= (r+1) P_t \left( 1 - \frac \right) \end. The new version looks a lot uglier. But, we can make it look pretty by defining new parameters. Let \begin R &= r+1\\ K &= M(r+1)/r. \end In terms of these new parameters, the dynamical system in function iteration form looks much better: \begin P_ &= R P_t \left( 1 - \frac \right). \label \end. Equation \eqref is a more traditional form for the logistic discrete dynamical system. Unfortunately, in this form, the parameter $K$ is not the carrying capacity. Instead, it is the carrying capacity multiplied by $(r+1)/r = R/(R-1)$. In equation \eqref , the equilibria are always $P_t=0$ and $P_t=M$. In the new form of \eqref , one has to do some algebra to determine that the equilibria are $P_t=0$ and $P_t=K(R-1)/R$. Normalized function iteration form. One can do the same trick to normalized the new form of the logistic equation. Divide by $K$ to obtain \begin \frac > &= R \frac \left( 1 - \frac \right). \end Define a new variable, the population size as a fraction of $K$. We might as well use the same letter as above even though it is a different quantity, and define \begin Q_t = \frac . \end In terms of $Q_t$, the function iteration form of the logistic dynamical system is \begin Q_ &= R Q_t ( 1 - Q_t). \end.