# Phenological model intercomparison for estimating grapevine budbreak date (Vitis vinifera L.) in Europe

Luisa Leolini, Sergi Costafreda-Aumedes, João A. Santos, Christoph Menz, Helder Fraga, Daniel Molitor, Paolo Merante, Jürgen Junk, Thomas Kartschall, Agnes Destrac-Irvine, et al.

► To cite this version:

Luisa Leolini, Sergi Costafreda-Aumedes, João A. Santos, Christoph Menz, Helder Fraga, et al. Phenological model intercomparison for estimating grapevine budbreak date (Vitis vinifera L.) in Europe. Applied Sciences, 2020, 10 (11), pp.1-22. 10.3390/app10113800. hal-02902577

HAL Id: hal-02902577

https://hal.inrae.fr/hal-02902577

Submitted on 20 Jul 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Phenological Model Intercomparison for Estimating Grapevine Budbreak Date (*Vitis vinifera* L.) in Europe

Luisa Leolini <fnref n="1" />,* ID, Sergi Costafreda-Aumedes <fnref n="1" /> ID, João A. Santos <fnref n="2" /> ID, Christoph Menz <fnref n="3" />, Helder Fraga <fnref n="2" /> ID, Daniel Molitor <fnref n="4" /> ID, Paolo Merante <fnref n="5" />, Jürgen Junk <fnref n="4" /> ID, Thomas Kartschall <fnref n="3" />, Agnès Destrac-Irvine <fnref n="6" />, Cornelis van Leeuwen <fnref n="6" /> ID, Aureliano C. Malheiro <fnref n="2" /> ID, José Eiras-Dias <fnref n="7" /> ID, José Silvestre <fnref n="7" /> ID, Camilla Dibari <fnref n="1" /> ID, Marco Bindi <fnref n="1" /> ID and Marco Moriondo <fnref n="8" /> ID

<figure><img src="image_2.png" /></figure>

check for
updates

Received: 16 April 2020; Accepted: 22 May 2020; Published: 29 May 2020

**Abstract:** Budbreak date in grapevine is strictly dependent on temperature, and the correct simulation of its occurrence is of great interest since it may have major consequences on the final yield and quality. In this study, we evaluated the reliability for budbreak simulation of two modeling approaches, the chilling-forcing (CF), which describes the entire dormancy period (endo- and eco-dormancy) and the forcing approach (F), which only describes the eco-dormancy. For this, we selected six phenological models that apply CF and F in different ways, which were tested on budbreak simulation of eight grapevine varieties cultivated at different latitudes in Europe. Although none of the compared models showed a clear supremacy over the others, models based on CF showed a generally higher estimation accuracy than F where fixed starting dates were adopted. In the latter models, the accurate simulation of budbreak was dependent on the selection of the starting date for forcing accumulation that changes according to the latitude, whereas CF models were independent. Indeed, distinct thermal requirements were found for the grapevine varieties cultivated in Northern and Southern Europe. This implies the need to improve modeling of the dormancy period to avoid under- or over-estimations of budbreak date under different environmental conditions.

**Keywords:** chilling-forcing models; dormancy; forcing models; phenology

<footnote n="1">¹ Department of Agriculture, Food, Environment and Forestry, DAGRI, University of Florence, UNIFI, 50144 Florence, Italy; sergi.costafredaaumedes@unifi.it (S.C.-A.); camilla.dibari@unifi.it (C.D.); marco.bindi@unifi.it (M.B.)</footnote>

<footnote n="2">² Centre for the Research and Technology of Agro-Environmental Sciences, CITAB, Universidade de Trás-os-Montes e Alto Douro, UTAD, 5000-801 Vila Real, Portugal; jsantos@utad.pt (J.A.S.); hfraga@utad.pt (H.F.); amalheir@utad.pt (A.C.M.)</footnote>

<footnote n="3">³ Potsdam Institute for Climate Impact Research, PIK, 14473 Potsdam, Germany; menz@pik-potsdam.de (C.M.); thomas@pik-potsdam.de (T.K.)</footnote>

<footnote n="4">⁴ Luxembourg Institute of Science and Technology, LIST, 4422 Luxembourg, Luxembourg; daniel.molitor@list.lu (D.M.); juergen.junk@list.lu (J.J.)</footnote>

<footnote n="5">5 Tropical Plant Production and Agricultural Systems Modelling—TROPAGS, University of Göttingen, 37077 Göttingen, Germany; paolo.merante@uni-goettingen.de</footnote>

<footnote n="6">6 EGFV, Bordeaux Sciences Agro, INRAE, University of Bordeaux, ISVV, F-33882 Villenave d'Ornon, France; agnes.destrac-irvine@inra.fr (A.D.-I.); vanleeuwen@agro-bordeaux.fr (C.v.L.)</footnote>

<footnote n="7">⁷ Instituto Nacional de Investigação Agrária e Veterinária, INIAV, Estação Vitivinícola Nacional, Quinta da Almoínha, 2565-191 Dois Portos, Portugal; eiras.dias@iniav.pt (J.E.-D.); jose.silvestre@iniav.pt (J.S.)</footnote>

<footnote n="8">⁸ Institute of Bioeconomy, National Research Council, CNR-IBE, 50019 Florence, Italy; marco.moriondo@cnr.it</footnote>

<footnote>* Correspondence: luisa.leolini@unifi.it; Tel.: +39-055-2755742</footnote>
# 1. Introduction

Grapevine phenological phases have a profound influence on the final grape and wine quality [1-3]. In particular, changes in the current climate conditions are strongly affecting the occurrence of the developmental stages and length of the grapevine growing cycle in the most renowned wine areas [4,5].

Several authors have reported the increasing temperature trend across the most traditional wine areas of the Mediterranean basin during the last decades. For example, Tomasi et al. [6] showed a 2.3 °C-temperature increase over the grape growing season in Veneto region (Italy) during the period 1964-2009, while van Leeuwen et al. [7] evidenced a 2.0 °C increase from 1951 to 2016 in the Bordeaux wine production area. Moreover, Ramos et al. [8] observed an overall warming from 1.0 to 2.2 °C during the grape growing season in north-eastern Spain. These changes are expected to greatly influence grapevine varietal suitability across the European region [9-12].

Since grapevine phenology is mainly temperature-driven [4], the current warming is affecting the phenological stages and shifting the annual growing cycle [13-16]. Among the most important developmental phases, budbreak represents the starting point of vegetative growth. A temperature increase can have a dual effect on this phase (anticipating or delaying its occurrence), with strong consequences on the following reproductive stages. Indeed, a shift in budbreak may determine poor shoot growth together with fewer clusters and an inhomogeneous grape development [17,18].

Under increasing temperature, a delay in budbreak phase can be determined by insufficient chilling during the dormancy period [19]. Indeed, the bud vegetative growth is limited by a first sub-period induced by physiological conditions (endo-dormancy), and then by a second sub-period in which bud opening is affected by the environmental conditions (eco-dormancy [17,20]). In this context, some studies showed that dormancy release is regulated by bud exposure period to temperatures between 0 to 10 °C [19,21], until fulfilment of the chilling requirement, and a subsequent period of temperatures (from 0 to 25 °C), until the reaching of the forcing requirement [22].

Budbreak occurrence being strongly influenced by the dormancy period, several studies attempted the implementation of phenological models for improving the prediction of this developmental stage [23-28]. Phenological models can be classified in two main categories: forcing and chilling-forcing models. Forcing models are simply based on the forcing unit accumulation, starting, in general, from a fixed date in the year. These models only focus on the description of the eco-dormancy period assuming that the chilling unit accumulation requirement was previously fulfilled and the endo-dormancy period had concluded. Instead, chilling-forcing models can explain both the endo- and eco-dormancy periods by taking into consideration the chilling unit and the forcing accumulation with respect to specific temperature thresholds [24,29,30].

Hence, the selection of one of these two categories for an accurate budbreak simulation often gave rise to a wide debate [29,30]. Chilling-forcing models can take into account the influences of temperature on the entire dormancy period, evidencing not only the shortening effect of warmer temperatures during eco-dormancy but also a possible delay of budbreak caused by a failed chilling requirement [13,31]. However, the limited improvement in the performances of these models and the easier implementation provided by forcing models, focus on the issue about the most suitable modeling approach for budbreak estimation. In modeling studies, the choice of phenological model plays an important role in the use of crop simulation models and in determining potential risk of late spring frost events. In crop simulation models, phenology defines the plant growing cycle with important consequences on future biomass partitioning [32]. In this context, an accurate estimation of budbreak stage may avoid randomly selecting a starting point for vegetative biomass accumulation and over- or under-detection of late spring frost events under different climate scenarios [31,33].

Thus, the phenological intercomparison may be useful to highlight the performances of different models in simulating budbreak [29,30]. In order to take into account the influences of the grapevine precocity cycle and the effects of environmental conditions on budbreak estimation, model intercomparison experiments should consider multiple grapevine varieties and study sites. Indeed, grapevine varieties are characterized by various thermal requirements and focusing on a specific variety
does not allow the range of different physiological responses derived from the impact of temperature on dormancy period to be highlighted [30]. In addition, the use of a single date as a starting point for the chilling and forcing accumulation [29] does not allow the influence of temperature on the beginning of endo- or eco-dormancy periods to be appreciated, especially when the model intercomparison is performed with phenological data of grapevine varieties collected in different environments.

The objective of this study is to evaluate the performances of six phenological models (BRIN tested using hourly and daily temperature) at reproducing the budbreak date of eight grapevine varieties collected at different latitudes in Europe. In order to consider the influences of temperature during the dormancy period on model representations, two starting dates were set during chilling-forcing and forcing models' simulation, respectively. The analysis of budbreak date estimation is thus addressed to obtain responses of grapevine varieties' behavior under different environments.

## 2. Materials and Methods

### 2.1. *Observations*

