Insight diagram
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

Clone of Final Midterm Student version of A More Realistic Model of Isle Royale: Predator Prey Interactions
Insight diagram
This model is a modified version of the 'Very Simple Ecosystem Model' (VSEM; Hartig et al. 2019). Controls have been added to gross primary productivity (GPP) and heterotrophic respiration (Rhetero) based on evapotranspiration rates.

Reference:
Hartig, F., Minunno, F., and Paul, S. (2019). BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R package version 0.1.7. https://CRAN.R-project.org/package=BayesianTools
Clone of Very Simple Ecosystem Model with Evapotranspiration (VSEM-ET)
Insight diagram
This model is a modified version of the 'Very Simple Ecosystem Model' (VSEM; Hartig et al. 2019). Controls have been added to gross primary productivity (GPP) and heterotrophic respiration (Rhetero) based on evapotranspiration rates.

Reference:
Hartig, F., Minunno, F., and Paul, S. (2019). BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R package version 0.1.7. https://CRAN.R-project.org/package=BayesianTools
Clone of Very Simple Ecosystem Model with Evapotranspiration (VSEM-ET)
Insight diagram
This model is to be used with Mr. Roderick's AP biology activity on population growth. See steveroderick.net for a copy of the activity worksheet.

Use the sliders below to quickly change the initial values of components of the model.
Clone of wolf ~ boom and bust
Insight diagram

This model has two main components. First is modelling the change in population composition as non-First Nations immigration increases with the opening of new mines in the region. The second is modelling the increasing income disparity between First Nations and non-First Nations as mining jobs are disproportionately gained by non-First Nations workers.

Clone of Northern Ontario Demographic and Income Trend Model
Insight diagram

A simulation illustrating how the size of a fish population affects the number of fish caught and vice versa.

Fish Population and Fisheries
Insight diagram
Clone of wolf ~ delayed logistic growth
Insight diagram
This model is a modified version of the 'Very Simple Ecosystem Model' (VSEM; Hartig et al. 2019). Controls have been added to gross primary productivity (GPP) and heterotrophic respiration (Rhetero) based on evapotranspiration rates.

Reference:
Hartig, F., Minunno, F., and Paul, S. (2019). BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R package version 0.1.7. https://CRAN.R-project.org/package=BayesianTools
Clone of Very Simple Ecosystem Model with Evapotranspiration (VSEM-ET)
Insight diagram
STEM-SM combines a simple ecosystem model (modified version of VSEM; Hartig et al. 2019) with a soil moisture model (Guswa et al. (2002) leaky bucket model). Outputs from the soil moisture model influence ecosystem dynamics in three ways. 
(1) The ratio of actual transpiration to maximum evapotranspiration (T/ETmax) modifies gross primary productivity (GPP).
(2) Degree of saturation of the soil (Sd) modifies the rate of soil heterotrophic respiration.
(3) Water limitation of GPP (by T/ETmax) and of soil nutrient availability (approximated by Sd) combine with leaf area limitation (approximated by fraction of incident photosynthetically-active radiation that is absorbed) to modify the allocation of net primary productivity to aboveground and belowground parts of the vegetation.

Ecosystem dynamics in turn influence flows of water in to and out of the soil moisture stock. The size of the aboveground biomass stock determines fractional vegetation cover, which modifies interception, soil evaporation and transpiration by plants.

References:
Guswa, A.J., Celia, M.A., Rodriguez-Iturbe, I. (2002) Models of soil moisture dynamics in ecohydrology: a comparative study. Water Resources Research 38, 5-1 - 5-15.

Hartig, F., Minunno, F., and Paul, S. (2019). BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R package version 0.1.7. https://CRAN.R-project.org/package=BayesianTools

Clone of Simple Terrestrial Ecosystem Model - Soil Moisture (STEM-SM)
Insight diagram
This is a model which explains the difference between Mountain bikes riding compared to logging in the Tasmanian forests.
Simulation of Derby Mountain bikes riding versus logging
Insight diagram
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

