growing_degree_days: Degree days as the sum of capped average daily temperature...

Description Usage Arguments Details Value References See Also Examples

View source: R/growing_degree_days.R

Description

Calculate the sum of of degree days (average temperature - base temperature base_temp for each day) during specified dates for a tibble or data frame of daily weather data. Alternatively, a number of days before or after a specific date may be specified. Note that the maximum temperature is capped at maxt_cap when calculating the average temperature.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
growing_degree_days(
  data,
  varmax = NULL,
  varmin = NULL,
  datevar = NULL,
  maxt_cap = 30,
  base_temp = 5,
  na.rm = FALSE,
  ndays = 5,
  startdate = NULL,
  enddate = NULL,
  monitor = FALSE,
  warn.consecutive = TRUE
)

Arguments

data

Tibble or dataframe of daily weather data

varmax

Name of variable containing max temp (default 'maxt')

varmin

Name of variable containing min temp (default 'mint')

datevar

Date variable specifying day (Default: date_met)

maxt_cap

A numeric value set to the temperature considered to be the maximum necessary for plant growth. Maximum temperature is capped at this value for calculating average daily temp (Default: 30)

base_temp

A numeric value set to the temperature considered to be the minimum necessary for plant growth (Default: 5)

na.rm

Used for calculations (Default: FALSE)

ndays

Number of days after/before the start or end date, respectively. Ignored of both the start and end dates are specified (Default: 5)

startdate

Start date of data to be extracted

enddate

Final date of data to be extracted

monitor

For debugging. Prints data and dates. (Default: FALSE)

warn.consecutive

A logical indicating whether to check that dates are consecutive, that none are missing and provide a warning if not (Default:TRUE)

Details

The value returned is the sum of of degree days GDD = ∑_i (Tmax_i + Tmin_i) / 2 - Tbase during specified dates for a tibble or data frame of daily weather data. The maximum temperature Tmax is capped at maxt_cap degrees when calculating average temp (see https://farmwest.com/climate/calculator-information/gdd/ or McMaster, GS, & Wilhelm, WW (1997)). Baskerville, G & Emin, P (1969) provide variations on this method.

The sum of degree days is returned but if there are any missing values, then the value returned will be NA since the default na.rm is TRUE. Note that if there are any missing dates, then a warning is issued but the sum of non-missing values is returned.

If any values are missing, while the sum of degree days may prove useful, it will not include all the data and so will lead to biased underestimates. Hence, in these cases it is unlikely that the value returned is a good estimate but the appropriateness of the estimate will depend on the exact circumstances of the missing data and so this decision is left to the user.

Value

Numeric variable containing the sum of degree days during the period

References

See Also

cumulative, daily_mean, stress_days_over, and weather_extract

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
## Selected calculations
library(tidyverse)
growing_degree_days(boonah, enddate = crop$flower_date[4], ndays = 3,
                    varmax = maxt, varmin = mint,
                    monitor = TRUE)
growing_degree_days(boonah, enddate = crop$harvest_date[4], ndays = 3,
                    varmax = maxt, varmin = mint,
                    monitor = TRUE)
growing_degree_days(boonah, startdate = crop$flower_date[4],
                    varmax = maxt, varmin = mint,
                    enddate = crop$harvest_date[4], monitor = TRUE)

## Add selected growing degree days at 'boonah' to 'crop' tibble
## using 'map2_dbl' from the 'purrr' package
## Note: using equivalent 'furrr' functions can speed up calculations 
crop2 <- crop %>%
  mutate(gddays8_post_sow_7d =
          map_dbl(sowing_date, function(x)
            growing_degree_days(boonah, startdate = x, ndays = 7,
                                base_temp = 8)),
          gddays_flower_harvest =
            map2_dbl(flower_date, harvest_date, function(x, y)
              growing_degree_days(boonah, startdate = x, enddate = y)))
crop2

cropgrowdays documentation built on Dec. 11, 2021, 9:23 a.m.