The phenological data were collected from 15 experimental vineyards in three European countries (France, Luxembourg and Portugal, see Figure 1 and Table 1). The weather data (i.e., daily maximum and minimum air temperature, °C) were freely available from weather services. More specifically, for France, the land-based daily weather data were collected by the National Centers for Environmental Information and then archived at NOAA (https://www.ncdc.noaa.gov/ [34]) and from a weather station near the VitAdapt project [35]. Hourly temperatures, rainfall, wind speed and reference evapotranspiration were recorded and stored in Climatik, a database managed by INRA, US 1116 Agroclim research station, 84914 Avignon, France (https://intranet.inra.fr/climatik_v2; limited access). For Luxembourg, the weather data were collected by the Administration des Services Techniques de l'Agriculture at the Institut Viti-Vinicole (https://agriculture.public.lu/) while the weather data for Portugal derived from the Instituto Nacional de Investigação Agrária e Veterinária (INIAV http://www.iniav.pt/) and the Associação para o Desenvolvimento da Viticultura Duriense (ADVID https://www.advid.pt/). The weather data were obtained from the closest and more representative stations for the vineyard weather conditions (Figure 1 and Table 1).

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Map of the vineyards and weather stations situated in three European countries (Luxembourg, France and Portugal). The vineyards are displayed with triangles and points in relation to the use of phenological data for model fitting or application. The arrows indicate that both vineyards and weather stations are located in the same site. 10 weather stations and 10 vineyards are situated in Douro valley (Portugal).</figcaption></figure>
<table><caption>**Table 1.** Observed weather (daily maximum and minimum air temperature, °C) and budbreak data of eight grapevine varieties used for model fitting, validation and application in different locations across Europe. The geographical positions of all weather stations and vineyards are identified with the latitude (Lat) and longitude (Long) coordinates (degrees).</caption><thead><tr><th colspan="3">Variety × Site</th><th colspan="4">Weather Data</th><th colspan="2">Budbreak Data</th></tr><tr><th rowspan="2">Procedure</th><th rowspan="2">Variety</th><th rowspan="2">Country</th><th colspan="2">Weather Station</th><th colspan="2">Period</th><th colspan="2">Vineyard</th><th rowspan="2">Cases<br>Number</th></tr><tr><th>Lat</th><th>Long</th><th>Years</th><th>Lat</th><th>Long</th><th>Years</th></tr></thead><tbody><tr><td rowspan="24">Fitting & Validation</td><td rowspan="6">Cabernet Sauvignon<br>Chardonnay<br>Gewürztraminer<br>Grenache<br>Pinot Gris<br>Riesling</td><td rowspan="6">France<br>Luxembourg<br>France<br>Luxembourg<br>Luxembourg<br>France<br>Luxembourg</td><td>43.58</td><td>3.96</td><td>1950-2017</td><td>43.33</td><td>3.56</td><td>1951-2012</td><td>39</td></tr><tr><td>49.50</td><td>6.35</td><td>2016-2018</td><td>49.50</td><td>6.35</td><td>2017-2018</td><td>2</td></tr><tr><td>43.58</td><td>3.96</td><td>1950-2017</td><td>43.33</td><td>3.56</td><td>1951-2012</td><td>33</td></tr><tr><td>49.50</td><td>6.35</td><td>2010-2018</td><td>49.50</td><td>6.35</td><td>2011-2018</td><td>8</td></tr><tr><td>49.50</td><td>6.35</td><td>1970-2018</td><td>49.50</td><td>6.35</td><td>1972-2018</td><td>50</td></tr><tr><td>43.58</td><td>3.96</td><td>1950-2017</td><td>43.33</td><td>3.56</td><td>1951-2012</td><td>41</td></tr><tr><td rowspan="6">Touriga Franca</td><td rowspan="6">Portugal</td><td>49.50</td><td>6.35</td><td>1970-2018</td><td>49.50</td><td>6.35</td><td>1971-2018</td><td>51</td></tr><tr><td>49.50</td><td>6.35</td><td>1970-2018</td><td>49.50</td><td>6.35</td><td>1971-2018</td><td>50</td></tr><tr><td>41.81</td><td>-8.41</td><td>2004-2009</td><td>41.81</td><td>-8.41</td><td>2005-2009</td><td>5</td></tr><tr><td>41.25</td><td>-7.11</td><td>2013-2018</td><td>41.25</td><td>-7.11</td><td>2014-2018</td><td>5</td></tr><tr><td>41.19</td><td>-7.54</td><td>2013-2018</td><td>41.19</td><td>-7.54</td><td>2014-2018</td><td>5</td></tr><tr><td>41.04</td><td>-7.04</td><td>2013-2017</td><td>41.04</td><td>-7.04</td><td>2014-2017</td><td>4</td></tr><tr><td rowspan="6">Touriga Nacional</td><td rowspan="6">Luxembourg</td><td>41.17</td><td>-7.56</td><td>2013-2017</td><td>41.17</td><td>-7.56</td><td>2014-2017</td><td>4</td></tr><tr><td>41.15</td><td>-7.62</td><td>2013-2018</td><td>41.15</td><td>-7.62</td><td>2014-2018</td><td>5</td></tr><tr><td>41.17</td><td>-7.55</td><td>2014-2018</td><td>41.17</td><td>-7.55</td><td>2015-2018</td><td>4</td></tr><tr><td>49.50</td><td>6.35</td><td>2016-2018</td><td>49.50</td><td>6.35</td><td>2017-2018</td><td>2</td></tr><tr><td>39.04</td><td>-9.18</td><td>1989-2014</td><td>39.04</td><td>-9.18</td><td>1990-2014</td><td>19</td></tr><tr><td>41.81</td><td>-8.41</td><td>2004-2009</td><td>41.81</td><td>-8.41</td><td>2005-2009</td><td>5</td></tr><tr><td rowspan="6">Cabernet Sauvignon<br>Chardonnay<br>Gewürztraminer<br>Grenache<br>Pinot Gris<br>Riesling<br>Touriga Franca<br>Touriga Nacional</td><td rowspan="6">Portugal</td><td>41.21</td><td>-7.43</td><td>2013-2018</td><td>41.21</td><td>-7.43</td><td>2014-2018</td><td>5</td></tr><tr><td>41.24</td><td>-7.76</td><td>2014-2017</td><td>41.24</td><td>-7.76</td><td>2015-2017</td><td>3</td></tr><tr><td>41.22</td><td>-7.54</td><td>2013-2018</td><td>41.22</td><td>-7.54</td><td>2014-2018</td><td>5</td></tr><tr><td>41.15</td><td>-7.76</td><td>2015-2018</td><td>41.15</td><td>-7.76</td><td>2016-2018</td><td>3</td></tr><tr><td>44.79</td><td>-0.58</td><td>2012-2019</td><td>44.79</td><td>-0.58</td><td>2013-2019</td><td>7</td></tr><tr><td>48.55</td><td>7.64</td><td>1950-2017</td><td>48.22</td><td>7.35</td><td>1976-1990</td><td>15</td></tr><tr><td rowspan="6">Application</td><td rowspan="6">France<br>Luxembourg<br>France<br>Luxembourg<br>France<br>Luxembourg</td><td>48.55</td><td>7.64</td><td>1950-2017</td><td>48.22</td><td>7.35</td><td>2006-2009</td><td>4</td></tr><tr><td>48.55</td><td>7.64</td><td>1950-2017</td><td>48.22</td><td>7.35</td><td>1976-1990</td><td>13</td></tr><tr><td>48.55</td><td>7.64</td><td>1950-2017</td><td>48.22</td><td>7.35</td><td>1976-1990</td><td>15</td></tr><tr><td>48.55</td><td>7.64</td><td>1950-2017</td><td>48.22</td><td>7.35</td><td>2006-2009</td><td>4</td></tr><tr><td>44.79</td><td>-0.58</td><td>2012-2019</td><td>44.79</td><td>-0.58</td><td>2013-2019</td><td>7</td></tr><tr><td>44.79</td><td>-0.58</td><td>2012-2019</td><td>44.79</td><td>-0.58</td><td>2013-2019</td><td>7</td></tr></tbody></table>

The phenological data were collected for eight grapevine varieties (Cabernet Sauvignon, Chardonnay, Gewürztraminer, Grenache, Pinot Gris, Riesling, Touriga Franca and Touriga Nacional). For France, these data were retrieved from PERPHECLIM platform (https://www6.inrae.fr/projet-accaf-perpheclim_eng/The-PERPHECLIM-project [36]) and VitAdapt project [35]. For Luxembourg, the data were obtained from the Administration des Services Techniques de l'Agriculture at the Institut Viti-Vinicole. Finally, the data for the Lisbon wine region (Portugal) were obtained from the national ampelographic collection at the INIAV, located in Dois Portos. The data for the Douro wine region were collected by the associates from ADVID. Most of these grapevine varieties are for high-quality wine production derived from the world's benchmark regions and can be grouped on the basis of average growing season temperatures and ripening potential [37]. The phenological observations were collected by different research institutions and farmers at budbreak date which typically corresponded to BBCH 09 [38], except for Cabernet Sauvignon, Touriga Franca and Touriga Nacional collected in France (Bordeaux site; BBCH 07). We assumed that the distance between BBCH 07 and 09 is covered in 2 days as reported in [27].

## 2.2. *Phenological Models*

In total, we tested 6 models for estimating the budbreak date, divided into forcing (F) and chilling-forcing (CF) models. We used Growing Degree Days (GDD; [39]), WANG [40] and UNIFORC [24] for F models and BRIN [29] (Hourly and Daily approach), UNICHILL and UNIFIED [24] for CF models.
### 2.2.1. GDD

The GDD model [39] is an F model that describes the eco-dormancy period based on accumulation of the daily average temperature ($T_{avg}$; °C) above a specific threshold ($T_b$; °C). The GDD accumulation ($Forc_{GDD}$) starts at a fixed date of the year ($t_0$) and is accumulated until the GDD requirement is satisfied at budbreak date (Bud; Equation (1)).

$$Forc_{GDD} = \sum_{t_0}^{Bud} (T_{avg} - T_b), \qquad (1)$$

### 2.2.2. WANG

The WANG model [40] is a temperature response function used, in this study, to describe the eco-dormancy period from a fixed date of the year ($t_0$) and not to calculate the chilling unit accumulation as reported for other similar functions [30]. The temperature function ($Forc_{WANG}$) allows the phenological development of grapevine to be evaluated using a non-linear approach based on three cardinal temperatures: $T_{opt}$ (optimum temperature, °C), $T_{max}$ (maximum temperature, °C) and $T_{min}$ (minimum temperature, °C; Equations (2a) and (2b)).

$$Forc_{WANG} = \sum_{t_0}^{Bud} Max\left[\frac{2 \cdot (T_{avg} - T_{min})^{\alpha} \cdot (T_{opt} - T_{min})^{\alpha} - (T_{avg} - T_{min})^{2\alpha}}{(T_{opt} - T_{min})^{2\alpha}}\right], 0, \qquad (2a)$$

$$\text{with } \alpha = \frac{\ln(2)}{\ln\left(\frac{T_{\text{max}} - T_{\text{min}}}{T_{\text{opt}} - T_{\text{min}}}\right)}, \qquad (2b)$$

where $T_{avg}$ is the daily average temperature (°C) and Bud is the budbreak date.

### 2.2.3. UNIFORC

The UNIFORC model [24] is an F model based on a sigmoid curve for calculating the forcing unit accumulation ($Forc_{UNIFORC}$) from a fixed day of the year ($t_0$) until budbreak occurrence (Bud; Equation (3)):

$$Forc_{UNIFORC} = \sum_{t_0}^{Bud} \frac{1}{1 + e^{d \cdot (T_{avg} - e)}}, \qquad (3)$$

where $T_{avg}$ is daily average temperature (°C), while $d$ and $e$ ($d < 0$ and $e > 0$) are empirical parameters of the function.

### 2.2.4. BRIN

The BRIN model [29] is a CF model based on the Bidabe's Cold Action Model [41,42] in which the endo-dormancy period is evaluated through the accumulation of chilling units ($Chill_{BRIN}$; Equation (4)). The eco-dormancy period is based on forcing unit accumulation ($Forc_{BRIN}$) calculated with the Richardson equation [23]. The $Forc_{BRIN}$ is computed using both a daily (Equation (5)) and hourly (Equations (6a) and (6b)) temperature accumulation approach.

$$Chill_{BRIN} = \sum_{t_0}^{C_{endo}} (Q_{10c} \frac{-T_{x(n)}}{10} + Q_{10c} \frac{-T_{n(n)}}{10}), \qquad (4)$$

$$Forc_{BRINDaily} = \sum_{C_{endo}}^{Bud} Richardson_{Daily}(n), \qquad (5)$$
$$Forc_{BRINHourly} = \sum_{C_{endo}}^{Bud} Richardson_{Hourly}(n) \qquad (6a)$$

$$Richardson = Max(Min(T_{avg} - T_{l}, T_{h} - T_{l}), 0), \qquad (6b)$$

where $t_0$ is the simulation starting date, $T_{avg}$ is daily average temperature (°C), $C_{endo}$ is the date at which the chilling requirement is satisfied, $Bud$ is budbreak date, $Q_{10c}$ is the geometric progression rate of the thermal dormancy response, $T_n$ and $T_x$ are daily minimum and maximum temperatures (°C), respectively. $Richardson_{Daily}$ and $Richardson_{Hourly}$ are daily and hourly Richardson equations used for representing the eco-dormancy period [29]. The Richardson equations are limited by two cardinal temperatures: a lower ($T_l$, °C) and an upper limit temperature ($T_h$, °C).

### 2.2.5. UNICHILL

The UNICHILL model [24] is a CF two-phase sequential model constituted by two functions describing chilling accumulation ($Chill_{UNICHILL}$) during endo-dormancy (Equation (7)) and forcing unit accumulation ($Forc_{UNICHILL}$) during eco-dormancy period (Equation (8)). After chilling requirement is satisfied at a specific date ($C_{endo}$), the forcing units start accumulating until the forcing requirement is reached, leading to the occurrence of the budbreak phase ($Bud$). The two functions are defined as followed:

$$Chill_{UNICHILL} = \sum_{t_0}^{C_{endo}} \frac{1}{1 + e^{a(T_{avg} - c) + b(T_{avg} - c)}}, \qquad (7)$$

$$Forc_{UNICHILL} = \sum_{C_{endo}}^{Bud} \frac{1}{1 + e^{d(T_{avg} - e)}}, \qquad (8)$$

where $t_0$ is the simulation starting date, $T_{avg}$ is daily average temperature (°C) and $a, b, c, d$ and $e$ are empirical parameters of the function.

### 2.2.6. UNIFIED

The UNIFIED model [24] is a CF overlapping model based on the same equations for chilling and forcing unit accumulation used in the UNICHILL model. However, UNIFIED model includes the effect of chilling units on forcing accumulation through implementation of the exponential decrease curve (Equation (9)) in which $ForcReq_{UNIFIED}$ (Fcrit) is the critical amount of forcing and $Chill_{UNIFIED}$ is the total amount of chilling at day $t_d$. The budbreak date is estimated when both chilling and forcing requirements are satisfied:

$$ForcReq_{UNIFIED} = we^{-zChill_{UNIFIED}(t_d)}, \qquad (9)$$

where $w$ and $z$ are empirical parameters of the function.

## 2.3. Model Fitting and Validation

### 2.3.1. Outliers Removal

The model fitting, validation and application were performed after removing outliers from the observed dataset. This procedure consists in the identification of possible anomalies by comparing the observed and simulated results obtained by a preliminary models' fitting against the original datasets for each variety. The outliers were identified by using the Tukey method [43], based on the non-parametric measure of the statistical dispersion of a given distribution quantified using the interquartile range (IQR). For each model and starting date, the outliers represent the values in which the difference between the observed and simulated budbreak Days of Year (DOYs) are outside the whiskers of the distribution (below Q1 - 1.5 IQR or above Q3 + 1.5 IQR values). An outlier of a specific variety is removed from the dataset when it is observed in at least 4 combinations (models x starting
date). The percentage of outliers removed from the fitting, validation and application datasets is: 4.17% (Cabernet Sauvignon), 5.36% (Chardonnay), 1.85% (Gewurztraminer and Riesling), 3.70% (Grenache), 1.52% (Pinot Gris), 10% (Touriga Franca and Nacional). These values were generally observed after extreme weather conditions (e.g., frosts of 1956, 1963 in France) or due to other influences (e.g., agronomic practices or errors in data collection) which cannot be reproduced during model simulation.

### 2.3.2. Parameters' Range and Starting Date

The parameters' range of variation in the 6 phenological models was defined according to the literature [29,30,44–46]. For the GDD model, the fitted parameter was $T_b$ which ranges from 0 to 10 °C. In this study, 0 °C was considered the limit temperature for bud growth while 10 °C was the maximum base temperature adopted for describing the eco-dormancy period as found in other studies [19,22,29,44]. For the WANG model, this function being generally used for wheat simulation [47], the parameter range of the cardinal temperatures was established considering the values obtained by other similar functions (Fbeta) implemented for simulating budbreak date in grapevine [30]. Thus, $T_{min}$ and $T_{max}$ were set at 0 and 45 °C respectively, while $T_{opt}$ ranged in a ±30% of the value reported in Fila et al. [30]. For UNIFORC model, the *d* and *e* parameter values varied in a ±30% interval with respect to the average value reported for different grapevine varieties (Glera, Cabernet Sauvignon, Chardonnay and Merlot [46]). For both daily and hourly approaches of BRIN model, the $Q_{10}$ parameter range was defined in a ±30% interval of the value estimated for several grapevine varieties [29]. Moreover, the lower limit temperature ($T_l$) was fitted in the range 0–10 °C, as reported for GDD model while the upper limit temperature ($T_h$) for post dormancy was fixed at 25 °C as indicated in Garcia de Cortazar-Atauri et al. [29]. For the UNICHILL model, the parameters *a*, *b*, *c*, *d* and *e* were in a range of ±30% of the overall minimum and maximum value reported by Fila [46] for Glera, Cabernet Sauvignon, Chardonnay and Merlot varieties (field calibration). However, in the specific case of Cabernet Sauvignon and Chardonnay, the parameter range of *a*, *b*, *c*, *d* and *e* varied in ±30% interval with respect to the parameter values found for these varieties. Similarly, the range adopted for *a*, *b*, *c*, *d* and *e* parameters in UNIFIED model was identified by considering ±30% of the minimum and maximum values reported by Fila et al. [46] (field calibration). Based on the same study, the parameters *w* and *z* were fixed by considering the average values reported for field calibration of Sangiovese and Montepulciano.

Concerning the use of fixed starting dates $t_0$ during model simulation, we set the 1st of January [30] and 1st of March [27] of the current year for the F models and the 1st September [24] and 1st of August [29] of the previous year for the CF models.

### 2.3.3. Model Fitting, Validation and Application

The methodology used for model fitting, validation and application is displayed in Figure 2. After outliers removal (Section 2.3.1), the first step was to randomly split the observed data (budbreak DOY per year) into 10 random subsamples (RS) in which 60% of the observed data was used for model fitting whereas the remaining 40% was retained for validation. Once the 10 random subsamples were prepared, the model fitting was performed for any combination of varieties, models and starting dates using the Phenological Modeling Platform (PMP) software (version 5.5) [48] and setting the parameter ranges as described in Section 2.3.2. The best fitting for any combination was found by exploring the range of each parameter through the Metropolis-Hastings algorithm [49–51] of PMP. Each fitted model was validated on the remaining 40% of the data and then applied on the validation datasets of the other subsamples by using the R software environment (version 3.5.0) [52]. Thus, the parameter set of each model combination (model x variety x starting date), which provided the highest-performance (HP) during fitting and validation procedures, was selected. Moreover, the average parameter sets (AV) derived from all fitted subsamples was also tested on the entire distribution of validation datasets. Finally, the HP and AV parameter sets were applied to the independent datasets of other study sites to evaluate the behavior of different parameterizations under different environmental conditions.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Workflow of the methodology adopted in this study. Random subsamples = Random distribution of the phenological observed data in subsample; Fitting and Validation = Model fitting on the selected 60% of the observed data and model validation on the remaining 40%; Applications = Model application on all validation and independent datasets. Fit = model fitting; Val = model validation; HP = highest-performance models; AV = average models; RS(n) = random subsamples; ID = independent datasets.</figcaption></figure>

## 2.4. Statistical Analysis

The statistical indices used for evaluating the phenological model performances during fitting, validation and application were: R-squared (R², Equation (10), [53]); Akaike information criterion (AIC, Equation 11, [54]); Root Mean Squared Error (RMSE, Equation (12), [55]) and modeling efficiency coefficient (EF, Equation (13), [56]), each defined as follows:

$$R^2 = 1 - \frac{\sum (O_i - P_i)^2}{\sum (O_i - \bar{O})^2}, \qquad (10)$$

$$AIC = 2k - 2 \ln(\hat{L}), \qquad (11)$$

$$RMSE = \left[ \sum_{i=1}^{n} \frac{(P_i - O_i)^2}{n} \right]^{0.5}, \qquad (12)$$

$$EF = 1 - \frac{\sum_{i=1}^{n} (P_i - O_i)^2}{\sum_{i=1}^{n} (O_i - \bar{O})^2}, \qquad (13)$$

where $O_i$ is the observed value, $\bar{O}$ is the average of the observed values, $P_i$ is the predicted value, $n$ is the number of observations, $k$ is the number of parameters in the model, and $\hat{L}$ is the maximized value of the likelihood function of the predicted model.

Finally, the minimum, 10% quantile, median, mean, 90% quantile and maximum values of the statistics were extracted for describing the simulation results on all validation datasets.

# 3. Results

## 3.1. Overall Model Behavior in Fitting and Validation

The results of the phenological model intercomparison were described comparing the performances of the different modeling approaches (e.g., F vs. CF) and starting dates (1st January/March vs. 1st September/August) during model fitting and validation. Initially, the model fitting and validation results evidenced higher average performances for CF (Fit, R² = 0.74; AIC = 150.42; RMSE = 4.22, Val,
$R^2 = 0.61$; AIC = 111.64; RMSE = 5.79) compared to F models (Fit, $R^2 = 0.61$; AIC = 160.62; RMSE = 5.51; Val, $R^2 = 0.59$; AIC = 112.01; RMSE = 6.10). Among all models (Figure 3 and Table S1), the highest performances were achieved by UNIFIED model (Fit, $R^2 = 0.77$; AIC = 147.60; RMSE = 3.96; Val, $R^2 = 0.64$; AIC = 110.74; RMSE = 5.42), while the lowest were obtained by the WANG model (Fit, $R^2 = 0.61$; AIC = 161.15; RMSE = 5.60) and GDD model (Val, $R^2 = 0.58$; AIC = 112.59; RMSE = 6.25). Moreover, higher performances were found for those models in which the starting date was fixed at 1st September in CF models (Fit, $R^2 = 0.74$; AIC = 150.24; RMSE = 4.21; Val, $R^2 = 0.61$; AIC = 111.59; RMSE = 5.75) and at 1st January in F models (Fit, $R^2 = 0.66$; AIC = 159.20; RMSE = 5.53; Val, $R^2 = 0.63$; AIC = 111.10; RMSE = 6.15). Finally, concerning grapevine varieties, Chardonnay was best fitted by the models (Fit, $R^2 = 0.90$; AIC = 138.66; RMSE = 4.90; Val, $R^2 = 0.89$; AIC = 101.54; RMSE = 5.97), while the poorest results were obtained for Touriga Nacional (Fit, $R^2 = 0.42$; AIC = 147.09; RMSE = 6.23; Val, $R^2 = 0.31$; AIC = 105.40; RMSE = 7.74) and Touriga Franca in terms of AIC values (Fit = 180.19; Val = 127.35).

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Statistical performances ($R^2$, AIC, RMSE, EF) of all combinations of the six phenological models, eight grapevine varieties and different starting dates during fitting and validation. The results of forcing (F) and chilling-forcing (CF) models (different bar colours) are displayed at the top and bottom of the figure, respectively. The two bars identify the model fitting (no bar texture) and validation (line bar texture). CS = Cabernet Sauvignon, CH = Chardonnay, GE = Gewürztraminer, GR = Grenache, PG = Pinot Gris, RS = Riesling, TF = Touriga Franca, TN = Touriga Nacional.</figcaption></figure>
3.2. *HP versus AV*

In order to compare the results of different models’ parameterizations, we selected the highest
performance parameter set (HP, Table 2) from the 10 fitting and validation subsamples for each
combination of model, starting date and variety. We then selected the average parameter set (AV, Table 3)
from the 10 fitting subsamples for each combination model, starting date and variety. On these bases,
Figures S1–S7 show the performances (R² and AIC) of the model fitting and validation (symbols) and
the applications of the different parameterizations on all validation subsamples (boxplot) obtained by
considering all grapevine varieties and starting dates. In detail, the selected parameter sets with highest
performances (HP; filled symbols with different shapes according to the starting dates) were the result of
the highest rank value of R² and AIC and the minimum discrepancy between fitting (red) and validation
(blue) models. By contrast, the AV parameter sets (boxplots in shadow area) were the application of the
average parameterization obtained for all subsamples. Despite that similar performances for HP and
AV parameter sets were evidenced in Figures S1–S7, higher discrepancy was found in the behavior of
CF and F models. In Figures S1–S4 (CF models), no differences were observed when different starting
dates were set during model fitting and application. On the other hand, high discrepancy was found
in model performances when different starting dates were adopted in Figures S5–S7 (F models). This
trend is also evidenced by the correlations between observed and simulated DOYs in Figure 4, Figures
S8 and S9. The results of model fitting showed accurate performances for all selected HP, with more
discrepancy between starting dates of F models (yellow and green correlations) compared to those used
in CF models (blue and red correlations). The values of R² also highlighted the higher performances
obtained with some varieties compared to others (Chardonnay vs. Touriga Nacional).

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Correlations between observed and simulated budbreak Days of Year (DOY) of the GDD and UNIFIED models. The results were obtained for all grapevine varieties and starting dates using the HP parameter set.</figcaption></figure>
<table><caption>**Table 2.** HP (highest-performance) parameter sets for each combination of models, grapevine varieties and starting dates. $C_{crit}$ = critical amount of chilling and $F_{crit}$ = critical amount of forcing.</caption><thead><tr><th rowspan="2">Variety</th><th rowspan="2">t<sub>0</sub></th><th colspan="3">GDD</th><th colspan="3">WANG</th><th colspan="8">UNIFORC/UNICHILL/UNIFIED</th><th colspan="4">BRIN h/BRIN d</th></tr><tr><th>F<sub>crit</sub></th><th>T<sub>b</sub></th><th>F<sub>crit</sub></th><th>T<sub>opt</sub></th><th>T<sub>min</sub></th><th>T<sub>max</sub></th><th>C<sub>crit</sub></th><th>F<sub>crit</sub></th><th>a</th><th>b</th><th>c</th><th>d</th><th>e</th><th>w</th><th>z</th><th>C<sub>crit</sub></th><th>F<sub>crit</sub></th><th>Q<sub>10c</sub></th><th>T<sub>1</sub></th><th>T<sub>h</sub></th></tr></thead><tbody><tr><td rowspan="3">Cabernet Sauvignon</td><td>1</td><td rowspan="3">763.47</td><td rowspan="3">0</td><td rowspan="3">67.32</td><td rowspan="3">18.15</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>16.39</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.47</td><td>-</td><td>-</td><td>190.14</td><td>6636.2</td><td>1.54</td><td>5.07</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>142.36</td><td>26.54</td><td>0.15</td><td>-23.61</td><td>-20.41</td><td>-0.19</td><td>11.80</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>184.86</td><td>295.3</td><td>1.54</td><td>4.45</td></tr><tr><td></td><td>151.14</td><td>20.62</td><td>1.31</td><td>-31.06</td><td>-1.45</td><td>-0.20</td><td>12.44</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Chardonnay</td><td>60</td><td rowspan="3">286.89</td><td rowspan="3">1.18</td><td rowspan="3">19.49</td><td rowspan="3">23.96</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>7.74</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.45</td><td>-</td><td>-</td><td>217.53</td><td>6597.1</td><td>1.52</td><td>5.07</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>178.96</td><td>23.19</td><td>0.16</td><td>-15.33</td><td>-24.24</td><td>-0.25</td><td>11.83</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>206.55</td><td>272.1</td><td>1.59</td><td>4.84</td></tr><tr><td></td><td>175.24</td><td>20.40</td><td>0.91</td><td>-26.60</td><td>-1.63</td><td>-0.14</td><td>14.69</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Grenache</td><td>1</td><td rowspan="3">432.23</td><td rowspan="3">2.34</td><td rowspan="3">14.63</td><td rowspan="3">30.42</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>8.58</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>18.46</td><td>-</td><td>-</td><td>148.48</td><td>5194.0</td><td>1.54</td><td>5.86</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>132.02</td><td>17.55</td><td>0.72</td><td>-29.03</td><td>-17.19</td><td>-0.24</td><td>12.56</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>148.34</td><td>238.3</td><td>1.54</td><td>4.99</td></tr><tr><td></td><td>123.96</td><td>19.71</td><td>0.85</td><td>-32.05</td><td>-2.20</td><td>-0.28</td><td>11.79</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Gewürrtztraminer</td><td>60</td><td rowspan="3">49.18</td><td rowspan="3">7.21</td><td rowspan="3">4.00</td><td rowspan="3">31.26</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>1.33</td><td>-</td><td>-</td><td>-</td><td>-0.30</td><td>18.36</td><td>-</td><td>-</td><td>173.38</td><td>5276.9</td><td>1.54</td><td>5.81</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>165.85</td><td>19.05</td><td>0.46</td><td>-19.82</td><td>-13.84</td><td>-0.24</td><td>11.93</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>172.90</td><td>211.9</td><td>1.54</td><td>5.44</td></tr><tr><td></td><td>172.71</td><td>19.63</td><td>1.19</td><td>-34.66</td><td>0.22</td><td>-0.25</td><td>10.46</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Pinot Gris</td><td>1</td><td rowspan="3">645.68</td><td rowspan="3">0</td><td rowspan="3">47.27</td><td rowspan="3">21.91</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>13.40</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.48</td><td>-</td><td>-</td><td>177.23</td><td>3769.3</td><td>1.67</td><td>6.63</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>156.67</td><td>8.05</td><td>0.61</td><td>-20.19</td><td>-11.39</td><td>-0.19</td><td>16.94</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>158.40</td><td>312.9</td><td>1.55</td><td>4.02</td></tr><tr><td></td><td>157.54</td><td>20.94</td><td>0.55</td><td>-30.18</td><td>-10.31</td><td>-0.26</td><td>9.98</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Riesling</td><td>60</td><td rowspan="3">85.11</td><td rowspan="3">5.99</td><td rowspan="3">5.17</td><td rowspan="3">30.40</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>4.17</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.49</td><td>-</td><td>-</td><td>215.25</td><td>5029.4</td><td>1.53</td><td>5.27</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>186.21</td><td>4.73</td><td>0.79</td><td>-26.99</td><td>-7.35</td><td>-0.18</td><td>21.33</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>210.27</td><td>319.2</td><td>1.54</td><td>2.80</td></tr><tr><td></td><td>186.89</td><td>20.41</td><td>1.04</td><td>-30.11</td><td>-1.81</td><td>-0.28</td><td>9.72</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Franca</td><td>1</td><td rowspan="3">137.17</td><td rowspan="3">7.10</td><td rowspan="3">12.11</td><td rowspan="3">31.26</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>9.37</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>18.50</td><td>-</td><td>-</td><td>283.90</td><td>3961.1</td><td>1.53</td><td>6.19</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>186.01</td><td>17.05</td><td>0.17</td><td>-32.51</td><td>-12.74</td><td>-0.35</td><td>10.79</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>280.72</td><td>181.3</td><td>1.54</td><td>5.39</td></tr><tr><td></td><td>172.62</td><td>20.42</td><td>1.18</td><td>-34.32</td><td>-10.15</td><td>-0.23</td><td>11.48</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">255.89</td><td rowspan="3">3.61</td><td rowspan="3">9.17</td><td rowspan="3">31.26</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>8.43</td><td>-</td><td>-</td><td>-</td><td>-0.23</td><td>16.45</td><td>-</td><td>-</td><td>312.24</td><td>4159.1</td><td>1.53</td><td>5.98</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>161.74</td><td>12.51</td><td>1.62</td><td>-34.01</td><td>-6.44</td><td>-0.33</td><td>12.71</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>311.34</td><td>182.5</td><td>1.54</td><td>5.36</td></tr><tr><td></td><td>214.68</td><td>20.29</td><td>0.50</td><td>-31.89</td><td>-10.95</td><td>-0.28</td><td>10.28</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Riesling</td><td>1</td><td rowspan="3">122.16</td><td rowspan="3">7.39</td><td rowspan="3">12.15</td><td rowspan="3">31.26</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>9.63</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>18.49</td><td>-</td><td>-</td><td>285.01</td><td>4633.1</td><td>1.53</td><td>5.22</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>183.87</td><td>14.70</td><td>0.56</td><td>-23.01</td><td>-19.53</td><td>-0.25</td><td>12.67</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>284.32</td><td>225.9</td><td>1.55</td><td>4.11</td></tr><tr><td></td><td>183.25</td><td>20.45</td><td>0.95</td><td>-30.94</td><td>-10.68</td><td>-0.28</td><td>10.26</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Franca</td><td>60</td><td rowspan="3">250.48</td><td rowspan="3">3.91</td><td rowspan="3">9.36</td><td rowspan="3">31.23</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>9.44</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.54</td><td>-</td><td>-</td><td>314.40</td><td>6374.7</td><td>1.53</td><td>3.52</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>174.41</td><td>12.60</td><td>1.13</td><td>-25.39</td><td>-8.01</td><td>-0.33</td><td>12.02</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>314.75</td><td>218.9</td><td>1.54</td><td>4.29</td></tr><tr><td></td><td>211.57</td><td>20.26</td><td>0.36</td><td>-11.66</td><td>-7.14</td><td>-0.28</td><td>10.31</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>1</td><td rowspan="3">156.14</td><td rowspan="3">6.74</td><td rowspan="3">12.60</td><td rowspan="3">31.26</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>9.90</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>18.46</td><td>-</td><td>-</td><td>267.03</td><td>4980.8</td><td>1.53</td><td>5.44</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>150.99</td><td>5.33</td><td>1.39</td><td>-27.15</td><td>-5.31</td><td>-0.18</td><td>21.68</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>280.95</td><td>208.2</td><td>1.56</td><td>4.93</td></tr><tr><td></td><td>173.34</td><td>20.47</td><td>1.09</td><td>-28.42</td><td>-10.79</td><td>-0.21</td><td>11.23</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">279.34</td><td rowspan="3">3.46</td><td rowspan="3">11.80</td><td rowspan="3">30.29</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>10.14</td><td>-</td><td>-</td><td>-</td><td>-0.22</td><td>16.45</td><td>-</td><td>-</td><td>297.05</td><td>5022.2</td><td>1.53</td><td>5.41</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>163.61</td><td>20.58</td><td>1.00</td><td>-19.62</td><td>-5.65</td><td>-0.26</td><td>10.30</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>313.12</td><td>221.8</td><td>1.53</td><td>4.68</td></tr><tr><td></td><td>16.61</td><td>24.35</td><td>0.83</td><td>-21.47</td><td>-1.62</td><td>-0.21</td><td>16.45</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Franca</td><td>1</td><td rowspan="3">743.93</td><td rowspan="3">0</td><td rowspan="3">36.53</td><td rowspan="3">26.70</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>142.03</td><td>24.35</td><td>0.42</td><td>-14.45</td><td>-10.02</td><td>-0.27</td><td>12.24</td><td>-0.00025</td><td>21.05</td><td>-</td><td>158.30</td><td>7933.9</td><td>1.53</td><td>4.80</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>141.22</td><td>20.61</td><td>0.42</td><td>-14.45</td><td>-10.02</td><td>-0.27</td><td>12.24</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>144.69</td><td>493.9</td><td>1.52</td><td>2.89</td></tr><tr><td></td><td>4.84</td><td>24.24</td><td>0.33</td><td>-18.24</td><td>-25.30</td><td>-0.23</td><td>11.70</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>180.41</td><td>6998.5</td><td>1.52</td><td>5.79</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">149.08</td><td rowspan="3">2.86</td><td rowspan="3">7.61</td><td rowspan="3">28.57</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>171.91</td><td>20.53</td><td>0.38</td><td>-13.35</td><td>-5.10</td><td>-0.19</td><td>14.67</td><td>-0.00025</td><td>21.05</td><td>-</td><td>183.23</td><td>293.1</td><td>1.52</td><td>5.28</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>159.86</td><td>20.53</td><td>0.38</td><td>-13.35</td><td>-5.10</td><td>-0.21</td><td>15.44</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td></td><td>16.11</td><td>21.11</td><td>0.19</td><td>-22.31</td><td>-18.25</td><td>-0.21</td><td>15.44</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>1</td><td rowspan="3">688.00</td><td rowspan="3">1.06</td><td rowspan="3">35.57</td><td rowspan="3">26.38</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>117.39</td><td>19.76</td><td>0.47</td><td>-27.59</td><td>-0.002</td><td>-0.14</td><td>17.19</td><td>-0.00025</td><td>21.05</td><td>-</td><td>133.81</td><td>13,660.2</td><td>1.59</td><td>2.15</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>124.23</td><td>5.20</td><td>0.47</td><td>-27.59</td><td>-0.002</td><td>-0.14</td><td>17.19</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>123.57</td><td>774.2</td><td>1.54</td><td>0.92</td></tr><tr><td></td><td>150.90</td><td>20.54</td><td>0.26</td><td>-8.36</td><td>-0.72</td><td>-0.14</td><td>17.87</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">199.94</td><td rowspan="3">0.80</td><td rowspan="3">8.19</td><td rowspan="3">28.71</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>148.45</td><td>20.97</td><td>0.61</td><td>-22.96</td><td>-6.80</td><td>-0.21</td><td>16.45</td><td>-</td><td>-</td><td>165.37</td><td>12,985.3</td><td>1.57</td><td>2.40</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>150.90</td><td>20.54</td><td>0.26</td><td>-8.36</td><td>-0.72</td><td>-0.14</td><td>17.87</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>
<table><caption>**Table 3.** AV (average) parameter sets for each combination of models, grapevine varieties and starting dates. $C_{crit}$ = critical amount of chilling and $F_{crit}$ = critical amount of forcing.</caption><thead><tr><th rowspan="2">Variety</th><th rowspan="2">t<sub>0</sub></th><th colspan="3">GDD</th><th colspan="3">WANG</th><th colspan="7">UNIFORC/UNICHILL/UNIFIED</th><th colspan="4">BRIN h/BRIN d</th></tr><tr><th>F<sub>crit</sub></th><th>T<sub>b</sub></th><th>F<sub>crit</sub></th><th>T<sub>opt</sub></th><th>T<sub>min</sub></th><th>T<sub>max</sub></th><th>C<sub>crit</sub></th><th>F<sub>crit</sub></th><th>a</th><th>b</th><th>c</th><th>d</th><th>e</th><th>w</th><th>z</th><th>C<sub>crit</sub></th><th>F<sub>crit</sub></th><th>Q<sub>10c</sub></th><th>T<sub>1</sub></th><th>T<sub>h</sub></th></tr></thead><tbody><tr><td rowspan="3">Cabernet Sauvignon</td><td>1</td><td rowspan="3">765.16</td><td rowspan="3">0</td><td rowspan="3">68.54</td><td rowspan="3">17.99</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>16.42</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.45</td><td>-</td><td>-</td><td>182.31</td><td>11,075.3</td><td>1.54</td><td>2.07</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>24.79</td><td>0.18</td><td>-20.62</td><td>-20.10</td><td>-0.19</td><td>11.91</td><td>-</td><td>-</td><td>177.95</td><td>508.7</td><td>1.55</td><td>1.49</td></tr><tr><td>-152</td><td>152.62</td><td>0.94</td><td>-23.60</td><td>-2.86</td><td>-0.22</td><td>12.08</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Chardonnay</td><td>60</td><td rowspan="3">225.25</td><td rowspan="3">3.16</td><td rowspan="3">11.83</td><td rowspan="3">28.64</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>7.31</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.46</td><td>-</td><td>-</td><td>196.01</td><td>11,942.8</td><td>1.56</td><td>1.91</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>176.29</td><td>0.18</td><td>-21.03</td><td>-18.95</td><td>-0.19</td><td>11.89</td><td>-</td><td>-</td><td>193.98</td><td>496.5</td><td>1.61</td><td>1.80</td></tr><tr><td>-121</td><td>180.78</td><td>0.87</td><td>-27.91</td><td>-4.13</td><td>-0.19</td><td>12.69</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Grenache</td><td>1</td><td rowspan="3">527.48</td><td rowspan="3">0.91</td><td rowspan="3">22.21</td><td rowspan="3">28.11</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>11.02</td><td>-</td><td>-</td><td>-</td><td>-0.22</td><td>16.91</td><td>-</td><td>-</td><td>150.00</td><td>6316.0</td><td>1.61</td><td>4.97</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>127.15</td><td>14.36</td><td>0.57</td><td>-22.95</td><td>-14.48</td><td>-0.25</td><td>14.42</td><td>-</td><td>-</td><td>167.46</td><td>254.1</td><td>1.54</td><td>4.63</td></tr><tr><td>-152</td><td>130.52</td><td>20.32</td><td>0.88</td><td>-24.70</td><td>-2.52</td><td>-0.27</td><td>11.40</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Gewürrtztraminer</td><td>60</td><td rowspan="3">47.34</td><td rowspan="3">7.24</td><td rowspan="3">4.10</td><td rowspan="3">21.26</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>1.50</td><td>-</td><td>-</td><td>-</td><td>-0.30</td><td>18.30</td><td>-</td><td>-</td><td>171.33</td><td>5931.7</td><td>1.61</td><td>5.35</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>160.41</td><td>17.18</td><td>0.65</td><td>-27.13</td><td>-13.07</td><td>-0.25</td><td>13.11</td><td>-</td><td>-</td><td>184.56</td><td>266.4</td><td>1.55</td><td>4.44</td></tr><tr><td>-121</td><td>163.81</td><td>20.45</td><td>0.85</td><td>-28.65</td><td>-3.25</td><td>-0.26</td><td>11.35</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Pinot Gris</td><td>1</td><td rowspan="3">653.93</td><td rowspan="3">0</td><td rowspan="3">54.66</td><td rowspan="3">19.74</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>13.61</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.50</td><td>-</td><td>-</td><td>184.09</td><td>6316.6</td><td>1.58</td><td>4.31</td><td rowspan="3">25</td></tr><tr><td>-121</td><td>126.06</td><td>8.07</td><td>1.04</td><td>-26.44</td><td>-8.24</td><td>-0.26</td><td>16.68</td><td>-</td><td>-</td><td>172.80</td><td>286.7</td><td>1.66</td><td>3.87</td></tr><tr><td>-152</td><td>151.30</td><td>20.59</td><td>0.78</td><td>-24.27</td><td>-5.82</td><td>-0.24</td><td>10.49</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Riesling</td><td>60</td><td rowspan="3">113.43</td><td rowspan="3">4.85</td><td rowspan="3">6.61</td><td rowspan="3">29.66</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>3.96</td><td>-</td><td>-</td><td>-</td><td>-0.23</td><td>17.10</td><td>-</td><td>-</td><td>182.40</td><td>8882.5</td><td>1.67</td><td>2.97</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>150.54</td><td>10.30</td><td>0.90</td><td>-28.01</td><td>-10.48</td><td>-0.24</td><td>16.05</td><td>-</td><td>-</td><td>198.81</td><td>324.5</td><td>1.60</td><td>3.11</td></tr><tr><td>-121</td><td>173.87</td><td>20.58</td><td>0.80</td><td>-25.81</td><td>-5.45</td><td>-0.24</td><td>11.04</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Franca</td><td>60</td><td rowspan="3">314.55</td><td rowspan="3">2.55</td><td rowspan="3">12.43</td><td rowspan="3">29.91</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>8.95</td><td>-</td><td>-</td><td>-</td><td>-0.23</td><td>18.17</td><td>-</td><td>-</td><td>251.39</td><td>5324.3</td><td>1.71</td><td>5.01</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>157.90</td><td>12.92</td><td>0.93</td><td>-26.85</td><td>-13.41</td><td>-0.25</td><td>14.65</td><td>-</td><td>-</td><td>264.22</td><td>228.9</td><td>1.60</td><td>4.53</td></tr><tr><td>-121</td><td>176.90</td><td>20.45</td><td>0.71</td><td>-25.17</td><td>-8.95</td><td>-0.26</td><td>10.95</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">416.44</td><td rowspan="3">0.93</td><td rowspan="3">16.95</td><td rowspan="3">28.35</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>7.61</td><td>-</td><td>-</td><td>-</td><td>-0.25</td><td>18.35</td><td>-</td><td>-</td><td>271.42</td><td>6580.6</td><td>1.56</td><td>3.75</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>169.54</td><td>14.04</td><td>0.89</td><td>-27.31</td><td>-13.36</td><td>-0.23</td><td>14.11</td><td>-</td><td>-</td><td>269.51</td><td>287.7</td><td>1.59</td><td>3.29</td></tr><tr><td>-121</td><td>173.20</td><td>20.38</td><td>0.69</td><td>-22.47</td><td>-9.57</td><td>-0.25</td><td>10.82</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Franca</td><td>60</td><td rowspan="3">359.67</td><td rowspan="3">1.77</td><td rowspan="3">15.09</td><td rowspan="3">28.93</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>9.48</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.47</td><td>-</td><td>-</td><td>292.53</td><td>6507.2</td><td>1.61</td><td>3.86</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>188.72</td><td>14.01</td><td>0.85</td><td>-25.77</td><td>-12.46</td><td>-0.24</td><td>14.19</td><td>-</td><td>-</td><td>299.62</td><td>290.8</td><td>1.57</td><td>3.24</td></tr><tr><td>-121</td><td>184.87</td><td>20.47</td><td>0.77</td><td>-21.17</td><td>-8.63</td><td>-0.27</td><td>10.65</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">143.12</td><td rowspan="3">3.71</td><td rowspan="3">7.49</td><td rowspan="3">28.97</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>8.51</td><td>-</td><td>-</td><td>-</td><td>-0.23</td><td>18.47</td><td>-</td><td>-</td><td>267.92</td><td>6755.6</td><td>1.65</td><td>3.75</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>167.21</td><td>13.93</td><td>0.94</td><td>-25.59</td><td>-11.82</td><td>-0.24</td><td>14.39</td><td>-</td><td>-</td><td>270.05</td><td>317.7</td><td>1.63</td><td>2.84</td></tr><tr><td>-121</td><td>166.51</td><td>20.56</td><td>0.95</td><td>-24.12</td><td>-8.63</td><td>-0.26</td><td>10.92</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">759.41</td><td rowspan="3">0.07</td><td rowspan="3">36.55</td><td rowspan="3">26.40</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>9.76</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.65</td><td>-</td><td>-</td><td>294.27</td><td>6967.0</td><td>1.66</td><td>3.58</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>176.49</td><td>13.62</td><td>0.99</td><td>-27.07</td><td>-12.29</td><td>-0.27</td><td>13.61</td><td>-</td><td>-</td><td>296.64</td><td>297.4</td><td>1.66</td><td>3.20</td></tr><tr><td>-121</td><td>174.36</td><td>20.30</td><td>1.03</td><td>-25.01</td><td>-7.57</td><td>-0.26</td><td>10.71</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Franca</td><td>60</td><td rowspan="3">710.02</td><td rowspan="3">0.77</td><td rowspan="3">37.39</td><td rowspan="3">26.29</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>15.97</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.87</td><td>-</td><td>-</td><td>154.21</td><td>8897.4</td><td>1.54</td><td>4.39</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>137.46</td><td>19.81</td><td>0.69</td><td>-22.91</td><td>-8.82</td><td>-0.24</td><td>14.33</td><td>-</td><td>-</td><td>149.90</td><td>363.8</td><td>1.57</td><td>4.42</td></tr><tr><td>-121</td><td>137.41</td><td>20.48</td><td>0.73</td><td>-23.67</td><td>-6.05</td><td>-0.23</td><td>13.26</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">108.60</td><td rowspan="3">0.77</td><td rowspan="3">37.39</td><td rowspan="3">26.29</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>4.96</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.46</td><td>-</td><td>-</td><td>178.92</td><td>8466.3</td><td>1.54</td><td>4.63</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>148.81</td><td>24.27</td><td>0.63</td><td>-25.76</td><td>-11.79</td><td>-0.22</td><td>14.28</td><td>-</td><td>-</td><td>178.72</td><td>387.1</td><td>1.56</td><td>4.20</td></tr><tr><td>-121</td><td>148.85</td><td>20.42</td><td>0.64</td><td>-22.47</td><td>-5.83</td><td>-0.21</td><td>13.90</td><td>-0.00025</td><td>21.05</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td rowspan="3">Touriga Nacional</td><td>60</td><td rowspan="3">176.18</td><td rowspan="3">1.89</td><td rowspan="3">9.42</td><td rowspan="3">26.87</td><td rowspan="3">0</td><td rowspan="3">45</td><td>-</td><td>4.86</td><td>-</td><td>-</td><td>-</td><td>-0.21</td><td>16.49</td><td>-</td><td>-</td><td>106.82</td><td>15,757.2</td><td>1.81</td><td>2.14</td><td rowspan="3">25</td></tr><tr><td>-152</td><td>123.73</td><td>20.27</td><td>0.62</td><td>-23.51</td><td>-12.11</td><td>-0.22</td><td>15.56</td><td>-</td><td>-</td><td>118.99</td><td>694.1</td><td>1.58</td><td>1.98</td></tr><tr><td>-121</td><td>123.73</td><td>20.27</td><td>0.62</td><td>-23.51</td><td>-12.11</td><td>-0.22</td><td>15.53</td><td>-</td><td>-</td><td>140.60</td><td>691.0</td><td>1.65</td><td>1.92</td></tr></tbody></table>
A deeper analysis of the model validation showed that the normalized difference between HP and AV, (HP-AV/HP)*100, was higher in CF ($R^2$ = 11.99%, AIC = 3.36%, RMSE = 22.20%) compared to F models ($R^2$ = 5.73%, AIC = 1.51%, RMSE = 8.16%). Among CF models, UNICHILL model showed the lowest difference between HP and AV in terms of $R^2$ and AIC ($R^2$ = 6.20%, AIC = 2.70%), while it evidenced the highest difference in terms of RMSE (40.02%). Considering F models, WANG showed the highest difference ($R^2$ = 7.07%, AIC = 2.25%, RMSE = 9.38%) and UNIFORC, the lowest ($R^2$ = 4.21%, AIC = 0.84%, RMSE = 7.13%). Following the same trend, the starting dates 1st September and 1st August ($R^2$ = 3.89%, AIC = 0.86%, RMSE = 15.62%) showed bigger differences compared to 1st January and 1st March ($R^2$ = 2.96%, AIC = 0.72%, RMSE = 7.62%). In particular, 1st September showed a lower difference in terms of AIC and RMSE (AIC = 0.86%, RMSE = 13.53%) compared to 1st August in CF models while 1st January showed the highest difference in F models ($R^2$ = 2.88%, AIC = 0.65%, RMSE = 7.23%). Finally, regarding grapevine varieties, the highest difference was found for Touriga Nacional ($R^2$ = 23.96%, RMSE = 24.20%) while the lowest was shown by Chardonnay in terms of $R^2$ (1.89%) and by Touriga Franca in terms of AIC and RMSE (1.10% and 9.55%, respectively).

### 3.3. *Model Applications on Other Independent Datasets*

The statistical performances of the phenological models, applied on the respective independent datasets, are reported in Figure 5 and Figures S10-S16 while the comparison between CF and F models during fitting and applications is shown in Figure 6 and Figure S17. In general, higher performances were found for CF compared to F models during fitting and applications (Figure 5), except for those varieties characterized by few observed data in the model application (Gewürztraminer, Riesling; Figure 6 and Figure S17).

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Highest-performance (HP) and average (AV) parameters' set applied on the independent datasets by considering all phenological models and starting dates. The example of Chardonnay variety.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Statistical results ($R^2$, AIC and RMSE) of the comparison between CF (black bars) and F (red bars) approaches in model fitting and application for all grapevine varieties using HP parameter set. The differences between CF and F models (CF-F) in both fitting and applications are displayed using green and yellow bars, respectively. CS = Cabernet Sauvignon, CH = Chardonnay, GE = Gewürztraminer, GR = Grenache, PG = Pinot Gris, RS = Riesling, TF = Touriga Franca, TN = Touriga Nacional.</figcaption></figure>

In particular, the model application showed similar performances for both HP ($R^2$ = 0.55, AIC = 57.97, RMSE = 9.41) and AV ($R^2$ = 0.54, AIC = 58.37, RMSE = 9.26) parameter sets. However, this difference slightly increased when CF ($R^2$ = 0.56, AIC = 59.05, RMSE = 8.75) and F models ($R^2$ = 0.52, AIC = 57.01, RMSE = 10.01) were compared without distinguishing between HP and AV parameterizations. Therefore, CF models showed better results in terms of $R^2$ (+8%) and RMSE (-15%) compared to F models while lower results of AIC value were obtained for CF compared to F models (+3.44%). A detailed analysis revealed that the CF models with highest performances ($R^2$ and AIC) were UNICHILL ($R^2$ = 0.64, AIC = 56.74) and UNIFIED ($R^2$ = 0.58, AIC = 56.61) models, while, in F models, the better results were provided by UNIFORC ($R^2$ = 0.54, AIC = 55.98). However, in terms of RMSE, highest performances were found for BRIN hourly (8.37, CF) and GDD (9.67, F). Concerning the starting dates, similar results were found when the 1st September ($R^2$ = 0.58, AIC = 58.92, RMSE = 8.55) and 1st of August ($R^2$ = 0.57, AIC = 59.17, RMSE = 8.95) were analyzed, with the 1st September showing better performance ($R^2$ = +1.83%, AIC = -0.44% and RMSE = -4.68%). Instead, higher differences were found when the 1st January ($R^2$ = 0.55, AIC = 60.47, RMSE = 7.74) and 1st of March ($R^2$ = 0.49, AIC = 52.57, RMSE = 12.47) were compared, evidencing the most relevant performances for 1st January ($R^2$ = +10.83% and RMSE = -61.23%, with the exception of AIC = +14.15%).

Moreover, the behavior of grapevine varieties in response to specific climate conditions was noticed during model application. The higher performances of CF compared to F, as well the relevant
differences between starting dates when F models were used, may be explained by the influence of temperatures on the dormancy period. The end of the first phase of dormancy (endo-dormancy) in CF models was generally simulated between the end of December (e.g., Touriga Nacional: DOY 356 for t₀ = -121 and -152, respectively) and end of February (e.g., Riesling: DOY 56 and 57 for t₀ = -121 and -152, respectively). More specifically, the period of endo-dormancy fulfilment was similarly reached regardless of starting date (1st September/August), allowing the next forcing unit accumulation during the eco-dormancy period to be correctly satisfied (differences in days from 0 to 3 days). The analysis further highlighted that the endo-dormancy release occurs later at northern than at southern latitudes and this is consistent with results obtained using F approach. These demonstrate that a late starting date (1st of March) for forcing unit accumulation is more suitable for budbreak simulation of varieties collected in Northern Europe (e.g., Gewürztraminer, Riesling and Pinot Gris), whose endo-dormancy release is simulated by CF late in the season (54 DOY on average). Conversely, an early starting date (1 January) provided the best results in simulating budbreak for grapevine varieties collected in Southern Europe (e.g., Touriga Franca and Nacional, Cabernet Sauvignon, Grenache and Chardonnay), which are characterized by an early end of endo-dormancy as simulated by CF models (13 DOY on average).

## 4. Discussion

This study compares the performance of six phenological models in simulating the budbreak date of eight grapevine varieties across different European sites. We found that the model performances characterized by a Chilling-Forcing approach were generally higher compared to those obtained for Forcing approach when different starting dates were set (Figure 6 and Figure S17).

Although model intercomparison studies on the estimation of budbreak date in grapevine have already been conducted by other authors [15,29,30], this study focused on evaluation of the models' reliability to reproduce budbreak in different environmental conditions, by considering the impact of temperature during the chilling and forcing accumulation period. In contrast to other works [27,30,57,58], this model intercomparison considers the use of a long-term series of observed phenology, which include grapevine varieties with different precocity levels collected at various latitudes in three European countries (France, Luxembourg and Portugal, Figure 1 and Table 1). The extensive observed dataset coupled to the methodology applied in this study allowed to evaluate the influence of study site-temperature variability on the grapevine development. Moreover, the cross-validation scheme adopted in this model intercomparison experiment provided an exhaustive analysis of the bias and variance of the CF and F models while the comparison between two different types of parameterization (HP or AV) pointed out the issue to select a specific parameter set for model applications (Results section, Figure 5). This in-depth investigation of model fitting, validation and application was particularly useful for evaluating the influence of temperature during the dormancy period. In this context, the different model responses during budbreak predictions are related to the physiological requirements needed for dormancy fulfilment of the grapevine varieties. Differently from the CF models, which describe the entire dormancy period (endo- and eco-dormancy), the F models only describe the eco-dormancy period through forcing units accumulation, under the assumption that the chilling requirement was already satisfied [24]. Although the debate on which model type provides the highest reliability is still open [29,30], our analysis shows that CF models, considering the influences of winter and spring temperatures on the accumulation of chilling and forcing units, may avoid the under- or over-estimation at the beginning of eco-dormancy, which occur when F models are set on fixed dates. In this study, we estimated a longer endo-dormancy for the vines from the Northern European sites (e.g., Gewürztraminer), which occurs later than the endo-dormancy of vines from the Southern European sites (e.g., Touriga Nacional), irrespective of the precocity level of the different varieties. In other terms, grapevine varieties such as Gewürztraminer (early budbreak) and Riesling (late budbreak) collected in Northern Europe showed a late endo-dormancy release (higher chilling requirement) with respect to Chardonnay (early budbreak) and Cabernet Sauvignon (late budbreak)
mainly collected in Southern Europe (lower chilling requirement). Although this trend was examined in this study only from a modeling point of view without any physiological evidence on the dormancy period, our results suggested that plants cultivated at different latitudes might show distinct dormancy period lengths. The higher accuracy of CF compared to F models is likely linked to the ability of these models to estimate the endo-dormancy release and to better account for the effect of temperature in this phase, avoiding the use of a fixed a priori value for the beginning of the eco-dormancy period as in F models [27,29]. Indeed, the goodness of F models is shown to significantly increase when the starting date is optimized against observed experimental data [30,45], demonstrating that the beginning of eco-dormancy period is dependent to the origin of the observed dataset and, thus, on the fulfilment of specific physiological requirements that change in relation to the climate conditions of the different environments. In fact, although no relevant differences were found for budbreak date estimation of grapevine varieties using CF models set on 1st of September/August starting dates, the discrepancy of the F model simulations was evident when 1st of January/March dates were used. Higher accuracy of F models to predict the budbreak date of grapevine varieties collected in Northern European countries was observed using the 1st of March starting date as also confirmed by [27], while the 1st January seems to be a more suitable starting date for representing budbreak date of the grapevine varieties from Southern Europe (Figures S5–S7). In this case, a compromise between both starting dates may be represented by the beginning of February [59]. However, the adoption of different starting dates in F models highlights the role of temperature on bud growth and the need to satisfy specific dormancy requirements for varieties cultivated at different latitudes before budbreak occurrence. These results were supported by the work of Reis-Pereira et al. [60] on the grapevine flowering stage. The authors pointed out that the selection of a priori value for the starting of thermal unit accumulation in F models may lead to an under- or over-estimation of the flowering phase, demonstrating the importance of a correct evaluation of the previous phenological stage in respect to the subsequent one. These findings were also observed in relation to the flowering stage of other species such as olive tree e.g., [61]). In this context, we infer that the dormancy period of grapevine may play a relevant role for predicting budbreak date especially under future scenarios. Indeed, more accurate estimations of budbreak date using CF models are provided when the endo-dormancy release is calibrated against observed data of perennial crops. This behavior is particularly evident under warmer temperatures (after 2050) where CF without endo-dormancy calibration and F models showed poor performances of phenology prediction [62]. The accurate identification of the endo-dormancy release avoids selecting a fixed a priori date, especially when this may have strong consequences on the estimation of budbreak date (as in F models).

