Skip to content

Parameter in the model

ParameterValueDescription
ϕ_TRSA0.0683575 m^2 g^-1Reference root surface area per total biomass, used in nutrient stress function and maintenance costs for roots function, set to mean of community: ϕ_TRSA=mean((1abp)rsa)
ϕ_TAMC0.108293Reference arbuscular mycorriza colonisation rate per total biomass, used in nutrient stress function and maintenance costs for mycorrhizae function, set to mean of community: ϕ_TAMC=mean((1abp)amc)
ϕ_sla0.008808 m^2 g^-1Reference specific leaf area, used in senescence function, set to mean of community: ϕ_sla=mean(sla)
γ_RUEmax0.003 kg MJ^-1Maximum radiation use efficiency
γ_RUE_k0.6Light extinction coefficient
α_RUE_cwmH0.998592Reduction factor of radiation use efficiency at a height of 0.2 m ∈ [0, 1]
α_WAT_rsa050.822698Water stress growth reduction factor for species with mean trait: TRSA=ϕ_TRSA, when the plant available water equals: W_p,t=0.5
β_WAT_rsa8.12609Slope of the logistic function that relates the plant available water to the water stress growth reduction factor
δ_WAT_rsa1.60378 g m^-2Controls how strongly species differ in their water stress growth reduction from the mean response
ω_NUT_totalN0.1 kg g^-1Controls the influence of the total soil nitrogen on the nutrient index
ω_NUT_fertilization0.001 ha kg^-1Controls the influence of the fertilization rate on the nutrient index
β_TS1.0Scales the similarity for nutrient uptake traits (root surface area and arbuscular and mycorriza colonisation rate), which is used for competition for nutrients (very dissimilar species compete less strongly for nutrients). If β_TS=0, all trait similiarities are set to 1 and the trait similarity has no effect on the competition for nutrients, with increasing β_TS the effect of the trait similarity on the nutrient competition becomes stronger.
α_NUT_TSB5001.26 kg ha^-1Reference value, if the sum of the product of trait similarity and biomass of all species equals: TSB<1, TSB=1, TSB>1 the nutrient adjustment factor NUT_adj,ts is higher than one, one and lower than one, respectively
α_NUT_maxadj10.0Maximum of the nutrient adjustment factor, fixed for calibration
α_NUT_amc050.700429Nutrient stress based on arbuscular mycorriza colonisation growth reduction factor for species with mean trait: TAMC=ϕ_TAMC, when the plant available nutrients equal: N_p,ts=0.5
α_NUT_rsa050.959764Nutrient stress based on root surface area growth reduction factor for species with mean trait: TRSA=ϕ_TRSA, when the plant available nutrients equal: N_p,ts=0.5
β_NUT_rsa8.13689Slope of the logistic function that relates the plant available nutrients to the nutrient stress growth reduction factor based on root surface area & calibrated
β_NUT_amc15.9892Slope of the logistic function that relates the plant available nutrients to the nutrient stress growth reduction factor based on arbuscular mycorriza colonisation
δ_NUT_rsa10.3304 g m^-2Controls how strongly species differ in their nutrient stress growth reduction based on root surface area from the mean response
δ_NUT_amc15.0Controls how strongly species differ in their nutrients stress growth reduction based on arbuscular mycorriza colonisation from the mean response & calibrated
κ_ROOT_amc0.0723991Maximum growth reduction due to maintenance costs for mycorrhizae based on arbuscular mycorriza colonisation rate
κ_ROOT_rsa0.00761945Maximum growth reduction due to maintenance costs for fine roots based on root surface area
γ_RAD14.45e-6 ha MJ^-1Controls the steepness of the linear decrease in radiation use efficiency for high PAR_t values
γ_RAD250000.0 MJ ha^-1Threshold value of PAR_t from which starts a linear decrease in radiation use efficiency
ω_TEMP_T14.0 °CMinimum temperature for growth
ω_TEMP_T210.0 °CLower limit of optimum temperature for growth
ω_TEMP_T320.0 °CUpper limit of optimum temperature for growth
ω_TEMP_T435.0 °CMaximum temperature for growth
ζ_SEA_ST1787.414 °CThreshold of the cumulative temperate since the beginning of the current year, the seasonality factor starts to decrease from ζ_SEAmax to ζ_SEAmin above ζ_SEA,ST_1100°C
ζ_SEA_ST21800.0 °CThreshold of the cumulative temperate since the beginning of the current year, above which the seasonality factor is set to ζ_SEAmin
ζ_SEAmin0.987514Minimum value of the seasonal growth effect
ζ_SEAmax2.44419Maximum value of the seasonal growth effect
α_SEN0.0452193Basic senescence rate
β_SEN_sla1.41955Controls the influence of the specific leaf area on the senescence rate
ψ_SEN_ST11765.52 °CThreshold of the cumulative temperate since the beginning of the current year above which the senescence begins to increase
ψ_SEN_ST23000.0 °CThreshold of the cumulative temperate since the beginning of the current year above which the senescence reaches the maximum senescence rate
ψ_SENmax1.55784Maximum senescence rate
β_GRZ_lnc0.578292Controls the influence of leaf nitrogen per leaf mass on grazer preference
β_GRZ_H0.0563242Controls the influence of height on grazer preference
η_GRZ2.0Scaling factor that controls at which biomass density additional feed is supplied by farmers, fixed for calibration
κ_GRZ22.0 kgConsumption of dry biomass per livestock and day
ϵ_GRZ_minH0.05 mMinimum height that is reachable by grazers
β_SND_WHC0.5678Slope parameter relating the sand content to the soil water content at the water holding capacity
β_SLT_WHC0.9228Slope parameter relating the silt content to the soil water content at the water holding capacity
β_CLY_WHC0.9135Slope parameter relating the clay content to the soil water content at the water holding capacity
β_OM_WHC0.6103Slope parameter relating the organic matter content to the soil water content at the water holding capacity
β_BLK_WHC-0.2696 cm^3 g^-1Slope parameter relating the bulk density to the soil water content at the water holding capacity
β_SND_PWP-0.0059Slope parameter relating the sand content to the soil water content at the permanent wilting point
β_SLT_PWP0.1142Slope parameter relating the silt content to the soil water content at the permanent wilting point
β_CLY_PWP0.5766Slope parameter relating the clay content to the soil water content at the permanent wilting point
β_OM_PWP0.2228Slope parameter relating the organic matter content to the soil water content at the permanent wilting point
β_BLK_PWP0.02671 cm^3 g^-1Slope parameter relating the bulk density to the soil water content at the permanent wilting point

