Problem sheet 4

Thomas E. Woolley |

Last edited on: |

April 2, 2018

Data on the number of Canadian lynx and snowshoe hair pelts traded where collected by the Hudson Bay trading
company^{1} .
The data are presented in Figure 1(b) and has been used as a proxy for population data.

In this question we consider a mathematical model that has been suggested to describe the features seen in the data. Specifically, it is a predator-prey interaction model called the Lotka-Volterra model. Let L be the lynx population and H be the hare population. The interaction equations are

| (1) |

- 1.
- Describe two features seen in the population data of Figure 1(b).
- 2.
- Name two troubling assumptions behind the Lotka-Volterra interaction equations and suggest how they could be fixed.
- 3.
- Write down the ODEs representing the interaction system.
- 4.
- Let u and v be non-dimensional variables of the hare and lynx, respectively. Non-dimensionalise the system to produce a
system of the form:
= u - uv, (2) = α(uv - v), (3) where the time derivative is with respect to some non-dimensionalised time.

Provide the dimensional scales of the population (i.e. [L], [H] and [t]) as well as the definition of α. Note, that you are not required to show that [L] and [H] have the right dimension and α is dimensionless, but it is a good way to check your working.

- 5.
- What are the steady states?
- 6.
- What is the linear stability of the states?
- 7.
- Consider
(4) Show that equation (4) can be directly integrated to show that the populations must satisfy the constraint

(5) where C is a constant of integration that depends on the initial conditions, which you should specify explicitly. (Hint: rearrange du∕ dv such that one side contains all u terms and the other contains all v terms).

Let α = 1∕2 and simulate equations (2) and (3). Plot the results in the (u,v) plane along with equation (5). Note you will have to choose appropriate initial conditions, specify C in terms of these initial conditions and use an implicit plotting algorithm such as fimplicit in MatLab.

- 1.
- Vary the initial conditions, what do you notice?
- 2.
- Do your discoveries accord with the results from question 1?
- 3.
- Now consider the plot of (t,u) and (t,v). Do the simulated curves match the data seen in Figure 1(b)?

Consider the following set of equations which model the interactions of two populations N_{1} and N_{2}:

Ṅ_{1} | = r_{1}N_{1}, | (6) |

Ṅ_{2} | = r_{2}N_{2}. | (7) |

- 1.
- What dynamics are occurring between the species N
_{1}and N_{2}?Hint 1: by the symmetry of N

_{1}and N_{2}in the equations whatever N_{1}is doing to N_{2}, N_{2}is doing to N_{1}. This means that the dynamics could be mutual creation (mutualism) or mutual destruction (competition).Hint 2: Compare the above equations to regular logistic curve = ru(1 - u∕K). What happens if we increase or decrease K? Thus, what influence does increasing or decreasing N

_{1}have on N_{2}? - 2.
- Use N
_{1}= K_{1}u_{1}, N_{2}= K_{2}u_{2}, t = τ∕r_{1}to non-dimensionalise the equations. Define ρ = r_{2}∕r_{1}, α_{12}= b_{12}K_{2}∕K_{1}and α_{21}= b_{21}K_{1}∕K_{2}. Rewrite the system parameters in terms of α_{12}, α_{21}and ρ. - 3.
- Show that the system has four steady states:
(8) where

(9) What restrictions (if any) do we need to place of the steady states?

- 4.
- Determine the linear stability of the steady states. Taking care to note any bifurcation conditions and how they relate to the restrictions on the steady state existence.

One of the assumptions in the Lotka-Volterra equation is that the predation effect is proportional to both the predator and prey population. However, as the number of prey increases the competition between predators will increase, thus, we consider the adapted equations

= u - uv, | (10) | |

= b(uv - v) - bv^{2}. | (11) |

- 1.
- What are the steady states of the system?
- 2.
- Characterise the stability of the valid steady states, noting any dependences on the parameter b.
- 3.
- How does predator competition change the outcome of the situation, compared to the basic Lotka-Volterra equation shown in question 1?

In 1963, Edward Lorenz developed a simplified mathematical model for atmospheric convection. The model is a system of three ordinary differential equations now known as the Lorenz equations:

ẋ | = σ(y - x), | (12) |

ẏ | = x(ρ - z) - y, | (13) |

ż | = xy - βz. | (14) |

The equations relate the properties of a two-dimensional fluid layer uniformly warmed from below and cooled from above. In particular, the equations describe the rate of change of three quantities with respect to time: x is proportional to the rate of convection; y is proportional to the horizontal temperature variation; z is proportional to the vertical temperature variation. The constants σ, ρ, and β are system parameters proportional to the Prandtl number, Rayleigh number, and certain physical dimensions of the layer itself.

For simplicity, let σ = β = 1.

- 1.
- What are the steady states, noting dependencies of ρ?
- 2.
- Characterise the stability of the zero steady state only, noting dependencies of ρ. Note that you will have to find the eigenvalues of a 3 × 3 matrix. Substituting in the values of the steady state will help you. The non-zero steady states are always either stable nodes or stable spirals when they exist.
- 3.
- Simulations for ρ > 1 and ρ < 1 are illustrated in Figure 2. Do these accord with your findings?

EXTENSION 1: If you are brave enough calculate the eigenvalues corresponding to the non-zero steady states and show that they always have negative real part when ρ > 1, but they may be complex.

EXTENSION 2: If you are even braver rerun the analysis with variables σ = 10, ρ = 28 and β = 3. This causes the system to act chaotically. Categorise the steady states in this case. What happens? An image of the chaotic trajectory can be seen in Figure 3.