The impact of warmer temperatures after the endo-dormancy release showed higher effect on the budbreak phase compared to the influence of temperatures during chilling accumulation, as reported by Martinez-Lüscher et al. [63] in a two-year ad hoc experiment with different grapevine varieties and temperature regimes. However, the starting date for thermal unit accumulation depends on the physiological features of the grapevine variety and their interaction with the environment. Indeed, we found that the effect of temperature during the last part of winter and beginning of spring was also evident when models were tested on external sites (independent datasets, Figure 5). Although most of the models showed satisfactory results during application (R² and AIC), the average error (RMSE) between observed and simulated data was relatively high (from 4 to 19 days). This situation was exacerbated in F models where the reliability in simulating budbreak is dependent to the use of a specific starting date (1st of January/1st of March) and thus the original observed dataset. The discrepancy between the two starting dates raises the question of which initial date to choose for forcing unit accumulation when simulations are conducted on a European scale [31,64].
The implications of this study are particularly relevant for growers especially when considering the impact of climate change. Our study, in line with a large part of literature, confirmed that chill units accumulations, the amount of which varies according to variety and location, is required for exiting from endo-dormancy phase. In a warmer climate, as expected in the future, this amount may not fully satisfy, leading to developmental issues during the growing season as outlined in the introduction of this paper. Some authors already highlighted that, on a global scale, wine areas may shift to match their specific climate niches in case of no adaptation [11,65], while new varietal selections and breeding programmes are needed for improving the adaptability of the grapevine varieties and maintaining the high quality production [12,59,66]. In this context, our parameterizations may be used as metric, in combination with climate simulations under future scenarios, to highlight which regions and varieties are expected to be most exposed to the risk of budbreak failure. Accordingly, the use of alternative varieties, better adapted to the expected future climate, may be proposed or the direction of specific breeding programmes can be indicated to develop varieties with a lower chill units requirement as requested in a warmer climate. Moreover, we should also take into account that, besides thermal requirements, dormancy phase may be influenced by the photoperiod [67]. Since each grapevine variety requires critical thresholds of day length for the onset of dormancy [68,69], this regulating factor should be considered when CF models are adopted for simulations at different latitudes. As such, comprehensive datasets of phenological observations and ad hoc experiments for evaluating the influence of more factors on dormancy [30,46,69] may be useful to take into account for improving model reliability under future climate scenarios and different environments [31,64,70].

