Problem 8 - Biocalculus Calculus for the Life Sciences

Question

8\. Metapopulations Consider a simple metapopulation in which subpopulation \(\mathrm{A}\) grows at a per capita rate of \(r_{\mathrm{A}}=1\) and subpopulation \(\mathrm{B}\) declines at a per capita rate of \(r_{\mathrm{B}}=-1 .\) Suppose the per capita rate of movement between subpopulation patches is \(m\) in both directions. This gives $$\begin{aligned} \frac{d x_{\mathrm{A}}}{d t} &=(1-m) x_{\mathrm{A}}+m x_{\mathrm{B}} \\ \frac{d x_{\mathrm{B}}}{d t} &=-(1+m) x_{\mathrm{B}}+m x_{\mathrm{A}} \end{aligned}$$ where \(x_{\text { A }}\) and \(x_{\mathrm{B}}\) are the numbers of individuals in patches \(\mathrm{A}\)and \(\mathrm{B},\) respectively. $$\begin{array}{l}{\text { (a) Classify the equilibrium at the origin. }} \\\ {\text { (b) Find the general solution. }} \\ {\text { (c) What is the solution to the initial-value problem if }} \\ {x_{\mathrm{A}}(0)=1 \text { and } x_{\mathrm{B}}(0)=0 ?}\end{array}$$

Step-by-Step Solution

Answer

The origin is a saddle point for some values of \( m \), and the solutions to the differential equations involve exponential functions derived from characteristic roots with initial conditions providing specific constants.

Step 1: Finding the Equilibrium at the Origin

To classify the equilibrium at the origin, we analyze the system of differential equations: \( \frac{d x_{\mathrm{A}}}{d t} =(1-m) x_{\mathrm{A}}+m x_{\mathrm{B}} \) and \( \frac{d x_{\mathrm{B}}}{d t} =-(1+m) x_{\mathrm{B}}+m x_{\mathrm{A}} \). The equilibrium points occur when both derivatives are zero. Setting \( \frac{d x_{\mathrm{A}}}{d t} = 0 \) and \( \frac{d x_{\mathrm{B}}}{d t} = 0 \), and substituting \( x_{\mathrm{A}} = 0 \) and \( x_{\mathrm{B}} = 0 \), we find that the origin (0,0) is an equilibrium point. The nature of the equilibrium can be classified by examining the eigenvalues of the Jacobian matrix evaluated at this point.

Step 2: Calculating the Jacobian Matrix

The Jacobian matrix \( J \) of the system is computed from the partial derivatives of the functions with respect to \( x_{\mathrm{A}} \) and \( x_{\mathrm{B}} \). This gives: \[ J = \begin{pmatrix} 1-m & m \ m & -(1+m) \end{pmatrix} \] Now, we calculate the eigenvalues of this matrix to classify the equilibrium.

Step 3: Determining Eigenvalues

The eigenvalues of \( J \) are found by solving the characteristic equation \( \det(J - \lambda I) = 0 \), where \( I \) is the identity matrix. For our matrix, this becomes: \[ \det \begin{pmatrix} 1-m-\lambda & m \ m & -(1+m)-\lambda \end{pmatrix} = (1-m-\lambda)(-(1+m)-\lambda) - m^2 = 0 \] This quadratic equation simplifies to: \( \lambda^2 + (2-m)\lambda + 1 = 0 \). Solving for \( \lambda \), we find: \( \lambda = \frac{-(2-m) \pm \sqrt{(2-m)^2 - 4}}{2} \). The solutions provide insight into the stability of the origin.

Step 4: Classifying Stability

The nature of the equilibrium depends on the eigenvalues: if they both have negative real parts, the equilibrium is stable; if at least one has a positive real part, it's unstable. Solving the quadratic equation, we find that the eigenvalues can have positive, negative, or complex parts based on \( m \), determining the stability at the origin. For smaller \( m \), the equilibrium tends to be a saddle point or unstable.

Step 5: Finding the General Solution

