Your browser (Internet Explorer 8 or lower) is out of date. It has known security flaws and may not display all features of this and other websites. Learn how to update your browser.

X

Menu

Populations

Day 22: Isle Royale: Predator/Prey Model for Moose and Wolves

Jacob Englert
This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale. It was "cloned" from a model that InsightMaker provides to its users, at
https://insightmaker.com/insight/2068/Isle-Royale-Predator-Prey-Interactions
Thanks Scott Fortmann-Roe.

I've created a Mathematica file that replicates the model, at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/Moose-n-Wolf-InsightMaker.nb

It allows one to experiment with adjusting the initial number of moose and wolves on the island.

I used steepest descent in Mathematica to optimize the parameters, with my objective data being the ratio of wolves to moose. You can try my (admittedly) kludgy code, at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/Moose-n-Wolf-InsightMaker-BestFit.nb

{WolfBirthRateFactorStart,
WolfDeathRateStart,
MooseBirthRateStart,
MooseDeathRateFactorStart,
moStart,
woStart} =
{0.000267409,
0.239821,
0.269755,
0.0113679,
591,
23.};

Environment Ecology Populations Math Modeling Mat375

  • 1 year 8 months ago

Clone of Final Midterm Student version of A More Realistic Model of Isle Royale: Predator Prey Interactions

Allison Zembrodt
Allison Zembrodt's Model
This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale.

We incorporate logistic growth into the moose dynamics, and we replace the death flow of the moose with a kill rate modeled from the kill rate data found on the Isle Royale website.

I start with these parameters:
Wolf Death Rate = 0.15
Wolf Birth Rate = 0.0187963
Moose Birth Rate = 0.4
Carrying Capacity = 2000
Initial Moose: 563
Initial Wolves: 20

I used RK-4 with step-size 0.1, from 1959 for 60 years.

The moose birth flow is logistic, MBR*M*(1-M/K)
Moose death flow is Kill Rate (in Moose/Year)
Wolf birth flow is WBR*Kill Rate (in Wolves/Year)
Wolf death flow is WDR*W

equations I used in kill rate :

power model - 12*0.1251361120909615*([Moose]/[Wolves])^.44491970277839954*[Wolves]


Kill rate sqrt = 12*(0.0933207+.0873463*([Moose]/[Wolves])^.5)*[Wolves]


Holling Type III - ((0.986198*([Moose]/[Wolves])^2)/ (601.468 +([Moose]/[Wolves])^2))*[Wolves]*12


linear - 12*[Wolves]*(.400271+.00560299([Moose]/[Wolves]))


Environment Ecology Populations Midterm

  • 1 year 8 months ago

A More Realistic Model of Isle Royale: Predator Prey Interactions

Andrew E Long
This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale.

We incorporate logistic growth into the moose dynamics, and we replace the death flow of the moose with a kill rate modeled from the kill rate data found on the Isle Royale website.

A decent match to the data is made with
Wolf Death Rate = 0.15
Wolf Birth Rate Factor = 0.0203
Moose Death Rate Factor = 1.08
Moose Birth Rate = 0.4
Carrying Capacity = 2000
Initial Moose: 563
Initial Wolves: 20

I used RK-4 with step-size 0.1, from 1959 for 60 years.

The moose birth flow is MBR*M*(1-M/K)
Moose death flow is MDRF*Sqrt(M*W)
Wolf birth flow is WBRF*Sqrt(M*W)
Wolf death flow is WDR*W

Environment Ecology Populations Midterm

  • 1 year 8 months ago

Isle Royale: Predator/Prey Model for Moose and Wolves, with Total Population

Andrew E Long
This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale. It was "cloned" from a model that InsightMaker provides to its users, at
https://insightmaker.com/insight/2068/Isle-Royale-Predator-Prey-Interactions
Thanks Scott Fortmann-Roe.

I've added in an adjustment to handle population.

I've created a Mathematica file that replicates the model, at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/Moose-n-Wolf-InsightMaker.nb

It allows one to experiment with adjusting the initial number of moose and wolves on the island.

I used steepest descent in Mathematica to optimize the parameters, with my objective data being the ratio of wolves to moose. You can try my (admittedly) kludgy code, at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/Moose-n-Wolf-InsightMaker-BestFit.nb

{WolfBirthRateFactorStart,
WolfDeathRateStart,
MooseBirthRateStart,
MooseDeathRateFactorStart,
moStart,
woStart} =
{0.000267409,
0.239821,
0.269755,
0.0113679,
591,
23.};

Environment Ecology Populations Math Modeling

  • 1 year 8 months ago

Sleek, non-dimensionalized Isle Royale: Predator/Prey Model for Moose and Wolves

Andrew E Long
This non-dimensionalized, sleekest most neatest model illustrates predator prey interactions using logistic growth for the moose population, for the wolf and moose populations on Isle Royale.

Thanks Scott Fortmann-Roe for the original model.

I've added in an adjustment to handle population sizes, by dividing by moose carrying capacity.

Time is scaled by the moose birth parameter:
tau=bm*t

There are therefore only three parameters left to account for any dynamics:

beta = bw/bm (relative wolf to moose births)
delta = dm/bm (relative death to birth ratio for moose)
gamma = dw/bm (wolf deaths to moose births)

The equations are thus

dM/dtau = M [ (1-M) - delta W ]
dW/dtau = W [beta M - gamma ]

There is a stable equilibrium pair of population values, relative to the carrying capacity:

M^* = gamma / beta
W^* = (1-gamma / beta) / delta

I have a sleek version with a logistical growth term for the moose, at

http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/Moose-n-Wolf-InsightMaker-sleek.nb

Environment Ecology Populations Math Modeling

  • 1 year 8 months ago

Pages