**5. Conclusions**

This model intercomparison evidenced the different performance of the phenological models in reproducing the budbreak date of eight grapevine varieties collected in different European regions. The reliability of two main phenological modeling approaches (CF and F) was evaluated considering different starting dates (1st January/March or 1st September/August) and the selection methodologies of the parameter set (HP and AV) for model application. The results showed, in general, higher performances of CF compared to F models where fixed starting dates were adopted, even if neither modeling approach highlighted a clear supremacy over the other. For this reason, CF and F models' behavior should be further investigated through ad hoc experiments for evaluating the role of dormancy on budbreak occurrence. Indeed, the findings of our model simulations evidenced distinct thermal requirements for grapevine varieties cultivated at different latitudes (i.e., Northern Europe vs. Southern Europe). The use of different starting dates highlighted that, unlike CF models, the accurate estimation of budbreak in F models is affected by the selection of a starting date for beginning of the eco-dormancy period. This implies that, if this phase is not accurately defined, F models may lead to over- or under-estimations of the budbreak date under different environmental conditions. In this context, comprehensive phenological datasets (e.g., grapevine phenological data collected at different latitudes) are needed to thoroughly investigate the effect of temperature variability on dormancy period and budbreak estimation. Finally, it is proposed the development and implementation of improved phenological models, which consider the effect of other climate (e.g., radiation) and soil (e.g., soil temperature) variables as well as some regulating factors of the endo-dormancy phase (e.g., day length).
**Supplementary Materials:** The following are available online at http://www.mdpi.com/2076-3417/10/11/3800/s1.
Figure S1: R² and AIC values of the BRIN Daily model obtained for all fitted and average models, grapevine varieties and starting dates. The symbols (points and triangles) show the model results at two starting dates (in this example $t_0 = -121$ and $-152$ DOY), while the two colours represent fitting (red) and validation (blue) procedure, respectively. The filled symbols constitute the R² and AIC values of the models obtained by the selected HP parameter sets, while the dashed (fitting) and dotted (validation) lines depict the average R² and AIC values for both starting dates (in this example orange line = $-121$ DOY and green line = $-152$ DOY). The boxplots represent the R² and AIC values distribution for HP (orange and light blue: HP model application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively) and AV (boxplot in the shaded area, orange and light blue: AV application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively). The outliers represent all values out of the 1.5\*InterQuartile Range (IQR) of the boxplot. Figure S2: R² and AIC values of the BRIN Hourly model obtained for all fitted and average models, grapevine varieties and starting dates. The symbols (points and triangles) show the model results at two starting dates (in this example $t_0 = -121$ and $-152$ DOY), while the two colours represent fitting (red) and validation (blue) procedure, respectively. The filled symbols constitute the R² and AIC values of the models obtained by the selected HP parameter sets, while the dashed (fitting) and dotted (validation) lines depict the average R² and AIC values for both starting dates (in this example orange line = $-121$ DOY and green line = $-152$ DOY). The boxplots represent the R² and AIC values distribution for HP (orange and light blue: HP model application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively) and AV (boxplot in the shaded area, orange and light blue: AV application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively). The outliers represent all values out of the 1.5\*InterQuartile Range (IQR) of the boxplot. Figure S3: R² and AIC values of the UNICHILL model obtained for all fitted and average models, grapevine varieties and starting dates. The symbols (points and triangles) show the model results at two starting dates (in this example $t_0 = -121$ and $-152$ DOY), while the two colours represent fitting (red) and validation (blue) procedure, respectively. The filled symbols constitute the R² and AIC values of the models obtained by the selected HP parameter sets, while the dashed (fitting) and dotted (validation) lines depict the average R² and AIC values for both starting dates (in this example orange line = $-121$ DOY and green line = $-152$ DOY). The boxplots represent the R² and AIC values distribution for HP (orange and light blue: HP model application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively) and AV (boxplot in the shaded area, orange and light blue: AV application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively). The outliers represent all values out of the 1.5\*InterQuartile Range (IQR) of the boxplot. Figure S4: R² and AIC values of the UNIFIED model obtained for all fitted and average models, grapevine varieties and starting dates. The symbols (points and triangles) show the model results at two starting dates (in this example $t_0 = -121$ and $-152$ DOY), while the two colours represent fitting (red) and validation (blue) procedure, respectively. The filled symbols constitute the R² and AIC values of the models obtained by the selected HP parameter sets, while the dashed (fitting) and dotted (validation) lines depict the average R² and AIC values for both starting dates (in this example orange line = $-121$ DOY and green line = $-152$ DOY). The boxplots represent the R² and AIC values distribution for HP (orange and light blue: HP model application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively) and AV (boxplot in the shaded area, orange and light blue: AV application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively). The outliers represent all values out of the 1.5\*InterQuartile Range (IQR) of the boxplot. Figure S5: R² and AIC values of the UNIFORC model obtained for all fitted and average models, grapevine varieties and starting dates. The symbols (points and triangles) show the model results at two starting dates (in this example $t_0 = 1$ and $60$ DOY), while the two colours represent fitting (red) and validation (blue) procedure, respectively. The filled symbols constitute the R² and AIC values of the models obtained by the selected HP parameter sets, while the dashed (fitting) and dotted (validation) lines depict the average R² and AIC values for both starting dates (in this example orange line = $-121$ DOY and green line = $-152$ DOY). The boxplots represent the R² and AIC values distribution for HP (orange and light blue: HP model application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively) and AV (boxplot in the shaded area, orange and light blue: AV application on all validation subsamples for $t_0 = -121$ and $t_0 = -152$ DOY, respectively). The outliers represent all values out of the 1.5\*InterQuartile Range (IQR) of the boxplot. Figure S6: R² and AIC values of the GDD model obtained for all fitted and average models, grapevine varieties and starting dates. The symbols (points and triangles) show the model results at two starting dates (in this example $t_0 = 1$ and $60$ DOY), while the two colours represent fitting (red) and validation (blue) procedure, respectively. The filled symbols constitute the R² and AIC values of the models obtained by the selected HP parameter sets, while the dashed (fitting) and dotted (validation) lines depict the average R² and AIC values for both starting dates (in this example orange line = $1$ DOY and green line = $60$ DOY). The boxplots represent the R² and AIC values distribution for HP (orange and light blue: HP model application on all validation subsamples for $t_0 = 1$ and $t_0 = 60$ DOY, respectively) and AV (boxplot in the shaded area, orange and light blue: AV application on all validation subsamples for $t_0 = 1$ and $t_0 = 60$ DOY, respectively). The outliers represent all values out of the 1.5\*InterQuartile Range (IQR) of the boxplot. Figure S7: R² and AIC values of the WANG model obtained for all fitted and average models, grapevine varieties and starting dates. The symbols (points and triangles) show the model results at two starting dates (in this example $t_0 = 1$ and $60$ DOY), while the two colours represent fitting (red) and validation (blue) procedure, respectively. The filled symbols constitute the R² and AIC values of the models obtained by the selected HP parameter sets, while the dashed (fitting) and dotted (validation) lines depict the average R² and
AIC values for both starting dates (in this example orange line = 1 DOY and green line = 60 DOY). The boxplots represent the R² and AIC values distribution for HP (orange and light blue: HP model application on all validation subsamples for t₀ = 1 and t₀ = 60 DOY, respectively) and AV (boxplot in the shaded area, orange and light blue: AV application on all validation subsamples for t₀ = 1 and t₀ = 60 DOY, respectively). The outliers represent all values out of the 1.5*InterQuartile Range (IQR) of the boxplot. Figure S8: Correlations between observed and simulated budbreak Days of Year (DOY) of the UNIFORC and UNICHILL models. The results were obtained for all grapevine varieties and starting dates using the HP parameter set. Figure S9: Correlations between observed and simulated budbreak Days of Year (DOY) of the WANG and BRIN Daily models. The results were obtained for all grapevine varieties and starting dates using the HP parameter set. Figure S10: Highest-Performance (HP) and average (AV) parameters' set applied on the independent datasets by considering all phenological models and starting dates. The example of Cabernet sauvignon variety. Figure S11: Highest-Performance (HP) and average (AV) parameters' set applied on the independent datasets by considering all phenological models and starting dates. The example of Gewürztraminer variety. Figure S12: Highest-Performance (HP) and average (AV) parameters' set applied on the independent datasets by considering all phenological models and starting dates. The example of Grenache variety. Figure S13: Highest-Performance (HP) and average (AV) parameters' set applied on the independent datasets by considering all phenological models and starting dates. The example of Pinot Gris variety. Figure S14: Highest-Performance (HP) and average (AV) parameters' set applied on the independent datasets by considering all phenological models and starting dates. The example of Riesling variety. Figure S15: Highest-Performance (HP) and average (AV) parameters' set applied on the independent datasets by considering all phenological models and starting dates. The example of Touriga Franca variety. Figure S16: Highest-Performance (HP) and average (AV) parameters' set applied on the independent datasets by considering all phenological models and starting dates. The example of Touriga Nacional variety. Figure S17: Statistical results (R², AIC and RMSE) of the comparison between CF (black bars) and F (red bars) models in model fitting and application for all grapevine varieties using AV parameters' set. The differences between CF and F models in model fitting and applications were displayed using green and yellow bars, respectively. CS = Cabernet Sauvignon, CH = Chardonnay, GE = Gewürztraminer, GR = Grenache, PG = Pinot Gris, RS = Riesling, TF = Touriga Franca, TN = Touriga Nacional. Table S1: Statistical indices (RMSE, EF, R², AIC) obtained for all models, varieties and starting dates after model fitting and validation on all subsamples.

