calc_vegperiod {LWFBrook90R}

Calculate the dates of budburst and beginning of leaf fall

Description

Wrapper for vegperiod

Usage

calc_vegperiod(
  budburst_method,
  leaffall_method,
  dates = NULL,
  tavg = NULL,
  out_yrs = NULL,
  budburstdoy.fixed = 121,
  leaffalldoy.fixed = 279,
  ...
)

Arguments

budburst_method

name of model for estimating budburst day of year. Either 'fixed' or one of the values accepted by the 'start.method'-argument of the function vegperiod.

leaffall_method

name of model for estimating day of year when leaffall begin. Either 'fixed' or one of the values accepted by the 'end.method'-argument of the function vegperiod.

dates

date vector passed to vegperiod, ignored if both leaffall_method and budburst_method = 'fixed'

tavg

vector of daily mean air temperature (deg C) passed to vegperiod, ignored if leaffall_method = 'fixed' and budburst_method = 'fixed'.

out_yrs

integer vector of the years to be returned. If not specified, values for the years in dates will be returned.

budburstdoy.fixed

vector of values to be returned if budburst_method = 'fixed'.

leaffalldoy.fixed

vector of values to be returned if leaffall_method = 'fixed'.

...

additional argument passed to vegperiod.

Value

a data.frame with columns year, start, end. If budburst_method = 'fixed' or leaffall_method = 'fixed', start and end contain the values specified in budburstdoy.fixed and leaffalldoy.fixed respectively.

Examples

# fixed budburst and leaffall doy
calc_vegperiod(out_yrs = 2001:2010,
               budburst_method = "fixed",
               leaffall_method = "fixed",
               budburstdoy.fixed = floor(runif(10, 120,130)),
               leaffalldoy.fixed = floor(runif(2, 260,280)))

# dynamic budburst and leaffall using air temperature
data(slb1_meteo)

calc_vegperiod(budburst_method = "Menzel",
               leaffall_method = "fixed",
               leaffalldoy.fixed = 280,
               dates = slb1_meteo$dates,
               tavg = slb1_meteo$tmean,
               species = "Fagus sylvatica",
               est.prev = 3)

calc_vegperiod(budburst_method = "Menzel",
               leaffall_method = "ETCCDI",
               dates = slb1_meteo$dates,
               tavg = slb1_meteo$tmean,
               species = "Quercus robur",
               est.prev = 3)