mcyst_chemostat: Microcystin Production Model (Chemostat Version)

In simecolModels: Model collection for the simecol package

Description

This is a chemostat version of the MCYST model with substrate (phosphorus) dependend growth.

Usage

mcyst_chemostat()

Details

To see all details, please have a look into the implementation and the original publications.

Value

S4 object according to the odeModel specification. The object contains the following slots:

main	The differential equations for cell abundance and Microcystin concentration: cells abundance of cells (1/L), phos phosphorus concentration (nmol/L). mcyst microcystin concentration (fg/L).
parms	Vector with the named parameters of the model: mu.m max. growth rate (1/d), km half sat. constant (mumol/L), Y Yield coefficient, p MCYST production coefficient (fg/cell), dM MCYST decay rate (1/d), S0 phosphorus in import (umol/L), D dilution rate of chemostat (1/d).
times	Simulation time and integration interval.
init	Vector with start values for cells, phos and mcyst.
solver	Character string with the integration method.

References

Jähnichen, S., Petzoldt, T. and Benndorf, J. (2001). Evidence for Control of Microcystin Dynamics in Bautzen Reservoir (Germany) by Cyanobacterial Population Growth Rates and Dissolved Inorganic Carbon. Archiv für Hydrobiologie, 2, 150, 177-196.

Jähnichen, S., Ihle, T. and Petzoldt, T. (2008). Variability of microcystin cell quota: A small model explains dynamics and equilibria. Limnologica - Ecology and Management of Inland Waters, 38, 339-349. URL http://dx.doi.org/10.1016/j.limno.2008.05.003

See Also

sim, parms, init, times.

Examples

## create model object
m <- mcyst_chemostat()
## inspect parameters, start values, time steps
parms(m)
init(m)
times(m)
## simulate the model
m <- sim(m)
o <- out(m)

par(mfrow=c(2,2))
plot(o$time, o$cells, col="green", type="l",
  xlab="time", ylab="algae", ylim=c(0, max(o$cells)))
plot(o$time, o$phos,  col="blue",  type="l",
  xlab="time", ylab="substrate")
plot(o$time, o$mcyst/1000, col="cyan",  type="l",
  xlab="time", ylab="mcyst mug/l")
plot(o$time, o$mcyst/o$cells, col="red",  type="l",
  xlab="time", ylab="spec. mcyst")

simecolModels documentation built on May 2, 2019, 4:59 p.m.