Clone of Final Midterm Student version of A More Realistic Model of Isle Royale: Predator Prey Interactions
Insight diagram
​Physical meaning of the equations
The Lotka–Volterra model makes a number of assumptions about the environment and evolution of the predator and prey populations:

1. The prey population finds ample food at all times.
2. The food supply of the predator population depends entirely on the size of the prey population.
3. The rate of change of population is proportional to its size.
4. During the process, the environment does not change in favour of one species and genetic adaptation is inconsequential.
5. Predators have limitless appetite.
As differential equations are used, the solution is deterministic and continuous. This, in turn, implies that the generations of both the predator and prey are continually overlapping.[23]

Prey
When multiplied out, the prey equation becomes
dx/dtαx - βxy
 The prey are assumed to have an unlimited food supply, and to reproduce exponentially unless subject to predation; this exponential growth is represented in the equation above by the term αx. The rate of predation upon the prey is assumed to be proportional to the rate at which the predators and the prey meet; this is represented above by βxy. If either x or y is zero then there can be no predation.

With these two terms the equation above can be interpreted as: the change in the prey's numbers is given by its own growth minus the rate at which it is preyed upon.

Predators

The predator equation becomes

dy/dt =  - 

In this equation, {\displaystyle \displaystyle \delta xy} represents the growth of the predator population. (Note the similarity to the predation rate; however, a different constant is used as the rate at which the predator population grows is not necessarily equal to the rate at which it consumes the prey). {\displaystyle \displaystyle \gamma y} represents the loss rate of the predators due to either natural death or emigration; it leads to an exponential decay in the absence of prey.

Hence the equation expresses the change in the predator population as growth fueled by the food supply, minus natural death.


Clone of Prey&Predator
Insight diagram
This model is a modified version of the 'Very Simple Ecosystem Model' (VSEM; Hartig et al. 2019). Controls have been added to gross primary productivity (GPP) and heterotrophic respiration (Rhetero) based on evapotranspiration rates.

Reference:
Hartig, F., Minunno, F., and Paul, S. (2019). BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R package version 0.1.7. https://CRAN.R-project.org/package=BayesianTools
Clone of Very Simple Ecosystem Model with Evapotranspiration (VSEM-ET)
Insight diagram
This model is a modified version of the 'Very Simple Ecosystem Model' (VSEM; Hartig et al. 2019). Controls have been added to gross primary productivity (GPP) and heterotrophic respiration (Rhetero) based on evapotranspiration rates.

Reference:
Hartig, F., Minunno, F., and Paul, S. (2019). BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R package version 0.1.7. https://CRAN.R-project.org/package=BayesianTools
Clone of Very Simple Ecosystem Model with Evapotranspiration (VSEM-ET)
Insight diagram
​Physical meaning of the equations
The Lotka–Volterra model makes a number of assumptions about the environment and evolution of the predator and prey populations:

1. The prey population finds ample food at all times.
2. The food supply of the predator population depends entirely on the size of the prey population.
3. The rate of change of population is proportional to its size.
4. During the process, the environment does not change in favour of one species and genetic adaptation is inconsequential.
5. Predators have limitless appetite.
As differential equations are used, the solution is deterministic and continuous. This, in turn, implies that the generations of both the predator and prey are continually overlapping.[23]

Prey
When multiplied out, the prey equation becomes
dx/dtαx - βxy
 The prey are assumed to have an unlimited food supply, and to reproduce exponentially unless subject to predation; this exponential growth is represented in the equation above by the term αx. The rate of predation upon the prey is assumed to be proportional to the rate at which the predators and the prey meet; this is represented above by βxy. If either x or y is zero then there can be no predation.

With these two terms the equation above can be interpreted as: the change in the prey's numbers is given by its own growth minus the rate at which it is preyed upon.

Predators

The predator equation becomes

dy/dt =  - 

