Document Detail

 A global analysis of the comparability of winter chill models for fruit and nut trees. Jump to Full Text MedLine Citation: PMID:  20730614     Owner:  NLM     Status:  In-Data-Review Abstract/OtherAbstract: Many fruit and nut trees must fulfill a chilling requirement to break their winter dormancy and resume normal growth in spring. Several models exist for quantifying winter chill, and growers and researchers often tacitly assume that the choice of model is not important and estimates of species chilling requirements are valid across growing regions. To test this assumption, Safe Winter Chill (the amount of winter chill that is exceeded in 90% of years) was calculated for 5,078 weather stations around the world, using the Dynamic Model [in Chill Portions (CP)], the Chilling Hours (CH) Model and the Utah Model [Utah Chill Units (UCU)]. Distributions of the ratios between different winter chill metrics were mapped on a global scale. These ratios should be constant if the models were strictly proportional. Ratios between winter chill metrics varied substantially, with the CH/CP ratio ranging between 0 and 34, the UCU/CP ratio between -155 and +20 and the UCU/CH ratio between -10 and +5. The models are thus not proportional, and chilling requirements determined in a given location may not be valid elsewhere. The Utah Model produced negative winter chill totals in many Subtropical regions, where it does not seem to be useful. Mean annual temperature and daily temperature range influenced all winter chill ratios, but explained only between 12 and 27% of the variation. Data on chilling requirements should always be amended with information on the location and experimental conditions of the study in which they were determined, ideally including site-specific conversion factors between winter chill models. This would greatly facilitate the transfer of such information across growing regions, and help prepare growers for the impact of climate change. Authors: Eike Luedeling; Patrick H Brown Related Documents : 17851524 - Spatial vegetation patterns and imminent desertification in mediterranean arid ecosystems.14967654 - Frost injury prediction model for douglas-fir seedlings in the pacific northwest.12859004 - Incorporating natural capital into economy-wide impact analysis: a case study from albe...12420964 - Environmetric approaches to estimate pollution impacts on a coastal area by sediment an...17968114 - Scalable hybrid unstructured and structured grid raycasting.18191524 - Investigating the neural basis for fmri-based functional connectivity in a blocked desi... Publication Detail: Type:  Journal Article     Date:  2010-08-22 Journal Detail: Title:  International journal of biometeorology     Volume:  55     ISSN:  1432-1254     ISO Abbreviation:  Int J Biometeorol     Publication Date:  2011 May Date Detail: Created Date:  2011-04-15     Completed Date:  -     Revised Date:  - Medline Journal Info: Nlm Unique ID:  0374716     Medline TA:  Int J Biometeorol     Country:  United States Other Details: Languages:  eng     Pagination:  411-21     Citation Subset:  IM Affiliation: Department of Plant Sciences, University of California Davis, Mail Stop Two, One Shields Ave, Davis, CA, 95616, USA, e.luedeling@cgiar.org. Export Citation: APA/MLA Format     Download EndNote     Download BibTex MeSH Terms Descriptor/Qualifier:

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

 Full Text Journal Information Journal ID (nlm-ta): Int J Biometeorol ISSN: 0020-7128 ISSN: 1432-1254 Publisher: Springer-Verlag, Berlin/Heidelberg Article Information Download PDF © The Author(s) 2010 Received Day: 9 Month: 2 Year: 2010 Revision Received Day: 25 Month: 7 Year: 2010 Accepted Day: 25 Month: 7 Year: 2010 Electronic publication date: Day: 22 Month: 8 Year: 2010 pmc-release publication date: Day: 22 Month: 8 Year: 2010 Print publication date: Month: 5 Year: 2011 Volume: 55 Issue: 3 First Page: 411 Last Page: 421 ID: 3077742 PubMed Id: 20730614 Publisher Id: 352 DOI: 10.1007/s00484-010-0352-y
 A global analysis of the comparability of winter chill models for fruit and nut trees Eike Luedeling12 Address: +254-20-7224120 +254-20-7224001 e.luedeling@cgiar.org Patrick H. Brown1 1Department of Plant Sciences, University of California Davis, Mail Stop Two, One Shields Ave, Davis, CA 95616 USA 2World Agroforestry Centre (ICRAF), United Nations Ave, Gigiri, P.O. Box 30677-00100, Nairobi, Kenya

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:

[Formula ID: Equa]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$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.$$\end{document}]

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:

[Formula ID: Equb]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$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.$$\end{document}]

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:

[Formula ID: Equc]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\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}$$\end{document}]
[Formula ID: Equd]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$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.$$\end{document}]
[Formula ID: Eque]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$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.$$\end{document}]
[Formula ID: Equf]
[\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$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.$$\end{document}]

The experimentally derived constants slp, tetmlt, a0, a1, e0 and e1 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).

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.

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.

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 (Tmean) and mean daily temperature range (Trange; 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).

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 Trange was much less clearly defined, but was characterized by decreasing ratios with increasing Trange.

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.

We greatly appreciate the help of Mikhail Semenov at Rothamsted Research, UK, for improving his LARS-WG weather generator to suit our data processing needs. We also acknowledge support by Heiner Lieth at UC Davis, who provided computing power, and by Minghua Zhang, UC Davis, and the California Department for Food and Agriculture for supporting earlier work on winter chill models. We acknowledge the journal field editor and three anonymous reviewers for providing invaluable comments that greatly improved the manuscript.

References
 Alburquerque N,García-Montiel F,Carrillo A,Burgos L. Chilling and heat requirements of sweet cherry cultivars and the relationship between altitude and the probability of satisfying the chill requirementsEnviron Exp BotYear: 200864216217010.1016/j.envexpbot.2008.01.003 Allan P,et al. Kiwifruit research in a subtropical areaActa HorticYear: 19974443742 Allan P,Rufus G,Linsley-Noakes GC,Matthee GW. Winter chill models in a mild subtropical area and effects of constant 6°C chilling on peach budbreakActa HorticYear: 1995409917 Almorox J,Hontoria C,Benito M. Statistical validation of daylength definitions for estimation of global solar radiation in Toledo, SpainEnergy Convers ManageYear: 2005469–101465147110.1016/j.enconman.2004.07.007 Balandier P,Bonhomme M,Rageau R,Capitan F,Parisot E. Leaf bud endodormancy release in peach-trees—evaluation of temperature models in temperate and tropical climatesAgric For MeteorolYear: 1993671–29511310.1016/0168-1923(93)90052-J Baldocchi D,Wong S. Accumulated winter chill is decreasing in the fruit growing regions of CaliforniaClim ChangeYear: 200887S153S16610.1007/s10584-007-9367-8 Chandler WH. Deciduous orchardsYear: 1942PhiladelphiaLea & Febiger Chmielewski FM,Rötzer T. Response of tree phenology to climate change across EuropeAgric For MeteorolYear: 2001108210111210.1016/S0168-1923(01)00233-7 Chuine I. A unified model for budburst of treesJ Theor BiolYear: 2000207333734710.1006/jtbi.2000.217811082304 Chuine I,Morin X,Bugmann H. Warming, photoperiods, and tree phenologyScienceYear: 2010329598927727810.1126/science.329.5989.277-e20647448 Erez A. Erez ABud dormancy; phenomenon, problems and solutions in the tropics and subtropicstemperate fruit crops in warm climatesYear: 2000DordrechtKluwer1748 Erez A,Fishman S,Linsley-Noakes GC,Allan P. The dynamic model for rest completion in peach budsActa HorticYear: 1990276165174 Fishman S,Erez A,Couvillon GA. The temperature-dependence of dormancy breaking in plants - computer-simulation of processes studied under controlled temperaturesJ Theor BiolYear: 1987126330932110.1016/S0022-5193(87)80237-0 Fishman S,Erez A,Couvillon GA. The temperature dependence of dormancy breaking in plants: mathematical analysis of a two-step model involving a cooperative transitionJ Theor BiolYear: 1987124447348310.1016/S0022-5193(87)80221-7 Garcia-Mozo H,Mestre A,Galan C,Data C. Phenological trends in southern Spain: a response to climate changeAgric For MeteorolYear: 2010150457558010.1016/j.agrformet.2010.01.023 Gardner RAW,Bertling I. Effect of winter chilling and paclobutrazol on floral bud production in Eucalyptus nitensS Afr J BotYear: 2005712238249 Hudson LI. Interdisciplinary approaches: towards new statistical methods for phenological studiesClim ChangeYear: 2010100114317110.1007/s10584-010-9859-9 IPCC (2007) Climate Change 2007 - Synthesis Report. Contributions of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Intergovernmental Panel on Climate Change, Geneva, Switzerland. Knight TA. Account of some experiments on the ascent of the sap in treesPhilos Trans R Soc LondYear: 18019133335310.1098/rstl.1801.0017 Körner C,Basler D. Phenology under global warmingScienceYear: 201032759721461146210.1126/science.118647320299580 Linkosalo T,Lappalainen HK,Hari P. A comparison of phenological models of leaf bud burst and flowering of boreal trees using independent observationsTree PhysiolYear: 200828121873188219193570 Linsley-Noakes GC,Allan P. Comparison of 2 models for the prediction of rest completion in peachesSci HorticYear: 199459210711310.1016/0304-4238(94)90077-9 Linvill DE. Calculating chilling hours and chill units from daily maximum and minimum temperature observationsHortscienceYear: 19902511416 Luedeling E,Blanke M,Gebauer J. Auswirkungen des Klimawandels auf die Verfügbarkeit von Kältewirkung (Chilling) für Obstgehölze in DeutschlandErwerbs-ObstbauYear: 200951819410.1007/s10341-009-0085-4 Luedeling E,Gebauer J,Buerkert A. Climate change effects on winter chill for tree crops with chilling requirements on the Arabian PeninsulaClim ChangeYear: 20099621923710.1007/s10584-009-9581-7 Luedeling E,Zhang M,Girvetz EH. Climatic changes lead to declining winter chill for fruit and nut trees in California during 1950–2099PLoS ONEYear: 200947e616610.1371/journal.pone.000616619606220 Luedeling E,Zhang M,Luedeling V,Girvetz EH. Sensitivity of winter chill models for fruit and nut trees to climate changeAgric Ecosyst EnvironYear: 2009133233110.1016/j.agee.2009.04.016 Luedeling E,Zhang M,McGranahan G,Leslie C. Validation of winter chill models using historic records of walnut phenologyAgric For MeteorolYear: 20091491854186410.1016/j.agrformet.2009.06.013 Menzel A,et al. European phenological response to climate change matches the warming patternGlob Chang BiolYear: 200612101969197610.1111/j.1365-2486.2006.01193.x NCDC (2008) Database of the National Climatic Data Center. National Oceanic and Atmospheric Administration. http://www.noaa.gov/ Perez FJ,Ormeno NJ,Reynaert B,Rubio S. Use of the dynamic model for the assessment of winter chilling in a temperate and a subtropical climatic zone of ChileChil J Agric ResYear: 200868198206 Richardson EA,Seeley SD,Walker DR. A model for estimating the completion of rest for Redhaven and Elberta peach treesHortscienceYear: 197494331332 Ruiz D,Campoy JA,Egea J. Chilling requirements of apricot varietiesActa HorticYear: 20067176770 Ruiz D,Campoy JA,Egea J. Chilling and heat requirements of apricot cultivars for floweringEnviron Exp BotYear: 20076125426310.1016/j.envexpbot.2007.06.008 Samish RM. Dormancy in woody plantsAnnu Rev Plant Physiol Plant Mol BiolYear: 19545183204 Saure MC. Dormancy release in deciduous fruit treesHortic RevYear: 19857239300 Schwartz MD. Advancing to full bloom: planning phenological research for the 21st centuryInt J BiometeorolYear: 199942311311810.1007/s004840050093 Semenov MA. Simulation of extreme weather events by a stochastic weather generatorClim ResYear: 200835320321210.3354/cr00731 Shaltout AD,Unrath CR. Rest completion prediction model for Starkrimson Delicious applesJ Am Soc Hortic SciYear: 19831086957961 Spencer JW. Fourier series representation of the position of the SunSearchYear: 197125172 Vegis A (1961) Samenkeimung und vegetative Entwicklung der Knospen. In: Ruhland W (ed) Handbuch der Pflanzenphysiologie, vol 16. Springer, Berlin, pp 168–298 Viti R,et al. Effect of climatic conditions on the overcoming of dormancy in apricot flower buds in two Mediterranean areas: Murcia (Spain) and Tuscany (Italy)Sci HorticYear: 2010124221722410.1016/j.scienta.2010.01.001 Weinberger JH (1950) Chilling requirements of peach varieties. Proc Am Soc Hortic Sci 56:122–128

