Introduction

Winter chill is essential for most perennial plants from cold climates (Erez 2000; Knight 1801; Samish 1954; Saure 1985; Vegis 1961). In order to avoid frost damage in the winter, such deciduous plants fall dormant in the fall and do not resume growth until a certain amount of winter chill has accumulated, fulfilling their chilling requirement. In spite of more than two centuries of research on the breaking of dormancy (starting with Knight 1801), the physiological and genetic processes occurring in trees during winter chill accumulation are poorly understood. Models of winter chill accumulation are thus purely empirical and based on either field observations (e.g., Linkosalo et al. 2008) or controlled temperature experiments (e.g., Fishman et al. 1987a) rather than on a functional understanding of tree physiology.

The study of plant phenology and how it responds to temperature is becoming increasingly important in the context of observed and future projected temperature increases (Chuine et al. 2010; Hudson 2010; Körner and Basler 2010). Many authors have reported on spring phenology changes in a wide range of species (e.g. Chmielewski and Rötzer 2001; Garcia-Mozo et al. 2010; Menzel et al. 2006; Schwartz 1999), and such changes are likely to continue, since all major climate models project progressively rising temperatures (IPCC 2007). Insufficient knowledge of species chilling requirements and a lack of accurate winter chill models are major limiting factors for the accuracy of projections of climate change effects on deciduous vegetation.

The desire to expand the production of deciduous fruits and nuts into warmer regions has motivated a lot of horticultural research on winter chill (e.g., Erez et al. 1990; Richardson et al. 1974; Weinberger 1950). Many commercially important fruit and nut trees, such as apples, pears, cherries, peaches, walnuts and almonds, require winter chill, and growers around the world select their tree cultivars according to the amount of winter chill they expect at their production sites. This selection process requires quantification of available winter chill and species chilling requirements, and a range of models have been developed for this purpose (e.g., Chandler 1942; Erez et al. 1990; Richardson et al. 1974). However, it is often historic precedence rather than systematic research that determines which model is used in a given growing region. Studies on winter chill and chilling requirements often tacitly assume that the choice of the model is not important, and that all models can be used interchangeably. Information about the chilling requirement of a given tree cultivar is thus often assumed to be valid in different locations, and summary tables of chilling requirements rarely specify where the data was obtained.

Recent research suggests that models may not be interchangeable (Luedeling et al. 2009d, e; Ruiz et al. 2007; Viti et al. 2010). Subjecting four common horticultural winter chill models to the same projected temperature changes for California, Luedeling et al. (2009d) showed that, depending on which model was used, winter chill changes between 1950 and 2050 ranged from −22 to −46%. This study showed substantial variation in the sensitivity of the tested models to climate change. Consequently, projections of winter chill that use the most sensitive Chilling Hours Model have projected rather dramatic declines in future winter chill in warm regions (Baldocchi and Wong 2008; Luedeling et al. 2009b), whereas changes indicated by the less sensitive Dynamic Model tend to be much less severe (Luedeling et al. 2009c, d). Compared to projected changes in warm regions, historical winter chill changes in relatively cool Germany were moderate, and differences between the Chilling Hours Model and the Dynamic Model were less pronounced (Luedeling et al. 2009a). Luedeling et al. (2009e) analyzed ratios between total accumulated winter chill during entire winter seasons, as calculated with different winter chill models, for eight locations in California and for all years between 1952 and 2008. Ratios varied substantially between locations and between years. Particularly strong differences were detected between temperature regimes prevailing in a cold chamber (constant temperatures of 6°C) and ambient orchard temperatures, indicating that chilling requirements determined by exposing trees to constant low temperatures may not be valid under field conditions (Luedeling et al. 2009e).

It thus appears that the choice of the model that is used for quantifying winter chill is important, and that assuming that chilling requirements determined in a given location are valid elsewhere is problematic. Nevertheless, converting between units of different models using site-specific conversion factors may be possible. At present, however, site-specific conversion factors have been determined only for a few locations and there is no information on the variability of these factors around the world. The objective of this paper is to assess the comparability of several winter chill models and to determine site-specific conversion factors, based on a global temperature dataset.