In this equation, {\displaystyle \displaystyle \delta xy} represents the growth of the predator population. (Note the similarity to the predation rate; however, a different constant is used as the rate at which the predator population grows is not necessarily equal to the rate at which it consumes the prey). {\displaystyle \displaystyle \gamma y} represents the loss rate of the predators due to either natural death or emigration; it leads to an exponential decay in the absence of prey.

Hence the equation expresses the change in the predator population as growth fueled by the food supply, minus natural death.


Clone of Prey&Predator
Insight diagram
This model is a modified version of the 'Very Simple Ecosystem Model' (VSEM; Hartig et al. 2019). Controls have been added to gross primary productivity (GPP) and heterotrophic respiration (Rhetero) based on evapotranspiration rates.

Reference:
Hartig, F., Minunno, F., and Paul, S. (2019). BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R package version 0.1.7. https://CRAN.R-project.org/package=BayesianTools
Clone of Very Simple Ecosystem Model with Evapotranspiration (VSEM-ET)
Insight diagram
This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale.

Experiment with adjusting the moose birth-rate to simulate Over-shoot followed by environmental recovery
Clone of Royal Island- Resilience
Insight diagram
This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale.

Experiment with adjusting the moose birth-rate to simulate Over-shoot followed by environmental recovery
Clone of Clone of Royal Island- Resilience
Insight diagram
Overview
This model is a working simulation of the competition between the mountain biking tourism industry versus the forestry logging within Derby Tasmania.

How the model works
The left side of the model highlights the mountain bike flow beginning with demand for the forest that leads to increased visitors using the forest of mountain biking. Accompanying variables effect the tourism income that flows from use of the bike trails.
On the right side, the forest flow begins with tree growth then a demand for timber leading to the logging production. The sales from the logging then lead to the forestry income.
The model works by identifying how the different variables interact with both mountain biking and logging. As illustrated there are variables that have a shared effect such as scenery and adventure and entertainment.

Variables
The variables are essential in understanding what drives the flow within the model. For example mountain biking demand is dependent on positive word mouth which in turn is dependent on scenery. This is an important factor as logging has a negative impact on how the scenery changes as logging deteriorates the landscape and therefore effects positive word of mouth.
By establishing variables and their relationships with each other, the model highlights exactly how mountain biking and forestry logging effect each other and the income it supports.

Interesting Insights
The model suggests that though there is some impact from logging, tourism still prospers in spite of negative impacts to the scenery with tourism increasing substantially over forestry income. There is also a point at which the visitor population increases exponentially at which most other variables including adventure and entertainment also increase in result. The model suggests that it may be possible for logging and mountain biking to happen simultaneously without negatively impacting on the tourism income.
Clone of Simulation of Derby Mountain biking versus logging
Insight diagram
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.};

Clone of Isle Royale: Predator/Prey Model for Moose and Wolves
Insight diagram
​Physical meaning of the equations
The Lotka–Volterra model makes a number of assumptions about the environment and evolution of the predator and prey populations:

1. The prey population finds ample food at all times.
2. The food supply of the predator population depends entirely on the size of the prey population.
3. The rate of change of population is proportional to its size.
4. During the process, the environment does not change in favour of one species and genetic adaptation is inconsequential.
5. Predators have limitless appetite.
As differential equations are used, the solution is deterministic and continuous. This, in turn, implies that the generations of both the predator and prey are continually overlapping.[23]

Prey
When multiplied out, the prey equation becomes
dx/dtαx - βxy
 The prey are assumed to have an unlimited food supply, and to reproduce exponentially unless subject to predation; this exponential growth is represented in the equation above by the term αx. The rate of predation upon the prey is assumed to be proportional to the rate at which the predators and the prey meet; this is represented above by βxy. If either x or y is zero then there can be no predation.

With these two terms the equation above can be interpreted as: the change in the prey's numbers is given by its own growth minus the rate at which it is preyed upon.

