This model implements the one-dimensional version of the advection-dispersion equation for an estuary.  It develops the basic 1D model by simulating the mass of salt rather than the concentration. This makes it straightforward to deal with multiple boxes of different volume.  The equation is:  dS/dt
This model implements the one-dimensional version of the advection-dispersion equation for an estuary.

It develops the basic 1D model by simulating the mass of salt rather than the concentration. This makes it straightforward to deal with multiple boxes of different volume.

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.
Simple model of the global economy, the global carbon cycle, and planetary energy balance.    The planetary energy balance model is a two-box model, with shallow and deep ocean heat reservoirs. The carbon cycle model is a 4-box model, with the atmosphere, shallow ocean, deep ocean, and terrestrial c
Simple model of the global economy, the global carbon cycle, and planetary energy balance.

The planetary energy balance model is a two-box model, with shallow and deep ocean heat reservoirs. The carbon cycle model is a 4-box model, with the atmosphere, shallow ocean, deep ocean, and terrestrial carbon. 

The economic model is based on the Kaya identity, which decomposes CO2 emissions into population, GDP/capita, energy intensity of GDP, and carbon intensity of energy. It allows for temperature-related climate damages to both GDP and the growth rate of GDP.

This model was originally created by Bob Kopp (Rutgers University) in support of the SESYNC Climate Learning Project.
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 System Dynamics Causal Loop Diagram developed herein shows an initial study o
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 System Dynamics Causal Loop Diagram developed herein shows an initial study of the dynamics among the variables that influences the accumulation of PM in the air, in particular the case of Temuco, in the South of Chile. In Temuco, 97% of the PM inventories comes from the combustion of low quality firewood, which in turns is being burned due to its low price and cultural habits/tradition.
 Overview:   This simulation will show the relationship between tree logging forestry and how this can affect mountain biking tourism in Derby Park Tasmania. The main goal of this simulation is to show these two industries can co-exist in the same environment, or increase in demand or production in
Overview: 
This simulation will show the relationship between tree logging forestry and how this can affect mountain biking tourism in Derby Park Tasmania. The main goal of this simulation is to show these two industries can co-exist in the same environment, or increase in demand or production in one sector will affect the result of another.  

Function of the model:
In comparison there are both pros and cons for both sectors working correspondently. Demand for derby park is caused by individual past experience when visiting the park or friends recommendation which increase in the number of demands. Increase in demands will increase in the number of visitors. When visitors visits the park they require make a purchase a bike and pay the park for using the park facilities. All this will adds up to bikers total spending when visiting Derby. When consumer spend it is booting the economy especially in the tourism sector. Similarly tree logging will also contribute financially towards the Tasmania economy. The regeneration stage is relatively low compare to the logging rate. The growth will not cover the loss which can cause some level of damage in the scenery of the park, affecting tourist to view when mountain biking. Visitors overall experience will have the impact towards the demand for mountain biking in derby park, if visitors experience is satisfied they will come back to visit again or visit with group of friends, even words of mouth recommendation will also increase the level of demand of visiting Derby. 

Some key insights base on the simulation:
Based on the simulation of the two models we can see there are some key changes.
Tree logging increase will cause the disturbance of the natural scenery, thus change the overall experience of the visitors, decrease in the level of demand. Tree logging will also have negative impact towards the overall tourist experience thus affect the park facility and track. The natural scenery and the overall experience can affect their experience and if they would continue to recommend this area to friends to increase the demand. 

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.
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.
This is an edited copy of the original simple building heat flow. I addressed the insulation aspect of this problem by reducing the conductivity values of both the walls and windows. To address a chiller system, I integrated mechanical cooling into the process. To address solar panels, I calculated
This is an edited copy of the original simple building heat flow. I addressed the insulation aspect of this problem by reducing the conductivity values of both the walls and windows. To address a chiller system, I integrated mechanical cooling into the process. To address solar panels, I calculated the effect they may have on electricity prices by adding Solar Panel variables to the equation.
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-se
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.
Interactions between wolves and livestock depending on abundance
Interactions between wolves and livestock depending on abundance
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-se
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.
 Harvested fishery with endogenous investment and ship deployment policy. Ch 9 p345-360 John Morecroft (2007) Strategic Modelling and Business Dynamics. See simpler models at  IM-2990  and  IM-2991

Harvested fishery with endogenous investment and ship deployment policy. Ch 9 p345-360 John Morecroft (2007) Strategic Modelling and Business Dynamics. See simpler models at IM-2990 and IM-2991

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 websi
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 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 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.
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 System Dynamics Causal Loop Diagram developed herein shows an initial study o
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 System Dynamics Causal Loop Diagram developed herein shows an initial study of the dynamics among the variables that influences the accumulation of PM in the air, in particular the case of Temuco, in the South of Chile. In Temuco, 97% of the PM inventories comes from the combustion of low quality firewood, which in turns is being burned due to its low price and cultural habits/tradition.
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 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 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.
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.
   THE 2018 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 R

THE 2018 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 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.
BY SETTING THE NEW ECOLOGICAL POLICIES TO 2015 WE CAN SEE THAT SOME POPULATIONS CAN BE SAVED BUT CITIES WILL SUFFER MOST. 
CURRENT MARKET SATURATION PLATEAU OF SOLID PRODUCTS AND BEHAVIORAL SINK FACTORS ARE ALSO ADDED

Use the sliders to experiment with the initial amount of non-renewable resources to see how these affect the simulation. Does increasing the amount of non-renewable resources (which could occur through the development of better exploration technologies) improve our future? Also, experiment with the start date of a low birth-rate, environmentally focused policy.

 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.

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.

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


 The purpose of this deer management model is to explore the capacity of wildlife management actions to help us adapt to the effects of climate change.

The purpose of this deer management model is to explore the capacity of wildlife management actions to help us adapt to the effects of climate change.

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 System Dynamics Causal Loop Diagram developed herein shows an initial study o
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 System Dynamics Causal Loop Diagram developed herein shows an initial study of the dynamics among the variables that influences the accumulation of PM in the air, in particular the case of Temuco, in the South of Chile. In Temuco, 97% of the PM inventories comes from the combustion of low quality firewood, which in turns is being burned due to its low price and cultural habits/tradition.