This is a basic BIDE (birth, immigration, death, emigration) model.  Not all parts are implemented, however Birth and Death are.

This is a basic BIDE (birth, immigration, death, emigration) model.  Not all parts are implemented, however Birth and Death are.

    Clone of Bio103 Predator-Prey Model ("Lotka'Volterra")  Tags:  Education ,  Chaos ,  Ecology ,  Biology ,  Population   Thanks to Insight Author:  John Petersen       Edits by Andy Long     Everything that follows the dashes was created by John Petersen (or at least came from his Insight model).

Clone of Bio103 Predator-Prey Model ("Lotka'Volterra")
Thanks to Insight Author: John Petersen

Edits by Andy Long

Everything that follows the dashes was created by John Petersen (or at least came from his Insight model). I just wanted to make a few comments.

We are looking at Hare and Lynx, of course. Clone this insight, and change the names.

Then read the text below, to get acquainted with one of the most important and well-known examples of a simple system of differential equations in all of mathematics.

http://www.nku.edu/~longa/classes/mat375/mathematica/Lotka-Volterra.nb
------------------------------------------------------------

Dynamic simulation modelers are particularly interested in understanding and being able to distinguish between the behavior of stocks and flows that result from internal interactions and those that result from external forces acting on a system. 

For some time modelers have been particularly interested in internal interactions that result in stable oscillations in the absence of any external forces acting on a system. 

The model in this last scenario was independently developed by Alfred Lotka (1924) and Vito Volterra (1926).  Lotka was interested in understanding internal dynamics that might explain oscillations in moth and butterfly populations and the parasitoids that attack them.  Volterra was interested in explaining an increase in coastal populations of predatory fish and a decrease in their prey that was observed during World War I when human fishing pressures on the predator species declined. 

Both discovered that a relatively simple model is capable of producing the cyclical behaviors they observed. 

Since that time, several researchers have been able to reproduce the modeling dynamics in simple experimental systems consisting of only predators and prey.  It is now generally recognized that the model world that Lotka and Volterra produced is too simple to explain the complexity of most predator-prey dynamics in nature.  And yet, the model significantly advanced our understanding of the critical role of feedback in predator-prey interactions and in feeding relationships that result in community dynamics.

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.


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 websi
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

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 websi
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

 This is a basic model for use with our lab section.  The full BIDE options.

This is a basic model for use with our lab section.  The full BIDE options.

Westley, F. R., O. Tjornbo, L. Schultz, P. Olsson, C. Folke, B. Crona and Ö. Bodin. 2013. A theory of transformative agency in linked social-ecological systems.  Ecology and Society   18 (3): 27.  link
Westley, F. R., O. Tjornbo, L. Schultz, P. Olsson, C. Folke, B. Crona and Ö. Bodin. 2013. A theory of transformative agency in linked social-ecological systems. Ecology and Society 18(3): 27. link

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 websi
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

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 websi
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

 This is a basic BIDE (birth, immigration, death, emigration) model.  Not all parts are implemented, however Birth and Death are.

This is a basic BIDE (birth, immigration, death, emigration) model.  Not all parts are implemented, however Birth and Death are.

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 websi
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

This model shows the changing happened in forest industry and mountain tourism in Derby Tasmania. Logging will degrade mountain tourism while benefit the forestry industry. Simulation borrowed from the Easter Island simulation.    According to the analysis, logging does not reduce tourism income. Wi
This model shows the changing happened in forest industry and mountain tourism in Derby Tasmania. Logging will degrade mountain tourism while benefit the forestry industry. Simulation borrowed from the Easter Island simulation.

According to the analysis, logging does not reduce tourism income. With the increase of number of bike guide, tourism income will increase as well. Also, in forest industry, timber income is higher than the harvest spending which means the industry always gain profits from logging. Therefore, the main concern is that the logging should be balanced between the Mountain Tourism and the forest industry.
 Addition of Allee effect to Logistic Growth Insight 1540. From the Appendix of Marten Scheffer's 2009 Book Critical Transitions in Nature and Society p332 http://bit.ly/yrd3GN. The Allee affect describes a threshold density (smaller than the carrying capacity K) below which populations go into free

