The simulation integrates or sums (INTEG) the Nj population, with a change of Delta N in each generation, starting with an initial value of 5.Nj+1 = Nj  + mu (1- Nj / Nmax ) Nj

Delta Nj: is the “change in number of items” as we go from the present generation into the next generation. This is just the number of items born minus the number of items who have died.

mu: is the growth or birth rate parameter, similar to that in the exponential growth and decay model. However, as we extend our model it will no longer be the actual growth rate, but rather just a constant that tends to control the actual growth rate without being directly proportional to it.

F(Nj) = mu(1‐Nj/Nmax): is our model for the effective “growth rate”, a rate that decreases as the number of items approaches the maximum allowed by external factors such as food supply, disease or predation. (You can think of mu as the growth or birth rate in the absence of population pressure from other items.) We write this rate as F(Nj), which is a mathematical way of saying F is affected by the number of items, i.e., “F is a function of Nj”. It combines both growth and all the various environmental constraints on growth into a single function. This is a good approach to modeling; start with something that works (exponential growth) and then modify it incrementally, while still incorporating the working model.

Nj+1 = Nj + Delta Nj : This is a mathematical way to say, “The new number of items equals the old number of items plus the change in number of items”.

Nj/Nmax: is what fraction a population has reached of the maximum "carrying capacity" allowed by the external environment. We use this fraction to change the overall growth rate of the population. In the real world, as well as in our model, it is possible for a population to be greater than the maximum population (which is usually an average of many years), at least for a short period of time. This means that we can expect fluctuations in which Nj/Nmax is greater than 1.

Simple box-model of the global carbon cycle

A storytelling of the nitrogen cycle.

Simulation of MTBF with controls


The inverse curve is the trust time
On the right the increase in failures brings its inverse which is loss of trust and move into suspicion and lack of confidence.
This can be seen in strategic social applications with those who put economy before providing the priorities of the basic living infrastructures for all.

This applies to policies and strategic decisions as well as physical equipment.
A) Equipment wears out through friction and preventive maintenance can increase the useful lifetime, 
B) Policies/working practices/guidelines have to be updated to reflect changes in the external environment and eventually be replaced when for instance a population rises too large (constitutional changes are required to keep pace with evolution, e.g. the concepts of the ancient Greeks, 3000 years ago, who based their thoughts on a small population cannot be applied in 2013 except where populations can be contained into productive working communities with balanced profit and loss centers to ensure sustainability)

Early Life
Useful Life

Wearout

A draft model of the techonomy

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.

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.

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.

Interplay between wolves eating sheep and farmers killing wolves.

Forcings and feedbacks based on Tom Fiddaman, James Hansen and other feedback and cycle diagrams

A simulation illustrating simple predator prey dynamics. You have two populations.

OVERSHOOT GROWTH GOES INTO TURBULENT CHAOTIC DESTRUCTION

The existing global capitalistic growth paradigm is totally flawed

The chaotic turbulence is the result of the concept of infinite bigness this has been the destructive influence on all empires and now shown up by Feigenbaum numbers and Dunbar numbers for neural netwoirks

See Guy Lakeman Bubble Theory for more details on keeping systems within finite limited size working capacity containers (villages communities)

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.

This model implements the equations proposed by Ketchum in 1954. The rationale behind the concept is that only phytoplankton that grows above a certain rate will not be flushed out of an estuary.

For biological processes:

Pt  =  Po exp(kt)

Where Pt is the phytoplankton biomass at time t, Po is the initial biomass, and k is the growth rate.

For physical processes:

Pm  =  Po (1-r)^m

Where Pm is the phytoplankton biomass after m tidal cycles, and r is the exchange ratio (proportion of estuary water which does not return each tidal cycle).

By substitution, and replacing t by m in the first equation, we get:

Pm = Poexp(km).(1-r)^m

For phytoplankton to exist in an estuary, Pm = Po (at least), i.e. 1 / (1-r)^m = exp(km)
ln(1) - m.ln(1-r) = km
-m.ln(1-r) = km
k = -ln(1-r)

Ketchum (1954) Relation between circulation and planktonic populations in estuaries. Ecology 35: 191-200.

In 2005, Ferreira and co-workers showed that this balance has direct implications on biodiversity of estuarine phytoplankton, and discussed how this could be relevant for water management, in particular for the EU Water Framework Directive 60/2000/EC (Ecological Modelling, 187(4) 513-523).

THE 2017 MODEL (BY GUY LAKEMAN) EMPHASIZES THE PEAK IN POLLUTION BEING CREATED BY OVERPOPULATION WITH THE CARRYING CAPACITY OF ARABLE LAND NOW BEING 1.5 TIMES OVER A SUSTAINABLE FUTURE (PASSED IN 1990) AND NOW INCREASING IN LOSS OF HUMAN SUSTAINABILITY DUE TO SEA RISE AND EXTREME GLOBAL WATER RELOCATION IN WEATHER CHANGES IN FLOODS AND DROUGHTS AND EXTENDED TROPICAL AND HORSE LATTITUDE CYCLONE ACTIVITY AROUND HADLEY CELLS

THE MODEL IS ZONE SPECIFIC AS GLOBAL WEATHER IS NOT HOMOGENEOUS BUT A COLLECTION OF HEAT BUMBPS DEPENDENT ON POPULATION SIZE OF URBAN HEAT ISLANDS AND MASSED CONURBATIONS AND AGGLOMERATIONS 

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.

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

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

Experiment with adjusting the moose birth-rate to simulate Over-shoot followed by environmental recovery

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

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.

Combining electromobility and renewable energies since 2014.

http://www.amsterdamvehicle2grid.nl/

This model shows the cycling of Mercury within a coastal wetland system. This cycling shows Elemental Mercury, Hg 2+, and Methylmercury within the soil, water, and air, and also interaction with the plants in the system.

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

The converter in this file contains precipitation for Phoenix only.