Challenges in sustainability are multilevel.
This diagram attempts to summarize levels of self reinforcing destructive dynamics, authors that deal with them, and point of leverage for change.

The base of the crisis is a mechanistic rather than ecological worldview. This mechanistic worldview is based on outdated science that assumed the universe to be a large machine. In a machine there is an inside and an outside. The health of the inside is important for the machine, the outside not. In an ecological view everything is interconnected, there is no clear separation in the future of self and other. All parts influence the health of other parts. To retain health sensitivity and democracy are inherent. The sense of separation from other that keeps the mechanistic worldview dominant is duality. Being cut off from spiritual traditions due to a mechanistic view of science people need access to inter-spirituality to reconnect with the human traditions and tools around connectedness, inner discovery, and compassion. Many books on modern physics and biology deal with the system view implications. "The coming interspiritual age" deals with the need to connect spiritual traditions and science.

At the bottom for the dynamic is an individual a sense of disconnectedness leads to a dependency on spending and having rather than connecting. The connecting has become too painful and dealing with it unpopular in our culture. Joanna Macy deals with this in Active Hope. 

This affluenza and disconnection is worsened by a market that floods one with advertisements aimed at creating needs and a sense of dissatisfaction with that one has.

National economies are structured around maximising GDP which means maximising consumption and financial capital movement. This is at the cost of local economies. These same local economies are needed for balanced happiness as well as for sustainability.

Generally institutions focus on maximising consumption rather than sustaining life support systems. David Korten covers this well.

Power and wealth is confused in this worldview. In striving for wealth only power is striven for in the form of money and monopoly.

Those at the head of large banks and corporations tend to be there because they exemplify this approach. They have few scruples about enforcing this approach onto everyone through wars and disaster capitalism. Naomi Klein and David Estulin documented this.

Power has become so centralized that we need this understanding to be widespread and include many of those in power. Progress of all of these levels are needed to show them and all that another way is possible.

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.

For Sustainability & Eco Innovation class

WIP Stock Flow representation of Panarchy Adaptive Cycles

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

Marine plastic is rapidly increasing due to increasing production and use of plastic in all economic activities, short use times and long life times of plastic, and large mismanagement of plastic waste. With this, the threat plastic poses to the marine biosphere is also increasing and will continue to increase over a long time into the future. Risk knowledge is limited and risk perception and awareness are not resulting in significant mitigation efforts. The case study will aim at modeling the use and life cycles of plastic and the transport paths that lead to plastic entering the ocean. The models will be used to simulate possible futures based on a scenario approach. The results of these efforts will be visualized with the goal to increase risk awareness.

This model describes the flow of energy from generation to consumption for neighborhoods in the metro Atlanta area. It also calculates the cost of energy production and the number of years it will take to recover that cost.

Interplay between wolves eating sheep and farmers killing wolves who kill deer that eat crops that feed sheep.

Model of how different features impact water supply and how water access disparity can influence conflict.

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.

This model explains the primary production of phytoplankton, forced by light and nutrients over a year period.

A clone of the first model with the addition of a converter to describe the competition between rabbits for available vegetation based on the relationship between rabbit density and rabbit birth rate

This stock and flow diagram is an updated working draft of a conceptual model of a dune-lake system in the Northland region of New Zealand.

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.

A simple model for breeding plants from generation to generation in 3 different locations, with one "yield" variable (e.g. height) and 4 combinations of plants from the parents. Simulation tracks the frequencies of each combination in each generation as well as the overall average height by generation.

Simulate an impact of an asteroid of any Diameter at any given Speed!

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

Here is the Covid 19 Statistics model based on the Philippines.

This is step 5 in making a climate model based on our insights of how trees actively contribute to the cooling capcacity of the Earth.​

In the prevoious step we added the reflection of sun energy by Clouds