Organisms as resources and consumers
Part 2

Overarching questions in Consumer–resource dynamics

  • What allows persistence of consumers and resources in the long term?
    • Are there conditions in which consumers “over-consume” and engineer their own demise?
    • How does perturbation of the system (e.g. removal of a “top predator”) affect the network as a whole?

Review from Wednesday

  • Simple models as “building blocks” for more complicated systems

Consider a system with two species (\(R_1\) and \(R_2\)) that compete with one another for resources, but that are also predated on by the same predator (\(C\)).

Write a model that captures the interactions in this system.

For competing species \(R_1\) and \(R_2\):

\[\frac{dR_1}{dt} = r_1R_1(1-\alpha_{11}R_1 - \alpha_{12}R_2)\] \[\frac{dR_2}{dt} = r_2R_2(1-\alpha_{22}R_2 - \alpha_{21}R_1)\]

Modify these to account for consumer effects:

\[\frac{dR_1}{dt} = r_1R_1(1-\alpha_{11}R_1 - \alpha_{12}R_2)-a_1R_1C\] \[\frac{dR_2}{dt} = r_2R_2(1-\alpha_{22}R_2 - \alpha_{21}R_1)-a_2R_2C\]

The consumer equation also needs to reflect both sources of food:

\[\frac{dC}{dt} = e_1 a_1 R_1 C + e_2 a_2 R_2 C - mC\]

  • In the reductionist tradition, we can use these building blocks to logically think through the types of communities that might exist in nature

Figure from Bell and Fortier-Dubois 2017, Proc B

Such simplification enables predictions, e.g. what should we expect to happen if one of the herbivores is eliminated?

Figure from Bell and Fortier-Dubois 2017, Proc B

Of course, there’s some limit to how complicated a system we can reasonably model

  • We will begin discussing less reductionist approaches next week.

Three-trophic system

  • One producer, an herbivore, and a carnivore

  • e.g. Kelp, sea urchins, and otters

  • Watch the following two videos if you missed Wednesday

How would we model this?

  • Write a simple model of the kelp–urchin–otter system

Assumptions:

  • Kelp experiences exponential growth but gets consumed by urchins
  • Urchins specialize on kelp as their food source, and are consumer by otters
  • Otters primarily eat otters, but they also consume other prey species (e.g. various fish species)

General approach

  • Identify the state variables (in our case, \(K\)elp, \(U\)rchin, \(O\)tter)
  • Determine what assumptions you are willing to make about the system
  • Use our building blocks to start putting together a model

Applying this logic to build a model of infectious disease

Let’s start with some history

In 1926-27, Volterra and Lotka developed the model of consumer–resource dynamics

Kermack and McKendrick’s model

Assumptions

  • In our focal population, individuals are either \(S\)usceptible to a given disease, are currently \(I\)nfected with the disease, or have recently \(R\)ecovered from the disease.
    • (The model is often called the \(SIR\) model for this reason)
  • \(S\)usceptible individuals get infected at a rate \(\beta\) when the encounter individuals that are currently \(I\)nfected
  • \(I\)nfected individuals recover from infection at a rate \(\gamma\)
    • This is equivalent to saying that infected individuals stay infected for \(\frac{1}{\gamma}\) amount of time
  • \(R\)ecovered individuals have immunity from the disease for life
  • This could make sense for a non-fatal disease that moves through quickly
  • E.g. the disease moves quickly enough that we don’t need to think about any births happening in the population
  • We can modify the model for diseases that behave differently (e.g. what if recovered individuals lose their immunity over some period of time?) . . .

flowchart LR
  Susceptible --> Infected --> Recovered

Assumptions

  • In our focal population, individuals are either \(S\)usceptible to a given disease, are currently \(I\)nfected with the disease, or have recently \(R\)ecovered from the disease.
    • (The model is often called the \(SIR\) model for this reason)
  • \(S\)usceptible individuals get infected at a rate \(\beta\) when the encounter individuals that are currently \(I\)nfected
  • \(I\)nfected individuals recover from infection at a rate \(\gamma\)
    • This is equivalent to saying that infected individuals stay infected for \(\frac{1}{\gamma}\) amount of time
  • \(R\)ecovered individuals have immunity from the disease for life

To model this system we need three equations: \(\frac{dS}{dt}\), \(\frac{dI}{dt}\), \(\frac{dR}{dt}\)

\[\frac{dS}{dt} = \beta SI\]

\[\frac{dI}{dt} = \beta SI - \gamma I\]

\[\frac{dR}{dt} = \gamma I\]

What do we expect to happen in this system?

\[\frac{dS}{dt} = \beta SI\]

\[\frac{dI}{dt} = \beta SI - \gamma I\]

\[\frac{dR}{dt} = \gamma I\]

  • The number of susceptible individuals is always declining
  • Eventually, every individual will get infected, and will eventually recover
  • At equilibrium, all individuals will be in the Recovered category

What if the system is more complicated?

  • E.g. individuals lose immunity over time

flowchart LR
  Susceptible --> Infected --> Recovered
  Recovered --> Susceptible

  • E.g. what if there is a vaccine that allows people to gain immunity without becoming infected?

flowchart LR
  Susceptible --> Infected --> Recovered
  Susceptible --> Vaccinated

  • E.g. what if the disease is fatal, such that not all infected individuals recover?

flowchart LR
  Susceptible --> Infected --> Recovered
  Infected --> Dead

  • Some disease move through populations slowly enough that vital rates (birth and daeth rates) become important

    • Are newborns susceptible or immune (innate/maternal immunity)?

  • Some diseases exist in which there are individuals who are carriers but not infected. These carriers don’t get sick but can get others sick.

  • Some diseases have an ‘Incubation’ period after infection (i.e. infected individual doesn’t immediately show signs of infection)

  • Stage/age structure: some diseases cause more effects in older or younger patients – combining matrix model approaches with SIR approach

Next class

  • Some lessons from compartment models for understanding disease dynamics