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 + pMs - Mout  The model was first used in the 1960s to determine the phosphorus concentration in lakes and reservoirs for eutrophication assessment.  This version considers mercury, and adds diagenesis, using
Simple mass balance model for lakes based on the Vollenweider equation:

dMw/dt = Min - sMw + pMs - Mout

The model was first used in the 1960s to determine the phosphorus concentration in lakes and reservoirs for eutrophication assessment.

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

 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

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. 

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

 WIP Cloned from WIP  Africa Just Transition insight  including (Fig 3.1 from Jorgen Randers  book  2052 a Global Forecast for the Next Forty Years) with Fadhel Kaboub MMT Perspective to continue top down work on my slides of clds and macroeconomics

WIP Cloned from WIP Africa Just Transition insight including (Fig 3.1 from Jorgen Randers book 2052 a Global Forecast for the Next Forty Years) with Fadhel Kaboub MMT Perspective to continue top down work on my slides of clds and macroeconomics

10 months ago
Modeling the change in concentration of O2 in a lake with a continuous loading of BOD, modeled as a CSTR
Modeling the change in concentration of O2 in a lake with a continuous loading of BOD, modeled as a CSTR
My AP Environmental Homework for the Cats Over Borneo Assignment
My AP Environmental Homework for the Cats Over Borneo Assignment
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
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.
Clone of Pesticide Use in Central America for Lab work        This model is an attempt to simulate what is commonly referred to as the “pesticide treadmill” in agriculture and how it played out in the cotton industry in Central America after the Second World War until around the 1990s.     The cotto
Clone of Pesticide Use in Central America for Lab work


This model is an attempt to simulate what is commonly referred to as the “pesticide treadmill” in agriculture and how it played out in the cotton industry in Central America after the Second World War until around the 1990s.

The cotton industry expanded dramatically in Central America after WW2, increasing from 20,000 hectares to 463,000 in the late 1970s. This expansion was accompanied by a huge increase in industrial pesticide application which would eventually become the downfall of the industry.

The primary pest for cotton production, bol weevil, became increasingly resistant to chemical pesticides as they were applied each year. The application of pesticides also caused new pests to appear, such as leafworms, cotton aphids and whitefly, which in turn further fuelled increased application of pesticides. 

The treadmill resulted in massive increases in pesticide applications: in the early years they were only applied a few times per season, but this application rose to up to 40 applications per season by the 1970s; accounting for over 50% of the costs of production in some regions. 

The skyrocketing costs associated with increasing pesticide use were one of the key factors that led to the dramatic decline of the cotton industry in Central America: decreasing from its peak in the 1970s to less than 100,000 hectares in the 1990s. “In its wake, economic ruin and environmental devastation were left” as once thriving towns became ghost towns, and once fertile soils were wasted, eroded and abandoned (Lappe, 1998). 

Sources: Douglas L. Murray (1994), Cultivating Crisis: The Human Cost of Pesticides in Latin America, pp35-41; Francis Moore Lappe et al (1998), World Hunger: 12 Myths, 2nd Edition, pp54-55.

This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale. It was "cloned" from a model that InsightMaker provides to its users, at  https://insightmaker.com/insight/2068/Isle-Royale-Predator-Prey-Interactions  Thanks Scott Fortmann-Roe.
This model illustrates predator prey interactions using real-life data of wolf and moose populations on the Isle Royale. It was "cloned" from a model that InsightMaker provides to its users, at
https://insightmaker.com/insight/2068/Isle-Royale-Predator-Prey-Interactions
Thanks Scott Fortmann-Roe.

I've created a Mathematica file that replicates the model, at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/Moose-n-Wolf-InsightMaker.nb

It allows one to experiment with adjusting the initial number of moose and wolves on the island.

I used steepest descent in Mathematica to optimize the parameters, with my objective data being the ratio of wolves to moose. You can try my (admittedly) kludgy code, at
http://www.nku.edu/~longa/classes/2018spring/mat375/mathematica/Moose-n-Wolf-InsightMaker-BestFit.nb

{WolfBirthRateFactorStart,
WolfDeathRateStart,
MooseBirthRateStart,
MooseDeathRateFactorStart,
moStart,
woStart} =
{0.000267409,
0.239821,
0.269755,
0.0113679,
591,
23.};

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 uses nitrogen and adds phytoplankton growth b
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 uses nitrogen and adds phytoplankton growth based on a Michaelis-Menten equation
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 model to illustrate an annual cycle for phytoplankton biomass in temperate waters. Potential primary production uses Steele's equation and a Michaelis-Menten (or Monod) function for nutrient limitation. Respiratory losses are only a function of biomass.
Simple model to illustrate an annual cycle for phytoplankton biomass in temperate waters.
Potential primary production uses Steele's equation and a Michaelis-Menten (or Monod) function for nutrient limitation. Respiratory losses are only a function of biomass.