Simple mass balance model for lakes, based on the Vollenweider equation:  dMw/dt = Min - sMw - Mout   The model was first used in the 1960s to determine the phosphorus concentration in lakes and reservoirs, for eutrophication assessment.     This version uses nitrogen and adds phytoplankton growth b
Simple mass balance model for lakes, based on the Vollenweider equation:

dMw/dt = Min - sMw - Mout

The model was first used in the 1960s to determine the phosphorus concentration in lakes and reservoirs, for eutrophication assessment.

This version uses nitrogen and adds phytoplankton growth based on a Michaelis-Menten equation
The time-variable solution to a step-function change in inflow concentration for an ideal, completely mixed lake.
The time-variable solution to a step-function change in inflow concentration for an ideal, completely mixed lake.
UDB 100 Assignment 1D   This simulation shows the drastic decline in the SEQ koala population due to various issues including human urban development, environmental issues, disease, habitat loss and domestic animals.
UDB 100 Assignment 1D


This simulation shows the drastic decline in the SEQ koala population due to various issues including human urban development, environmental issues, disease, habitat loss and domestic animals.

In Chile,  60% of its population are exposed to levels of Particulate Matter (PM) above international standards . Air Pollution is causing  4,000 premature deaths per year , including health costs over US$8 billion.    The System Dynamics Causal Loop Diagram developed herein shows an initial study o
In Chile, 60% of its population are exposed to levels of Particulate Matter (PM) above international standards. Air Pollution is causing 4,000 premature deaths per year, including health costs over US$8 billion.

The System Dynamics Causal Loop Diagram developed herein shows an initial study of the dynamics among the variables that influences the accumulation of PM in the air, in particular the case of Temuco, in the South of Chile. In Temuco, 97% of the PM inventories comes from the combustion of low quality firewood, which in turns is being burned due to its low price and cultural habits/tradition.
This model is a classic simulation of the production cycle in the ocean, including the effects of the thermocline in switching off advection of dissolved nutrients and detritus to the surface layer.  It illustrates a number of interesting features including the coupling of three state variables in a
This model is a classic simulation of the production cycle in the ocean, including the effects of the thermocline in switching off advection of dissolved nutrients and detritus to the surface layer.

It illustrates a number of interesting features including the coupling of three state variables in a closed cycle, the use of time to control the duration of advection, and the modulus function for cycling annual temperature data over multiple years.

The model state variables are expressed in nitrogen units (mg N m-3), and the calibration is based on:

Baliño, B.M. 1996. Eutrophication of the North Sea, 1980-1990: An evaluation of anthropogenic nutrient inputs using a 2D phytoplankton production model. Dr. scient. thesis, University of Bergen.
 
Fransz, H.G. & Verhagen, J.H.G. 1985. Modelling Research on the Production Cycle of Phytoplankton in the Southern Bight of the Northn Sea in Relation to Riverborne Nutrient Loads. Netherlands Journal of Sea Research 19 (3/4): 241-250.

This model was first implemented in PowerSim some years ago by one of my M.Sc. students, who then went on to become a Buddhist monk. Although this is a very Zen model, as far as I'm aware, the two facts are unrelated.
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

 A simulation illustrating simple predator prey dynamics. You have two populations.

A simulation illustrating simple predator prey dynamics. You have two populations.

Simple model to illustrate oyster growth based on primary production of Phytoplankton as a state variable, forced by light and nutrients, running for a yearly period.  Phytoplankton growth based on on Steele's and Michaelis-Menten equations), where:   Primary Production=(([Pmax]*[I]/[Iopt]*exp(1-[I]
Simple model to illustrate oyster growth based on primary production of Phytoplankton as a state variable, forced by light and nutrients, running for a yearly period.

Phytoplankton growth based on on Steele's and Michaelis-Menten equations), where: 

Primary Production=(([Pmax]*[I]/[Iopt]*exp(1-[I]/[Iopt])*[S])/([Ks]+[S]))

Pmax: Maximum production (d-1)
I: Light energy at depth of interest (uE m-2 s-1)
Iopt: Light energy at which Pmax occurs (uE m-2 s-1)
S: Nutrient concentration (umol N L-1)
Ks: Half saturation constant for nutrient (umol N L-1).

Further developments:
- Nutrients as state variable in cycle with detritus from phytoplankton and oyster biomass.
- Light limited by the concentration of phytoplankton.
- Temperature effect on phytoplankton and Oyster growth.


 Interplay between wolves eating sheep and farmers killing wolves who kill deer that eat crops that feed sheep.

Interplay between wolves eating sheep and farmers killing wolves who kill deer that eat crops that feed sheep.

 Interplay between wolves eating sheep and farmers killing wolves.

Interplay between wolves eating sheep and farmers killing wolves.

 Description:    A simple model for breeding plants from generation to generation in 3 different locations, with one "yield" variable (e.g. height) and 4 combinations of plants from the parents. Simulation tracks the frequencies of each combination in each generation as well as the overall average h
Description:

A simple model for breeding plants from generation to generation in 3 different locations, with one "yield" variable (e.g. height) and 4 combinations of plants from the parents. Simulation tracks the frequencies of each combination in each generation as well as the overall average height by generation.

The slider will select from 1 of 5 presets that changes the characteristics of each location's plants.

The graph of A1A2 Proportion represents both A1A2 and A2A1 since they are interchangeable.

Simple (Kind of) food web of the Cane Toad Species. Includes different levels of consumers including predators.
Simple (Kind of) food web of the Cane Toad Species. Includes different levels of consumers including predators.
Simulation of MTBF with controls   F(t) = 1 - e ^ -λt   Where    • F(t) is the probability of failure    • λ is the failure rate in 1/time unit (1/h, for example)   • t is the observed service life (h, for example)  The inverse curve is the trust time On the right the increase in failures brings its
Simulation of MTBF with controls

F(t) = 1 - e ^ -λt 
Where  
• F(t) is the probability of failure  
• λ is the failure rate in 1/time unit (1/h, for example) 
• t is the observed service life (h, for example)

The inverse curve is the trust time
On the right the increase in failures brings its inverse which is loss of trust and move into suspicion and lack of confidence.
This can be seen in strategic social applications with those who put economy before providing the priorities of the basic living infrastructures for all.

This applies to policies and strategic decisions as well as physical equipment.
A) Equipment wears out through friction and preventive maintenance can increase the useful lifetime, 
B) Policies/working practices/guidelines have to be updated to reflect changes in the external environment and eventually be replaced when for instance a population rises too large (constitutional changes are required to keep pace with evolution, e.g. the concepts of the ancient Greeks, 3000 years ago, who based their thoughts on a small population cannot be applied in 2013 except where populations can be contained into productive working communities with balanced profit and loss centers to ensure sustainability)

Early Life
If we follow the slope from the leftmost start to where it begins to flatten out this can be considered the first period. The first period is characterized by a decreasing failure rate. It is what occurs during the “early life” of a population of units. The weaker units fail leaving a population that is more rigorous.

Useful Life
The next period is the flat bottom portion of the graph. It is called the “useful life” period. Failures occur more in a random sequence during this time. It is difficult to predict which failure mode will occur, but the rate of failures is predictable. Notice the constant slope.  

Wearout
The third period begins at the point where the slope begins to increase and extends to the rightmost end of the graph. This is what happens when units become old and begin to fail at an increasing rate. It is called the “wearout” period.