This story presents a conceptual model of nitrogen cycling in a dune-lake system in the Northland region of New Zealand. It is based on the concept of a stock and flow diagram. Each orange ellipse represents an input, while each blue box represents a stock. Each arrow represents a flow. A flow involves a loss from the stock at which it starts and an addition to the stock at which it ends.

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.

Dissolved oxygen mass balance in a tide pool, forced by tides and light.

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

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

This model illustrates predator prey interactions using real-life data of rabbit and fox populations in Chile
Experiment with adjusting the initial number of moose and wolves on the island.

This model provides a dynamic simulation of the Sverdrup (1953) paper on the vernal blooming of phytoplankton.

The model simulates the dynamics of the mixed layer over the year, and illustrates how it's depth variation leads to conditions that trigger the spring bloom. In order for the bloom to occur, production of algae in the water column must exceed respiration.

This can only occur if vertical mixing cannot transport algae into deeper, darker water, for long periods, where they are unable to grow.

Sverdrup, H.U., 1953. On conditions for the vernal blooming of phytoplankton. J. Cons. Perm. Int. Exp. Mer, 18: 287-295

THE BROKEN LINK BETWEEN SUPPLY AND DEMAND CREATES TURBULENT CHAOTIC DESTRUCTION

The existing global capitalistic growth paradigm is totally flawed

Growth in supply and productivity is a summation of variables as is demand ... when the link between them is broken by catastrophic failure in a component the creation of unpredictable chaotic turbulence puts the controls ito a situation that will never return the system to its initial conditions as it is STIC system (Lorenz)

The chaotic turbulence is the result of the concept of infinite bigness this has been the destructive influence on all empires and now shown up by Feigenbaum numbers and Dunbar numbers for neural netwoirks

See Guy Lakeman Bubble Theory for more details on keeping systems within finite working containers (villages communities)

March 12, 2019

A basic model of the short-term carbon cycle.

This model describes phosphorus cycling in a dune-lake system in the Northland region of New Zealand. It is based on stock and flow diagrams where each orange oval represents an input, while each blue box represents a stock. Each arrow represents a flow. Flows involve a loss from the stock at which they start and add to the stock at which they end.

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.

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

EIGENE MODIFIKATIONEN

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

 

Use the sliders to experiment with the initial amount of non-renewable resources to see how these affect the simulation. Does increasing the amount of non-renewable resources (which could occur through the development of better exploration technologies) improve our future? Also, experiment with the start date of a low birth-rate, environmentally focused policy.

This model describes the flow of energy from generation to consumption for neighborhoods in the metro Atlanta area. It also calculates the cost of energy production and the number of years it will take to recover that cost.

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.

Here is my systems flow chart for a moose population!

European Masters in System Dynamics 2016
New University of Lisbon, Portugal

Simple model to represent oyster individual growth by simulating feeding and metabolism.

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.

Modèle simple pour illustrer la croissance des huîtres sur la base de la production primaire de phytoplancton comme une variable d'état, forcé par la lumière et les éléments nutritifs, en cours d'exécution pour une période annuelle.

La croissance du phytoplancton sur la base de Steele et équations de Michaelis-Menten), où:

Production primaire = (([Pmax] * [I] / [Iopt] * exp (1 - [I] / [Iopt]) * [S]) / ([K] + [S]))

Pmax: production maximale (d-1)
I: L'énergie lumineuse en profondeur de l'intérêt (Ue m-2 s-1)
Iopt: L'énergie lumineuse à laquelle se produit Pmax (Ue m-2 s-1)
S: concentration des éléments nutritifs (N umol L-1)
KS: Demi constants de saturation en nutriments (N umol L-1).

D'autres développements:
- Les éléments nutritifs comme variable d'état dans le cycle de détritus de phytoplancton et d'huîtres de la biomasse.
- Lumière limitée par la concentration de phytoplancton.
- Effet de la température sur le phytoplancton et la croissance des huîtres.

This model shows how a persistent pollutant such as mercury or DDT can be bioamplified along a trophic chain to levels that result in reduction of top predator populations.

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