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. The equation for DeltaN is a version of  Nj+1 = Nj  + mu (1- Nj / Nmax ) Nj  the maximum population is set to be one million, and the growth rate constant mu
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.
The equation for DeltaN is a version of 
Nj+1 = Nj  + mu (1- Nj / Nmax ) Nj
the maximum population is set to be one million, and the growth rate constant mu = 3.
 
Nj: is the “number of items” in our current generation.

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.

This equation is a form of what is known as the logistic map or equation. It is a map because it "maps'' the population in one year into the population of the next year. It is "logistic'' in the military sense of supplying a population with its needs. It a nonlinear equation because it contains a term proportional to Nj^2 and not just Nj. The logistic map equation is also an example of discrete mathematics. It is discrete because the time variable j assumes just integer values, and consequently the variables Nj+1 and Nj do not change continuously into each other, as would a function N(t). In addition to the variables Nj and j, the equation also contains the two parameters mu, the growth rate, and Nmax, the maximum population. You can think of these as "constants'' whose values are determined from external sources and remain fixed as one year of items gets mapped into the next year. However, as part of viewing the computer as a laboratory in which to experiment, and as part of the scientific process, you should vary the parameters in order to explore how the model reacts to changes in them.
This model shows how a persistent pollutant such as mercury or DDT can be bioamplified along a trophic chain to levels that result in reduction of top predator populations.
This model shows how a persistent pollutant such as mercury or DDT can be bioamplified along a trophic chain to levels that result in reduction of top predator populations.
This stock and flow diagram provides a broad description of the key nutrient pathways (N and P) that exist in a dune-lake system subject to external loadings emanating from intensive agriculture.
This stock and flow diagram provides a broad description of the key nutrient pathways (N and P) that exist in a dune-lake system subject to external loadings emanating from intensive agriculture.
Simulation of MTBF with controls   F(t) = 1 - e ^ -λt   Where    • F(t) is the probability of failure    • λ is the failure rate in 1/time unit (1/h, for example)   • t is the observed service life (h, for example)  The inverse curve is the trust time On the right the increase in failures brings its
Simulation of MTBF with controls

F(t) = 1 - e ^ -λt 
Where  
• F(t) is the probability of failure  
• λ is the failure rate in 1/time unit (1/h, for example) 
• t is the observed service life (h, for example)

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
If we follow the slope from the leftmost start to where it begins to flatten out this can be considered the first period. The first period is characterized by a decreasing failure rate. It is what occurs during the “early life” of a population of units. The weaker units fail leaving a population that is more rigorous.

Useful Life
The next period is the flat bottom portion of the graph. It is called the “useful life” period. Failures occur more in a random sequence during this time. It is difficult to predict which failure mode will occur, but the rate of failures is predictable. Notice the constant slope.  

Wearout
The third period begins at the point where the slope begins to increase and extends to the rightmost end of the graph. This is what happens when units become old and begin to fail at an increasing rate. It is called the “wearout” period. 
 In 2012, the City of Vancouver created a sustainability strategy for staying on the leading edge of urban development called the “Greenest City: 2020 Action Plan (GCAP)” [ 1 — Open in Pop-up ]. In the report, the GCAP noted that its highest priority action was to encourage the use of electric vehic

In 2012, the City of Vancouver created a sustainability strategy for staying on the leading edge of urban development called the “Greenest City: 2020 Action Plan (GCAP)” [1Open in Pop-up]. In the report, the GCAP noted that its highest priority action was to encourage the use of electric vehicle transport in both public and private sectors. Since then, programs such as the Clean Energy Vehicle (CEV) program have been revamped to encourage consumers to choose the greener choice, often rewarding owners with up to $5000 in incentives for battery-powered vehicles and plug-in hybrids. However, the benefits of choosing electric cars are not all clear as several reports have found that hybrid electric vehicles (HEV), plug-in electric hybrid vehicles (PHEV), and battery electric cars (BEV) generate more carbon emissions during their production than current conventional vehicles [2]. I thought it would be interesting to study this sustainability issue through a systems model to determine how much impact it has on the environment compared to conventional vehicles. 

https://insightmaker.com/insight/159243/CO2-Emissions-by-Vehicle-Type-Gasoline-vs-Electric

Our model explores both carbon emissions of standard gasoline vehicles and electric vehicles from production to distribution in Canada specifically. Unfortunately, we were unable to find any statistics regarding the number of electric vehicles in production in Canada, so we have used the sales number as our production number estimate. For CO2 emission statistics, we made sure to carefully separate different types of electric vehicles as the production of the battery in battery electric vehicles have significantly more carbon emissions during production.

