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.
-
BMStanding biomass in kg DM per ha. -
BMGStanding green biomass (kg DM / ha). -
cBMCumulativeley grown biomass (kg DM / ha). -
dBMDaily grown biomass (kg DM / ha). -
hvBMCumulative harvested biomass (kg DM / ha). -
OMDOrganic matter digestibility (kg / kg). -
OMDGOMD of green matter (kg / kg). -
STTemperature sum in degree Celsius days. -
REPReproductive function. Gives the fraction of growth that is assigned to reproductive growth. The remainder goes into vegetative growth. Dimensionless. -
PGROPotential growth in kg DM / ha. -
GROEffective growth in kg DM / ha. -
LAILeaf area index, accounting for the proportion of light intercepted by the sward. Dimensionless. -
LAIGVLAI of green vegetative biomass. Dimensionless. -
AETActual evapotranspiration in mm. -
WRWater reserves in mm. -
ENVFunction representing environmental effects on growth. Acts as a multiplicative factor. Dimensionless. -
ENVfPARPart of ENV due to strength of incident radiation. Dimensionless. -
ENVfTPart of ENV due to temperature. Dimensionless. -
ENVfWPart of ENV due to water limitation. Dimensionless.
Initial conditions
AgeGVAge of green vegetative matter in degree Celsius days.AgeGRAge of green reproductive matter in degree Celsius days.AgeDVAge of dead vegetative matter in degree Celsius days.AgeDRAge of dead reproductive matter in degree Celsius days.BMGVbiomass of GV (kg DM per ha).BMGRbiomass of GR (kg DM per ha).BMDVbiomass of DV (kg DM per ha).BMDRbiomass of DR (kg DM per ha).BMDRbiomass of DR (kg DM per ha).SENGsenescence of GV (kg DM per ha).SENGsenescence of GR (kg DM per ha).ABSGabscission of DV (kg DM per ha).ABSGabscission of DR (kg DM per ha).STthermal time (degree days).cBMcumulative total biomass (kg per ha).
Model parameters
LONgeographic longitude of site in degree.LATgeographic latitude of site in degree.ELVgeographic elevation of site in m.a.s.l.WHCwater-holding capacity of site in mm.NIsite nutritional index (dimensionless).RUEmaxmaximum radiuation use efficiency in g DM per MJ.w_FGArelative weight of functional group A.w_FGBrelative weight of functional group B.w_FGCrelative weight of functional group C.w_FGDrelative weight of functional group D.sigmaGVrate of GV respirative biomass loss (dimensionless).sigmaGRrate of GR respirative biomass loss (dimensionless).T0photosynthesis activation temperature (degree C).T1photosynthesis plateau temperature (degree C).T2photosynthesis max temperature (degree C).KGVbasic senescence rate GV (dimensionless).KGRbasic senescence rate GR (dimensionless).KlGVbasic abscission rate GV (dimensionless).KlGRbasic abscission rate GR (dimensionless).maxOMDDVorganic matter digestibility in gram per gram DV.minOMDDRorganic matter digestibility in gram per gram DR.CO2_growth_factorstrength of effect of CO2 concentration on growth. See parameter b in [fCO2_growth_mod()].crop_coefficientmultiplicative factor Kc by which reference evapotranspiration ET0 has to be multiplied to get the crop evapotranspiration ETc: ETc = Kc ET0senescence_capfraction 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_heightfloat. Minimum height the grass can assume. The biomass will not fall below that height. This effectively presents a simple model of plant reserves.SGS_methodstring. 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
SLASpecific Leaf Area in m2 per g.pcLAMPercentage of laminae (number between 0 and 1).ST1Temperature sum in degree Celsiues days after which the seasonality functionSEAstarts to decrease from its maximum plateau. See also [SEA()].ST2Temperature sum in degree Celsiues days after which the seasonality functionSEAhas decreased back to its minimum value. See also [SEA()].maxSEAMaximum value of the seasonality function [SEA()]minSEAMinimum value of the seasonality function [SEA()]. Usually,minSEA = 1 - (maxSEA - 1).maxOMDGVMaximum organic matter digestability for green vegetative matter in arbitrary units.minOMDGVMinimum organic matter digestability for green vegetative matter in arbitrary units.maxOMDGRMaximum organic matter digestability for green reproductive matter in arbitrary units.minOMDGRMinimum organic matter digestability for green reproductive matter in arbitrary units.BDGVBulk density of green vegetative dry matter in g per m3.BDDVBulk density of dead vegetative dry matter in g per m3.BDGRBulk density of green reproductive dry matter in g per m3.BDDRBulk density of dead reproductive dry matter in g per m3.fg_parameter_namesVector of strings of the variable names of all vegetation parameters governed by functional group composition.
Weather variables
DOYday of year in given yearTaaverage temperature of given day (Celsius).precipprecipitation in millimeter per day.PARphotosynthetically active radiation in MJ/m2. Can be calculated from average sunlight irradianceSRadin J/s/m2 as:PAR = SRad * 0.47 * 24 * 60 * 60 / 1e6ET0evapotranspiration in mm.