**Author Contributions:** Conceptualization, L.L., S.C.-A., J.A.S. and M.M.; methodology, L.L., S.C.-A., C.M., M.M., H.F., C.v.L., P.M.; data provision, D.M., J.J., A.D.-I., C.v.L., J.A.S., J.E.-D. and J.S.; writing—original draft preparation, L.L. and S.C.-A.; writing—review and editing, L.L., S.C.-A., J.A.S., C.M., H.F., D.M., P.M., J.J., T.K., A.D.-I., C.v.L., A.C.M., C.D., M.B. and M.M.; supervision, M.M.; funding acquisition, M.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the European Union's Horizon 2020 Research and Innovation Programme, under the Clim4Vitis project: "Climate change impact mitigation for European viticulture: knowledge transfer for an integrated approach", grant agreement no. 810176. It was also supported by FCT-Portuguese Foundation for Science and Technology, under the project UIDB/04033/2020 and the French National Research Agency (ANR) in the frame of the Investments for the Future Program, within the cluster of excellence COTE (ANR-10-LABX-45).

**Acknowledgments:** We wish to acknowledge the Clim4Vitis project "Climate change impact mitigation for European viticulture: knowledge transfer for an integrated approach", funded by European Union's Horizon 2020 Research and Innovation Program, under grant agreement no. 810176 that supported this study. H.F. thanks the FCT for contract CEECIND/00447/2017 and project COA/CAC/0030/2019.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Mullins, M.G.; Bouquet, A.; Williams, L.E. *Biology of the Grapevine*; Cambridge University Press: Cambridge, UK, 1992.

