russet | R Documentation |
This function assesses the risk of russet in fruit skins. The risk is defined by the number of hours with the relative humidity (RH) above a threshold during a given period. For reference, in 'Conference' pear the risk is defined by the number of hours with RH> 75% from 12 to 30 days after full bloom (Alegre, 2013). In 'Golden' apple, the risk is defined by the number of hours with RH> 55% from 30 to 34 days after full bloom (Barcelo-Vidal et al., 2013). The function requires hourly temperatures and humidity, if only daily data is available, the function hourly_RH can be used to estimate them.
russet(climdata, fendata, RH_crit, init_d, end_d)
climdata |
a dataframe with hourly temperature and RH data. Required columns are Date, Year, Month, Day, DOY (julian day), Hour and RH. |
fendata |
a dataframe with julian day of occurrence of the full bloom (F2) phenological stage. Must contain the columns Year and Fday in that order. |
RH_crit |
the relative humidity threshold |
init_d |
the initial date (as days after full bloom) of the sensitive period |
end_d |
the end date (as days after full bloom) of the sensitive period |
data frame with the number of risk hours (Russet_hours) in the sensitive period for each year in the series.
Carlos Miranda, carlos.miranda@unavarra.es
Alegre S. 2013. Tecnicas de cultivo. In. VII Foro INIA "adaptacion a cambio climatico en la produccion fruticola de hueso y pepita". Madrid, Spain, pp 1-18 Barcelo-Vidal C, Bonany J, Martin-Fernandez JA and Carbo J. 2013. Modelling of weather parameters to predict russet on 'Golden Delicious' apple. J. Hort. Sci. Biotech. 88: 624-630.
# Select the appropiate columns from the example dataset
# Dates_BT and rename column names to make the file compatible
# with the function
library(magrittr)
library(dplyr)
library(lubridate)
Bloom <- Dates_BT %>%
select(Year, sbloom) %>%
rename(Fday=sbloom) %>%
filter(Year==2003)
# Obtain estimated hourly RH from the example dataset Tudela_DW
Weather <- Tudela_DW %>%
filter (Tudela_DW$Year==2003)
RH_h <- hourly_RH(Weather, 42.13132)
# Estimate the number of russet-inducing days for a RH>55\%
# between 30 to 34 days after full bloom for each season
Russet_Risk <-russet(RH_h,Bloom,55,30,34)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.