slight: SALMO-Light: an absolutely simplistic lake model

In simecolModels: Model collection for the simecol package

Description

This model is an extremely cut-down version of the ecological lake model SALMO of Benndorf (1979). The small demo was inspired by the simplicity of the AQUAMOD-I model by Straskraba and Gnauck (1983) and is intended solely for teaching purposes.

Usage

slight()

Value

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

main	The differential equations for two phytoplankton groups, zooplankton and phosporus: phyto biovolume of phytoplankton (Chlorophyll, mu g/L), zoo biovolume of zooplankton (Carbon, mu g/L), phos phosphorus concentration (P, mu g/L),
parms	Vector with the named parameters of the model, see parms

,

times	Simulation time and integration interval.
init	Vector with start values for the state variables.
solver	Character string with the integration method.
parms	zin = 0, zooplankton-import (carbon, mu g/L) kxg = 0.25, half sat. ingest. zooplankton (SALMO, converted to carbon) gmax = 0.8, max. ingestion of zooplankton mean value between min (0.26) and max (1.3) (SALMO), temperature dependence ignored ae = 0.6, assimilation efficiency ep = 0.2, extinction of plankton free water (1/m) (lake specific) epx = 0.425, specific extinktion of phytoplankton carbon (m^2/mg) Lampert & Sommer, converted from chlorophyll to carbon (* 500/20) ki = 30, half saturation const light (J/(cm^2*d)) (SALMO, rounded) kp = 1.7, half saturation constant for phosphorus (mg/m^3) (SALMO, Diatoms) kx = 0.125, strongly simplified from (SALMO/C) (value of kxmin) photxmax = 1.8, maximum photosynthesis rate (1/d) (SALMO) photxmin = 0.17, minimum photosynthesis rate (1/d) (SALMO) cp = 0.04, C:P ratio mg/micro g, (Redfield, SALMO Diatoms) toptx = 20, optimal temperature for photosynthesis (degrees C) (SALMO Diatoms) vs = 0.1, sedimentation velocity (m/d) (SALMO Diaomst) rxtopt = 0.06, resp rate 1/d at opt temp (SALMO, Diatoms) rxtmin = 0.02, T dependence of resp (1/degrees C) (SALMO, Diatoms) rxmf = 0.3, light dpendend part of phytoplankton respiration (SALMO) kmo = 0.0175, half saturation constant of zooplankton mortality (carbon, mg/L) (SALMO, converted to carbon) momin = 0.015, Zooplankton mortality near 0 degrees C (1/d) (SALMO) mot = 0.006, temperature dependence of zooplankton mortality (degrees C/d) (SALMO) rz = 0.2, respiration of zooplankton (20% of grazing), assumption dtb = 0.3720601, inverse of egg development time at 20 degrees C (d) (BOTTRELL) dtmin = 5, minimum egg development time (d) (SALMO) tref = 293.15, reference temperature 20 degrees C (K) taegg = 11500, Arrhenius temperature of egg development time (Kelvin temperature), approximated after the curve of Bottrell (BOTTRELL, KOOIJMAN) iz = TRUE, i-depth version (TRUE) or i-mean version xcomp = TRUE, phytoplankton self limiting (TRUE) or not

References

Benndorf, J. (1979). Kausalanalyse, theoretische Synthese und Simulation des Eutrophierungsprozesses in stehenden und gestauten Gewässern. Habilitationsschrift. Technische Universität Dresden, Fakultät Bau-, Wasser- und Forstwesen, 165 S.

Benndorf, J. and Recknagel, F. (1982). Problems of application of the ecological model SALMO to lakes and reservoirs having various trophic states. Ecol. Model. 17: 129-145.

Straskraba, M. and Gnauck, A. (1983). Modellierung limnischer Ökosysteme. Fischer, Jena.

See Also

sim, parms, init, times.

Examples

slight <- slight()
## Not run: 
plot(sim(slight))

## show input matrix
View(inputs(slight))

## reduce phosphorus
inputs(slight)[,"pin"] <- 10
plot(sim(slight))


## End(Not run)

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