2. Van Leeuwen, C.; Seguin, G. The concept of terroir in viticulture. *J. Wine Res.* **2006**, *17*, 1–10. [CrossRef]

3. Jones, G.V. Winegrape phenology. In *Phenology: An Integrative Environmental Science*; Schwartz, M.D., Ed.; Springer: Dordrecht, The Netherlands, 2013; pp. 563–584.

4. Jones, G.V.; Davis, R.E. Climate influences on grapevine phenology, grape composition, and wine production and quality for Bordeaux, France. *Am. J. Enol. Vitic.* **2000**, *51*, 249–261.

5. Menzel, A. A 500 year pheno-climatological view on the 2003 heatwave in Europe assessed by grape harvest dates. *Meteorol. Zeitschrift* **2005**, *14*, 75–77.

6. Tomasi, D.; Jones, G.V.; Giust, M.; Lovat, L.; Gaiotti, F. Grapevine Phenology and Climate Change: Relationships and Trends in the Veneto Region of Italy for 1964-2009. *Am. J. Enol. Vitic.* **2011**, *62*, 329–339.

7. Van Leeuwen, C.; Destrac-Irvine, A.; Dubernet, M.; Duchêne, E.; Gowdy, M.; Marguerit, E.; Pieri, P.; Parker, A.; de Rességuier, L.; Ollat, N. An Update on the Impact of Climate Change in Viticulture and Potential Adaptations. *Agronomy* **2019**, *9*, 514. [CrossRef]
8. Ramos, M.C.; Jones, G.V.; Martínez-Casasnovas, J.A. Structure and trends in climate parameters affecting winegrape production in northeast Spain. *Clim. Res.* **2008**, *38*, 1–15. [CrossRef]

