Economic growth cannot go on forever, although politicians and most economist seem to think so. The activity involved in economic growth necessarily  generates entropy (disorder and environmental degradation). Entorpy in turn generates powerful negative feedback loops which will, as a response from nature, ensure that economic activity will eventually grind to a complete halt.  In these circumstances organised society cannot persist and will collapse. The negative feedback loops shown in this graph have already started to operate. The longer economic growth continues unabated, the more powerful these negative feedback loops will become. How long can economic growth continue before it is overwhelmed? It may not be very far in the future.

My AP Environmental Homework for the Cats Over Borneo Assignment

this is the Australian food web of the water buffalo

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

A storytelling of the nitrogen cycle.

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.

Clone of Pesticide Use in Central America for Lab work

Simple Model of the Food Chain

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 shows how a persistent pollutant such as mercury or DDT can be bioamplified along a trophic chain to levels that result in reduction of top predator populations.

Simple model to illustrate   algal  ,   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.

  Biogas, model  as well birefineray option to seperate c02 , chp from bogas model are proposed

My model is on global population and its impact on the availability of natural resources. The stocks in my system include food availability, soil resources and water resource availability. One question I believe my model can address is, what are the connections between food availability, soil resources and water resource availability; or in other words, are these stocks influenced equally by variables?  I hope to show a direct correlation between all three of these stocks. Food availability as stated in The Impact of Population Growth on Food Supplies and the Environment stated that, “The continued production of an adequate food supply is directly dependent on ample fertile land, fresh water, energy, plus the maintenance of biodiversity.” As population continues to grow so will the inputs to natural resources including water, fertilizer, and the need to have more available land.  What is more astonishing is that if these natural resources are never completely tapped dry, on a per capita perspective availability these resources will decline on astronomical levels since it has to be split amongst people (Pimental et al, 1996).

The flows in my system include food production,drought, water pollution, and greenhouse gases. I picked drought as a flow since it directly impacts the level of water available. Take for instance in California, the five year drought has caused scarcity and triggered state-wide executive orders to conserve water (California Department of Water Resources, 2017). Drought and water pollution can be affected by the number of people living in a country, which is why I picked these elements as flows. Furthermore food production, water pollution and greenhouse gases have strong influences on the availability of natural resources.

I picked mortality rates, birth rates, water scarcity, and industrial development as my variables. Since birth rates and mortality rates vary depending on the country I picked these as variables on my system since population growth is influenced by these variables.   Impact of Population Growth describes how the U. S. is already being affected by population growth, as stated here, “In populous industrial nations such as the United States, most economies of scale are already being exploited; we are on the diminishing returns part of most of the important curves.”

I have decided to change “developed countries” and undeveloped countries” as stocks to variables since these factors actually act more like variables. One question I hope to address with my model is how developed countries can  reduce their impact on resources? Furthermore, Population growth rate does depend on whether a country is developed versus undeveloped, so a country's level of economic development is more of a variable. I have decided to change food production from a stock to a flow, since it seems to be more of a flow that might affect the level of a stock of available food. I have also changed water scarcity from a stock to a variable because it actually affects the flow of water into an overall stock of fresh drinking water

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 version adds diagenesis, using an extra state variable (phosphorus in the sediment) and incorporates desorption processes that release phosphorus 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 P concentration in lake
2. High loading, increasing P concentration in lake
3. Desorption rate is low, P in sediment increases
4. Measures implemented for source control, loading reduces
5. P in lake gradually decreases, but below a certain point, desorption increases, and lake P concentration does not improve
6. Recovery only occurs when the secondary load in the sediment is strongly reduced.

Simple (Kind of) food web of the Cane Toad Species. Includes different levels of consumers including predators.

General global ocean box model, including isotopes.

A food web for Africa. :)

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 top 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 second 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 third slider allows for adjustment of different aquaculture densities.

Wild filter-feeding species are included in the model, using an identical clearance rate to the cultivated oysters. Wild species can be turned on or off in the model using the fourth slider.

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

A system diagram for the Mojave Desert for an assignment at OSU- RNG 341.

Harvested fishery with stepwise changes in fleet size. Ch 9 p337-339 John Morecroft (2007) Strategic Modelling and Business Dynamics

Clone of Pesticide Use in Central America for Lab work