This model describes the N isotope dynamics in a marine system where the shallowness of the chemocline prevents significant nitrification near the base of the photic zone.

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

This model implements the equations proposed by Ketchum in 1954. The rationale behind the concept is that only phytoplankton that grows above a certain rate will not be flushed out of an estuary.

For biological processes:

Pt  =  Po exp(kt)

Where Pt is the phytoplankton biomass at time t, Po is the initial biomass, and k is the growth rate.

For physical processes:

Pm  =  Po (1-r)^m

Where Pm is the phytoplankton biomass after m tidal cycles, and r is the exchange ratio (proportion of estuary water which does not return each tidal cycle).

By substitution, and replacing t by m in the first equation, we get:

Pm = Poexp(km).(1-r)^m

For phytoplankton to exist in an estuary, Pm = Po (at least), i.e. 1 / (1-r)^m = exp(km)
ln(1) - m.ln(1-r) = km
-m.ln(1-r) = km
k = -ln(1-r)

Ketchum (1954) Relation between circulation and planktonic populations in estuaries. Ecology 35: 191-200.

In 2005, Ferreira and co-workers showed that this balance has direct implications on biodiversity of estuarine phytoplankton, and discussed how this could be relevant for water management, in particular for the EU Water Framework Directive 60/2000/EC (Ecological Modelling, 187(4) 513-523).

HANDY Model of Societal Collapse from Ecological Economics Paper 

This diagram provides a stylised description of important feedbacks within a shallow-lake system.

Simple model to illustrate   algal  ,   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.

  Biogas, model  as well birefineray option to seperate c02 , chp from bogas model are proposed

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 system dynamics model of a predator-prey lifecycle relationship

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

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.

Z212 from System Zoo 1 p142-148

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 used in a world studies extended essay research. The research question is: In what ways would the water desalination method used in Singapore benefit   water-stressed and economically less developed countries, using Palestine as a case study.

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.

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 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 basic model of the short-term carbon cycle.

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.

Primitives for Watershed modeling project. Click Clone Insight at the top right to make a copy that you can edit.

The converter in this file contains precipitation for Phoenix only.

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

In the prevoious step we added the reflection of sun energy by Clouds 

Harvested fishery with endogenous investment and ship deployment policy. Ch 9 p345-360 John Morecroft (2007) Strategic Modelling and Business Dynamics. See simpler models at IM-2990 and IM-2991