9. Moriondo, M.; Bindi, M.; Fagarazzi, C.; Ferrise, R.; Trombi, G. Framework for high-resolution climate change impact assessment on grapevines at a regional scale. *Reg. Environ. Chang.* **2011**, *11*, 553–567.

10. Gaál, M.; Moriondo, M.; Bindi, M. Modelling the impact of climate change on the Hungarian wine regions using Random Forest. *Appl. Ecol. Environ. Res.* **2012**, *10*, 121–140.

11. Moriondo, M.; Jones, G.V.; Bois, B.; Dibari, C.; Ferrise, R.; Trombi, G.; Bindi, M. Projected shifts of wine regions in response to climate change. *Clim. Change* **2013**, *119*, 825–839. [CrossRef]

12. Wolkovich, E.M.; García de Cortázar-Atauri, I.; Morales-Castilla, I.; Nicholas, K.A.; Lacombe, T. From Pinot to Xinomavro in the world's future wine-growing regions. *Nat. Clim. Chang.* **2018**, *8*, 29–37. [CrossRef]

13. Fraga, H.; Pinto, J.G.; Santos, J.A. Climate change projections for chilling and heat forcing conditions in European vineyards and olive orchards: A multi-model assessment. *Clim. Chang.* **2019**, *152*, 179–193. [CrossRef]

14. Malheiro, A.C.; Campos, R.; Fraga, H.; Eiras-Dias, J.; Silvestre, J.; Santos, J.A. Winegrape phenology and temperature relationships in the Lisbon wine region, Portugal. *J. Int. Sci. Vigne Vin* **2013**, *47*, 287–299.

15. Costa, R.; Fraga, H.; Fonseca, A.; De Cortázar-Atauri, I.G.; Val, M.C.; Carlos, C.; Reis, S.; Santos, J.A. Grapevine phenology of cv. Touriga Franca and Touriga Nacional in the Douro wine region: Modelling and climate change projections. *Agronomy* **2019**, *9*, 1–20.

16. Molitor, D.; Junk, J. Climate change is implicating a two-fold impact on air temperature increase in the ripening period under the conditions of the Luxembourgish grapegrowing region. *OENO One* **2019**, *53*, 409–422. [CrossRef]

17. Shulman, Y.; Nir, G.; Fanberstein, L.; Lavee, S. The effect of cyanamide on the release from dormancy of grapevine buds. *Sci. Hortic.* (Amsterdam) **1983**, *19*, 97–104.

18. Wicks, A.S.; Johnson, J.O.; Bracho, E.; Jensen, F.L.; Neja, R.A.; Lider, L.A.; Weaver, R.J. Induction of early and more uniform budbreak in *Vitis vinifera* L. vs. Perlette, Thompson Seedless, and Flame Seedless. In Proceedings of the Bud Dormancy in Grapevine: Potential and Practical Uses of Hydrogen Cyanamide on Grapevine, Fresno, CA, USA, 20 August 1984; Weaver, R.J., Ed.; University of California: Fresno, CA, USA, 1984; pp. 48–58.

