The time-variable solution to a step-function change in inflow concentration for an ideal, completely mixed lake.

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.

​Summary of Hermans Scale dynamics of grassroots innovations through parallel pathways of  transformative change Ecological Economics 2016 article (paywalled) This is applied to health in a subsequent insight

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.

Examining the ecosystem of the sea turtle and how that influences its population as an endangered species.

Simple (Kind of) food web of the Cane Toad Species. Includes different levels of consumers including predators.

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.

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.

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 Tucson only. Tucson watersheds are Arroyo Chico, Canada Agua, and Lower Canada del Oro.

Simple model to illustrate Steele's equation for primary production of phytoplankton.

The equation is:

Ppot = Pmax I/Iopt exp(1-I/Iopt)

Where:

Ppot: Potential production (e.g. d-1, or mg C m-2 d-1)
Pmax: Maximum production (same units as Ppot)
I: Light energy at depth of interest (e.g. uE m-2 s-1)
Iopt: Light energy at which Pmax occurs (same units as I)

The model contains no state variables, just illustrates the rate of production, by making the value of I equal to the timestep (in days). Move the slider to the left for more pronounced photoinhibition, to the right for photosaturation.

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

Our final beer game!

Simple model to illustrate a simple simulation of the microalgae biomass production, focusing on the dependent variables such as light, nutrients and other factor that is running for a yearly period.

The biomass model uses an example, Phytoplankton growth based 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).

Once this is understood, it looks upon the viability of biogas production from the microalgae biomass.

This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale.

Experiment with adjusting the moose birth-rate to simulate Over-shoot followed by environmental recovery

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.

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

Eastern oyster growth model calibrated for Great Bay.

Developed and implemented by Joao G. Ferreira and Camille Saurel; growth data, driver data, and culture practice from Phil Trowbridge, Ray Grizzle, and Suzanne Bricker.

This is a one box model for an idealized farm with one million oysters seeded (one hectare @ a stocking density of 100 oysters per square meter)

1. Run WinShell individual growth model for one year with Great Bay growth drivers;

2. Determine the scope for growth (in dry tissue weight per day) for oysters centered on the five weight classes)
 
3. Apply a classic population dynamics equation:

dn(s,t)/dt = -d[n(s,t)g(s,t)]/ds - u(s)n(s,t)

s: Weight (g)
t: Time
n: Number of individuals of weight s
g: Scope for growth (g day-1)
u: Mortality rate (day-1)

4. Set mortality at 30% per year, slider allows scenarios from 30% to 80% per year

5. Determine harvestable biomass, i.e. weight class 5, 40-50 g (roughly three inches length)