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 diagram provides a stylised description of important feedbacks within a shallow-lake 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.

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.

Testing for MSc

Very simple model demonstrating growth of phytoplankton using Steele's equation for potential production and Michaelis-Menten equation for nutrient limitation.

Both light and nutrients (e.g. nitrogen) are modelled as forcing functions, and the model is "over-calibrated" for stability.

The phytoplankton model approximately reproduces the spring-summer diatom bloom and the (smaller) late summer dinoflagellate bloom.
 
Oyster growth is modelled only as a throughput from algae. Further developments would include filtration as a function of oyster biomass, oyster mortality, and other adjustments.

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

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 implements the one-dimensional version of the advection-dispersion equation for an estuary.

It develops the basic 1D model by simulating the mass of salt rather than the concentration. This makes it straightforward to deal with multiple boxes of different volume.

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 upper 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 lower slider allows you to simulate a variable river flow, and understand how dispersion compensates for changes in freshwater input.

This model implements a very simple shellfish carrying capacity simulation for tidal creeks with freshwater input.

Physics

The 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 upper slider to turn off dispersion (set to zero). If the variable being simulated is (a) salinity, you will see that if the tidal wave did not mix with the estuary water due to turbulence, the estuary would quickly become a freshwater system; (b) POM, then the ocean (which typically has less POM) will not contribute a flushing effect and the concentration of POM in the tidal creek or estuary will be higher.

The middle slider allows you to simulate a variable river flow, and understand how dispersion compensates for changes in freshwater input.

Biology

Two biological functions are implemented in CREEK, both extremely simplified.

1. Primary production - a constant primary production rate is considered in gC m-3 d-1

2. Oyster filtration - a constant clearance rate (CR) is considered in L ind- 1 h-1, scaled to a certain stocking density S (ind m-3)

Units are normalized, and food depletion is CR * S * POM, in g POM m-3 d-1

The lower slider allows for adjustment of different densities.

The model provides three outputs:
1. POM concentration in mg L-1
2. Equivalent in chlorophyll (ug L-1)
3. Total oyster biomass in kg for the whole system

The body of research and studies generated on the Fryingpan River between the 1940s and the present supports the development of a conceptual model of ecosystem responses to hydrological regime behavior and streamflow management activities. This conceptual model should encourage conversations about system behavior and collective understanding among stakeholders regarding connections between specific hydrological regime characteristics affected by management of Ruedi Reservoir and the ecological or biological variables important to local communities. For the sake of simplicity, the model includes mostly unidirectional relationships—feedback loops are exploded to reveal intermediate connections between variables. This approach increases the number of variables represented in the system, perhaps increasing its complexity at first glance. However, the primary benefit to the end user is that the model becomes more readable and explicit in its representation of system behavior. 

 

The conceptual model presented here likely differs by degrees from those held by the various investigators who considered Fryingpan River processes over the previous 80 years. However, it affectively aggregates the ideas main presented by each of those individuals. This model focuses on hydrological and biological variables and does not incorporate the entire diversity of human uses and needs for water from the Fryingpan River (e.g. hydropower production for the City of Aspen, revenue generated in the Town of Basalt by angling activities, etc.).  Rather it attempts to illustrate how the conditional state of important ecosystem characteristics might respond to reservoir management activities that impact typical spring flows, peak flow timing and magnitude, summer flows, fall flows, and winter flows. 

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)

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.

A system diagram for the Mojave Desert including example socio-economic factors for an assignment at OSU- RNG 341.

This diagram provides a stylised description of important feedbacks within a shallow-lake system.

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

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

We incorporate logistic growth into the moose dynamics, and we replace the death flow of the moose with a kill rate modeled from the kill rate data found on the Isle Royale website.

I start with these parameters:
Wolf Death Rate = 0.15
Wolf Birth Rate = 0.0187963
Moose Birth Rate = 0.4
Carrying Capacity = 2000
Initial Moose: 563
Initial Wolves: 20

I used RK-4 with step-size 0.1, from 1959 for 60 years.

The moose birth flow is logistic, MBR*M*(1-M/K)
Moose death flow is Kill Rate (in Moose/Year)
Wolf birth flow is WBR*Kill Rate (in Wolves/Year)
Wolf death flow is WDR*W

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.

World4 is a predictive model for world population. Population has grown hyper-exponentially in the last millenium, with the doubling time decreasing from 900 years  in 1000 CE to a minimum of ~35 years in 1963 CE. Technology is defined as that which decreases the death rate and/or increases the effective birth rate (i.e. by decreasing infant mortality). Technology grows exponentially, therefore population fits a hyper-exponential (exponent within an exponent). Models for the end of growth are explored using equations that express the ways humans are depleting Earth's biocapacity, the nature of resource depletion, and the relationship between natural resources and human carrying capacity. This simple model, containing just two closed systems, captures the subtle shifts in the population trajectory of the last 50 years. Specifically, hyperexponential growth has given way to subexponential growth. A peak is predicted for the time around 2028.  [Bystroff, C. (2021). Footprints to singularity: A global population model explains late 20th century slow-down and predicts peak within ten years. PloS one, 16(5), e0247214.]

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 clearance rate per gram is used if we multiply by the oyster biomass

​ This model is used in a world studies extended essay research. The research question is: In what ways would the water desalination method used in Singapore benefit   water-stressed and economically less developed countries, using Palestine as a case study.

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

We incorporate logistic growth into the moose dynamics, and we replace the death flow of the moose with a kill rate modeled from the kill rate data found on the Isle Royale website.

Thanks to Jacob Englert for the model if-then-else structure.

I start with these parameters:
Wolf Death Rate = 0.15
Wolf Birth Rate = 0.0187963
Moose Birth Rate = 0.4
Carrying Capacity = 2000
Initial Moose: 563
Initial Wolves: 20

I used RK-4 with step-size 0.1, from 1959 for 60 years.

The moose birth flow is logistic, MBR*M*(1-M/K)
Moose death flow is Kill Rate (in Moose/Year)
Wolf birth flow is WBR*Kill Rate (in Wolves/Year)
Wolf death flow is WDR*W

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)