This non-dimensionalized, sleekest most neatest model illustrates predator prey interactions using logistic growth for the moose population, for the wolf and moose populations on Isle Royale.   Thanks Scott Fortmann-Roe for the original model.  I've added in an adjustment to handle population sizes,
This non-dimensionalized, sleekest most neatest model illustrates predator prey interactions using logistic growth for the moose population, for the wolf and moose populations on Isle Royale.

Thanks Scott Fortmann-Roe for the original model.

I've added in an adjustment to handle population sizes, by dividing by moose carrying capacity.

Time is scaled by the moose birth parameter:
tau=bm*t

There are therefore only three parameters left to account for any dynamics:

beta = bw/bm (relative wolf to moose births)
delta = dm/bm (relative death to birth ratio for moose)
gamma = dw/bm (wolf deaths to moose births)

The equations are thus

dM/dtau = M [ (1-M) - delta W ]
dW/dtau = W [beta M - gamma ]

There is a stable equilibrium pair of population values, relative to the carrying capacity:

M^* = gamma / beta
W^* = (1-gamma / beta) / delta

I have a sleek version with a logistical growth term for the moose, at

http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/Moose-n-Wolf-InsightMaker-sleek.nb
Simple mass balance model for lakes based on the Vollenweider equation:  dMw/dt = Min - sMw + pMs - Mout  The model was first used in the 1960s to determine the phosphorus concentration in lakes and reservoirs for eutrophication assessment.  This version considers mercury, and adds diagenesis, using
Simple mass balance model for lakes based on the Vollenweider equation:

dMw/dt = Min - sMw + pMs - Mout

The model was first used in the 1960s to determine the phosphorus concentration in lakes and reservoirs for eutrophication assessment.

This version considers mercury, and adds diagenesis, using an extra state variable (mercury in the sediment), and incorporates desorption processes that release mercury 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 Hg concentration in lake
2. High loading, increasing Hg concentration in lake
3. Desorption rate is low, Hg in sediment increases
4. Measures implemented for source control, loading reduces
5. Hg in lake gradually decreases, but below a certain point, desorption increases, and lake Hg concentration does not improve
6. Recovery only occurs when the secondary load in the sediment is strongly reduced.
 Harvested fishery with endogenous investment. Ch 9 p340-345 John Morecroft (2007) Strategic Modelling and Business Dynamics

Harvested fishery with endogenous investment. Ch 9 p340-345 John Morecroft (2007) Strategic Modelling and Business Dynamics

4 months ago
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
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 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
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.
M.Sc. in Environmental Engineering SIMA 2018 New University of Lisbon, Portugal   Model to represent oyster individual growth by simulating feeding and metabolism. Model (i) partitions metabolic costs into feeding and fasting catabolism; (ii) adds allometry to clearance rate; (iii) adds temperature
M.Sc. in Environmental Engineering SIMA 2018
New University of Lisbon, Portugal

 Model to represent oyster individual growth by simulating feeding and metabolism. Model (i) partitions metabolic costs into feeding and fasting catabolism; (ii) adds allometry to clearance rate; (iii) adds temperature dependence to clearance rate; (iv) illustrates how coupled model requires a substantial volume of water (a single oyster typically clears 20-30 m3 of water in one growth cycle)
Find the steady state completely mixed model with reaction decay and the three-compartment steady state model with reaction decay of a non-conservative tracer.
Find the steady state completely mixed model with reaction decay and the three-compartment steady state model with reaction decay of a non-conservative tracer.
Examining the population growth of a species and how reproductive technology influences its population as an endangered species.
Examining the population growth of a species and how reproductive technology influences its population as an endangered species.
Model prognozuje wpływ dystrybucji subwencji w ramach WPR na strukturę obszarową i intensywność produkcji gospodarstw rolnych w Polsce. Struktura obszarowa gospodarstw decyduje o stopniu mozaikowatości krajobrazu wiejskiego a intensywność produkcji przekłada się na bilans węgla i azotu użytków rolny
Model prognozuje wpływ dystrybucji subwencji w ramach WPR na strukturę obszarową i intensywność produkcji gospodarstw rolnych w Polsce. Struktura obszarowa gospodarstw decyduje o stopniu mozaikowatości krajobrazu wiejskiego a intensywność produkcji przekłada się na bilans węgla i azotu użytków rolnych. W oparciu o te dane szacowane są popularne wskaźniki jakości środowiska: zmiany zawartości materii organicznej w glebach, odpływ azotu do wód gruntowych i liczebność pospolitych ptaków krajobrazu wiejskiego.

Model opracowano w ramach Programu Wieloletniego IUNG-PIB w zadaniu 2.1 i jest obecnie rozwijany w zadaniu 1.9.
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 sat
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.
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.
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.
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
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 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
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 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
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.
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 System Dynamics Causal Loop Diagram developed herein shows an initial study o
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 System Dynamics Causal Loop Diagram developed herein shows an initial study of the dynamics among the variables that influences the accumulation of PM in the air, in particular the case of Temuco, in the South of Chile. In Temuco, 97% of the PM inventories comes from the combustion of low quality firewood, which in turns is being burned due to its low price and cultural habits/tradition.
This model uses simple functions (converters, cosine) to simulate the water balance inside a reservoir.
This model uses simple functions (converters, cosine) to simulate the water balance inside a reservoir.
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]
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 simulation illustrating how the size of a fish population affects the number of fish caught and vice versa.

A simulation illustrating how the size of a fish population affects the number of fish caught and vice versa.

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