19. Dokoozlian, N.K. Chilling temperature and duration interact on the budbreak of "Perlette" grapevine cuttings. *HortScience* **1999**, *34*, 1054–1056. [CrossRef]

20. Sarvas, R. Investigations on the annual cycle of development of forest trees. Autumn dormancy and winter dormancy. *Commun. Inst. For. Fenn.* **1974**, *84*, 1–101.

21. Pouget, R. Considérations Générales sur le rythme végétatif et la dormance des bourgeons de la Vigne. *Vitis* **1972**, *217*, 198–217.

22. Pouget, R. Nouvelle conception du seuil de croissance chez la vigne. *Vitis* **1968**, *7*, 201–205.

23. Richardson, E.A.; Seeley, S.D.; Walker, D.R.; Walker, D.I. A model for estimating the completion of rest for "Redhaven" and "Elberta" peach trees. *Hortic. Sci.* **1974**, *9*, 331–332.

24. Chuine, I. A unified model for budburst of trees. *J. Theor. Biol.* **2000**, *207*, 337–347. [CrossRef] [PubMed]

25. Cesaraccio, C.; Spano, D.; Snyder, R.L.; Duce, P. Chilling and forcing model to predict bud-burst of crop and forest species. *Agric. For. Meteorol.* **2004**, *126*, 1–13. [CrossRef]

26. García de Cortazar-Atauri, I.; Brisson, N.; Seguin, B.; Gaudillere, J.; Baculat, E. Simulation of budbreak date for vine. The bring model. Some application in climate change study. In Proceedings of the XIV International GESCO Viticulture Congress, Geisenheim, Germany, 23–27 August 2005; pp. 485–490.

27. Nendel, C. Grapevine bud break prediction for cool winter climates. *Int. J. Biometeorol.* **2010**, *54*, 231–241. [CrossRef]

28. Caffarra, A.; Eccel, E. Projecting the impacts of climate change on the phenology of grapevine in a mountain area. *Aust. J. Grape Wine Res.* **2011**, *17*, 52–61. [CrossRef]

29. García de Cortázar-Atauri, I.; Brisson, N.; Gaudillere, J.P. Performance of several models for predicting budburst date of grapevine (*Vitis vinifera* L.). *Int. J. Biometeorol.* **2009**, *53*, 317–326. [CrossRef]

30. Fila, G.; Gardiman, M.; Belvini, P.; Meggio, F.; Pitacco, A. A comparison of different modelling solutions for studying grapevine phenology under present and future climate scenarios. *Agric. For. Meteorol.* **2014**, *195–196*, 192–205. [CrossRef]
31. Leolini, L.; Moriondo, M.; Fila, G.; Costafreda-Aumedes, S.; Ferrise, R.; Bindi, M. Late spring frost impacts on future grapevine distribution in Europe. *Field Crops Res.* **2018**, *222*, 197–208. [CrossRef]

32. Moriondo, M.; Ferrise, R.; Trombi, G.; Brilli, L.; Dibari, C.; Bindi, M. Environmental Modelling & Software Modelling olive trees and grapevines in a changing climate *. *Environ. Model. Softw.* **2015**, *72*, 387–401.

33. Molitor, D.; Junk, J.; Evers, D.; Hoffmann, L.; Beyer, M. A high-resolution cumulative degree day-based model to simulate phenological development of grapevine. *Am. J. Enol. Vitic.* **2014**, *65*, 72–80. [CrossRef]

34. Menne, M.J.; Durre, I.; Vose, R.S.; Gleason, B.E.; Houston, T.G. An overview of the global historical climatology network-daily database. *J. Atmos. Ocean. Technol.* **2012**, *29*, 897–910. [CrossRef]

35. Destrac-Irvine, A.; van Leeuwen, C. VitAdapt: An experimental program to study the behavior of a wide range of Vitis vinifera varieties in a context of climate change in the Bordeaux vineyards. In *Proceedings of the Sustainable Grape and Wine Production in the Context of Climate Change*. ClimWine 2016 International Symposium, Bordeaux, France, 10–13 April 2016.

36. Garcia de Cortazar-Atauri, I.; Audergon, J.M.; Delzon, S.; Anger, C.; Duchêne, E.; Bonhomme, M.; Legave, J.M.; Chuine, I.; Raynal, H.; Davi, H.; et al. PERPHECLIM ACCAF Project-perennial fruit crops and forest phenology evolution facing climatic changes. In *Proceedings of the XVI International Symposium on Apricot Breeding and Culture and XV Chinese National Symposium on Plum and Apricot, Shenyang, China*, 29 June–3 July 2015.

37. Jones, G.V. Climate and Terroir: Impacts of Climate Variability and Change on Wine. In *Fine Wine and Terroir—the Geoscience Perspective*; Meinert, L.D., Ed.; Geological Association of Canada: St. John's, NL, Canada, 2006; p. 247.

38. Eichhorn, V.K.W.; Lorenz, D.H. Phänologische Entwicklungsstadien der Rebe [Phenological development stages of the grapevine]. *Nachr. Dtsch. Pflanzenschutzd.* **1977**, *29*, 119–120.

39. Amerine, M.A.; Winkler, A.J. Composition and Quality of Musts and Wines of California Grapes. *Hilgardia* **1944**, *15*, 493–675. [CrossRef]

40. Wang, E.; Engel, T. Simulation of phenological development of wheat crops. *Agric. Syst.* **1998**, *58*, 1–24. [CrossRef]

41. Bidabe, B. Contrôle de l'époque de floraison du pommier par une nouvelle conception de l'action de températures. C. R. Séances L'Acad. D'Agric. Fr. **1963**, *49*, 934–945.

42. Bidabe, B. *L'action Des Températures Sur L'évolution Des Bourgeons de L'entrée en Dormance à la Floraison*; Congrès Pomologique: Paris, France, 1965.

43. Tukey, J.W. Box-and-whisker plots. In *Exploratory Data Analysis*; Addison-Wesley: Reading, MA, USA, 1977; pp. 39–43. ISBN 978-0201076165.

44. Shellie, K.C. Viticultural performance of red and white wine grape cultivars in southwestern Idaho. *Horttechnology* **2007**, *17*, 595–603. [CrossRef]

45. Parker, A.K.; Garcia de Cortazar-Atauri, I.; Van Leeuwen, C.; Chuine, I. General phenological model to characterise the timing of flowering and veraison of Vitis vinifera L. Aust. J. Grape Wine Res. **2011**, *17*, 206–216. [CrossRef]

46. Fila, G. Mathematical Models for the Analysis of the Spatio-Temporal Variability of Vine Phenology [Modelli Matematici per l'analisi Della Variabilità Spazio-Temporale Della Fenologia Della Vite]. Ph.D. Thesis, University of Padua, Padua, Italy, 2012.

47. Streck, N.A.; Weiss, A.; Xue, Q.; Stephen Baenziger, P. Improving predictions of developmental stages in winter wheat: A modified Wang and Engel model. *Agric. For. Meteorol.* **2003**, *115*, 139–150. [CrossRef]

48. Chuine, I.; Garcia de Cortazar-Atauri, I.; Kramer, K.; Hänninen, H. Plant development models. In *Phenology: An Integrative Environmental Science*; Schwartz, M.D., Ed.; Springer: Dordrecht, The Nederlands, 2013; pp. 275–293. ISBN 9789400769250.

49. Hastings, W.K. Monte carlo sampling methods using Markov chains and their applications. *Biometrika* **1970**, *57*, 97–109. [CrossRef]

50. Kirkpatrick, S.; Gelatt, C.D.; Vecchi, M.P. Optimization by simulated annealing. *Science* **1983**, *220*, 671–680. [CrossRef] [PubMed]

51. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.H.; Teller, E. Equation of state calculations by fast computing machines. *J. Chem. Phys.* **1953**, *21*, 1087–1093.

52. R Development Core Team. *R: A Language and Environment for Statistical Computing*; Version 3.5.1; R Foundation for Statisrical Computing: Vienna, Austria, 2018; ISBN 3900051070.
53. Wright, S. Correlation and causation. *J. Agric. Res.* **1921**, *20*, 557-585.

54. Akaike, H. Information Theory and an Extension of the Maximum Likelihood Principle. In Proceedings of the Second International Symposium on Information Theory, Tsahkadsor, Armenia, 2-8 September 1973; Petrov, B.N., Caski, F., Eds.; Akademiai Kiado: Budapest, Hungary, 1973; pp. 267-281.

55. Kenney, J.F. Root mean square. In *Mathematics of Statistics*; Chapman & Hall, LTD.: Princeton, NJ, USA, 1962; pp. 59-60.

56. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. *J. Hydrol.* **1970**, *10*, 282-290.

57. Cola, G.; Failla, O.; Maghradze, D.; Megrelidze, L.; Mariani, L. Grapevine phenology and climate change in Georgia. *Int. J. Biometeorol.* **2017**, *61*, 761-773. [CrossRef] [PubMed]

58. Fila, G.; Di, B.; Gardiman, M.; Storchi, P.; Tomasi, D.; Silvestroni, O.; Pitacco, A. Calibration and validation of grapevine budburst models using growth-room experiments as data source. *Agric. For. Meteorol.* **2012**, *160*, 69-79. [CrossRef]

59. Duchêne, E.; Huard, F.; Dumas, V.; Schneider, C.; Merdinoglu, D. The challenge of adapting grapevine varieties to climate change. *Clim. Res.* **2010**, *41*, 193-204. [CrossRef]

60. Reis Pereira, M.; Ribeiro, H.; Abreu, I.; Eiras-Dias, J.; Mota, T.; Cunha, M. Predicting the flowering date of Portuguese grapevine varieties using temperature-based phenological models: A multi-site approach. *J. Agric. Sci.* **2018**, *156*, 865-876. [CrossRef]

61. Moriondo, M.; Orlandini, S.; Nuntiis, P.D.; Mandrioli, P. Effect of agrometeorological parameters on the phenology of pollen emission and production of olive trees (*Olea europea* L. *Aerobiologia* **2001**, *17*, 225-232.

62. Chuine, I.; Bonhomme, M.; Legave, J.M.; García de Cortázar-Atauri, I.; Charrier, G.; Lacointe, A.; Améglio, T. Can phenological models predict tree phenology accurately in the future? The unrevealed hurdle of endodormancy break. *Glob. Chang. Biol.* **2016**, *22*, 3444-3460.

63. Martínez-Lüscher, J.; Kizildeniz, T.; Vučetić, V.; Dai, Z.; Luedeling, E.; van Leeuwen, C.; Gomès, E.; Pascual, I.; Irigoyen, J.J.; Morales, F.; et al. Sensitivity of Grapevine Phenology to Water Availability, Temperature and CO2 Concentration. *Front. Environ. Sci.* **2016**, *4*, 48. [CrossRef]

64. Fraga, H.; García de Cortázar Atauri, I.; Malheiro, A.C.; Santos, J.A. Modelling climate change impacts on viticultural yield, phenology and stress conditions in Europe. *Glob. Chang. Biol.* **2016**, *22*, 3774-3788. [CrossRef] [PubMed]

65. Hannah, L.; Roehrdanz, P.R.; Ikegami, M.; Shepard, A.V.; Shaw, M.R.; Tabor, G.; Zhi, L.; Marquet, P.A.; Hijmans, R.J. Climate change, wine, and conservation. *Proc. Natl. Acad. Sci. USA* **2013**, *110*, 6907-6912. [PubMed]

66. Fraga, H.; Santos, J.A.; Malheiro, A.C.; Oliveira, A.A.; Moutinho-Pereira, J.; Jones, G.V. Climatic suitability of Portuguese grapevine varieties and climate change adaptation. *Int. J. Climatol.* **2016**, *36*, 1-12.

67. Camargo-Alvarez, H.; Salazar-Gutiérrez, M.; Keller, M.; Hoogenboom, G. Modeling the effect of temperature on bud dormancy of grapevines. *Agric. For. Meteorol.* **2020**, *280*, 107782. [CrossRef]

68. Fennell, A.; Mathiason, K.; Luby, J. Genetic segregation for indicators of photoperiod control of dormancy induction in vitis species. *Acta Hortic.* **2005**, *689*, 533-540. [CrossRef]

69. Rubio, S.; Dantas, D.; Bressan-Smith, R.; Pérez, F.J. Relationship Between Endodormancy and Cold Hardiness in Grapevine Buds. *J. Plant Growth Regul.* **2016**, *35*, 266-275. [CrossRef]

70. Leolini, L.; Moriondo, M.; Romboli, Y.; Gardiman, M.; Costafreda-Aumedes, S.; Costafreda-Aumedes, S.; Bindi, M.; Granchi, L.; Brilli, L. Modelling sugar and acid content in Sangiovese grapes under future climates: An Italian case study. *Clim. Res.* **2019**, *78*, 211-224.

<figure><img src="image_9.png" /></figure>