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.

From Schluter et al 2017 article A framework for mapping and comparing behavioural theories in models of social-ecological systems COMSeS2017 video. See also Balke and Gilbert 2014 JASSS article How do agents make decisions? (recommended by Kurt Kreuger U of S)

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 is a causal diagram story I made as an introduction to a workshop on systems thinking.

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 stock and flow diagram is an updated working draft of a conceptual model of a dune-lake system in the Northland region of New Zealand.

THE 2018 MODEL (BY GUY LAKEMAN) EMPHASIZES THE PEAK IN POLLUTION BEING CREATED BY OVERPOPULATION.
WITH THE CARRYING CAPACITY OF ARABLE LAND NOW BEING 1.5 TIMES OVER A SUSTAINABLE FUTURE (PASSED IN 1990) AND NOW INCREASING IN LOSS OF HUMAN SUSTAINABILITY DUE TO SEA RISE AND EXTREME GLOBAL WATER RELOCATION IN WEATHER CHANGES IN FLOODS AND DROUGHTS AND EXTENDED TROPICAL AND HORSE LATTITUDE CYCLONE ACTIVITY AROUND HADLEY CELLS

The World3 model is a detailed simulation of human population growth from 1900 into the future. It includes many environmental and demographic factors.

THIS MODEL BY GUY LAKEMAN, FROM METRICS OBTAINED USING A MORE COMPREHENSIVE VENSIM SOFTWARE MODEL, SHOWS CURRENT CONDITIONS CREATED BY THE LATEST WEATHER EXTREMES AND LOSS OF ARABLE LAND BY THE  ALBEDO EFECT MELTING THE POLAR CAPS TOGETHER WITH NORTHERN JETSTREAM SHIFT NORTHWARDS, AND A NECESSITY TO ACT BEFORE THERE IS HUGE SUFFERING.

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

The Logistic Map is a polynomial mapping (equivalently, recurrence relation) of degree 2, often cited as an archetypal example of how complex, chaotic behaviour can arise from very simple non-linear dynamical equations. The map was popularized in a seminal 1976 paper by the biologist Robert May, in part as a discrete-time demographic model analogous to the logistic equation first created by Pierre François Verhulst. 

Mathematically, the logistic map is written

where:

 is a number between zero and one, and represents the ratio of existing population to the maximum possible population at year n, and hence x0 represents the initial ratio of population to max. population (at year 0)r is a positive number, and represents a combined rate for reproduction and starvation. To generate a bifurcation diagram, set 'r base' to 2 and 'r ramp' to 1
To demonstrate sensitivity to initial conditions, try two runs with 'r base' set to 3 and 'Initial X' of 0.5 and 0.501, then look at first ~20 time steps

A simple phosphorus model for a generic lake.

The following insight shows the level of crime in the town of Bourke in comparison to the levels of Police and Community Engagement

This model implements the one-dimensional version of the advection-dispersion equation for an estuary. The equation is:

dS/dt = (1/A)d(QS)/dx - (1/A)d(EA)/dx(dS/dx) (Eq. 1)

Where S: salinity (or any other constituent such as chlorophyll or dissolved oxygen), (e.g. kg m-3); t: time (s); A: cross-sectional area (m2); Q: river flow (m3 s-1); x: length of box (m); E: dispersion coefficient (m2 s-1).

For a given length delta x, Adx = V, the box volume. For a set value of Q, the equation becomes:

VdS/dt = QdS - (d(EA)/dx) dS (Eq. 2)

EA/x, i.e. (m2 X m2) / (m s) = E(b), the bulk dispersion coefficient, units in m3 s-1, i.e. a flow, equivalent to Q

At steady state, dS/dt = 0, therefore we can rewrite Eq. 2 for one estuarine box as:

Q(Sr-Se)=E(b)r,e(Sr-Se)-E(b)e,s(Se-Ss) (Eq. 3)

Where Sr: river salinity (=0), Se: mean estuary salinity; Ss: mean ocean salinity

E(b)r,e: dispersion coefficient between river and estuary, and E(b)e,s: dispersion coefficient between the estuary and ocean.

By definition the value of E(b)r,e is zero, otherwise we are not at the head (upstream limit of salt intrusion) of the estuary. Likewise Sr is zero, otherwise we're not in the river. Therefore:

QSe=E(b)e,s(Se-Ss) (Eq. 4)

At steady state

E(b)e,s = QSe/(Se-Ss) (Eq 5)

The longitudinal dispersion simulates the turbulent mixiing of water in the estuary during flood and ebb, which supplies salt water to the estuary on the flood tide, and make the sea a little more brackish on the ebb.

You can use the slider to turn off dispersion (set to zero), and see that if the tidal wave did not mix with the estuary water due to turbulence, the estuary would quickly become a freshwater system.

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.

Fooodwaste happens everywhere and in every part of the food cycle even if nobody wants it to happen. 

We created a local solution to reduce the waste. This solution is situated in Belgium (Kotrijk) where an exchange system (for services) already exists and it is called letsleie http://www.letsleie.be.  We did choose letstlei because their exchange system doesn’t work with money but with a fictive money system "vlasbloemen". In their system we want to integrate the exchange of food leftovers. After some years the system could become world wide. 

Our solution begins with an event in a neighbourhood or apartments. This event brings the neighbours together who don't know each anymore. It explains the existing system and the problems of the food waste. Every person had to take a leftover and chefs will create a delicious meal of it. The members will receive a food box who is biodegradable and contains a QR code that will simplify the food/ service exchange. 

 People will talk to each other after the event and more and more people will join without needing new publicity.

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.

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 implements the one-dimensional version of the advection-dispersion equation for an estuary. The equation is:

dS/dt = (1/A)d(QS)/dx - (1/A)d(EA)/dx(dS/dx) (Eq. 1)

Where S: salinity (or any other constituent such as chlorophyll or dissolved oxygen), (e.g. kg m-3); t: time (s); A: cross-sectional area (m2); Q: river flow (m3 s-1); x: length of box (m); E: dispersion coefficient (m2 s-1).

For a given length delta x, Adx = V, the box volume. For a set value of Q, the equation becomes:

VdS/dt = QdS - (d(EA)/dx) dS (Eq. 2)

EA/x, i.e. (m2 X m2) / (m s) = E(b), the bulk dispersion coefficient, units in m3 s-1, i.e. a flow, equivalent to Q

At steady state, dS/dt = 0, therefore we can rewrite Eq. 2 for one estuarine box as:

Q(Sr-Se)=E(b)r,e(Sr-Se)-E(b)e,s(Se-Ss) (Eq. 3)

Where Sr: river salinity (=0), Se: mean estuary salinity; Ss: mean ocean salinity

E(b)r,e: dispersion coefficient between river and estuary, and E(b)e,s: dispersion coefficient between the estuary and ocean.

By definition the value of E(b)r,e is zero, otherwise we are not at the head (upstream limit of salt intrusion) of the estuary. Likewise Sr is zero, otherwise we're not in the river. Therefore:

QSe=E(b)e,s(Se-Ss) (Eq. 4)

At steady state

E(b)e,s = QSe/(Se-Ss) (Eq 5)

The longitudinal dispersion simulates the turbulent mixiing of water in the estuary during flood and ebb, which supplies salt water to the estuary on the flood tide, and make the sea a little more brackish on the ebb.

You can use the slider to turn off dispersion (set to zero), and see that if the tidal wave did not mix with the estuary water due to turbulence, the estuary would quickly become a freshwater system.

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.