This page describes the parameters that can be used to customize PHENIPS. The model was developed by Baier2007;textualbarrks. Look here to find out how to apply the model.
The date, when the degree days start to sum up ('MM-DD').
Base temperature to calculate degree days to trigger the onset.
Degree days that are required to trigger the onset of
infestation. Additionally, the maximum temperature must exceed tfly.
Minimum temperature that beetles need to fly.
Base temperature to calculate degree days for calculating the beetles development.
Degree days that are required for a generation to fully develop
Share in total development when the egg development starts and the juvenile beetle's development ends respectively. Usable if the development below/above these thresholds should account for mating, oviposition etc.
Share in the total development when a sister brood will be established.
The beetles are considered to be in
white stages (egg, larva, pupa) if their development exceeds dev_mortal_min
and subceeds dev_mortal_max. During these stages, the beetles could die
due to a mortality event. NULL means that no lower/upper threshold is
defined.
Temperature for optimal development.
Temperature below/above which no development happens.
Functions to calculate the effective bark temperature (see Baier2007;nobracketsbarrks, equations A.3 to A5). Each parameter will be passed as SpatRaster:
tmean: mean air temperature
tmax: maximum air temperature
rad: radiation
btmax: maximum bark temperature
Date when the model ends (no further development will be modeled).
When the daylength falls below this threshold, diapause will be initiated.
Date when all white stages (egg, larva, pupa) die.
In barrks, model() is used to customize a model. The following code
illustrates which parameters are available for PHENIPS and specifies their
default values.
model("phenips", # ==== onset ====
dd_onset_start_date = '04-01',
dd_onset_base = 8.3,
dd_onset_threshold = 140,
# ==== onset + development ====
tfly = 16.5,
# ==== development ====
dd_development_base = 8.3,
dd_total_dev = 557,
dev_start = 0,
dev_end = 1,
dev_sister_brood = 0.5,
dev_mortal_min = NULL,
dev_mortal_max = 0.6,
topt = 30.4,
tlow = 8.3,
tup = 38.9,
func_btmean = function(tmean, rad) {
-0.173 + 0.0008518 * rad + 1.054 * tmean
},
func_btmax = function(tmax, rad) {
1.656 + 0.002955 * rad + 0.534 * tmax + 0.01884 * tmax ^ 2
},
func_btdiff = function(btmax) {
(-310.667 + 9.603 * btmax) / 24
},
model_end_date = '10-31',
# ==== diapause ====
daylength_dia = 14.5,
# ==== mortality ====
mortality_date = '10-31'
)
model(), phenology(), model.phenips.apply
Other model customizations:
model.bso.customize,
model.chapy.customize,
model.joensson.customize,
model.lange.customize,
model.phenips_clim.customize,
model.rity.customize