FW 662 Lecture 1 - Density-independent population models

Text: Gotelli, 2001, A Primer of Ecology

What is a population?

Krebs (1972). A group of organisms of the same species occupying a particular space
at a particular time

Cole (1957). A biological unit at the level of ecological integration where it is

meaningful to speak of birth rate, death rate, sex ratios, and age structure in
describing properties or parameters of the unit.

Gotelli (1998) A group of plants, animals, or other organisms, all of the same species,
that live together and reproduce.

NOAA website (http://coris.noaa.gov/glossary/glossary_l_z.html) - a group of

individuals of the same species living in the same area at the same time and
sharing a common gene pool; a group of potentially interbreeding organisms in
a geographic area

Definitions taken together point to some unique aspects of populations

1. Abundance or Density
2. Have boundaries
3. Change over time
4. Composition (sex and age)
5. Distribution (pattern, scale)

What is "population dynamics"? It’s an area of investigation between the fields of

population biology and population mathematics.

Population Biology Population Dynamics Population Mathematics

Little or no concern for Balance of biology and Little or no concern about
mathematical mathematics biological reality
representation

Some questions we would like to be able to answer after this course:

1. Maximum number of adult females that we can harvest sustainably from a mule deer
population.
2. Number of individuals needed in a re-introduction, e.g., wolves into Yellowstone, to be
fairly certain that the population will persist 100 years.
3. Given data on survival and reproduction, what are the chances that an endangered
population will persist 100 years?
4. How large of a marine reserve is required to prevent the extinction of a harvested fish
species?
FW 662 Lecture 1 - Density-independent population models

Some of the topics that we will cover include

density dependence/independence, exponential growth, logistic growth, age-structured

populations, stage-structured populations, additive vs. compensatory mortality,
stochasticity (demographic, environmental, individual heterogeneity, genetic, sampling),
chaos, metapopulations, predation, competition, herbivory, evolution, natural selection,
genetics, social systems, management of a population, population viability analysis

Population Dynamics is concerned with changes in the density or numbers of organisms

and the processes that cause these changes.

Population Dynamics has relied heavily on mathematical modeling to describe theoretical

relationships and to estimate the parameters needed to describe population change.

Some thoughts on modeling:

Hence, our truth is the intersection of independent lies. Levins, R. 1966. The strategy of model
building in population biology. Am. Sci. 54:421-431.
The bottom line is that models are abstractions of reality. Levins (1966) points out that
modeling is essentially a tradeoff between
1. Generality
2. Realism
3. Precision
FW 662 Lecture 1 - Density-independent population models

The usefulness of any particular model depends on the modeler’s goals. To describe
general ecological principles, it is usually necessary to sacrifice realism and precision.
To describe a particular population, it is usually necessary to sacrifice generality.

Williams et al. (2002, Chapter 7) Describe model usefulness with the following
conceptual model.

Classification of population models based on biological understanding and empirical

strength. I= models incorporating good biological understanding and supported by data.
II= models incorporating good biological understanding but poorly supported by data.
III= models incorporating little biological understanding but supported by data. IV =
model based primarily on speculation and poorly supported by data.

Advantages to mathematical models

1. Equations can be manipulated
This leads to flexibility beyond what is possible with a conceptual model
2. Equations require explicit statements of rules and variables
This leads to more precise statements or hypotheses
3. Equations can result in emergent properties
Such properties are difficult or impossible to predict using non-mathematical
reasoning
FW 662 Lecture 1 - Density-independent population models

The number of adults (N) is some function F.

N(t)=F(t-1) the number alive at time t is some function of the number alive at time t-1.
So what is F???
F=(pB)
Then N(t)=pB(N(t-1))
Then N t = pB( N t −1 ) or N t +1 = pB( N t )

A mathematical aside
Change, difference, and ratio equations
When time periods are discrete we use difference equations to describe change over time. Difference
equations can be expressed different ways to highlight various aspects of the math. Three important
ones:
Change Equation N t +1 = f ( N t )

Difference Equation N t +1 − N t = f ( N t )

N t +1
Ratio Equation = f (Nt )
Nt
FW 662 Lecture 1 - Density-independent population models

Geometric Population Growth

In population dynamics we are generally concerned with four fundamental rates: Births
(B), Deaths (D), Immigration (I), and Emigration (E). We refer to models using these
parameters as BIDE models. Change in abundance over time is modeled as a function of
how many individuals are born, die, leave, or enter the population

N t +1 = N t + Bt − Dt + I t − E t

Initially, we will consider closed populations (i.e. no immigration or emigration) and

non-overlapping generations (or change is measured in discrete intervals). Biologically,
this means pulsed reproduction and finite survival.

N t +1 = N t + Bt − Dt
we need to convert the raw Birth and Death rate to per capita rates (per individual)
because per capita rates are easier to estimate

Bt Dt
bt = and dt =
Nt Nt
Substituting
N t +1 = N t + bt N t − d t N t
Now we assume that per capita birth and death rates are constant over time
N t +1 = N t + bN t − dN t
Simplifying this equation
N t +1 = N t + (b − d )N t
The term b-d is the geometric rate of increase and is given its own symbol R
N t +1 = N t + RN t
N t +1 − N t = RN t (a difference equation)
∆N t = RN t
∆N t
=R
N
The parameter R represents the per capita rate of change in the size of the population

Lets go back to N t +1 = N t + RN t and factor out Nt

FW 662 Lecture 1 - Density-independent population models

N t +1 = (1 + R )N t
In population ecology 1+R is given its own symbol lamda

(1 + R ) = λ
N t +1 = λN t a change equation

N t = λN t −1 = λ2 N t − 2 = λ3 N t −3 ... = λt N 0 this is called a Taylor Series

N t = λt N 0

using this equation there are three possibilities for

population growth.

Geometric Population Growth

lam=1.2
lam=1
lam=0.8
N

Time

Lamda can also be expressed as ratio.

N t +1
λ= a ratio equation
Nt

Hence we can calculate lamda without knowing the per capita birth and death rates;
however, in actuality it will probably be important to estimate b and d.
FW 662 Lecture 1 - Density-independent population models

Exponential Population Growth

Now let’s suppose that we have continuous population growth (or change is measured
over a very small interval). Biologically, this means continuous reproduction and
instantaneous survival rates. This requires a differential equation.

dN
= bN − dN
dt
Let b-d=r
dN
= rN
dt
r is called the instantaneous or intrinsic rate of increase

With integration we can get a predictive equation for N

dN
Put in terms of N = rdt
N
dx
from calculus = d ln(x )
x

so d ln( N ) = rdt
t t
Integrate ∫ d ln( N ) = ∫ rdt
0 0

ln( N t ) − ln( N 0 ) = rt
take exponential of both sides

e ln( N ( t )) e − ln( N ( 0 )) = e rt recall that e ln( a ) = a and e −ln( a ) =

1
a
Nt
= e rt
N0
Solve for N(t) and we get a continuous model that predicts population size

N t = N 0 e rt
FW 662 Lecture 1 - Density-independent population models

Using this equation there are three possibilities for

population growth.

r=-0.3
r=0
N r=0.3

Time

We have derived two predictive equations for population growth.

Discrete Continuous

Density Independent N t = λt N 0 N t = N 0 e rt

Comparing the two models.

λ = e r or ln(λ ) = r

if r is small then we can use a Taylor series to find the approximate relationship

r2 r3
er = 1 + r + + + ...
2! 3!
λ ≅ 1+ r
The two models diverge as r gets large and as the number of generations (t) increases.
Consider the following example with R=0.2=r.
FW 662 Lecture 1 - Density-independent population models

R=r 0.2

N t = N 0 (λ ) t
N = N e rt
t 0
t
0 10.00 10.00
1 12.00 12.21
2 14.40 14.92
3 17.28 18.22
4 20.74 22.26
5 24.88 27.18
6 29.86 33.20
7 35.83 40.55
8 43.00 49.53
9 51.60 60.50
10 61.92 73.89

The distinction between R (or λ) and r is important biologically, i.e. pulse versus
continuous reproduction, instantaneous versus finite survival rates. There is confusion in
the literature between r, R, and λ, so be careful to understand the author’s definition.
FW 662 Lecture 1 - Density-independent population models

Some Useful References

Begon, M., J. L. Harper and C. R. Townsend. 1996. Ecology: Individuals, populations

and communities. 3rd ed. Blackwell Scientific Ltd., Cambridge, Mass. 1068pp.

Caughley, G. 1977. Analysis of vertebrate populations. John Wiley and Sons, New York.
pp. 1-7.

Cushing, D. H. 1981. Fisheries biology: a study of population dynamics. University of

Wisconsin Press, Madison. (on reserve at library)

Donovan, T.M. and Welden C.W. Spreadsheet exercises in Ecology and Evolution.
Sinaueur Assoc., MA.

Hastings, A. 1996. Population Biology: concepts and models. Springer-Verlag Inc.

New York

Williams, B.K., J.D. Nichols, M.J. Conroy. 2002. Analysis and management of animal
populations. Academic Press. New York.