Figures
 [Figure ID: Fig1] Fig. 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) [Figure ID: Fig2] Fig. 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 [Figure ID: Fig3] Fig. 3  Ratios between the three Safe Winter Chill metrics, normalized by the median of the distributions among all 5,078 weather stations [Figure ID: Fig4] Fig. 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 [Figure ID: Fig5] Fig. 5  Relationship of all calculated Safe Winter Chill ratios with mean temperature and mean daily temperature range

Tables
[TableWrap ID: Tab1] 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

Region CH / CP UCU / CP UCU / CH
California – Central Valley 13.0 ± 0.8 8.9 ± 5.3 0.6 ± 0.8
Northern Mexico 15.5 ± 1.6 −14.3 ± 18.8 −1.1 ± 1.3
Southeastern US (Carolinas, Georgia, Northern Florida) 14.4 ± 1.0 −10.4 ± 27.1 −0.7 ± 2.0
Chile 14.8 ± 0.8 5.5 ± 11.2 0.6 ± 0.6
South Africa – Western Cape 8.9 ± 1.1 −17.8 ± 0.7 −2.9 ± 0.3
Maghreb (Morocco, Tunisia, Northern Algeria) 9.5 ± 1.7 −14.9 ± 9.8 −1.6 ± 1.0
Israel 7.0 ± 1.2 −23.8 ± 5.2 −0.5 ± 0.5
Northwestern Mediterranean (Portugal, Spain, Southern France, Italy) 11.3 ± 2.9 11.9 ± 9.1 1.1 ± 0.6
South Australia 7.2 ± 1.9 −13.1 ± 14.7 −0.3 ± 1.2
New Zealand 8.3 ± 1.1 15.8 ± 2.7 3.2 ± 0.1

