Insight diagram
​Modelo retirado do link 
https://insightmaker.com/insight/71649/Fern-Population-Model
Clone of Modelo da populacao de samambaias
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 simulation shows how plant, deer and wolf populations impact each other in a deciduous forest ecosystem.
Plant, Deer and Wolf Population Dynamics
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
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
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 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

Overview

This model simulates logging and mountain biking competition in Derby, Tasmania. The Simulation is referenced to simulate Derby mountain biking with logging.

 

Model Work

The tourism industry is represented on the model's left side, and the logging industry is on the right side. Interactions between these two industries generate tax revenues. Logging and tourism have different growth rates regarding people working/consuming. The initial values of these two industries in the model are not fixed but increase yearly due to inflation or economic growth.

 

Detail Insights

From the perspective of tourism, as the number of tourists keeps growing, the number of people who choose to ride in Derby City also gradually increases. And the people who ride rate the ride. The negative feedback feeds back into the cycling population. Similarly, positive cycling reviews lead to more customer visits. And all the customers will create a revenue through tourism, and a certain proportion of the income will become tourism tax.

From a logging perspective, it is very similar to the tourism industry. As the number of people working in the industry is forecast to increase, the industry's overall size is predicted to grow. And as the industry's size continues to rise, the taxes on the logging industry will also continue to rise. Since logging is an industry, the tax contribution will be more significant than the tourism excise tax.

 

This model assumption is illustrated below:

1. The amount of tax reflects the level of industrial development.

2. The goal of reducing carbon emissions lets us always pay attention to the environmental damage caused by the logging industry.

3. The government's regulatory goal is to increase overall income while ensuring the environment.

4. Logging will lead to environmental damage, which will decrease the number of tourists.

 

This model is based on tourism tax revenue versus logging tax revenue. Tourism tax revenue is more incredible than logging tax revenue, indicating a better environment. As a result of government policy, the logging industry will be heavily developed in the short term. Growth in the logging industry will increase by 40%. A growth rate of 0.8 and 0.6 of the original is obtained when logging taxes are 2 and 4 times higher than tourism taxes.

 

Furthermore, tourism tax and logging tax also act on the positive rate, which is the probability that customers give a positive evaluation. The over-development of the logging industry will lead to the destruction of environmental resources and further affect the tourism industry. The logging tax will also affect the tourism Ride Rate, which is the probability that all tourism customers will choose Derby city.

 

This model more accurately reflects logging and tourism's natural growth and ties the two industries together environmentally. Two ways of development are evident in the two industries. Compared to tourism, logging shows an upward spiral influenced by government policies. Government attitudes also affect tourism revenue, but more by the logging industry. 

Simulating Derby Mountain Biking Versus Logging
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 model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale.

Experiment with adjusting the initial number of moose and wolves on the island.
Clone of Clone of Isle Royale: Predator Prey Interactions
Insight diagram
This is a simulation that represents the populations of lions in the world over the last 200 years.
Clone of Lion Population Over The Last 200 Years
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 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 illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale.

Experiment with adjusting the initial number of moose and wolves on the island.
Clone of Isle Royale: Predator Prey Interactions (Español)
Insight diagram
See basic Dystopia model first. Idea here is to add a second "clean" energy source, which is is not constrained by pollution. The expectation, of course, is that, as the pollution constraint kicks in, this can take over the drive from the demand that would otherwise remain unsatisfied...
Dystopia-2: simple energy system "clean transition" model
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
This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale.

Experiment with adjusting the initial number of moose and wolves on the island.
Isle Royale Revised: Predator Prey Interactions
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 is a simulation that represents the populations of lions in the world over the last 200 years.
Clone of Lion Population Over The Last 200 Years
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 initial number of moose and wolves on the island.
Clone of (3) Copy of "Isle Royale: Predator Prey Interactions"
Insight diagram
Overview
This model which simulates the competition of Logging with Mountain Tourism in Derby, Tasmania.  This main reason of this simulation is to find if logging will affect the mountain tourism and by any chance they can co-exist.

How the model works.
Both Timber harvesting and mountain tourism can bring the economic contribution to Tasmania. In the Logging industry, it helps increase the need of employment and at the same time logging generate the profit through selling those timbers. In the Mountain Tourism industry, it can get the revenue through couple of ways which include accommodation (approximately 3 days find in paper), Restaurant and parking fee. However, the low growth rate of the trees is not keeping up with the rate of logging, if the trees getting less in Derby mountain, it will affect the sights and the riding experience for tourists, which will affect the satisfaction and expectation as it depends on the sights and experience. The satisfaction and expectation will influence the number of visitors, if they satisfied, they can come again or tell others about the great experience, if not, more and more people will not come again.

Interesting insights
It seems like logging has no significant negative effect to the mountain tourism, compare the forestry income with the tourism income, tourism income gradually higher than the forestry income at last, which means tourism is in a very important position, as long as the visitors are stable, tourism industry can provide greater economic contribution, stakeholders and governments can find the balance by maintain the status or better slightly reduce logging in order to make them co-exist.
Simulation of Derby Mountain biking versus logging
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)