Which method uses a parameter?

ParameterUsed in...
ϕ_TRSAroot_investment!; nutrient_reduction!; water_reduction!
ϕ_TAMCroot_investment!; nutrient_reduction!
ϕ_slainitialize_senescence_rate!
γ_RUEmaxpotential_growth!
γ_RUE_klight_competition_height_layer!; potential_growth!
α_RUE_cwmHpotential_growth!
α_WAT_rsa05water_reduction!
β_WAT_rsawater_reduction!
δ_WAT_rsawater_reduction!
ω_NUT_totalNinput_nutrients!
ω_NUT_fertilizationinput_nutrients!
β_TSsimilarity_matrix!
α_NUT_TSBnutrient_competition!
α_NUT_maxadjnutrient_competition!
α_NUT_amc05nutrient_reduction!
α_NUT_rsa05nutrient_reduction!
β_NUT_rsanutrient_reduction!
β_NUT_amcnutrient_reduction!
δ_NUT_rsanutrient_reduction!
δ_NUT_amcnutrient_reduction!
κ_ROOT_amcroot_investment!
κ_ROOT_rsaroot_investment!
γ_RAD1radiation_reduction!
γ_RAD2radiation_reduction!
ω_TEMP_T1temperature_reduction!
ω_TEMP_T2temperature_reduction!
ω_TEMP_T3temperature_reduction!
ω_TEMP_T4temperature_reduction!
ζ_SEA_ST1seasonal_reduction!
ζ_SEA_ST2seasonal_reduction!
ζ_SEAminseasonal_reduction!
ζ_SEAmaxseasonal_reduction!
α_SENinitialize_senescence_rate!
β_SEN_slainitialize_senescence_rate!
ψ_SEN_ST1seasonal_component_senescence
ψ_SEN_ST2seasonal_component_senescence
ψ_SENmaxseasonal_component_senescence
β_GRZ_lncgrazing!
β_GRZ_Hgrazing!
η_GRZgrazing!
κ_GRZgrazing!
ϵ_GRZ_minHgrazing!
β_SND_WHCinput_WHC_PWP!
β_SLT_WHCinput_WHC_PWP!
β_CLY_WHCinput_WHC_PWP!
β_OM_WHCinput_WHC_PWP!
β_BLK_WHCinput_WHC_PWP!
β_SND_PWPinput_WHC_PWP!
β_SLT_PWPinput_WHC_PWP!
β_CLY_PWPinput_WHC_PWP!
β_OM_PWPinput_WHC_PWP!
β_BLK_PWPinput_WHC_PWP!

