Computer Experiments
in Population Ecology (XGROW)

Self-Study Exercises


BIO150Y Home Page ~

This document is a copy of Chapter 8: Computer Experiments in Population Ecology in Studies in Evolution, Ecology and Behaviour: A Laboratory Manual for BIO150Y (Eighth Edition), Departments of Botany and Zoology, University of Toronto, 1997, by Corey A. Goldman.

Objectives:

About the XGROW Computer Program

Schedule:

Preparation:


Contents


Introduction

As surely as individuals are born and die, populations grow and decrease. At any time in its history the size of a population is the product of births and deaths, and immigration and emigration. The rate at which individuals enter the population (births and immigration) and individuals leave (deaths and emigration) will determine whether the population is growing, shrinking, or is stable, and if it is changing, how fast.

Logistic Growth

Populations seldom undergo unrestricted growth under constant conditions for long periods. Given the capacity for increase of bacteria, this is just as well! More commonly, as populations grow individuals affect one another's access to vital resources such as food and space. This is an example of intraspecific competition (intra = within). As a consequence, birth rates tend to decrease and death rates increase as populations become more crowded -- so-called density-dependent effects on birth and death rates. At some point, the carrying capacity (K) of the environment is reached by the population and there is no more growth (i.e., births balance deaths). The carrying capacity is the number of individuals which can be supported by an environment.

Interspecific Competition

Populations do not live alone, just as individuals do not. When individuals of the same species, or of two different species, depend on a common important resource then competition occurs. Competition can be defined as "interactions between individuals brought about by a shared requirement for a resource in limited supply leading to a reduction in survivorship, growth, and reproduction of individuals." Interspecific competition (inter = between) is the competition between two or more different species for a resource.

Predation

Predation is an important interaction between species that affects population levels of species. Predation can be defined as occurring when individuals of one population eat living individuals of another. There are several models for predation, but all models have two parts: the birth rate of the predator increases as the number of prey increases, and the death rate of the prey increases as the number of predators increases.


Part A: About XGROW

XGROW is a program which allows you to simulate

The underlying mathematical models for XGROW are detailed in the appendix and in the online Help.

There are two different versions of XGROW available (there are minor differences between the versions, mostly in the look of the interface): XGROW for Windows is recommended because it will run faster on most computers; XGROW on the Web, a Java applet, requires an Internet connection, and can often be slow.

Using XGROW

The XGROW for Windows screen is divided into two parts: the graph window on the left and the program controls (buttons and sliders) on the right. From the function bar above the graph window you can: Print a graph or Export a data file, Exit the program, change a species in the Single Species model, and obtain online Help.

