This stock and flow diagram is a working draft of a conceptual model of a dune-lake system in the Northland region of New Zealand.

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 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 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 simulates the growth of carp in an aquaculture pond, both with respect to production and environmental effects.

Both the anabolism and fasting catabolism functions contain elements of allometry, through the m and n exponents that reduce the ration per unit body weight as the animal grows bigger.

The 'S' term provides a growth adjustment with respect to the number of fish, so implicitly adds competition (for food, oxygen, space, etc).

 Carp are mainly cultivated in Asia and Europe, and contribute to the world food supply.

Aquaculture currently produces sixty million tonnes of fish and shellfish every year. In 2011, aquaculture production overtook wild fisheries for human consumption.

This paradigm shift last occurred in the Neolithic period, ten thousand years ago, when agriculture displaced hunter-gatherers as a source of human food.

Aquaculture is here to stay, and wild fish capture (fishing) will never again exceed cultivation.

Recreational fishing will remain a human activity, just as hunting still is, after ten thousand years - but it won't be a major source of food from the seas.

The best way to preserve wild fish is not to fish them.

This model prototypes the working of an Smart Grid with Electric Vehicles

The objective is testing the theoretical advantages of batteries (also batteries in Electric Vehicles) in combination with renewable energies. The model considers two houses, that store energy both in Electric Vehicles (Vehicle to Grid), and in a communal battery.

Except when specified otherwise, the units of all variables are expressed in W/h.

Press "Story" in the lower bar for a guided tour over the model. Better seen at 50% zoom.

by Carlos Varela (cvarela@gmx.at)

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

Collapse of the economy, not just recession, is now very likely. To give just one possible cause, in the U.S. the fracking industry is in deep trouble. It is not only that most fracking companies have never achieved a free cash flow (made a profit) since the fracking boom started in 2008, but that  an already very weak  and unprofitable oil industry cannot cope with extremely low oil prices. The result will be the imminent collapse of the industry. However, when the fracking industry collapses in the US, so will the American economy – and by extension, probably, the rest of the world economy. To grasp a second and far more serious threat it is vital to understand the phenomenon of ‘Global Dimming’. Industrial activity not only produces greenhouse gases, but emits also sulphur dioxide which converts to reflective sulphate aerosols in the atmosphere. Sulphate aerosols act like little mirrors that reflect sunlight back into space, cooling the atmosphere. But when economic activity stops, these aerosols (unlike carbon dioxide) drop out of the atmosphere, adding perhaps as much as 1° C to global average temperatures. This can happen in a very short period time, and when it does mankind will be bereft of any means to mitigate the furious onslaught of an out-of-control and merciless climate. The data and the unrelenting dynamic of the viral pandemic paint bleak picture.  As events unfold in the next few months,  we may discover that it is too late to act,  that our reign on this planet has, indeed,  come to an abrupt end?  

This model demonstrates the various activities occurring across South East Queensland can both negatively and positively impact our endangered Koala population.

By Austin Wallace, Michelle Tomaszewski, and Ravish Sharma

Z212 from System Zoo 1 p142-148

My take on simulating the peak oil production

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

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

Diagrams of theories of control of destiny at multiple scales as fundamental causes of social determinants of health from Whitehead 2016 article in Health and Place

A draft model of the techonomy

This model simulates the growth of carp in an aquaculture pond, both with respect to production and environmental effects.

Both the anabolism and fasting catabolism functions contain elements of allometry, through the m and n exponents that reduce the ration per unit body weight as the animal grows bigger.

The 'S' term provides a growth adjustment with respect to the number of fish, so implicitly adds competition (for food, oxygen, space, etc).

 Carp are mainly cultivated in Asia and Europe, and contribute to the world food supply.

Aquaculture currently produces sixty million tonnes of fish and shellfish every year. In 2011, aquaculture production overtook wild fisheries for human consumption.

This paradigm shift last occurred in the Neolithic period, ten thousand years ago, when agriculture displaced hunter-gatherers as a source of human food.

Aquaculture is here to stay, and wild fish capture (fishing) will never again exceed cultivation.

Recreational fishing will remain a human activity, just as hunting still is, after ten thousand years - but it won't be a major source of food from the seas.

The best way to preserve wild fish is not to fish them.

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.

From science Dec 2106 article showing a CLD diagram