How to change a parameter value

julia
import GrasslandTraitSim as sim
using Unitful

# default parameter values
sim.SimulationParameter()

# optimized/calibrated parameter values
sim.optim_parameter()

# you can change parameter values with keyword arguments, when you create the parameter object
p = sim.SimulationParameter(γ_RUE_k  = 0.65,  ϕ_TRSA = 0.05u"m^2 / g")

# or you can change the parameter values after the object is created
p.ϕ_TAMC = 0.1
p.ϕ_sla = 0.01u"m^2 / g"

p
┌─────────────────────┬───────────────────┐
           Parameter │ Value             │
├─────────────────────┼───────────────────┤
│              ϕ_TRSA │ 0.05 m^2 g^-1     │
│              ϕ_TAMC │ 0.1               │
│               ϕ_sla │ 0.01 m^2 g^-1     │
│            γ_RUEmax │ 0.003 kg MJ^-1    │
│             γ_RUE_k │ 0.65              │
│          α_RUE_cwmH │ 0.95              │
│         α_WAT_rsa05 │ 0.9               │
│           β_WAT_rsa │ 7.0               │
│           δ_WAT_rsa │ 20.0 g m^-2       │
│        ω_NUT_totalN │ 0.1 kg g^-1       │
│ ω_NUT_fertilization │ 0.001 ha kg^-1    │
│                β_TS │ 1.0               │
│           α_NUT_TSB │ 15000.0 kg ha^-1  │
│        α_NUT_maxadj │ 10.0              │
│         α_NUT_amc05 │ 0.95              │
│         α_NUT_rsa05 │ 0.95              │
│           β_NUT_rsa │ 15.0              │
│           β_NUT_amc │ 15.0              │
│           δ_NUT_rsa │ 20.0 g m^-2       │
│           δ_NUT_amc │ 10.0              │
│          κ_ROOT_amc │ 0.02              │
│          κ_ROOT_rsa │ 0.01              │
│              γ_RAD1 │ 4.45e-6 ha MJ^-1  │
│              γ_RAD2 │ 50000.0 MJ ha^-1  │
│           ω_TEMP_T1 │ 4.0 °C            │
│           ω_TEMP_T2 │ 10.0 °C           │
│           ω_TEMP_T3 │ 20.0 °C           │
│           ω_TEMP_T4 │ 35.0 °C           │
│           ζ_SEA_ST1 │ 775.0 °C          │
│           ζ_SEA_ST2 │ 1450.0 °C         │
│            ζ_SEAmin │ 0.9               │
│            ζ_SEAmax │ 1.5               │
│               α_SEN │ 0.05              │
│           β_SEN_sla │ 1.5               │
│           ψ_SEN_ST1 │ 775.0 °C          │
│           ψ_SEN_ST2 │ 3000.0 °C         │
│            ψ_SENmax │ 1.5               │
│           β_GRZ_lnc │ 1.2               │
│             β_GRZ_H │ 2.0               │
│               η_GRZ │ 2.0               │
│               κ_GRZ │ 22.0 kg           │
│          ϵ_GRZ_minH │ 0.05 m            │
│           β_SND_WHC │ 0.5678            │
│           β_SLT_WHC │ 0.9228            │
│           β_CLY_WHC │ 0.9135            │
│            β_OM_WHC │ 0.6103            │
│           β_BLK_WHC │ -0.2696 cm^3 g^-1 │
│           β_SND_PWP │ -0.0059           │
│           β_SLT_PWP │ 0.1142            │
│           β_CLY_PWP │ 0.5766            │
│            β_OM_PWP │ 0.2228            │
│           β_BLK_PWP │ 0.02671 cm^3 g^-1 │
└─────────────────────┴───────────────────┘

API

GrasslandTraitSim.SimulationParameter Type

Parameter of the GrasslandTraitSim.jl model

source