Predators

The predator equation becomes

dy/dt =  - 

In this equation, {\displaystyle \displaystyle \delta xy} represents the growth of the predator population. (Note the similarity to the predation rate; however, a different constant is used as the rate at which the predator population grows is not necessarily equal to the rate at which it consumes the prey). {\displaystyle \displaystyle \gamma y} represents the loss rate of the predators due to either natural death or emigration; it leads to an exponential decay in the absence of prey.

Hence the equation expresses the change in the predator population as growth fueled by the food supply, minus natural death.


Clone of Prey&Predator
Insight diagram
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

Clone of Final Midterm Student version of A More Realistic Model of Isle Royale: Predator Prey Interactions
Insight diagram
​Physical meaning of the equations
The Lotka–Volterra model makes a number of assumptions about the environment and evolution of the predator and prey populations:

1. The prey population finds ample food at all times.
2. The food supply of the predator population depends entirely on the size of the prey population.
3. The rate of change of population is proportional to its size.
4. During the process, the environment does not change in favour of one species and genetic adaptation is inconsequential.
5. Predators have limitless appetite.
As differential equations are used, the solution is deterministic and continuous. This, in turn, implies that the generations of both the predator and prey are continually overlapping.[23]

Prey
When multiplied out, the prey equation becomes
dx/dtαx - βxy
 The prey are assumed to have an unlimited food supply, and to reproduce exponentially unless subject to predation; this exponential growth is represented in the equation above by the term αx. The rate of predation upon the prey is assumed to be proportional to the rate at which the predators and the prey meet; this is represented above by βxy. If either x or y is zero then there can be no predation.

With these two terms the equation above can be interpreted as: the change in the prey's numbers is given by its own growth minus the rate at which it is preyed upon.

Predators

The predator equation becomes

dy/dt =  - 

In this equation, {\displaystyle \displaystyle \delta xy} represents the growth of the predator population. (Note the similarity to the predation rate; however, a different constant is used as the rate at which the predator population grows is not necessarily equal to the rate at which it consumes the prey). {\displaystyle \displaystyle \gamma y} represents the loss rate of the predators due to either natural death or emigration; it leads to an exponential decay in the absence of prey.

Hence the equation expresses the change in the predator population as growth fueled by the food supply, minus natural death.


Clone of Prey&Predator
Insight diagram
​Physical meaning of the equations
The Lotka–Volterra model makes a number of assumptions about the environment and evolution of the predator and prey populations:

1. The prey population finds ample food at all times.
2. The food supply of the predator population depends entirely on the size of the prey population.
3. The rate of change of population is proportional to its size.
4. During the process, the environment does not change in favour of one species and genetic adaptation is inconsequential.
5. Predators have limitless appetite.
As differential equations are used, the solution is deterministic and continuous. This, in turn, implies that the generations of both the predator and prey are continually overlapping.[23]

Prey
When multiplied out, the prey equation becomes
dx/dtαx - βxy
 The prey are assumed to have an unlimited food supply, and to reproduce exponentially unless subject to predation; this exponential growth is represented in the equation above by the term αx. The rate of predation upon the prey is assumed to be proportional to the rate at which the predators and the prey meet; this is represented above by βxy. If either x or y is zero then there can be no predation.

With these two terms the equation above can be interpreted as: the change in the prey's numbers is given by its own growth minus the rate at which it is preyed upon.

Predators

The predator equation becomes

dy/dt =  - 

In this equation, {\displaystyle \displaystyle \delta xy} represents the growth of the predator population. (Note the similarity to the predation rate; however, a different constant is used as the rate at which the predator population grows is not necessarily equal to the rate at which it consumes the prey). {\displaystyle \displaystyle \gamma y} represents the loss rate of the predators due to either natural death or emigration; it leads to an exponential decay in the absence of prey.

Hence the equation expresses the change in the predator population as growth fueled by the food supply, minus natural death.


Clone of Prey&Predator