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.

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?  

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

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.

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. 

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.

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.

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.

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

The causal-effect of an illegal squatting in a mini agroecosystem in a private residence

From Schluter et al 2017 article A framework for mapping and comparing behavioural theories in models of social-ecological systems COMSeS2017 video. See also Balke and Gilbert 2014 JASSS article How do agents make decisions? (recommended by Kurt Kreuger U of S)

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

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

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

THIS MODEL BY GUY LAKEMAN, FROM METRICS OBTAINED USING A MORE COMPREHENSIVE VENSIM SOFTWARE MODEL, SHOWS CURRENT CONDITIONS CREATED BY THE LATEST WEATHER EXTREMES AND LOSS OF ARABLE LAND BY THE  ALBEDO EFECT MELTING THE POLAR CAPS TOGETHER WITH NORTHERN JETSTREAM SHIFT NORTHWARDS, AND A NECESSITY TO ACT BEFORE THERE IS HUGE SUFFERING.

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.

March 12, 2019

A model of an infectious disease and control

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