Graphs: The graphs show population size (N or lnN) as a function of time. In the Single Species model you can choose to graph population size as arithmetic or natural logarithmic (ln) values of population size by selecting the appropriate button in the Graph Type box; on the arithmetic graphs, the magnitude to which population size is expressed is given at the top of the y-axis (e.g., 000's refers to thousands). The Competition and Predation models display only natural logarithmic values of population size.

Models: Three models can be selected: Single Species, Competition, and Predation. The program starts each time with the Single Species model, with the starting species selected at random. Select Species from the pull-down menu above the graph and click+drag to a new species.

Controls: The Control panel has three functions: Start, Clear, and Reset. "Start" begins your experiment (simulation). "Clear" erases lines from the graphs. "Reset" returns the parameter sliders to the default values (in XGROW, an experiment is set up for you by default). Five simulations can be displayed on a graph at one time, each with a different coloured line. If a sixth simulation is run then the line generated by the first simulation is erased, if a seventh simulation is run then the line generated by the second simulation is erased, etc.

Time: For each species in the Single Species model the time scale along the x-axis (i.e., duration of the experiment) can be varied from the default values. This enables you to "zoom in" on the graph line.

Noise: To more closely simulate population growth, the Noise function randomly fluctuates birth and death rates. The default setting is "ON"; to turn "OFF" select Options from the function bar above the graph window.

Time: To more closely simulate population growth, the Time Delay function builds each arithmetic graph over time. The default setting is "ON"; to turn "OFF" select Options from the function bar above the graph window. The Time Delay function is "OFF" for logarithmic graphs.

Grid: To assist in determining the x,y-coordinates of a data point directly from a graph line, the Grid function superimposes grid lines on the graph.

Population Characteristics

In each model, particular characteristics of the population (parameters) can be varied, such as birth rate, death rate and the intensity of competition. In each model, a default experiment is set up initially. To change a parameter, click+drag the appropriate slider until the desired value is obtained. Alternatively, click on the variable box to the right of the slider to activate the box and type in the desired variable. Note that each parameter has maximum and minimum values that cannot be exceeded. The parameters are explained below:

Single Species | Interspecific Competition | Predation

Single Species


Interspecific Competition


Predation

Prey:

Predator:


Part B: Self-Study Exercises

You are to conduct computer experiments (simulations) with each of the three models:

Your understanding of these exercises with be evaluated on Quiz 6 (in Lab 8) and Test 4; you will not be evaluated on the equations in the appendix.

1 - Single Species Population Growth

XGROW always starts with a simulation in the Single Species Model already set up for you, and ready to run. It automatically selects one of the five species which you can study (human, rat, fruit fly, Paramecium caudatum, Paramecium aurelia, and T-phage virus). You can of course choose the species you want (from the Menu option bar, select Species and click+drag to select the species you want to study).

The simplest form of population growth results if individuals reproduce and die at a constant rate. As you will learn, these conditions seldom hold for long, if ever, in reality, but they are a good way to introduce basic population parameters: individual birth and death rates (b0 and d0), and the intrinsic rate of growth, r. In the Single Species model, b0 and d0 (not r per se) are given; these values cover the realistic range observed in nature for the different species.

  1. A default species has been assigned to you and a default simulation has also been set up -- a population of individuals which give birth at the maximum rate for the species and die at the minimum rate for the species. Record the simulation settings in Table 10-1. Record the species name, b0, and d0.

    Table 10-1.
    Simulating unrestricted population growth (r* is estimated from your graph).
    Species                 
    b0  
    d0  
    r    
    r*    
    ______________________ _____ _____ _____ _____
    ______________________ _____ _____ _____ _____

  2. Point+click on Start to run the default simulation. Observe the results: arithmetic values of population size as a function of time.

What is r?

One way to summarize the population growth is the intrinsic rate of growth, r, which is the per capita growth rate.

  1. View the logarithmic graph of the simulation by selecting the Logarithmic Graph Type. Note how the curve is now a straight line; this enables you to estimate the rate of population growth, r, as the slope of the line. From the graph, estimate r and enter your estimate in Table 10-1. (The units of r must also be reported. Since r is a rate, it must be reported with real time units -- the specific time units vary with the species. For humans, r is normally reported per year.)

  2. How does your observed estimate of r from the graph compare with the theoretically expected growth rate obtained from the mathematical relationship (i.e., r = b0 - d0 )? To confirm your estimate, use the actual birth and death rates to calculate the exact r for your species. Complete Table 10-1.

    How similar are the two r values?

How does r depend on b0 and d0?

  1. How does variation in basic birth and death rates influence population growth rate and population size? For example, in historic times the basic death rate of human individuals has decreased as medicine and technology have improved health. Simulate what consequences this has for growth rates and population size.

  2. For comparison, leave the results of the default simulation on the screen.

  3. Decrease d0. Record your new settings in Table 10-1. Run the new simulation.

  4. Compute the new r and record it in Table 10-1. Experiment some more until you are clear about the relationship between birth and death rates and population growth rates.

Intraspecific Competition

In the previous experiments, population growth was unchecked and as you observed, population grew at the theoretical maximum rate for the species. These conditions rarely hold for long for any species. Individuals normally interact with each other, often in a negative way -- intraspecific competition. As population density increases, competition for resources is likely to increase. A consequence of competition, is that individuals can no longer reproduce at the maximum rate or die at the minimum rate. What are the consequences of such density-dependence for the growth of a population?

What is the carrying capacity, K?

  1. Point+click on Reset to restore the default values for the species. Run the default. (We recommend that a species other than human be used here.)

  2. Now, increase the density-dependent birth and death rates (kb and kd) from zero. Record the new settings in Table 10-2. Run the simulation.

    What type of population growth is there now (i.e., what is the shape of the population growth curve)?

    Table 10-2.
    Simulating density-dependent population growth (K* is the observed value estimated from your graph).

    Species

    b0

    d0

    kb

    kd

    K

    K*

    ______________________ _____ _____ _____ _____ _____ _____
    ______________________ _____ _____ _____ _____ _____ _____

  3. Identify the carrying capacity, K, which is the point at which the population stops growing and levels off. From the graph, read off your estimate of K from the y-axis and record it in the above table. The K should be reported in real numbers.

  4. Since you know the values of kb and kd, you can compute the expected carrying capacity for your simulated population using equation 7 (in the appendix). Compute the expected K value and record it in Table 10-2.

How does K depend on the intensity of intraspecific competition?

What happens if you increase the intensity of intraspecific competition? For instance, suppose that habitat conditions change so that food levels decreased for a population.


2 - Interspecific Competition

In the previous simulations you examined how density-dependent effects (intraspecific competition) result in logistic population growth in a single species. But in nature, a species never lives alone. Thus, in addition to intraspecific competition, a population can expect to experience interspecific competition which will affect its growth.

  1. Select the Competition Model. XGROW simulates competition between any two species (Paramecium caudatum and P. aurelia are simply given as examples). The underlying model of how and affect individual growth is given in equations 10 and 11 in the appendix.

    Note the new population characteristics: (1) Birth and death rates have been replaced by r. (2) The density coefficients have been replaced by K. (3) The intensity of competition from the other species is given as for P. aurelia and for P. caudatum.

  2. Point+click on Start to run the default simulation which simulates growth of each Paramecium species alone (i.e., no interspecific competition). Confirm for yourself that both populations reach their maximum population size (K= 1000, or 6.9 on the ln-scale). Record the carrying capacity (K alone) in Table 10-3.

  3. Now simulate a situation where both species are coexisting in the same habitat and where there is interspecific competition. Let us assume that the resources used by an individual of P. caudatum are the same as the resources used by an individual of P. aurelia. To simulate this, you must increase the competition coefficients ( and ) for both species to 1.0.

  4. Run the new simulation. Your new results will be displayed on the same graph. Estimate the stable, equilibrium population size (N) for each species when they live in the same habitat. Record your results in Table 10-3.

    Table 10-3.
    Population growth of two species with and without interspecific competition.

    Simulation 1
    Simulation 1
    Simulation 2
    Simulation 2
    Simulation 3
    Simulation 3
    No interspecific competition
    K alone
    Interspecific competition
    N mixed
    Interspecific competition
    N mixed
    P. aurelia
    = 0
    _____
    = 1.0
    _____
    =
    _____
    P. caudatum
    = 0
    _____
    = 1.0
    _____
    =
    _____

  5. Now, simulate a situation in which one species experiences strong interspecific competition relative to the other species. Leave the carrying capacities for each species set at K = 1000. (This simplifies the interpretation of the different outcomes of interspecific competition outlined in the appendix, since the ratio of the Ks becomes equal to 1.)

  6. Set the competition coefficient for P. aurelia at > 1.0 (i.e., aurelia experiences relatively more inter- than intra-specific competition), and the coefficient for P. caudatum at < 1.0. Record your settings in the table above.

  7. Run the new simulation. Record final population sizes in the above table.

    You have just simulated an experiment like that which Gause conducted with real protozoans and which led to his formalization of the Competitive Exclusion Principle.

  8. Reverse the competitive effects so that now P. caudatum experiences relatively stronger interspecific competition.

    Probably many interspecific interactions are unequal, but relatively weak. And the environment is more favourable for one species than the other, in which case the carrying capacity, K, for each species differs.

  9. Adjust the settings to simulate an environment which can theoretically support about 1.5 times as many individuals of one species than the other. For example, you could set aurelia at K = 1500 and leave caudatum at K = 1000. Run the simulation without interspecific competition to check that each species attains a different K.

  10. Now add in weak competition; for example, = = 0.5 for both species. Run the new simulation.

    Table 10-4.
    Population growth of two species which differ in both the intensity of intra- and interspecific competition.

    P. aurelia

    P. aurelia

    P. caudatum

    P. caudatum

    Simulation

    K

    K

    Result

    ___________ _________ _________ _________ _________ _________
    ___________ _________ _________ _________ _________ _________
    ___________ _________ _________ _________ _________ _________
    ___________ _________ _________ _________ _________ _________

  11. Now, increase the strength of interspecific competition experienced by one species. For instance, increase to 1.0. Run a simulation and check the results. Now what is the outcome?

  12. Next, run a new simulation in which aurelia is experiencing more interspecific competition. Reset = 0.5, and set = 1.0. What is the outcome this time? Evidently, in some circumstances, strength of interspecific competition alone is not the only factor controlling the final outcome.

  13. Experiment with to find a value beyond which aurelia is excluded by caudatum.

  14. Confirm that the conditions for coexistence when the Ks differ holds for your simulations.

You have just seen with these simulations that the outcome of a two-species interaction depends on the interplay between interspecific competition and intraspecific competition for environmental resources (i.e., carrying capacities).


3 - Predation

In a situation where a predator population is feeding on a prey population, the two most important factors influencing the population sizes of each species are the intrinsic growth rate of the prey population and the predator death rate.

  1. Select the Predation Model. XGROW simulates lynx preying on snowshoe hares, but these organisms are just examples and the model could apply to any two-species predator-prey interaction; the simulation is based on equation 14 in the appendix.

    Note some new population characteristics for the predator: (1) Satiation -- predators can get full eating prey. (2) Efficiency -- you can vary the efficiency rate at which prey are converted to predator births. (3) Prey capture rate -- you can vary how many prey a predator can catch per day.

  2. Run the default simulation in which both predator and prey populations are able to coexist in the same habitat for a long time. Estimate the average population size each species has over the long term and record these values in Table 10-5.

    Table 10-5.
    Simulating population growth of predator and prey (N = population size attained).
    Prey
    Prey
    Prey
    Predator Predator
    Simulation

    r

    K
    N

    do

    N
    Default _____ _____ _____

    _____

    _____

    _______ _____ _____ _____

    _____

    _____

    _______ _____ _____ _____

    _____

    _____

    What happens to a stable predator-prey relationship such as you have just simulated if the prey's growth rate is perturbed? Suppose, for example, that the hare's basic birth rate is increased.

  3. Increase the prey growth rate, r, a small amount and run the new simulation. The new results (to simulate increased birth rate) will be displayed on the same graph, so that you can compare both simulations. Record the new long-term average population sizes reached by predator and by prey.

    Now simulate a different scenario. Assume for example, that there is increased hunting pressure by humans on the predator species. As a consequence, predator death rate goes up.

  4. Point+click on Reset and run the default simulation again.

  5. Now increase predator death rate. (Do not increase it by too much to begin with.) Keep a record of your simulations in the above table.

    You have just completed some simulations which demonstrate the Volterra Principle, named after the scientist who first described it, and which says that "if predator and prey are in stable coexistence, change in prey growth rate will affect predator population size more and changes in predator death rate will affect prey population size more."

    What happens to the stable predator-prey relationship if the quality of the prey's environment is improved. Imagine, for example, that a group of well-meaning citizens decide to supplement the natural diet of wild rabbits with commercial rabbit chow. In effect, they are increasing the K for the prey. Simulate what could happen to prey and predators.

  6. Run the default simulation for predator-prey.

  7. Now, conduct a series of simulations in which you increase prey K by two times, three times, etc.

    Does food supplementation have any effect on the long-term average population size for the prey?

  8. Look at the fluctuations in population size of each species. How does this change with increasing prey K?

  9. If prey K is increased dramatically, what is the outcome?

    These simulations have demonstrated a phenomenon known as the "paradox of enrichment." If the carrying capacity of the prey population is increased, what was previously a stable coexistence with its predator species can become unstable and actually result in extinction of either species. In our hypothetical example, the best intentions of the citizens could have back-fired. You will find more details about this paradox in the appendix.


Review Questions

  1. What does an r-value of 0.5/individual/month mean (explain what the number means)?

  2. If the birth rate increases but death rate remains the same in a population what happens to r? What is the outcome in terms of population growth?

  3. If birth and death rates are constant through time, what sort of population growth do we get?

  4. If there are density-dependent factors acting in a population, what sort of population growth do we get?

  5. At the carrying capacity, population growth stops. Does this mean that the birth and death rates equal zero?

  6. Give some examples of types of density-dependent factors.

  7. Can two species coexist if there is interspecific competition between them?

  8. What is the Competitive Exclusion Principle?

  9. Species A has an = 1.4, Species B has an = 0.7. Which species is experiencing the most interspecific competition?

  10. Biologists want to predict the population size of Canada in the year 2050. What data would they need to obtain about the population in order to use population growth models? How would they collect such data?


Glossary

Complete this glossary by adding the definitions for these terms:

Birth rate:

Carrying capacity:

Competition coefficients:

Death rate:

Density dependent effects:

Exponential growth:

Interspecific competition:

Intraspecific competition:

Intrinsic growth rate:

Logistic growth:

Predation:


Appendix: Mathematical Models of Population Growth


[ BIO150Y Home Page ]

Copyright © 1997 University of Toronto. All rights reserved.
Comments to BIO150Y staff.