Lemna_Schmitt {cvasi} | R Documentation |
Lemna model (Schmitt et al. 2013)
Description
The model is a mechanistic combined toxicokinetic-toxicodynamic (TK/TD) and growth model for the aquatic macrophytes Lemna spp. The model simulates the development of Lemna biomass under laboratory and environmental conditions and was developed by Schmitt et al. (2013). Growth of the Lemna population is simulated on basis of photosynthesis and respiration rates which are functions of environmental conditions. The toxicodynamic sub-model describes the effects of growth-inhibiting substances by a respective reduction in the photosynthesis rate based on internal concentrations. This is the historical version of the Lemna model. For current uses, we recommend the Lemna (SETAC) model, which is a more recent version of the Schmitt model.
Usage
Lemna_Schmitt(param, init)
Lemna_SchmittThold(param, init)
Arguments
param |
optional named |
init |
optional named numeric |
Details
Constructors to ease creation of scenarios based on the Lemna model by
Schmitt et al. (2013).
A variant of this Lemna model, Lemna_SchmittThold()
, provides an additional
cumulative exposure threshold parameter. The Lemna biomass stops growing
if the integral of exposure over time exceeds the threshold. The integral
of exposure is internally accounted for by an additional state variable
AUC
(Area Under Curve).
Value
an S4 object of type LemnaSchmitt
Functions
-
Lemna_SchmittThold()
: model variant with cumulative exposure threshold
State variables
The following list describes the default names and standard units of the model's state variables:
BM, g_dw/m2, dry weight biomass per square meter
E, -, effect [0,1]
M_int, ug, internal toxicant mass
AUC, ug/L, cumulative exposure (only for
LemnaThreshold
model)
Biomass (BM) and internal toxicant mass (M_int) are initialized to zero by
default. See set_init()
on how to set the initial states.
Model parameters
The following model parameters are required:
Fate and biomass
-
k_phot_fix
, logical, TRUE then k_phot_max is not changed by environmental factors, else FALSE -
k_phot_max
, 1/d, maximum photosynthesis rate -
k_resp
, 1/d, respiration rate -
k_loss
, 1/d, rate of loss (e.g. flow rate) -
mass_per_frond
, g_dw/frond, dry weight per frond -
BMw2BMd
, g_fw/g_dw, Fresh weight/dry weight
-
Effect
-
Emax
, -, maximum effect [0,1] -
EC50
, ug/L, midpoint of effect curve -
b
, -, slope of effect curve
-
Toxicokinetics
-
P_up
, cm/d, Permeability for uptake -
AperBM
, cm2/g_dw, A_leaf / d_leaf = 1/d_leaf (for circular disc, d=0.05 cm) -
Kbm
, -, Biomass(fw) : water partition coefficient -
P_Temp
, logical, TRUE to enable temperature dependence of cuticle permeability, else FALSE -
MolWeight
, g/mol, Molmass of molecule (determines Q10_permeability)
-
Temperature dependence
-
Tmin
, deg C, minimum temperature for growth -
Tmax
, deg C, maximum temperature for growth -
Topt
, deg C, optimal temperature for growth -
t_ref
, deg C, reference temperature for respiration rate -
Q10
, -, temperature dependence factor for respiration rate
-
Light dependence
-
k_0
, 1/d, light dependence: intercept of linear part -
a_k
, (1/d)/(kJ/m2.d), light dependence: slope of linear part
-
Phosphorus dependence (Hill like dep.)
-
C_P
, mg/L, phosphorus concentration in water -
CP50
, mg/L, phosphorus conc. where growth rate is halfed -
a_p
, -, Hill coefficient -
KiP
, mg/L, p-inhibition constant for very high p-conc.
-
Nitrogen dependence (Hill like dep.)
-
C_N
, mg/L, nitrogen concentration in water -
CN50
, mg/L, n-conc. where growth rate is halfed -
a_N
, -, Hill coefficient -
KiN
, mg/L, n-inhibition constant for very high p-conc.
-
Density dependence
-
BM50
, g_dw/m2, cut off BM
-
The Lemna_SchmittThold
model requires the following additional parameter:
-
threshold
, ug/L, cumulative exposure threshold
Forcings
Besides exposure, the Lemna model requires two environmental properties as
time-series input: global radiation (rad
, kJ/m2.d) and temperature (temp
, deg C).
Forcings time-series are represented by data.frame
objects consisting of two
columns. The first for time and the second for the environmental factor in question.
Entries of the data.frame
need to be ordered chronologically. A time-series
can consist of only a single row; in this case it will represent constant
environmental conditions. See scenarios for more details.
Effects
Supported effect endpoints include BM (biomass) and r (average growth rate during simulation). The effect on biomass is calculated from the last state of a simulation. Be aware that endpoint r is incompatible with frond transfers.
Parameter boundaries
Default values for parameter boundaries are set for all parameters by expert judgement, for calibration purposes. Values can be access from the object, and defaults overwritten.
Simulation output
Simulation results will contain two additional columns besides state variables:
-
C_int
, ug/L, internal concentration of toxicant -
FrondNo
, -, number of fronds
It is possible to amend the output of simulate()
with additional model
quantities that are not state variables, for e.g. debugging purposes or to
analyze model behavior. To enable or disable additional outputs, use the
optional argument nout
of simulate()
, see examples below. nout=1
enables reporting of internal concentration (C_int), nout=14
enables all
additional outputs, and nout=0
will disable additional outputs.
The available output levels are as follows:
-
nout >= 1
:C_int
, internal concentration (ug/L) -
nout >= 2
:FrondNo
, number of fronds (-) -
nout >= 3
:C_int_u
, unbound internal concentration (ug/l) Growth and TK/TD
-
nout >= 4
:BM_fresh
, fresh weight biomass (g_fw/m2) -
nout >= 5
:k_photo_eff
, current photosynthesis rate (1/d) -
nout >= 6
:k_resp_eff
, current respiration rate (1/d) -
nout >= 7
:f_Eff
, toxic effect factor (-) -
nout >= 8
:P_up_eff
, current permeability for uptake (cm/d)
-
Environmental variables
-
nout >= 9
:actConc
, current toxicant concentration in surrounding medium (ug/L) -
nout >= 10
:actTemp
, current environmental temperature (deg C) -
nout >= 11
:actRad
, current environmental radiation (kJ/m2.d)
-
Derivatives
-
nout >= 12
:d BM/dt
, current change in state variable BM -
nout >= 13
:d E/dt
, current change in effect -
nout >= 14
:d M_int/dt
, current change in internal toxicant mass
-
Solver settings
The arguments to ODE solver deSolve::ode()
control how model equations
are numerically integrated. The settings influence stability of the numerical
integration scheme as well as numerical precision of model outputs. Generally, the
default settings as defined by deSolve are used, but all deSolve settings
can be modified in cvasi workflows by the user, if needed. Please refer
to e.g. simulate()
on how to pass arguments to deSolve in cvasi
workflows.
Some default settings of deSolve were adapted for this model by expert judgement to enable precise, but also computationally efficient, simulations for most model parameters. These settings can be modified by the user, if needed:
-
hmax = 0.1
Maximum step length in time suitable for most simulations.
Biomass transfer
Models supporting biomass transfer can be instructed to move a fixed amount of biomass to a new medium after a period of time. This feature replicates a procedure occurring in e.g. Lemna effect studies and may be necessary to recreate study results.
The biomass transfer feature assumes that always a fixed amount of
biomass is transferred. Transfers can occur at any fixed point in time or
in regular intervals. During a transfer, the biomass is reset to the
transferred amount and additional compartments can be scaled 1:1 accordingly,
to e.g. reflect the change in internal toxicant mass when biomass is modified.
Transfer settings can be modified using set_transfer()
.
If a transfer occurs, simulation results of that time point will report the model state
before the transfer. Be aware that if transfers are defined using the
interval
argument, the transfers will always occur relative to time point
zero (t = 0
). As an example, setting a regular transfer of seven days,
interval = 7
, will result at transfers occurring at time points which are
integer multiplicates of seven, such as t=0
, t=7
, t=14
and so forth.
The starting and end times of a scenario do not influece when a regular
transfer occurs, only if it occurs.
References
Schmitt W., Bruns E., Dollinger M., and Sowig P., 2013: Mechanistic TK/TD-model simulating the effect of growth inhibitors on Lemna populations. Ecol Model 255, pp. 1-10. doi:10.1016/j.ecolmodel.2013.01.017
See Also
Lemna-models, Macrophyte-models, Transferable, Scenarios
Other Lemna models:
Lemna-models
,
Lemna_SETAC()
Other macrophyte models:
Lemna_SETAC()
,
Macrophyte-models
,
Myrio()
,
Myrio_log()