[TableWrap ID: Tab2] Table 2

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

CH / CP
Summary of fit Parameter estimates
R2 0.27 Intercept 17.32 ± 0.23
Mean of response 13.79 Tmean −0.05 ± 0.01
RMSE 3.89 Trange −0.30 ± 0.02
N 4,170 (Tmean-8.61) × (Trange-10.05) 0.10 ± 0.00
ANOVA P value <0.0001
UCU / CP
Summary of fit Parameter estimates
R2 0.12 Intercept 30.90 ± 1.37
Mean of response 7.74 Tmean −1.11 ± 0.05
RMSE 23.47 Trange −1.36 ± 0.12
N 4,170 (Tmean-8.61) × (Trange-10.05) −0.02 ± 0.01
ANOVA P value <0.0001
UCU / CH
Summary of fit Parameter estimates
R2 0.17 Intercept 2.53 ± 0.11
Mean of response 0.57 Tmean −0.07 ± 0.00
RMSE 1.82 Trange −0.14 ± 0.01
N 4,066 (Tmean-8.56) × (Trange-10.14) −0.02 ± 0.00
ANOVA P value <0.0001
UCU / CP
Summary of fit Parameter estimates
R2 0.12 Intercept 30.90 ± 1.37
Mean of response 7.74 Tmean −1.11 ± 0.05
RMSE 23.47 Trange −1.36 ± 0.12
N 4,170 (Tmean-8.61) × (Trange-10.05) −0.02 ± 0.01
ANOVA P value <0.0001
UCU / CH
Summary of fit Parameter estimates
R2 0.17 Intercept 2.53 ± 0.11
Mean of response 0.57 Tmean −0.07 ± 0.00
RMSE 1.82 Trange −0.14 ± 0.01
N 4,066 (Tmean-8.56) × (Trange-10.14) −0.02 ± 0.00
ANOVA P value <0.0001

 Article Categories:Original Paper Keywords: Keywords Chilling Hours, Chilling requirement, Chill Portions, Dynamic Model, Utah Model, Winter chill.

Previous Document:  CT fluoroscopy-guided lung biopsy versus conventional CT-guided lung biopsy: a prospective controlle...
Next Document:  A rapid genotyping test for the simultaneous detection and subtyping of group B streptococci: the fr...