As expected, the carbon emissions from electric vehicles are much lower than those of gasoline vehicles after taking into account the lifecycle emissions from an average lifespan of 8 years on the road (which is the standard warranty length offered from most car companies). Some interesting things to note are that with our current rise in electric vehicle adoption, electric vehicles will dominate the roads in about 100 years. This transformation may be further accelerated by the large-scale initiatives offered by governmental organizations and increased awareness for sustainable practices. Furthermore, it was very surprising to find that electric vehicle carbon emissions will exceed that of gasoline vehicles after nearly 1000 years, but after further analysis, this makes sense as by then electric vehicles will greatly outnumber gasoline vehicles. This means that electric vehicles are not only the greener choice -- electric vehicles are by far the greenest choice as it will take nearly a thousand years before its emissions will be equal to that of its gasoline counterpart. In fact, it may even take longer than 1000 years for electric vehicles to emit more carbon emissions than gasoline vehicles if we continue looking for more sustainable methods for producing electricity and proactively choose renewable energy over fossil fuels.

Sources:

[1] https://vancouver.ca/files/cov/Greenest-city-action-plan.pdfOpen in Pop-up

[2] http://www.ccsenet.org/journal/index.php/jsd/article/view/64183

Statistics for number of gasoline and electric vehicle sales:

Gasoline Vehicles: https://www150.statcan.gc.ca/t1/tbl1/en/tv.action?pid=2010000201

Electric Vehicles: https://www.fleetcarma.com/electric-vehicle-sales-canada-2017/
 Economic growth cannot go on forever, although politicians and most economist
seem to think so. The activity involved in economic growth necessarily  generates entropy (disorder and environmental degradation). Entorpy in turn generates powerful negative feedback loops which will, as
a response from

Economic growth cannot go on forever, although politicians and most economist seem to think so. The activity involved in economic growth necessarily  generates entropy (disorder and environmental degradation). Entorpy in turn generates powerful negative feedback loops which will, as a response from nature, ensure that economic activity will eventually grind to a complete halt.  In these circumstances organised society cannot persist and will collapse. The negative feedback loops shown in this graph have already started to operate. The longer economic growth continues unabated, the more powerful these negative feedback loops will become. How long can economic growth continue before it is overwhelmed? It may not be very far in the future.

 STEM-SM combines a simple ecosystem model (modified version of VSEM; Hartig et al. 2019) with a soil moisture model (Guswa et al. (2002) leaky bucket model). Outputs from the soil moisture model influence ecosystem dynamics in three ways.   (1) The ratio of actual transpiration to maximum evapotran
STEM-SM combines a simple ecosystem model (modified version of VSEM; Hartig et al. 2019) with a soil moisture model (Guswa et al. (2002) leaky bucket model). Outputs from the soil moisture model influence ecosystem dynamics in three ways. 
(1) The ratio of actual transpiration to maximum evapotranspiration (T/ETmax) modifies gross primary productivity (GPP).
(2) Degree of saturation of the soil (Sd) modifies the rate of soil heterotrophic respiration.
(3) Water limitation of GPP (by T/ETmax) and of soil nutrient availability (approximated by Sd) combine with leaf area limitation (approximated by fraction of incident photosynthetically-active radiation that is absorbed) to modify the allocation of net primary productivity to aboveground and belowground parts of the vegetation.

Ecosystem dynamics in turn influence flows of water in to and out of the soil moisture stock. The size of the aboveground biomass stock determines fractional vegetation cover, which modifies interception, soil evaporation and transpiration by plants.

References:
Guswa, A.J., Celia, M.A., Rodriguez-Iturbe, I. (2002) Models of soil moisture dynamics in ecohydrology: a comparative study. Water Resources Research 38, 5-1 - 5-15.

Hartig, F., Minunno, F., and Paul, S. (2019). BayesianTools: General-Purpose MCMC and SMC Samplers and Tools for Bayesian Statistics. R package version 0.1.7. https://CRAN.R-project.org/package=BayesianTools

