Parameters
ModVege makes use of a sizable number of different variables. Some of them are used as inputs to the simulation, some tune model behaviour and others are the actual state variables of the simulated grassland system.
The following section attempt to provide an as comprehensive as possible list of all of these. For a thorough understanding, the reader is referred to the original publication by [@jouven2006ModelPredictingDynamics].
State Variables
These values evolve during a model run and are stored in the respective fields in the [ModvegeSite] instance that was run.
-
BM
Standing biomass in kg DM per ha. -
BMG
Standing green biomass (kg DM / ha). -
cBM
Cumulativeley grown biomass (kg DM / ha). -
dBM
Daily grown biomass (kg DM / ha). -
hvBM
Cumulative harvested biomass (kg DM / ha). -
OMD
Organic matter digestibility (kg / kg). -
OMDG
OMD of green matter (kg / kg). -
ST
Temperature sum in degree Celsius days. -
REP
Reproductive function. Gives the fraction of growth that is assigned to reproductive growth. The remainder goes into vegetative growth. Dimensionless. -
PGRO
Potential growth in kg DM / ha. -
GRO
Effective growth in kg DM / ha. -
LAI
Leaf area index, accounting for the proportion of light intercepted by the sward. Dimensionless. -
LAIGV
LAI of green vegetative biomass. Dimensionless. -
AET
Actual evapotranspiration in mm. -
WR
Water reserves in mm. -
ENV
Function representing environmental effects on growth. Acts as a multiplicative factor. Dimensionless. -
ENVfPAR
Part of ENV due to strength of incident radiation. Dimensionless. -
ENVfT
Part of ENV due to temperature. Dimensionless. -
ENVfW
Part of ENV due to water limitation. Dimensionless.
Initial conditions
AgeGV
Age of green vegetative matter in degree Celsius days.AgeGR
Age of green reproductive matter in degree Celsius days.AgeDV
Age of dead vegetative matter in degree Celsius days.AgeDR
Age of dead reproductive matter in degree Celsius days.BMGV
biomass of GV (kg DM per ha).BMGR
biomass of GR (kg DM per ha).BMDV
biomass of DV (kg DM per ha).BMDR
biomass of DR (kg DM per ha).BMDR
biomass of DR (kg DM per ha).SENG
senescence of GV (kg DM per ha).SENG
senescence of GR (kg DM per ha).ABSG
abscission of DV (kg DM per ha).ABSG
abscission of DR (kg DM per ha).ST
thermal time (degree days).cBM
cumulative total biomass (kg per ha).
Model parameters
LON
geographic longitude of site in degree.LAT
geographic latitude of site in degree.ELV
geographic elevation of site in m.a.s.l.WHC
water-holding capacity of site in mm.NI
site nutritional index (dimensionless).RUEmax
maximum radiuation use efficiency in g DM per MJ.w_FGA
relative weight of functional group A.w_FGB
relative weight of functional group B.w_FGC
relative weight of functional group C.w_FGD
relative weight of functional group D.sigmaGV
rate of GV respirative biomass loss (dimensionless).sigmaGR
rate of GR respirative biomass loss (dimensionless).T0
photosynthesis activation temperature (degree C).T1
photosynthesis plateau temperature (degree C).T2
photosynthesis max temperature (degree C).KGV
basic senescence rate GV (dimensionless).KGR
basic senescence rate GR (dimensionless).KlGV
basic abscission rate GV (dimensionless).KlGR
basic abscission rate GR (dimensionless).maxOMDDV
organic matter digestibility in gram per gram DV.minOMDDR
organic matter digestibility in gram per gram DR.CO2_growth_factor
strength of effect of CO2 concentration on growth. See parameter b in [fCO2_growth_mod()].crop_coefficient
multiplicative factor Kc by which reference evapotranspiration ET0 has to be multiplied to get the crop evapotranspiration ETc: ETc = Kc ET0senescence_cap
fraction cs of GRO to which SEN is limited: SENimax = cs GROi for i in GV, GR. Makes it less likely for grass population to die out. Can be set to large values in order to effectively disable senescence capping.stubble_height
float. Minimum height the grass can assume. The biomass will not fall below that height. This effectively presents a simple model of plant reserves.SGS_method
string. Choice of method to determine the start of the growing season. Can be either"MTD"
for the multicriterial thermal definition (see [start_of_growing_season_mtd()]) or"simple"
for a commonly used approach as described in [start_of_growing_season()]).
Functional group parameters
SLA
Specific Leaf Area in m2 per g.pcLAM
Percentage of laminae (number between 0 and 1).ST1
Temperature sum in degree Celsiues days after which the seasonality functionSEA
starts to decrease from its maximum plateau. See also [SEA()].ST2
Temperature sum in degree Celsiues days after which the seasonality functionSEA
has decreased back to its minimum value. See also [SEA()].maxSEA
Maximum value of the seasonality function [SEA()]minSEA
Minimum value of the seasonality function [SEA()]. Usually,minSEA = 1 - (maxSEA - 1)
.maxOMDGV
Maximum organic matter digestability for green vegetative matter in arbitrary units.minOMDGV
Minimum organic matter digestability for green vegetative matter in arbitrary units.maxOMDGR
Maximum organic matter digestability for green reproductive matter in arbitrary units.minOMDGR
Minimum organic matter digestability for green reproductive matter in arbitrary units.BDGV
Bulk density of green vegetative dry matter in g per m3.BDDV
Bulk density of dead vegetative dry matter in g per m3.BDGR
Bulk density of green reproductive dry matter in g per m3.BDDR
Bulk density of dead reproductive dry matter in g per m3.fg_parameter_names
Vector of strings of the variable names of all vegetation parameters governed by functional group composition.
Weather variables
DOY
day of year in given yearTa
average temperature of given day (Celsius).precip
precipitation in millimeter per day.PAR
photosynthetically active radiation in MJ/m2. Can be calculated from average sunlight irradianceSRad
in J/s/m2 as:PAR = SRad * 0.47 * 24 * 60 * 60 / 1e6
ET0
evapotranspiration in mm.