Materials and methods

Weather data

Daily temperature and precipitation records were downloaded for all 11,361 available weather stations at the National Climatic Data Center (NCDC 2008). This dataset was filtered, removing all stations that had less than 5,000 daily records between 1973 and 2002, and excluding all stations with more than 25% of daily minimum or maximum temperatures or 50% of daily rainfall data missing. For the remaining 5,078 weather stations, all available temperature and precipitation records were used to calculate site parameters for use in the LARS-WG stochastic weather generator (Semenov 2008). This tool describes the weather at a given location in statistical terms, with separate statistics for observed dry and wet periods. This information can then be used to generate synthetic weather records for the site. For each station, we evaluated the entire weather record for all days between 1973 and 2002, calculating separate linear regression equations for minimum temperature, maximum temperature and precipitation for each month of the year. Based on these regression analyses, we developed a climate scenario representing typical climatic conditions for the year 2000. This scenario contains the typical deviation of monthly minimum and maximum temperatures and precipitation from the mean of the calibration period (1973–2002) for each station, as calculated by the regression equations. Using the weather generator, we then generated 101 years of synthetic weather data for the ‘2000’ scenario. This data set can be interpreted as 101 replicates of temperature and precipitation for the year 2000, with a random seed introducing variability that corresponds to variability observed during the calibration period. For each weather station and each day of the 101-year period for each scenario, hourly temperatures were then generated using the idealized temperature curve proposed by Linvill (Linvill 1990; Luedeling et al. 2009a). Linvill’s equations, which use a sine curve for daytime temperatures and a logarithmic decline curve for nighttime cooling, require sunset and sunrise hours, as well as daylength, as input parameters. These data were generated using the equations of Spencer (1971) and Almorox et al. (2005). Resulting from these processing steps were 101 years of continuous hourly weather records for each weather station, representing typical climatic conditions around the year 2000.

Winter chill models

We calculated winter chill using three important winter chill models: (1) the Chilling Hours Model, 2) the Utah Model and 3) the Dynamic Model. The Chilling Hours Model (Chandler 1942) is the oldest method to quantify winter chill that is still widely used. It considers all hours with temperatures between 0 and 7.2°C as equally effective for winter chill accumulation. The number of Chilling Hours at time t (CHt; t is measured in hours since the start of the dormancy season) can be calculated as:

$$ C{H_t} = \sum\limits_{i = 1}^t {{T_{7.2}}}, \,{\hbox{with}}\,{T_{7.2}} = \left\{ {\begin{array}{*{20}{c}} {0^\circ C < T < 7.2^\circ C} \hfill & {:1} \hfill \\{else} \hfill & {:0} \hfill \\\end{array} } \right. $$

The Utah Model (Richardson et al. 1974) contains a weight function assigning different chilling efficiencies to different temperature ranges, including negative contributions by high temperatures. While a few modified versions of the weight function exist, we only used the weights from the original publication (Richardson et al. 1974), which are most widely used. The number of Utah Chill Units at time t (UCUt) can be expressed as:

$$ UC{U_t} = \sum\limits_{i = 1}^t {{T_U}}, \,{\hbox{with}}\,{T_U} = \left\{ {\begin{array}{*{20}{c}} {T \leqslant 1.4^\circ C} \hfill & {:0} \hfill \\{1.4^\circ C < T \leqslant 2.4^\circ C} \hfill & {:0.5} \hfill \\{2.4^\circ C < T \leqslant 9.1^\circ C} \hfill & {:1} \hfill \\{9.1^\circ C < T \leqslant 12.4^\circ C} \hfill & {:0.5} \hfill \\{12.4^\circ C < T \leqslant 15.9^\circ C} \hfill & {:1} \hfill \\{15.9^\circ C < T \leqslant 18.0^\circ C} \hfill & {: - 0.5} \hfill \\{T \geqslant 18.0^\circ C} \hfill & {:1} \hfill \\\end{array} } \right. $$

The Dynamic Model (Erez et al. 1990; Fishman et al. 1987a, b), which was originally developed for the warm winters in Israel, takes a different approach to quantifying winter chill. It assumes that winter chill results from a two-step process, in which an intermediate product is first formed in a process promoted by cold temperatures. Warm temperatures can destroy this intermediate product. As soon as a certain quantity of intermediate has accumulated, it is irreversibly transformed into a Chill Portion, which can no longer be destroyed. The equations for the number of Chill Portions at time t (CPt) are:

$$ \begin{gathered} {x_i} = \frac{{{e^{slp \cdot tetmlt \cdot \frac{{{T_K} - tetmlt}}{{{T_K}}}}}}}{{1 + {e^{slp \cdot tetmlt \cdot \frac{{{T_K} - tetmlt}}{{{T_K}}}}}}} \hfill \\{x_s} = \frac{{{a_0}}}{{{a_1}}} \cdot {e^{\frac{{{e_1} - {e_0}}}{{{T_K}}}}} \hfill \\a{k_1} = {a_1} \cdot {e^{ - \frac{{{e_1}}}{{{T_K}}}}} \hfill \\inte{r_E} = {x_s} - \left( {{x_s} - a{k_1}} \right) \cdot {e^{ - a{k_1}}} \hfill \\\end{gathered} $$
$$ inte{r_S} = \left\{ {\begin{array}{*{20}{c}} {t = {t_0}} \hfill & {:0} \hfill \\{t > {t_0} \wedge inte{r_{{E_{t - 1}}}} < 1} \hfill & {:inte{r_{{E_{t - 1}}}}} \hfill \\{t > {t_0} \wedge inte{r_{{E_{t - 1}}}} \geqslant 1} \hfill & {:inte{r_{{E_{t - 1}}}} \cdot (1 - {x_i})} \hfill \\\end{array} } \right. $$
$$ delt = \left\{ {\begin{array}{*{20}{c}} {t = {t_0}} \hfill & {:0} \hfill \\{t > {t_0} \wedge inte{r_E} < 1} \hfill & {:0} \hfill \\{t > {t_0} \wedge inte{r_E} \geqslant 1} \hfill & {:{x_i} \cdot inte{r_E}} \hfill \\\end{array} } \right. $$
$$ C{P_t} = \left\{ {\begin{array}{*{20}{c}} {t = {t_0}} \hfill & {:delt} \hfill \\{t \geqslant {t_0}} \hfill & {:delt + C{P_{t - 1}}} \hfill \\\end{array} } \right. $$

The experimentally derived constants slp, tetmlt, a 0, a 1, e 0 and e 1 were set to 1.6, 277, 139,500, 2.567 × 1018, 4,153.5 and 12,888.8, respectively, according to standard practice in horticultural applications. While it was the original intention of the model developers that constants be adjusted for each species and cultivar (Fishman et al. 1987a, b), this does not happen often in practice, and we therefore used only the commonly used values for all parameters.

Various other models, such as the Positive Utah Model (Linsley-Noakes and Allan 1994) or the North Carolina Model (Shaltout and Unrath 1983) were developed for the climatic conditions of specific regions and are not often used outside those regions. We therefore did not include these models in this study.

The units of the various models are Chilling Hours (CH) for the Chilling Hours Model, Utah Chill Units (UCU) for the Utah Model, and Chill Portions (CP) for the Dynamic Model. Winter chill was calculated for each of these models for all 101 years of all 5,078 weather stations, resulting in data for 100 winter seasons per station. Winter chill was summarized for the time period between 1 October and 1 May for all stations in the Northern Hemisphere, and between 1 April and 1 November for stations in the South. For each station, we then calculated Safe Winter Chill, which is defined as the 10% quantile of the resulting distribution over all 100 winters on record (Luedeling et al. 2009c). This metric defines the minimum amount of winter chill that growers of fruit and nut trees can expect in 90% of all years. It can be interpreted as the maximum recommendable chilling requirement of tree species or cultivars grown in a given location. Species or cultivars with chilling requirements that are higher than locally available Safe Winter Chill are likely to not fulfill these requirements in more than 10% of years, which may render their production uneconomical. For all combinations of two winter chill models, we then calculated the ratios between Safe Winter Chills as calculated with the two models, for each weather station location. Ratios were calculated between Chilling Hours and Chill Portions (CH/CP), between Utah Chill Units and Chill Portions (UCU/CP) and between Utah Chill Units and Chilling Hours (UCU/CH).

For two winter chill models to be proportional, this ratio should be constant regardless of the location of the weather station. To avoid spurious outliers in calculated ratios, only denominators of >100 CH or >10 CP were accepted. All calculations were programmed in JSL, the scripting language of JMP 8 (SAS Institute, Cary, NC). Using ArcGIS (ESRI, Redlands, CA), all data were spatially interpolated to a 0.1 degree resolution using 12-neighbor Kriging with a spherical semivariogram. In all maps, areas that were more than 4° away from the closest weather station were shaded out, since interpolations for these areas are unreliable. For ten major growing regions of Mediterranean fruits and nuts, means and standard deviations of winter chill metric ratios were summarized using spatial statistics.

Factors influencing winter chill ratios

The two main factors that have led researchers and growers to prefer one winter chill model over another are mean temperature of the site and continentality. For each climate scenario, we therefore calculated the mean temperature for every weather station based on the 101-year synthetic hourly weather records. As a proxy for continentality, we used the mean daily temperature range, calculated as the difference between mean daily maximum and minimum temperature over the entire temperature record. Effects of mean temperature and mean daily temperature range on winter chill metric ratios were determined by two-way analysis of variance with interaction (JMP 8, SAS Institute).

Results

Spatial distribution of Safe Winter Chill

Safe Winter Chill varied substantially around the world, for all three winter chill metrics analyzed (Fig. 1). According to the Dynamic and Utah Models, Safe Winter Chill was highest in the maritime climates around the North Sea in Europe, along the Pacific Coast of Oregon (US), Washington (US) and British Columbia (Canada), and on New Zealand’s South Island, where it exceeded 130 CP or 2,500 UCU, respectively. Safe Winter Chill of more than 2,000 CH was registered around the North Sea in Europe, in Washington and British Columbia, in Southern Chile and Argentina, and in the northern half of Honshu Island (Japan).

Fig. 1
figure 1

Global distribution of Safe Winter Chill, calculated with the Dynamic Model, the Chilling Hours Model and the Utah Model. The spatial interpolation is based on 5,078 point estimates at weather station locations, with poor coverage in the shaded areas (more than 4° distance from closest weather station)

Intermediate Safe Winter Chill (70–130 CP; 700–2,000 CH; 700–2,000 UCU) in the Northern Hemisphere prevails in most of the continental United States and Southern Canada, as well in a long strip spanning most of Europe, and extending throughout Central Asia and into China and Japan. In the Southern Hemisphere, regions of intermediate winter chill are much scarcer, and exist only at the Southern tip of South America and in Southeastern Australia, as well as likely in the higher reaches of the Andes, where no useable weather stations were located. Regions with low Safe Winter Chill (30–70 CP; 300–700 CH; −1,000–700 UCU) were detected in most Subtropical regions, as well as at high latitudes adjacent to the intermediate Safe Winter Chill regions. In spite of sufficient Safe Winter Chill, the latter areas are unlikely to be suitable for fruit and nut production due to very cold winters and lack of heat in the summer. All three winter chill models identified virtually no Safe Winter Chill in the Tropics.

Differences among models

For some combinations of two winter chill models, correlations between Safe Winter Chills were fairly good, with coefficients of determination between 0.65 and 0.87, but ratios between different winter chill metrics varied substantially in all cases (Fig. 2). Normalizing all Safe Winter Chill ratio distributions by the median illustrates that there is no clearly defined and globally valid conversion factor between the models (Fig. 3). For example, while the median of the ratio distribution was 14.4, one Chill Portion could correspond to anywhere between 0 and 51 CH, depending on local climate. While these are the extreme cases, 50% of all weather stations had CH / CP ratios of >16.4 or <11.4. For half of all weather stations, a chilling requirement of 70 CP would thus correspond to between 800 and 1,150 CH. For 80% of all stations, the corresponding number of CH would be between 570 and 1,281 (CH / CP of 8.2 and 18.3, respectively). Even between growing regions with seemingly similar climate and comparable tree species distributions, and even within these regions, winter chill metric ratios varied substantially, as shown in Table 1 for some major producing regions of Mediterranean fruits and nuts.

Fig. 2
figure 2

Correlations between Safe Winter Chill estimates, based on calculations for 5,078 weather stations around the world with three winter chill models. Box plots on the right show the distribution of ratios between the same two winter chill metrics among all weather stations

Fig. 3
figure 3

Ratios between the three Safe Winter Chill metrics, normalized by the median of the distributions among all 5,078 weather stations

Table 1 Variation of Safe Winter Chill ratios within and between major growing regions of Mediterranean fruit and nut trees. For each region and ratio, means and standard deviations are given. CH Chilling Hours, UCU Utah Chill Units, CP Chill Portions

This large variation illustrates the problems inherent in using an inappropriate winter chill model combined with estimates of a tree’s chilling requirement determined in a location with a different climate. To allow meaningful conversions between winter chill metrics, the conversion factors need to be specified based on climatic attributes of the growing region where the models are applied.

Patterns in winter chill metric ratios

Analyzing the geographic distribution of Safe Winter Chill ratios reveals clear spatial patterns (Fig. 4). Generally, ratios between CH and CP, as well as between UCU and CP, were highest at temperate latitudes, and often substantially lower towards the poles and Subtropics. The ratio of CH per CP increased gradually from 0 in the Tropics and boreal regions to more than 25 in parts of Scandinavia and Greenland. All ratios that included the Utah Model were negative in the growing regions of the warm Subtropics (Fig. 4). This indicates that, when winter chill is quantified in UCU, fruit and nut trees in these regions would have negative chilling requirements.

Fig. 4
figure 4

Global distributions of ratios between two Safe Winter Chill metrics, including winter chill estimates by the Dynamic Model, the Chilling Hours Model and the Utah Model. The spatial interpolation is based on 5,078 point estimates at weather station locations, with poor coverage in the shaded areas

Influence of temperature and continentality on winter chill metric ratios

The ratios between different winter chill metrics clearly varied with mean daily temperature of the site (T mean) and mean daily temperature range (T range; Table 2). The relationships are not very well defined, as evidenced by low coefficients of determination, and visual interpretation of the relationships between Safe Winter Chill metric ratios and the individual factors shows the potentials and limitations of conversion between models more clearly (Fig. 5).

Table 2 Results of the two-way analysis of variance with interaction between the factors mean temperature (T mean) and mean daily temperature range (T range) and the Safe Winter Chill ratios (CH / CP, UCU / CP and UCU / CH). RMSE  Root mean square error
Fig. 5
figure 5

Relationship of all calculated Safe Winter Chill ratios with mean temperature and mean daily temperature range

None of the Safe Winter Chill metric ratios were constant over the whole range of mean temperatures or temperature ranges. All calculated ratios showed strong shifts between mean temperatures of −10 and about 8°C and erratic behavior at higher temperatures. The response to T range was much less clearly defined, but was characterized by decreasing ratios with increasing T range.

Discussion

All ratios between winter chill metrics varied substantially around the globe, showing clearly that the models are not proportional. There are large regions of relatively homogeneous Safe Winter Chill ratios, but in those areas where winter chill is most critical, in the warmer growing regions, ratios vary even over relatively short distances. This strongly suggests that a given tree cultivar may not require the same number of CH or UCU when grown in different climatic zones. This is particularly problematic when transferring tree species or cultivars from temperate to Subtropical growing regions. Especially in Mediterranean climates, conversion factors between models differed greatly, in particular for ratios that include the Utah Model, with its winter chill negation mechanism (Table 1). For areas with mean temperatures above approximately 12°C, chill negation is clearly overemphasized, resulting in negative conversion factors, even in major growing regions, such as the Mediterranean region or South Africa.

To our knowledge, all studies comparing the accuracy of winter chill models that included the Dynamic Model have found this model to be superior or equivalent to all other models. Evidence has been presented for Israel (Erez 2000; Erez et al. 1990; Fishman et al. 1987a, b), South Africa (Allan et al. 1995, 1997; Erez et al. 1990; Gardner and Bertling 2005), Chile (Perez et al. 2008), France (Balandier et al. 1993), Spain (Alburquerque et al. 2008; Ruiz et al. 2006, 2007) and California (Luedeling et al. 2009e). All horticultural winter chill models, including the Dynamic Model, fail to explain peach phenology in the Tropical climate of Réunion Island (Balandier et al. 1993), and none of the models allow reliable prediction of spring phenology (Luedeling et al. 2009e), indicating that no model perfectly describes tree physiology. Phenology modeling strategies from other disciplines, such as the Unified Model (Chuine 2000), have been successful in explaining budburst of temperate trees, after detailed calibration with observed phenology over many years or locations, but these models do not provide generally applicable chilling functions that can be used in warmer regions. Among the three winter chill models tested in this study, the Dynamic Model seems most accurate in warmer climates. We are not aware of any tests of the applicability of this model in cold or temperate regions. For such regions, however, the conversion factors between the winter chill metrics are much more homogeneous than in the Subtropics, indicating that conversion between winter chill models is possible there at relatively predictable conversion factors. These factors can be estimated based on mean annual temperature, but site-specific calibration is clearly preferable due to much variation that is not explained by mean temperature or temperature range (Fig. 5).

Assuming for the purpose of this discussion that the Dynamic Model is really the best existing model for most growing regions (as the previous discussion suggests but cannot prove without further testing), it would be preferable to specify the chilling requirements of all tree cultivars in CP rather than CH or UCU, especially when transferring cultivars from one growing region to another. Unfortunately, chilling requirements are often only available in CH or UCU and thus probably not valid outside the location where they were determined. Our results suggest a new strategy for dealing with this problem. Known chilling requirements determined at a specific location should first be converted to CP using the site-specific conversion factor for the location of the experiment. The conversion factor can be obtained by determining winter chill with different metrics for several winter seasons and calculating the ratio between the results. The chilling requirement resulting from the conversion (then in CP) should then be more reliable and more universally valid than the original CH or UCU estimate. If local industries in target growing regions are more comfortable with CH or UCU, converting the CP requirement back into these units can then be done using the site-specific factors of the target region.

In using this strategy, however, one should keep in mind that conversion factors are not necessarily valid over time. According to the relationship between mean temperature and winter chill metric ratios (Fig. 5), temperature changes can shift these ratios. For six weather stations in California, Luedeling et al. (2009e) already detected a gradual decline of the ratio between CH and CP from 15 to under 12 between 1951 and 2008, in response to increasing winter temperatures. Assuming again that CP are a better measure of winter chill than CH, such a shift would indicate that chilling requirements in CH that were valid in the 1950s are not necessarily valid today. For example, a chilling requirement of 50 CP would correspond to 750 CH in the 1950s and 600 CH today. Using the full range of Safe Winter Chill ratios encountered during this time span indicates CH-requirements ranging from 401 to 864 CH. It is thus necessary that trends in winter chill and in Safe Winter Chill ratios are observed and estimates adjusted accordingly, as the climate changes. The fruit and nut industry would also benefit from enhanced research efforts on winter chill and chilling requirements. At present, the physiological and genetic processes determining tree dormancy are inadequately understood, and their temperature responses have not been well characterized. Until much more functional understanding of these processes is gained, all winter chill models will only be proxies of the true biological processes. While all available models are useful for practical guidance of orchard management, they still leave a lot to be desired. In our opinion, more research on the effects of temperature on tree phenology, in particular during the dormancy season, is urgently needed, in order to prepare tree crop growers for the likely impacts of future climate change.

Conclusions

Our study showed that winter chill models are not equivalent, and that conversion factors between different winter chill metrics vary substantially around the globe. Data on chilling requirements should thus always be supplemented with information on the location, study duration and study conditions under which the requirements were determined. Only with this additional information will it be possible to compile data from various locations into a database that will be useful for matching tree cultivars with climate, a process that is likely to increase in importance as climate change accelerates.