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

 Model to represent oyster individual growth by simulating feeding and metabolism. Builds on the core model in three ways: (i) partitions metabolic costs into feeding and fasting catabolism; (ii) adds allometry to clearance rate; (iii) adds temperature dependence to clearance rate.

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 an accessible description of the key processes that influence the water quality within a lake.

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.

Westley, F. R., O. Tjornbo, L. Schultz, P. Olsson, C. Folke, B. Crona and Ö. Bodin. 2013. A theory of transformative agency in linked social-ecological systems. Ecology and Society 18(3): 27. link

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 May 2013, 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.

Primitives for Watershed modeling project. Click Clone Insight at the top right to make a copy that you can edit.

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.

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 time-variable solution to a step-function change in inflow concentration for an ideal, completely mixed lake.

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

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.

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.

For Sustainability & Eco Innovation class

This model illustrates the key processes that influence the water level within Lake Okeechobee.

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 describes the N isotope dynamics in a marine system where the shallowness of the chemocline prevents significant nitrification near the base of the photic zone.

For my school assignment, I created a net showing what the ideal energy system would be like.

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.

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.