Very simple model demonstrating growth of phytoplankton using Steele's equation for potential production and Michaelis-Menten equation for nutrient limitation.

Both light and nutrients (e.g. nitrogen) are modelled as forcing functions, and the model is "over-calibrated" for stability.

The phytoplankton model approximately reproduces the spring-summer diatom bloom and the (smaller) late summer dinoflagellate bloom.
 
Oyster growth is modelled only as a throughput from algae. Further developments would include filtration as a function of oyster biomass, oyster mortality, and other adjustments.

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.

Model of how different features impact water supply and how water access disparity can influence conflict.

This model explains the primary production of phytoplankton, forced by light and nutrients over a year period.

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

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

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.

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

HANDY Model of Societal Collapse from Ecological Economics Paper 

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.

Integrated Climate Model

An attempt to model some of the features in the national water market in Australia.

Challenges in sustainability are multilevel.
This diagram attempts to summarize levels of self reinforcing destructive dynamics, authors that deal with them, and point of leverage for change.

The base of the crisis is a mechanistic rather than ecological worldview. This mechanistic worldview is based on outdated science that assumed the universe to be a large machine. In a machine there is an inside and an outside. The health of the inside is important for the machine, the outside not. In an ecological view everything is interconnected, there is no clear separation in the future of self and other. All parts influence the health of other parts. To retain health sensitivity and democracy are inherent. The sense of separation from other that keeps the mechanistic worldview dominant is duality. Being cut off from spiritual traditions due to a mechanistic view of science people need access to inter-spirituality to reconnect with the human traditions and tools around connectedness, inner discovery, and compassion. Many books on modern physics and biology deal with the system view implications. "The coming interspiritual age" deals with the need to connect spiritual traditions and science.

At the bottom for the dynamic is an individual a sense of disconnectedness leads to a dependency on spending and having rather than connecting. The connecting has become too painful and dealing with it unpopular in our culture. Joanna Macy deals with this in Active Hope. 

This affluenza and disconnection is worsened by a market that floods one with advertisements aimed at creating needs and a sense of dissatisfaction with that one has.

National economies are structured around maximising GDP which means maximising consumption and financial capital movement. This is at the cost of local economies. These same local economies are needed for balanced happiness as well as for sustainability.

Generally institutions focus on maximising consumption rather than sustaining life support systems. David Korten covers this well.

Power and wealth is confused in this worldview. In striving for wealth only power is striven for in the form of money and monopoly.

Those at the head of large banks and corporations tend to be there because they exemplify this approach. They have few scruples about enforcing this approach onto everyone through wars and disaster capitalism. Naomi Klein and David Estulin documented this.

Power has become so centralized that we need this understanding to be widespread and include many of those in power. Progress of all of these levels are needed to show them and all that another way is possible.

My AP Environmental Homework for the Cats Over Borneo Assignment

Simple population dynamics examples based on ​Lotka-Volterra equations.

This is step 7 based on model 5 in making a climate model based on our insights of how trees actively contribute to the cooling capcacity of the Earth.​

In this step we added the PeTa-effect and the Atmosferic Window. 

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.

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.

A draft model of the techonomy

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

Simple model to illustrate an annual cycle for phytoplankton biomass in temperate waters.
Potential primary production uses Steele's equation and a Michaelis-Menten (or Monod) function for nutrient limitation. Respiratory losses are only a function of biomass.