150 2 days ago
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 th
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).
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 adds diagenesis, using an extra state variable (ph
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 adds diagenesis, using an extra state variable (phosphorus in the sediment) and incorporates desorption processes that release phosphorus 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 P concentration in lake
2. High loading, increasing P concentration in lake
3. Desorption rate is low, P in sediment increases
4. Measures implemented for source control, loading reduces
5. P in lake gradually decreases, but below a certain point, desorption increases, and lake P concentration does not improve
6. Recovery only occurs when the secondary load in the sediment is strongly reduced.
  My model is on global population and its impact on the availability of natural resources. The stocks in my system include food availability, soil resources and water resource availability. One question I believe my model can address is, what are the connections between food availability, soil reso

My model is on global population and its impact on the availability of natural resources. The stocks in my system include food availability, soil resources and water resource availability. One question I believe my model can address is, what are the connections between food availability, soil resources and water resource availability; or in other words, are these stocks influenced equally by variables?  I hope to show a direct correlation between all three of these stocks. Food availability as stated in The Impact of Population Growth on Food Supplies and the Environment stated that, “The continued production of an adequate food supply is directly dependent on ample fertile land, fresh water, energy, plus the maintenance of biodiversity.” As population continues to grow so will the inputs to natural resources including water, fertilizer, and the need to have more available land.  What is more astonishing is that if these natural resources are never completely tapped dry, on a per capita perspective availability these resources will decline on astronomical levels since it has to be split amongst people (Pimental et al, 1996).



The flows in my system include food production,drought, water pollution, and greenhouse gases. I picked drought as a flow since it directly impacts the level of water available. Take for instance in California, the five year drought has caused scarcity and triggered state-wide executive orders to conserve water (California Department of Water Resources, 2017). Drought and water pollution can be affected by the number of people living in a country, which is why I picked these elements as flows. Furthermore food production, water pollution and greenhouse gases have strong influences on the availability of natural resources.


I picked mortality rates, birth rates, water scarcity, and industrial development as my variables. Since birth rates and mortality rates vary depending on the country I picked these as variables on my system since population growth is influenced by these variables.   Impact of Population Growth describes how the U. S. is already being affected by population growth, as stated here, “In populous industrial nations such as the United States, most economies of scale are already being exploited; we are on the diminishing returns part of most of the important curves.”


I have decided to change “developed countries” and undeveloped countries” as stocks to variables since these factors actually act more like variables. One question I hope to address with my model is how developed countries can  reduce their impact on resources? Furthermore, Population growth rate does depend on whether a country is developed versus undeveloped, so a country's level of economic development is more of a variable. I have decided to change food production from a stock to a flow, since it seems to be more of a flow that might affect the level of a stock of available food. I have also changed water scarcity from a stock to a variable because it actually affects the flow of water into an overall stock of fresh drinking water


A system diagram for the Mojave Desert for an assignment at OSU- RNG 341.
A system diagram for the Mojave Desert for an assignment at OSU- RNG 341.
all pictures sourced from google images
all pictures sourced from google images

This model implements a very simple proxy for vertical dispersion of heat in a lake based on the equation:  dT/dt = 1/A d(EA)/dz (dT/dz)  where: T: temperature (oC); t: time (days); z: depth (m); A: cross-sectional area (m2); E: vertical dispersion coefficient (m2 d-1)  If we consider that E is cons
This model implements a very simple proxy for vertical dispersion of heat in a lake based on the equation:

dT/dt = 1/A d(EA)/dz (dT/dz)

where: T: temperature (oC); t: time (days); z: depth (m); A: cross-sectional area (m2); E: vertical dispersion coefficient (m2 d-1)

If we consider that E is constant (it is in this model), then the equation becomes dT/dt = (EA/A)(d^2T/dz^2) = E(d^2T/dz^2), the classic diffusion equation

The model is simplified by exchanging temperature as a state variable, rather than executing  the full heat balance. This would require a computation of fluxes of atmospheric longwave and shortwave radiation, water longwave radiation, water conduction and convection, and water evaporation and condensation.

The vertical dispersion coefficients are adjusted artificially so that mixing increases at lower temperatures, thus quickly homogenizing the water column in colder months of the year.
Simple model to illustrate   algal  ,   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
Simple model to illustrate   algal  ,   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.

  Biogas, model  as well birefineray option to seperate c02 , chp from bogas model are proposed
This model illustrates predator prey interactions using real-life data of bison and wolf populations at Yellowstone National Park.
This model illustrates predator prey interactions using real-life data of bison and wolf populations at Yellowstone National Park.


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.        Ecocity model , Joanna       http://www.divaportal.se/sm
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.


Ecocity model , Joanna