chillR (version 0.75)

tempResponse_daily_list: Calculation of climatic metrics from lists of daily temperature records

Description

Wrapper for the tempResponse function, to facilitate its use on lists of daily temperature records, e.g. those produced by the temperature_generation function. Daily temperature records are converted into hourly records using either the stack_hourly_temps function or an empirical relationship between observed hourly temperatures and daily temperature extremes (see Empirical_hourly_temperatures for details). These hourly records are then used as input into the tempResponse function, to which most parameters are passed. See the documentation of tempResponse for more details.

Usage

tempResponse_daily_list(
  temperature_list,
  latitude,
  Start_JDay = 1,
  End_JDay = 366,
  models = list(Chilling_Hours = Chilling_Hours, Utah_Chill_Units = Utah_Model,
    Chill_Portions = Dynamic_Model, GDH = GDH),
  misstolerance = 50,
  whole_record = FALSE,
  empirical = NULL,
  mean_out = FALSE
)

Value

data frame showing totals for all specified models for the respective periods for all seasons included in the temperature records. Columns are Season, End_year (the year when the period ended) and Days (the duration of the period), as well as one column per model, which receives the same name as the function in the models list. If the weather input consisted of a list with elements stack and QC, the output also contains columns from QC that indicate the completeness of the weather record that the calculations are based on.

Arguments

temperature_list

list of daily temperature records, as produced by temperature_generation.

latitude

latitude of the location of interest (used for generating hourly records).

Start_JDay

the start date (in Julian date, or day of the year) of the period, for which chill and heat should be quantified.

End_JDay

the end date (in Julian date, or day of the year) of the period, for which chill and heat should be quantified.

models

named list of models that should be applied to the hourly temperature data. These should be functions that take as input a vector of hourly temperatures. This defaults to the set of models provided by the chilling function.

misstolerance

maximum percentage of values for a given season that can be missing without the record being removed from the output. Defaults to 50.

whole_record

boolean parameter indicating whether the metrics should be summed over the entire temperature record. If set to TRUE (default is FALSE), then the function ignores the specified start and end dates and simply returns the totals of each metric that accumulated over the entire temperature record.

empirical

indicates whether hourly temperatures should be generated based on an idealized temperature curve (set to NULL, the default) or an empirically derived relationship between hourly temperatures and daily temperature extremes (see Empirical_hourly_temperatures and Empirical_daily_temperature_curve, also for the format of the empirical prediction coefficient data.frame). If the latter, this parameter needs to be a data.frame including columns Month, Hour and Prediction_coefficients. See Empirical_daily_temperature_curve for further details on the format.

mean_out

boolean parameter indicating whether the mean of the input metric (e.g. temperature) should be returned in a column named "Input_mean".

Author

Eike Luedeling

References

The chillR package:

Luedeling E, Kunz A and Blanke M, 2013. Identification of chilling and heat requirements of cherry trees - a statistical approach. International Journal of Biometeorology 57,679-689.

Examples

weather<-fix_weather(KA_weather[which(KA_weather$Year>2006),])
temperature_list<-list(weather,weather,weather)

tempResponse_daily_list(temperature_list,latitude=50.4)