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. 

THE BROKEN LINK BETWEEN SUPPLY AND DEMAND CREATES TURBULENT CHAOTIC DESTRUCTION

The existing global capitalistic growth paradigm is totally flawed

Growth in supply and productivity is a summation of variables as is demand ... when the link between them is broken by catastrophic failure in a component the creation of unpredictable chaotic turbulence puts the controls ito a situation that will never return the system to its initial conditions as it is STIC system (Lorenz)

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 working containers (villages communities)

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.

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

As initially proposed by Pr. William M White of Cornell University:

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 diagram provides a stylised description of important feedbacks within a shallow-lake system.

Two households with PV systems and Electric Vehicles, sharing a battery and connected to the grid. What are the advantages?

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.

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.

DRAFT conceptual model of climate change connections in Yamuna river project.

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.

The beginning of a systems dynamics model for teaching NRM 320.

サンプル

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

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

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

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.

Ce modèle est une simulation classique du cycle de productivité dans l'océan, incluant les effets de la thermocline pour désactiver l'advection d'éléments nutritifs dissous et de détritus à la couche superficielle.  

Eastern oyster growth model calibrated for Long Island Sound

This is a one box model for an idealized farm with one million oysters seeded (one hectare @ a stocking density of 100 oysters per square meter)

1. Run WinShell individual growth model for one year with Long Island Sound growth drivers;

2. Determine the scope for growth (in dry tissue weight per day) for oysters centered on the five weight classes)
 
3. Apply a classic population dynamics equation:

dn(s,t)/dt = -d[n(s,t)g(s,t)]/ds - u(s)n(s,t)

s: Weight (g)
t: Time
n: Number of individuals of weight s
g: Scope for growth (g day-1)
u: Mortality rate (day-1)

4. Set mortality at 30% per year, slider allows scenarios from 30% to 80% per year

5. Determine harvestable biomass, i.e. weight class 5, 40-50 g (roughly three inches length)

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.

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.