Addition of Allee effect to Logistic Growth Insight 1540. From the Appendix of Marten Scheffer's 2009 Book Critical Transitions in Nature and Society p332 http://bit.ly/yrd3GN. The Allee affect describes a threshold density (smaller than the carrying capacity K) below which populations go into free fall extinction.

 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 evapotran
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

9 months ago
A model situmalte the relationship between moutain bikes and logging industry in Derby, Tasmania, It explains more when the number of visitors increases or decreses.    How the model works  The left side shows when the number of travellers increase, the income from travellers rental of bike and stay
A model situmalte the relationship between moutain bikes and logging industry in Derby, Tasmania, It explains more when the number of visitors increases or decreses. 

How the model works
The left side shows when the number of travellers increase, the income from travellers rental of bike and stay of hotel increase simultaneously. However, there is a capacity for both parking lots and hotel venues, which means that the top ability of hospitality of Derby. The right side shows the logging industry of Derby and income from logging. It has a impact on how travellers would value Derby moutain.

Insights
As the number of travellers increase, it increases the total income of Derby, and in return, the local government will re-revest in Derby Moutain and will also maintain the forrestry logging industry. 
This is my first Insight from the BIDE model of environmental science for BIO 190
This is my first Insight from the BIDE model of environmental science for BIO 190
  Физический смысл уравнений    Модель Лотки-Вольтерры делает ряд предположений об окружающей среде и эволюции популяций хищников и жертв:         1. Хищная популяция всегда находит достаточно пищи.  2. Продовольственная обеспеченность популяции хищника полностью зависит от размера популяции жертвы.
Физический смысл уравнений
Модель Лотки-Вольтерры делает ряд предположений об окружающей среде и эволюции популяций хищников и жертв:

1. Хищная популяция всегда находит достаточно пищи.
2. Продовольственная обеспеченность популяции хищника полностью зависит от размера популяции жертвы.
3. Скорость изменения численности населения пропорциональна его численности.
4. В ходе этого процесса окружающая среда не меняется в пользу одного вида, и генетическая адаптация не имеет существенного значения.
5. Хищники обладают безграничным аппетитом.
Поскольку используются дифференциальные уравнения, решение является детерминированным и непрерывным. Это, в свою очередь, означает, что поколения как хищника, так и жертвы постоянно пересекаются.

Добыча
Когда умножается, уравнение добычи становится
dx/dt = αx - βxy
  Предполагается, что добыча имеет неограниченный запас пищи и размножается экспоненциально, если только она не подвержена хищничеству; этот экспоненциальный рост представлен в приведенном выше уравнении термином  αx. Предполагается, что скорость хищничества на добыче пропорциональна скорости, с которой встречаются хищники и добыча; это представлено выше в виде βxy.Если либо x, либо y равно нулю, то хищничества быть не может.
С помощью этих двух терминов приведенное выше уравнение можно интерпретировать следующим образом: изменение численности добычи определяется ее собственным ростом минус скорость, с которой она охотится.
ХищникиУравнение хищника становится

dy/dt =  - 

В этом уравнении,  представляет рост популяции хищника. (Обратите внимание на сходство со скоростью хищничества; однако используется другая константа, поскольку скорость роста популяции хищника не обязательно равна скорости, с которой он потребляет добычу).  представляет собой уровень потерь хищников вследствие естественной смерти или эмиграции; это приводит к экспоненциальному распаду в отсутствие добычи.


Следовательно, уравнение выражает изменение популяции хищников как рост, подпитываемый запасом пищи, минус естественная смерть.


This simulation shows how plant, deer and wolf populations impact each other in a deciduous forest ecosystem.
This simulation shows how plant, deer and wolf populations impact each other in a deciduous forest ecosystem.
 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). Baye
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