Sandusky Treatment Plant for Arsenic and Water purification

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.

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.

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.

Combining electromobility and renewable energies since 2014.

http://www.amsterdamvehicle2grid.nl/

This is a two-stock (ocean and atmosphere) climate model simulating the behavior of the earth climate from time zero. The initial conditions of the stocks are also set zero, so it demonstrates how long the earth takes to reach the temperature suitable for life.

Simple Model of the Food Chain

Simple model to illustrate Michaelis-Menten equation for nutrient uptake by phytoplankton.

The equation is:

P = Ppot S / (Ks + S)

Where:

P: Nutrient-limited production (e.g. d-1, or mg C m-2 d-1)
Ppot: Potential production (same units as P)
S: Nutrient concentation (e.g. umol N L-1)
Ks: Half saturation constant for nutrient (same units as S)

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

ENP 65000 assignment

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.

Here is my systems flow chart for a moose population!

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 adds diagenesis, using an extra state variable (phosphorus in the sediment) and incorporates desorption processes that release phosphorus trapped in the sediment back to the water column.

The temporal dynamics of the model simulate the typical development of pollution in time.

1. Low loading, low P concentration in lake
2. High loading, increasing P concentration in lake
3. Desorption rate is low, P in sediment increases
4. Measures implemented for source control, loading reduces
5. P in lake gradually decreases, but below a certain point, desorption increases, and lake P concentration does not improve
6. Recovery only occurs when the secondary load in the sediment is strongly reduced.

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.

A clone of the first model with the addition of a converter to describe the competition between rabbits for available vegetation based on the relationship between rabbit density and rabbit birth rate

The beginning of a systems dynamics model for teaching NRM 320.

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.

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.

The purpose of this deer management model is to explore the capacity of wildlife management actions to help us adapt to the effects of climate change.

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

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.