The general solution for \( x_{\mathrm{A}}(t) \) and \( x_{\mathrm{B}}(t) \) is expressed as a linear combination of eigenvectors and exponential terms derived from the eigenvalues found earlier. Assume a solution of the form \( x(t) = c_1 v_1 e^{\lambda_1 t} + c_2 v_2 e^{\lambda_2 t} \), where \( v_1 \) and \( v_2 \) are eigenvectors corresponding to \( \lambda_1 \) and \( \lambda_2 \). This gives full solutions for both populations.

Step 6: Solving the Initial-Value Problem

To solve the initial-value problem with \( x_{\mathrm{A}}(0) = 1 \) and \( x_{\mathrm{B}}(0) = 0 \), substitute these initial conditions into the general solution found in Step 5 to determine specific constants \( c_1 \) and \( c_2 \). You solve the system of equations resulting from substituting \( t = 0 \) in the general solution and equating them to initial conditions.

Key Concepts

Equilibrium Points

Understanding equilibrium points is essential to analyzing how different systems behave over time. In the context of metapopulations, equilibrium points represent the state where the populations of patches A and B do not change.
To find the equilibrium, we look at the system of differential equations provided and set the time derivatives to zero. This means finding values of \(x_A\) and \(x_B\) such that \(\frac{d x_{A}}{dt} = 0 \) and \(\frac{d x_{B}}{dt} = 0\). In the given exercise, this leads us to the origin (0,0) as an equilibrium point.
The origin being an equilibrium point means if both subpopulations start at zero, they will remain stable without external influences or changes. To further understand these points, we need to analyze their stability through Jacobian matrices and eigenvalues.

Jacobian Matrix

The Jacobian matrix plays a crucial role in determining the stability of equilibrium points in systems described by differential equations.
In our example of metapopulations, we construct the Jacobian matrix by finding the partial derivatives of the functions involved with respect to each variable.
The Jacobian for this system, calculated from the equations:\[ J = \begin{pmatrix} 1-m & m \m & -(1+m) \end{pmatrix} \]represents local behavior near the equilibrium point. It captures how the system responds to small changes in the populations of patches A and B.
To understand and classify this behavior, we then calculate the eigenvalues of the Jacobian matrix. These will help in conducting a stability analysis by showing whether the system diverges or returns to equilibrium when perturbed.

Stability Analysis

Stability analysis is key to predicting how systems behave when they're slightly disturbed. In our metapopulation model, we use eigenvalues of the Jacobian matrix to assess the stability of the equilibrium point at the origin.
We solve for the eigenvalues using the characteristic equation derived from the Jacobian:\( \det(J - \lambda I) = 0 \).This resolves into a quadratic equation, where solutions \( \lambda_1 \) and \( \lambda_2 \) inform us about system behavior: - If both \( \lambda_1 \) and \( \lambda_2 \) have negative real parts, the equilibrium is stable; the system returns to the origin after disturbances.- If at least one has a positive real part, the equilibrium is unstable; the system tends to diverge.- Complex eigenvalues indicate oscillating dynamics.
This analysis, crucial in systems biology and ecology, reveals how changes in parameters affect the balance between subpopulations.

Differential Equations

Differential equations are vital tools for modeling how systems change over time. In our example, they describe the dynamics of populations in subpopulations A and B.
These equations show rates of change based on per capita growth rates and migration between patches, represented as \(\frac{d x_{A}}{dt} = (1-m)x_A + mx_B\) and \(\frac{d x_{B}}{dt} = -(1+m)x_B + mx_A\). Differential equations enable us to formulate general solutions predicting population behavior. This involves methods such as finding eigenvectors and exponential terms that describe how populations evolve based on initial conditions.

• General solutions are combinations of terms derived from the system's eigenvalues and eigenvectors.
• These allow us to solve initial-value problems, determining future population dynamics given starting conditions.

By understanding these equations, we learn how multiple interacting populations might fluctuate over time under different conditions.