Quantifying the seasonal variations in grapevine yield components based on pre and postflowering weather conditions
Abstract
Seasonal differences in weather conditions cause marked variation in grapevine yield. However, quantitative relationships between various yield components and climatic factors at field scales are still lacking. By using a longterm field trial, we quantified the correlation between weather conditions during the key development stages and the yield components of Vitis vinifera L. Sauvignon blanc growing under coolclimate conditions. A longterm phenology and yield monitoring trial using both twocane and fourcane trained vertically shoot positioned (VSP) Sauvignon blanc vines was established in four vineyards in Marlborough, New Zealand in 2004. Phenology, bunch number, berry mass, yield and meteorology records were collated. A multivariable mixed linear model was used to assess the relationship between various yield components and weather conditions. The critical periods for each yield component and weather factor were optimised based on the maximum likelihood returned from the mixed linear model. The optimised critical periods of temperature for all yield components occurred mainly before 50 % flowering either in the previous season (during inflorescence initiation) and the current season, indicating the importance of the preflowering period on yield formation. Out of all weather factors, maximum daily temperature had the largest effect on bunch number and overall yield and strongly influenced berry number and bunch mass. Rainfall near flowering time had a negative effect on berry mass and bunch mass, but postflowering rainfall had a strong positive effect. The statistical model explained 60 to 85 percent of the seasonal variations in bunch number, berry number, berry and bunch mass and yield per vine.
Introduction
Seasonal or interseasonal factors, especially temperature, radiation and water status, have pronounced effects on the fruitfulness and overall yield of grapevine (Buttrose, 1969b; Buttrose, 1974). Maintaining stable yields from year to year is essential for achieving a consistent fruit quality and supply, which are crucial within the context of increasing competition in the international market (Kliewer and Dokoozlian, 2005; Trought, 2000).
Grapevine initiate inflorescence primordia in the summer of the year prior to that in which they flower (Figure 1). This enables us to anticipate potential yield early in the season (Perold, 1927), allowing adjustments of crop load to be made during winter pruning and after fruit set. In this study, we use longterm field experimental data to investigate the relationship between weather conditions and yield.
Figure 1. The reproductive cycle of a grapevine under southern hemisphere conditions. The reproductive cycle starts at the time of induction and ends at harvest, which occurs about 1518 months after induction.
The numbers before the phenology stage represent the Modified EL system code for growth stages (Coombe, 1995a). The idea of using a circle to represent the reproductive cycle was adapted from Wilson (1996). A figure detailing the yield formation processes is provided in Supplementary Figure S1.
Yield per hectare is the product of vines per hectare, shoots per vine, inflorescences per shoot, berry number per bunch and berry mass. Grapevines may generally form up to four inflorescences on a shoot (Watt, 2010). An undifferentiated mass of cells (anlagen, potential bunches for the following season) develop in the basal node approximately 15 days before bloom, with approximately 12 expanded leaves on the shoot (Swanepoel and Archer, 1988; Watt, 2010; Watt et al., 2008). This process progresses along the developing shoot as further leaves appear (Vasconcelos et al. 2009). Whether or not anlagen become inflorescence or tendril primordia will depend on temperature and radiation conditions around the bud at the time differentiation takes place (Buttrose, 1969a, 1969b, 1970). High radiation intensities or temperatures during initiation will encourage inflorescence primordia development, while shaded or cool conditions will lead to tendril formation (Buttrose, 1970; Sanchez and Dokoozlian, 2005; Trought, 2005). In addition to bud fruitfulness, temperature can influence the primary branching of inflorescence primordia, which can account for 51 to 81 % of the flower number per inflorescence the following season (Dunn and Martin, 2007). In this paper, we use ‘inflorescence’ to refer to the fruit bearing part of the vine that has not set berries and ‘bunch’ to refer to the fruit bearing part of the vine that has set berries.
Berry number per bunch is the product of flowers per inflorescence and percentage of flowers that set fruit. Cool temperatures shortly before budburst can increase the total flower number per shoot (Eltom et al., 2017). However, cold weather after budburst may prevent the formation of individual flowers (Buttrose and Hale, 1973). Fruit set (or the proportion of flowers that are retained as berries) represents a changeover from the static condition of the fully developed flower to the rapidly growing condition of the young fruit (Coombe, 1962). Both low and high temperature can affect fruit set in grapevines. Preflowering low temperatures can disrupt both the formation and function of ovules and pollen during flowering. Likewise, exposure to temperatures of 15 or 40 °C in Pinot noir and Carignane can significantly lower fruit set (Kliewer, 1977). Temperatures below 15°C and above 32 °C are considered detrimental. However, no strong relationships between daily assessments of temperature at cap fall and berry set or berry development were found by Friend (2005), who suggests that light intensity (due to its impact on current photoassimilate supply) may be a more valid environmental index than temperature.
Berry fresh mass is strongly linked with seed dry weight, while berry dry weight is more strongly influenced by leaf area to fruit number ratio (Friend et al., 2009; Petrie et al., 2000). The variation in total seed dry weight is related to changes in the extent of ovule fertilisation, which is related to the environmental and plant conditions before and after flowering (Ebadi et al., 1996; Friend et al., 2009). Preflowering or immediate post flowering low temperatures can slow pollen tube growth or disrupt ovule and seed development in Chardonnay (Ebadi et al., 1996; Staudt, 1982). However, no strong relationship between daily temperature at cap fall and berry weight in Merlot was found by Friend (2005), who suggested that radiation intensity (correlated with current carbohydrates supply) may be a more valid environmental index than temperature.
Despite the importance of grapevine yield on various viticultural and winemaking practices, little effort has been made to predict yield (Dunn and Martin, 2007; Trought, 2005). Santos et al. (2011) have developed a multivariate linear regression model of grapevine yield (19862008) in response to monthly mean temperatures and monthly precipitation totals from March to June for the Demarcated Region of Douro in northeast Portugal. Bock et al. (2013) have analysed the correlation between longterm (1805–2010) grapevine yield records in Lower Franconia, Germany at the mean MayAugust temperature and sugar content at the mean AprilAugust temperature. However, they did not consider the physiological basis of the grapevine yield formation, nor the effects of weather events during key developmental stages on yield. Furthermore, current plant growth models predict grapevine yield as a fraction of total dry matter produced during the berry growth period divided by the berry water content (Faluomi, 2017; Fraga et al., 2016), using a sourcesink carbon allocation approach. Thus, there is an urgent need to include the variations in shoot, bunch and berry number in plant models in order to make reliable yield predictions.
This study aimed to quantify the relationships between various yield components and climatic conditions during critical periods using a longterm field experiment of Vitis vinifera L. Sauvignon blanc. Sauvignon blanc was chosen as it is the main variety grown in New Zealand, representing 61 % of the total production area and 80 % of the yield (New Zealand Winegrowers, 2018). About 90 % of Sauvignon blanc in New Zealand is grown in Marlborough, which is considered to be a cool climate grape growing region (Jackson and Cherry, 1988; Jon et al., 2012) and which contributes to the production of grapes with distinctive flavours. We aimed to answer the following questions (Figure S1): 1) what are the correlations between each yield component and different weather factors? (e.g., maximum or minimum temperature and radiation) 2) Is there a dominant weather factor? 3) What are the critical periods that determine each yield component? 4) Can we explain the variations in different yield components based on weather conditions during critical periods?
Materials and methods
1. Experiment
Starting in 2004, data regarding phenology (budburst, flowering and fruit development) and yield per vine, bunch number (from which bunch mass was calculated) and berry mass (from which berry number per bunch was calculated) were collected from Sauvignon blanc vines in four commercial vineyards in Marlborough. Vines planted in rows 2.4 m apart and 1.8 m within each row were cane pruned and trained using vertical shoot positioning (VSP). The lowest fruiting wire was 90 cm from the ground and the top fruiting wire was 110 cm. Foliage wires were used to maintain a tight VSP canopy. Vines were either twocane or fourcane pruned by the authors, retaining 20 or 40 buds respectively. This helped separate out the seasonal differences in yield from the management effects mainly caused by the buds retained after pruning. Active canopy management was practised throughout the growing seasons, including wirelifting, canopy trimming and leaf plucking after flowering and before véraison. Vines were trickle irrigated, the timing and volumes determined by the vineyard manager. Pest and diseases were managed using industry protocols (www.nzwine.com/swnz/).
Data were collected from eight replicate plots of four vines (four in recent two seasons), planted between vineyard posts (bays), in each vineyard. Trunk diameters were measured on all vines in eight rows of the vineyards and plots were chosen to represent average size vines. The four chosen vineyards are spread across the predominant vineyard area of Marlborough, namely Upper Brancott (UB) [41.56569 S; 173.85154 E], Western Wairau (WW) [41.51113 S; 173.77012 E], Seaview Awatere (SA) [41.62863 S; 174.12950 E] and Central Rapaura (CR) [41.47617 S; 173.88981 E]. The current dataset contains 15 seasons (20042005 to 20182019) of phenology and 12 seasons of yield component records (except CR which has records of 15 seasons) for twocane pruned Sauvignon blanc and 5 seasons of phenology and 13 seasons of yield component records for fourcane pruned Sauvignon blanc (see detailed history records in Supplementary Method S1). One additional vineyard with four seasons (20142015 to 20172018) of phenology and yield component records for threecane pruned vines was added, namely Glasneven in the Canterbury region [43.10592 S; 172.74163 E].
Meteorological stations on or close to each site were used to examine the relationship between weather conditions and yield components. Four Marlborough stations were located between 0.1 and 2 km from the observation rows, namely UB [41.54248 S; 173.84736 E], WW [41.51343 S; 173.75983 E], SA [41.62959 S; 174.13096 E] and CR [41.49137 S, 173.8891 E]. Data were available from HortPlus NZ LTD. One National Institute of Water and Atmospheric Research (NIWA) station (Agent Number 26607, 43.06861 S, 172.65346 E) was used for the Waipara observation and was located 8 km from the observation rows. The meteorological stations complied with World Meteorological Organisation specifications, which require weather instruments to be sited in an open area away from buildings and shelter. The temperature and rainfall instruments were calibrated annually and data were downloaded and checked weekly. All stations recorded maximum and minimum daily temperature, rainfall and relative humidity. Radiation was not recorded in UB and CR stations. The radiation records from Blenheim central meteorology station (NIWA Agent Number 12430, 41.49722 S, 173.96292 E) were used for those two sites.
The recorded mean annual rainfall between 2002 and 2019 was 684 mm for UB, 910 mm for WW, 610 mm for SA and 723 mm for CR. The annual Penman evapotranspiration demand was 1022.6 mm according to records from the Blenheim meteorology station. The soil profile texture was loam originating from alluvium for UB and WW, silty loam originating from loess for SA and silty loam over sandy loam originating from alluvium for CR. The soil profile available water ranged from 80 mm to 140 mm per meter across sites and locations within the sites. Seasonal irrigation (mainly during summer) was between 200 mm to 450 mm. Predawn leaf water potential varied from 0.1 to 0.4 MPa during the growing season based on one of our irrigation trials (unpublished). The maximum leaf area per vine was about 7 m^{2 }for twocane pruned vines and 10 m^{2} for fourcane pruned vines. The yield per vine was about 5 kg for twocane pruned vines and 7 kg for fourcane pruned vines with significant variations among years and sites.
Flowering progression was estimated visually, twice a week, from late November through to late December (depending on the season) by assessing each of the inflorescences on all of the shoots arising from one cane in each bay (i.e., four canes per vineyard). The proportion of opened flowers per inflorescence was recorded in 5 % increments.
A random 32berry sample was collected weekly from eight different bunches across all canes in each of the eight or four monitored bays. Berry sample collection started shortly before véraison and continued until harvest, to determine the berry mass and total soluble solids. A threshold soluble solids concentration of 8°Brix was used as an alternate measure of the midpoint (50 %) of véraison, which was interpolated/extrapolated from soluble solids accumulation data (Parker et al., 2014).
At harvest, eight monitored bays (four in recent two seasons) were hand harvested. All bunches were counted and weighed from each bay. Bunches with severe botrytis infection (>10 % visually assessed) were counted and weighed separately. The average bunch number per vine was calculated from the total number of bunches harvested from the four vines in the bay. Average bunch number per vine included bunches from the shoots along the canes (count shoots), as well as bunches on shoots arising from quiescent buds (noncount shoots) from the vine’s head and trunk. Average bunch mass in this study was calculated based on bunches with less than 10 % botrytis infection. Berry numbers per bunch were estimated based on the average bunch mass and mean berry mass determined by the 32berry sample at harvest. In order to exclude the effects of botrytis infection on observed yield, final yield was calculated as the average bunch mass × total bunch number × mean berry mass. On average, the calculated yield excluding botrytis infection was 1.7 % higher than the direct yield measurement.
During the analyses, berry number per bunch and berry mass of different treatments were grouped. There was a tendency for the berry mass in twocane pruned vines to be slightly higher than that in fourcane pruned vines during the berry development. However, this difference was not consistent between sites and seasons and was also affected by the time of harvest. Fourcane pruned vines were normally harvested one or two weeks later than twocane pruned vines. The Glasneven vineyard and data were only used for assessing the berry number per bunch, berry mass and bunch mass; the data were excluded from the analysis of yield per vine as it was threecane pruned.
2. Analysis procedures
We hypothesised that bunch number per vine was determined by the weather conditions during the flowering periods of the previous season (inflorescence initiation) and berry number was determined by the weather conditions around flowering of the current season based on previous studies (Buttrose, 1970, 1974; May, 2000; Vasconcelos et al., 2009). We further hypothesised that potential berry mass is determined by the environmental and plant conditions during flowering and fruit set period of the current season, while potential berry size is affected by radiation, leaf area to fruit number ratio, water status, etc., between fruit set and harvest. Thus weather conditions during three periods (flowering periods in the previous season (inflorescence initiation), flowering periods of the current season and post flowering periods, but before véraison of the current season) were analysed in more detail than during other periods (see detailed descriptions below). Other periods (e.g., before and after budburst of the current season and post véraison) were also explored, but they showed little effect and were therefore not reported.
In summary, the effects of mean daily temperature (Tmean, °C), daily maximum temperature (Tmax, °C), daily minimum temperature (Tmin, °C), radiation intensity (Ra, MJ day^{1}), cumulative rainfall around flowering (RainTotFlow, mm) and number of rainfall days around flowering (RainDay) in the previous season (denoted by Ini after the factor; e.g., TmaxIni) and in the current season (denoted by adding Flow; e.g., TmeanFlow) on different yield components were tested. In addition, the effects of rainfall, vapour pressure deficit, potential transpiration and the difference between potential transpiration and rainfall after flowering  but before véraison in the current season (noted by Ver, e.g. RainTotVer)  were also tested. The effects of each weather factor were tested on all yield components: bunch number per vine, berry number per bunch, berry mass, bunch mass and yield per vine.
An optimisation procedure was developed to find the critical period which would give the maximum likelihood between a certain weather factor during that period and the yield component in question (see the overall analysis procedure in Figure S2). The procedure used the recorded 50 % flowering (or 50 % véraison) time as input and tried to optimise two parameters that defined the period: one parameter defined the time before flowering and the other one defined the time after flowering (Eq. 1). The concept of thermal day was used to standardise the periods in different years with different temperatures (Eq. 2 and 3).
${T}_{i}=f\left(Temp,{FL}_{season1,}{TD}_{backward},{TD}_{forward}\right)$
Eq. 1
${f}_{T}=\frac{2{({T}_{\alpha}{T}_{\mathrm{m}\mathrm{i}\mathrm{n}})}^{\alpha}{({T}_{\mathrm{o}\mathrm{p}\mathrm{t}}{T}_{\mathrm{m}\mathrm{i}\mathrm{n}})}^{\alpha}{({T}_{\alpha}{T}_{\mathrm{m}\mathrm{i}\mathrm{n}})}^{2\alpha}}{{({T}_{\mathrm{o}\mathrm{p}\mathrm{t}}{T}_{\mathrm{m}\mathrm{i}\mathrm{n}})}^{2\alpha}}\begin{array}{cc}& {T}_{\mathrm{m}\mathrm{i}\mathrm{n}}\le T\le {T}_{\mathrm{m}\mathrm{a}\mathrm{x}}\end{array}$
Eq. 2
$\alpha =\mathrm{l}\mathrm{n}2/\mathrm{l}\mathrm{n}\left[({T}_{\mathrm{m}\mathrm{a}\mathrm{x}}{T}_{\mathrm{m}\mathrm{i}\mathrm{n}})/({T}_{\mathrm{o}\mathrm{p}\mathrm{t}}{T}_{\mathrm{m}\mathrm{i}\mathrm{n}})\right]$
Eq. 3
T_{i} is the average max, min or mean daily temperature of the period defined by flowering time in the previous season (FL_{season1}) minus ${TD}_{backward}$ and FL_{season1} plus ${TD}_{forward}$. Thermal day (td) is the sum of thermal time units calculated using a temperature response curve (Eq. 2). f_{T} is the thermal time unit for a day with daily mean temperature T_{mean}. For seasonal accumulated thermal days td was summed from July 1^{st}. T_{min} is the minimum temperature for growth, which is assumed to be 0 °C for the period before budburst and 4 °C for the growth after budburst (García de CortázarAtauri et al., 2010). T_{opt} is the optimal temperature for growth, which is assumed to be 22 °C for the period before budburst and 28 °C for the growth after budburst. T_{max} is the maximum temperature beyond which no growth will occur, which is assumed to be 40 °C for the whole period. One thermal day in our calculation corresponds to 1.04 actual days if the mean daily temperature is 25 ^{o}C and to 1.56 actual days if the mean daily temperature is 18 ^{o}C.
A wide exploration of the parameter values that defined the periods was initially carried out to find the most plausible periods for fine optimisation (Figure S3). The package DEoptim (Ardia et al., 2011) in R (R Development Core Team, 2019) was used for conducting the optimisation in parallel, thus improving the speed of optimisation. A critical period was optimised for each factor. The critical periods for twocane and fourcane treatments were checked separately at first. As they did not practically differ from each other, a single critical period for each factor was optimised for the whole dataset simultaneously. The optimised periods for each factor on different yield components are shown in Supplementary Table S1.
The linear mixedeffects model (lmer) from the R package of ‘lme4’ (Bates et al., 2014) was used to assess the relationship between different weather factors and yield component, as well as to assess maximum likelihood during optimisation (Eq. 4). Vineyard was included as a random intercept factor in the regression. Random slopes were not included, as they rapidly decreased the degrees of freedom of the model.
$BerryNum{T}_{mean,flow}*{RainTot}_{flow}+{T}_{max,ini}+\left(1\vee Vineyard\right)$
Eq. 4
The ~ sign is the notation for the formula in R. The plus sign in the equation means that interaction between those two factors was not included; star (*) sign means interaction was included.
Maximum likelihood was returned by the basic R function logLik with restricted maximum likelihood (REML) equal to false; e.g., logLik(mod, REML=F). The REML was set to false, because we wanted to compare models with different fixed effects using likelihood ratio test and ANOVA (analysis of variance) test and the REML method was more used for estimating random effects (Hui et al., 2019). The package lmerTest developed by Kuznetsova et al. (2017) was used to obtain the pvalue of each factor.
The potential bias of parameter estimated by the mixed linear model caused by the year and site was evaluated by the bootstrap method using the function bootMer in the ‘lme4’ package (Bates et al., 2014). Bootstrap is a general approach to statistical inference based on building a sampling distribution for a statistic by resampling repeatedly from the data at hand. The bias and standard error of the parameter from the true value estimated by the boot method were reported (Supplementary Table S2).
3. Model selection
The correlation between each yield component and weather conditions was first analysed with the R package of PerformanceAnalytics (Peterson et al., 2018) to obtain an overview of the correlations. Afterwards, the combinations of different factors in the linear mixed model for improving prediction performance were tested, resulting in a large number of models. Model selection was based on loglikelihood, R squared (R^{2}) and Akaike information criterion AIC (Sakamoto et al., 1986). AIC rates model in terms of parsimony and efficiency, where the lowest value is associated with the best model.
$AIC=2*loglikelihood+k*{N}_{p}$
Eq. 5
Where N_{p} represents the number of parameters in the fitted model, k = 2 for the usual AIC, or k = log(n) (n being the number of observations) for the socalled BIC. A revised R^{2} for the linear mixedeffects model was calculated using the r.squaredGLMM function in the R package of MuMIn (Bartoń, 2018). This function returns a marginal R^{2} for the fixed effects and a conditional R^{2} of the entire model including both fixed and random effects.
When there were two or more factors in the regression, the contribution of all factors were first checked and only the factors with a significant contribution were retained. Afterwards, the interaction term between all the factors was checked. The criteria for including the interaction were: 1) the interaction term has significant contribution; we accepted that the main factor would become nonsignificant after introducing the interaction term and 2) the model with the interaction term improves significantly compared to the model without the interaction term (ANOVA test).
When the best model differed between treatments, we first checked the best model in each treatment to see whether all the factors were significant. If one factor was significant in one treatment and not in another, we tended to include this factor to increase the stability of the model performance under different conditions. A list of all tested models and their regression results for bunch number per vine is shown in Supplementary Table S3; berry number per bunch is shown in Supplementary Table S4; berry mass is shown in Supplementary Table S5; bunch mass is shown in Supplementary Table S6; yield per vine is shown in Supplementary Table S7. The relationships between each yield component with the highest correlation factors are shown in separate figures. It should be noted that the final selected model only represent the highest model parsimony and efficiency. It may not include all the factors that would affect the yield Nonlinear response functions (e.g., logistic responses) were also tested. However, we could not justify the nonlinear response in our dataset.
Results
1. Overview of the yield components
The mean yield across sites and years for twocane pruned vines was 5.1 ± 0.21 (standard error) kg per vine and for fourcane pruned vines it was 7.9 ± 0.30 kg per vine (Figure 2 a, b). UB and WW were the highest yielding sites of the four sites and SA was the lowest yielding site in both twocane and fourcane treatments. Of all the sites, UB had the highest variation in vine yield between years.
The mean bunch number per vine for twocane pruned vines was 39.3 ± 0.8 per vine and for fourcane pruned vines it was 64.7 ± 1.2 per vine (Figure 2 c, d). Mean berry number per bunch across sites and years was 63.4 ± 1.1 (Figure 2e) and mean berry mass was 1.99 ± 0.02 g (Figure 2f).
Figure 2. Sauvignon blanc grapevine yield (a and b), bunch number per vine (b and d) for twocane (a and c) and fourcane (b and d) pruned vines and mean berry number per bunch (e) and berry mass (f) for all treatments.
UB represents Upper Brancott, WW represents Western Wairau, SA represents Seaview Awatere and CR represents Central Rapaura in Marlborough region, GL represents Glasneven in Canterbury region. GL vineyard is trained with three canes per vine. The other four vineyards have both twocane and fourcane pruned vines. Only four years of data were obtained for the GL vineyards, while for the other four vineyards more than 12 years of observations were carried out for both treatments. Mean berry number per bunch and berry mass were combined for different treatments as no differences were found between treatments. Note: The bold black line in each box represents the median value for each vineyard. The middle “box” represents the middle 50 % of scores for the group, ranging from lower (25 percent, Q1) to upper (75 percent, Q3) quartile. Upper whisker represents the range to Q3 + 1.5 * IQR and lower whisker represents the range to Q1 – 1.5 * IQR where IQR equals to Q3 – Q1, the box length.
Figure 3. The correlation between mean bunch per vine of a fourcanepruned Sauvignon blanc vine with sum of rainfall, mean daily radiation (MJ day1), mean maximum daily temperature (°C), mean minimum daily temperature (°C) and mean daily temperature (°C) during the inflorescence initiation period.
The distribution of each variable is shown on the diagonal. To the right of the diagonal, the values of the correlation between each factor pair plus the significance level as stars are shown. Each significance level is associated with a symbol based on pvalues: 0.001 (***), 0.01 (**), 0.05 (*), 0.1 (.). To the left of the diagonal, the bivariate scatter plots with a fitted loess line are displayed. The first row and column show the correlation between bunch number per vine and climatic factor in question. Plots were made with the R package of PerformanceAnalytics (Peterson et al., 2018) .
2. Bunch number
When tested with a single factor for fourcanepruned Sauvignon blanc vines, mean Tmax during the inflorescence initiation period (TmaxIni) gave the highest correlation with bunch number per vine (correlation index R^{ }= 0.77, Figure 3), followed by TmeanIni (R = 0.52), RadIni (R = 0.30), TminIni (R = 0.26) and RainTotIni (R = 0.18). Adding TminIni into the regression between TmaxIni and bunch number per vine did not improve the regression, indicating TmaxIni had a dominant effect on inflorescence initiation. Adding RadIni into the regression improved the overall R^{2} and loglikelihood under both two and fourcane conditions (Table S3). The interaction between TmaxIni, RadIni and RainTotIni were not significant. Thus TmaxIni + RadIni was chosen for the prediction of bunch number per vine.
Mean bunch number per vine increased linearly with the TmaxIni (Figure 4). On average, a one degree increase in temperature in TmaxIni was associated with a 2.87 bunch increase in twocane pruned vines and a 4.6 bunch increase in fourcane pruned vines. The optimised period for TmaxIni was 15.9 td before 50 % flowering until 1.27 td after 50 % flowering and the optimised period for RadIni was 10.4 td before 50 % flowering until 0.14 td after 50 % flowering (Table 1), indicating that the critical period affecting bunch number per vine was mainly before 50 % flowering.
Figure 4. The correlation between mean maximum temperature during inflorescence initiation and mean bunch number per vine at harvest for twocane (a) and fourcane (b) pruned Sauvignon blanc vines.
Lines are the linear regression for each vineyard region without random factors. The actual relationship could be nonlinear. UB represents Upper Brancott, WW represents Western Wairau, SA represents Seaview Awatere and CR represents Central Rapaura in the Marlborough region. Note: SA is in the Awatere Valley, which is significantly cooler and windier than the other three sites.
Table 1. Selected prediction model for different yield components and the optimised periods for the factors that were included in the selected model
Yield component 
Factors^{1} 
TD backward^{2} 
TD forward 
Fixed effects^{3} 
Random effects 

Bunch number per vine (#) 
TmaxIni 
15.90 
1.27 
2Cane: 34.48+2.87*TmaxIni+0.58*RadIni 
0 
RadIni 
10.42 
0.14 
4Cane: 43.43+4.60*TmaxIni+0.42*RadIni 
UB = 1.06; WW = 0.38; SA = 0.26; CP = 0.42 

Berry number per bunch (#) 
TmeanFlow 
7.08 
0.02 
44.3+3.18*TmeanFlow0.83*RainTotFlow+2.58*TmaxIni+ 0.047*TmeanFlow*RainTotFlow 
UB = 3.80; WW = 0.26; SA = 1.61; CR =3.69; GL=1.95 
RainTotFlow 
10.50 
2.69 

TmaxIni 
7.92 
30.37 

Berry mass (gram) 
TmeanFlow 
10.02 
1.17 
1.4 1.65e2*TmeanFlow 1.67e2 *RainTotFlow + 3.8e2 * RadFlow+2.33e2 * RainTotVer+9.87e4 * TmeanFlow*RainTotFlow – 8.98e4 * RadFlow*RainTotVer 
UB = 1e2; WW = 2e2; SA = 1e2; CR = 3e3; GL = 1e3 
RainTotFlow 
13.42 
2.92 

RadFlow 
11.37 
11.42 

RainTotVer 
24.78 
13.10 

Bunch mass (gram) 
TmeanFlow 
8.49 
0.17 
47.23+6.42*TmeanFlow2.88 * RainTotFlow 0.09* Et0_RainTotVer + 0.17 * TmeanFlow * RainTotFlow 
UB = 9.6; WW = 4.36; SA = 5.84 ; CR = 6.19; GL = 1.92 
RainTotFlow 
10.50 
2.72 

Et0_RainTotVer 
27.60 
13.58 

Yield (kg/per vine) 
TmaxIni 
8.13 
12.78 
2Cane: 13.55 + 0.2*TmaxIni + 0.34* RadIni + 0.33 * TmaxFlow 9.6e3 * Et0_RainTotVer 
0 
RadIni 
8.12 
12.77 

TmaxFlow 
8.49 
2.46 
4Cane: 28.24 + 0.37*TmaxIni + 0.30* RadIni + 0.92 * TmaxFlow + 4.6e3* Et0_RainTotVer 
UB = 0.07; WW = 0.02; SA = 0.01; CR = 0.06; 

Et0_RainTotVer 
27.83 
13.57 
^{1}Ini in each factor denotes the period during inflorescence initiation. The calculation for Ini uses flowering time of the previous season as the reference point. Flow denotes the period during the flowering time of the current season, using flowering time as the reference point. Ver denotes the post flowering and prevéraison period of the current season, using time of véraison as the reference point.
^{2}TD backward refers to the thermal days before the reference point. TD forward refers to the thermal days after the reference point.
^{3}When applying the fitted equations, the biological limits of bunch number per bud, berry number per bunch and berry mass need to be considered.
3. Berry number per bunch
Tmean around the flowering period (TmeanFlow) gave the highest correlation with berry number per bunch (R = 0.74, Figure 5a and Figure S4) when only one factor was considered, followed by TmaxFlow (R^{ }= 0.71), TminFlow (R^{ }= 0.48), RainTotFlow (R = 0.47) and TmaxIni (R^{ }= 0.29) and RadFlow (R = 0.22). Combining TmaxIni or RainTotFlow with Tmean improved the overall regression (Table S4). The best model was TmeanFlow * RainTotFlow + TmaxIni, which had the lowest AIC value and an overall R^{2} of 0.75. Berry number per bunch decreased with the amounts of cumulative rainfall around the flowering period (Figure 5b). However, Tmean and RainTotFlow had a positive interaction on berry number.
The period that gave the highest correlation between berry number per bunch and Tmean was from 7.08 td before 50 % flowering to 0.02 td after 50 % flowering, while the highest correlation with RainTot was for the period from 10.5 td before 50 % flowering to 2.69 td after 50 % flowering (Table 1). A close check of the response of R^{2} to changes in the value of forward td showed that R^{2 }reached the peak when forward td was around 0 and then decreased when forward td further increased (Figure S3), indicating that berry number per bunch was more sensitive to temperature conditions before 50 % flowering.
Figure 5. The correlation between mean temperature (a) and total rainfall (b) around the flowering period with mean berry number per bunch of Sauvignon blanc at harvest.
Lines are the linear regression for each vineyard without random factors. The actual relationship could be nonlinear. UB represents Upper Brancott, WW represents Western Wairau, SA represents Seaview Awatere and CR represents Central Rapaura in the Marlborough region, GL represents Glasneven in the Canterbury region.
4. Berry mass
TmeanFlow and RainTotVer had strong effects on berry mass (Figure 6 and Figure S5), followed by TmaxFlow (R^{ }= 0.55) and RainTotFlow (R = 0.55), Et0_RainVer (R = 0.48), RadFlow (R = 0.44) and Tmin (R = 0.34). The best model for predicting berry mass was TmeanFlow * RainTotFlow + RadFlow * RainTotVer (R^{2} = 0.68) with all the interactions being significant.
RainTotFlow had a strong negative effect on berry mass, while RainTotVer had a strong positive effect, indicating the sensitivity of berry mass to changes in rainfall events during different development stages. The period that gave the highest negative correlation between RainTotFlow and berry mass was from 13.4 td before 50 % flowering until 2.9 td after 50 % flowering, while the highest positive correlation for RainTotVer was from 24.8 td to 13.1 before véraison (Table 1). On average, there were 40 td from 50 % flowering to véraison.
Figure 6. The correlation between mean daily temperature around flowering period (a) and total rainfall between flowering and véraison (b) with mean berry mass of Sauvignon blanc at harvest.
Lines are the linear regression for each vineyard without random factors. The actual relationship could be nonlinear. UB represents Upper Brancott, WW represents Western Wairau, SA represents Seaview Awatere, CR represents Central Rapaura in the Marlborough region and GL represents Glasneven in the Canterbury region.
Figure 7. The correlation between mean daily temperature around flowering period (a) and total rainfall around flowering (b) with mean bunch mass of Sauvignon blanc at harvest.
Lines are linear regression for each vineyard without random factors. The actual relationship could be nonlinear. UB represents Upper Brancott, WW represents Western Wairau, SA represents Seaview Awatere, CR represents Central Rapaura in Marlborough region and GL represents Glasneven in Canterbury region.
5. Bunch mass
Similar to berry mass, temperature, radiation and post flowering rainfall had positive effects on bunch mass, while rainfall around flowering showed negative effects (Figure 7 and Figure S6). However, the correlation between bunch mass and TmeanFlow (R = 0.78) was much stronger than for berry mass (R = 0.59). The best model for predicting bunch mass was TmeanFlow * RainTotFlow + ET0_RainTotVer with an overall R^{2} of 0.83.
6. Yield per vine
Among all the factors, TmaxFlow (R = 0.71) and TmaxIni (0.56) stood out with the highest correlation with yield per vine (Figure S7). Despite the fact that all the vines were irrigated, yield per vine still negatively correlated with potential water deficit (ET0 minus RainTotVer). RadIni also had a marginal positive contribution on yield per vine. These four factors constitute the final model for predicting yield per vine.
The final yield  calculated by multiplying the estimated bunch number per vine, berry number per bunch and berry mass determined by weather conditions at critical periods  corresponded well with the observed yield for both twocane and fourcane pruned vines (Figure 8). The slope between predicted yield and observed yield was 0.9 (estimated by linear mixed effects model) and explained 80 percent of the total variance. Similar results for R^{2} were obtained by the direct yield estimation (Figure S8, R^{2} = 0.81) and by bunch number times bunch mass (Figure S9, R^{2} = 0.79).
Figure 8. The relationship between seasonal observed yield per vine and predicted yields (calculated from the predicted bunch number per vine, berry number per bunch and berry mass). Individual vineyards
Orange symbols are twocane and green fourcane pruned vines. The slope of the linear mixed regression without intercept was 0.9 and R^{2} was 0.8.
Discussion
Using data from a longterm yield monitoring experiment with meteorology data, this study quantified the relationship between grapevine yield components (bunch number per vine, berry number per bunch, berry mass, bunch mass and yield per vine) and weather conditions during critical periods of grapevine development. Among all the weather factors, temperature was shown to have the strongest effects on all yield components. Rainfall near flowering time proved to have a negative effect, while post flowering rainfall had positive effects on berry mass, bunch mass and overall yield. Radiation had a moderate effect under our experimental conditions. We further show that weather conditions before 50 % flowering have stronger effects than post flowering weather conditions on berry number per bunch of the current season and bunch number per vine in the following season (Table 1). For instance, the optimised critical periods of temperature for berry number, berry mass and bunch mass all mainly occurred before 50 % flowering of the current season and for bunch number they mainly occurred before 50 % flowering of the previous season.
1. Maximum daily temperature
Our data indicated that maximum daily temperature had a dominant effect on bunch number and overall yield (Figure 1 and Table 2) and that it was one of the most influential factors regarding berry number and bunch mass, although it was sometimes surpassed by mean temperature. Buttrose (1969a) showed that the fruitfulness of Vitis vinifera L. Muscat Gordo Blanco buds was related to maximum temperature (rather than to temperature summation) when maintained for about four hours under growth chamber conditions. This was probably because high temperature stimulates cytokinin production and encourages greater inflorescence primordia development, which happens within a few hours (Buttrose, 1970). In contrast, low temperature led to a production of gibberellins and increased tendril formation (Jackson, 2008; Mullins et al., 1992).
Regarding berry number and berry mass, the correlation index of TmeanFlow slightly surpassed that of TmaxFlow, indicating that minimum temperature may also play a role (Figure S4 and S5). However, when considering the whole yield, the correlation of both TmaxFlow and TmaxIni was higher than that of TmeanFlow, revealing the importance of Tmax in the overall yield formation. A positive effect of TmaxIni was also found on berry number per bunch. This was likely due to the positive effects of temperature on primary branch initiation prior to buds entering dormancy in the previous season, while primary branching was strongly correlated with flower number per inflorescence (Dunn and Martin, 2007).
2. The critical periods
Weather events during critical developmental periods that affect bunch number and berry number per bunch have a strong influence on yield. For grapevine, flowering and inflorescence initiation are critical periods, as weather conditions during these periods not only affect the current season’s berry number and berry mass, but also greatly affect the following season’s bunch number per vine.
Flower development, which determines the number of bunches (inflorescences) and berries in grapevine, involves three main steps: (1) formation of anlagen or uncommitted primordia, (2) differentiation of anlagen when forming inflorescence or tendril primordia and (3) differentiation of flowers. For the number of bunches per vine, our estimated critical period for TmaxIni was 15.9 td before 50 % flowering until 1.27 td after 50 % flowering. The start of the critical period is in agreement with the findings of Bennett (2002) and Swanepoel and Archer (1988), who showed that induction and initiation of anlagen generally start approximately 20 days before 50 % flowering at the basal two nodes on Chardonnay and Chenin blanc. The initiation of the first anlage continues acropetally on all nodes and takes about two days at one node for Chenin blanc in South Africa (Swanepoel and Archer, 1988). The initiation and differentiation of second anlage starts when the differentiation of the first anlage is completed. This whole process continues on all nodes till véraison, which is about 65 days after 50 % flowering (Bennett, 2002; Swanepoel and Archer, 1988).
The end of our estimated critical period (1.27 td after 50 % flowering) roughly corresponds to the end of initiation of the first anlage at the 10th node according to the diagram in Vasconcelos et al. (2009). Ten being the number of nodes we laid down during winter pruning. The end of the estimated critical period was much earlier than the end of second anlage differentiation process which was ~45 days after 50 % flowering at the 10th node (Vasconcelos et al., 2009). We suspect that the period when the initiation of first anlage occurs is the most sensitive period for determining the bunch number per vine and weather conditions during this period will have more impact than the weather conditions which follow it. We further hypothesised that weather conditions during the initiation phase could greatly influence whether the anlage or uncommitted primordia will become inflorescence or tendril. The success of the initiation of the first anlage and its node position on the apical primordium could have carryon effects on the initiation and differentiation of subsequent anlage. However, the fact that the fact that the initiation of the second anlage happens in warmer conditions (middle of summer) in our climate could have reduced its sensitivity and thus was not picked up by the optimization procedure.
For berry number per bunch, we found the critical period for TmeanFlow was 7.08 td before 50 % flowering till 0.02 td after 50 % flowering. This supports the findings by Ebadi et al. (1996), who found fruit set was reduced by exposing the vine to a lower temperature regime for one week in Chardonnay and Shiraz at EL 15 (when inflorescence comprised single flowers in compact groups) or EL 17 (when inflorescence comprised separate single flowers) (Coombe, 1995b). They further examined the seed characteristics of Chardonnay and found total number of seeds per berry was not affected, but lower temperature resulted in greater proportions of floater seeds. They further concluded that low temperature before flowering would negatively affect the structure of the ovules and the function of the pollen, seed development and reduce the number of cells and mass of the pericarp (Ebadi et al., 1996; Ebadi et al., 1995).
3. Rainfall around flowering and after flowering
Rainfall near flowering time was found to have a negative effect on berry number and berry mass, while rainfall after flowering was found to have a positive effect on berry mass and overall yield. Rain during the flowering period can physically inhibit pollination and fertilisation (Mullins et al., 1992). Rain can also result in the failure of flower caps to be shed which will reduce fertilisation and thus fruit set (Keller, 2015). The negative effects of rain before flowering on berry number and berry mass were also confounded by reduction in radiation and temperature on rainy days (Figure 3); in overcast conditions , leaf photosynthesis will not proceed at optimal rates which could limit carbohydrate supply to the inflorescence (Friend, 2005). After flowering and fruit set, good soil moisture conditions can promote water uptake and berry development (Ojeda et al., 2001; Pagay et al., 2015). The positive effects of rainfall on berry and bunch mass after flowering for trickle irrigated vines as in this study (Table 1) suggest the vines could still experience some water stress. This could be because the inter row still has low soil water content under trickle irrigation. For nonirrigation practices with waterstressed plants, berry and bunch mass were generally more affected by rainfall.
A big variation in berry mass under conditions with low cumulative rainfall around flowering was found (Figure 6). This could be due to: 1) RainTotFlow not being the only factor determining berry mass; a strong interaction between TmeanFlow and RainTotFlow was also found and high temperature with high rainfall was certainly less harmful than the combination of low temperature and high rainfall, 2) a variation between seasons in the amount of irrigation applied and soil water conditions prior to the calculation of cumulated rainfall and 3) the distribution of rain events: continuous light rain would have more negative effects on berry mass than short and heavy bursts of rain.
4. Limitations of the current study
Radiation was found to only have a moderate effect on bunch number, berry number and final yield. This could be because 1) the variation of radiation intensity at the bud level may have been low due to the dense canopy, despite the changeable overall exterior radiation and 2) the overall radiation intensity under our field conditions was relatively high and was thus not the main limiting factor. Unfortunately, we did not record the pruning weight and cane diameter in each season and we were therefore unable to link the variation in radiation intensity during each season with the actual biomass production, which is an indication of plant carbon status.
Bunches from noncounted buds (quiescent buds) were not separated from the counted buds on canes at harvest. Bunches from quiescent buds account for approximately 15 % of total bunches (unpublished data). This may have added noise in the correlations between bunch number per vine and climate variables due to the fact that these bunches developed fully during the year of harvest. In addition, due to limited data on flower number, we could not quantify the effects of temperature before budburst on flower number and the effects of weather conditions on fruit set, although this could be inferred by berry number and berry mass.
We show that an increase of one degree in both TmaxIni and TmaxFlow is associated with an increase of 0.53 kg per vine in twocane pruned vines and an increase of 1.29 kg per vine in fourcane pruned vines, assuming other factors remain the same (Table 1). Extrapolation beyond the temperature range found in the current study is not warranted. For instance, no clear yield trends were found on Shiraz in Barossa Valley Australia by increasing daytime ambient temperature (1.8 to 4.1 ^{o}C) for 2 to 3 weeks during a single phenological window (Sadras and Soar, 2009; Victor et al., 2017). This was likely linked to the background temperature. The mean daily maximum temperature in Barossa Valley in January was 30.5 °C (source: weatherzone.com.au ), which was about 8 °C higher than Marlborough (Chappell, 2016).
Furthermore, direct application of the positive relationship between temperature and yield when evaluating the effects of global warming on yield is not encouraged. Warming will likely advance phenology in such a way that temperature conditions during the key periods for initiation and fruit set may not be much different to those being currently experienced. Such an advance in phenology could also occur with a reduction in radiation intensity at flowering, as flowering currently occurs about 10 days before the summer solstice in Marlborough conditions. Thus, a following step would be to integrate the effects of weather conditions on bunch number per vine, berry number per bunch and berry mass into a processedbased plant model for assessing any changes in phenology, as well as the effects of carbon assimilation on berry sugar accumulation, in order to evaluate the effects of climate change on grapevine yield.
Conclusion
We quantified the correlation between grapevine yield components and weather conditions during key developmental stages (e.g., flowering) by carrying out a longterm phenology and yield monitoring trial. We found daily maximum temperature played a critical role in inflorescence initiation, while both daily maximum and minimum temperature played essential roles in berry number and berry mass. Radiation and rainfall account for extra variation in yield components besides temperature. Incorporating the correlations between yield components and weather conditions into plant models will likely improve our yield prediction for grapevine.
Supplementary information
The following Additional Supplementary Data can be found in the online version of this article on the publisher’s website:
Method S1 Vineyard monitoring history
Table S1. The optimised critical periods of each weather factor for bunch number per plant, berry number per bunch, mean berry mass, bunch mass and yield per vine
Table S2. The potential bias and standard error of parameter estimated by the mixed linear model caused by the year and site as evaluated by the bootstrap method
Table S3. List of all models tested for bunch number per plant and weather conditions
Table S4. List of all models tested for berry number per bunch and weather conditions
Table S5. List of all models tested for berry mass and weather conditions
Table S6. List of all models tested between bunch mass and weather conditions
Table S7. List of all models tested between yield per vine and weather conditions
Figure S1. Reproductive sequence of different grapevine yield components and the potential influence of weather conditions on each yield component at different development stages
Figure S2 Illustration of the analysing procedure.
Figure S3. The change of R2 values of the linear regression between mean daily maximum temperature during flowering and berry number per bunch of a Sauvignon blanc vine in response to the parameters that defines the period for calculating the mean daily maximum temperature
Figure S4. The correlation between mean berry number per bunch with different weather factors
Figure S5. The correlation between mean berry mass with the different weather factors
Figure S6. The correlation between mean bunch mass with different weather factors
Figure S7. The correlation between mean yield per vine with different weather factors
Figure S8. The relationship between seasonal observed yield per vine and predicted yields (calculated directly by weather conditions around critical period).
Figure S9. The relationship between seasonal observed yield per vine and predicted yields yield (calculated from the predicted bunch number per plant × bunch mass).
Acknowledgements
The authors would like to thank all the staff at the Plant and Food Research (PFR) Marlborough site who helped maintain and conduct this longterm trial. Special thanks go to Rafidah Horner, Franziska Grab, Trevor Skilton and Rachel Bishell who conducted most of the phenology monitoring and berry sampling in Marlborough and to Tim Parker and Julian Smith in Waipara. The authors also thank Pernod Ricard New Zealand Limited and Delegat Wine Estate for allowing the monitoring to take place on their vineyards and Warrick Nelson and the PFR internal science publication office for their help in revising the manuscript. Over the past 15 years, this work was funded by the New Zealand Ministry of Business, Innovation and Employment (MBIE) Quality of New Zealand Wines programme contracted to the University of Auckland UOAX0404 (20042010), Marlborough Research Centre – Phenology Monitoring Project (20102014), New Zealand Winegrowers/Sustainable Farming Fund (20142017), New Zealand Winegrowers/Sustainable Winegrowing New Zealand (July 2017 onwards). This work was also part of the PFR Grape and Wine Research programme, funded by the MBIE Strategic Science Investment Fund.
References
 Ardia D., Boudt K., Carl P., Mullen K. and Peterson B.G., 2011. Differential evolution with DEoptim: an application to nonconvex portfolio optimization. The R Journal, 3(1), 2734. doi:10.32614/RJ2011005
 Bartoń K., 2018. MuMIn: MultiModel Inference R Package Version 0.12.2/r18. In. http://RForge.Rproject.org/projects/mumin/.
 Bennett J.S., 2002. Relationships between carbohydrate supply and reserves and the reproductive growth of grapevines (Vitis vinifera L.). Lincoln University,
 Bock A., Sparks, T. H., Estrella, N. and Menzel A., 2013. Climateinduced changes in grapevine yield and must sugar content in Franconia (Germany) between 1805 and 2010. PLoS ONE, 8(7), e69015. doi:10.1371/journal.pone.0069015
 Buttrose M.S., 1969a. Fruitfulness in grapevines: Effects of change in temperature and light regimes. Botanical Gazette, 130, 173179. doi:10.1086/336487
 Buttrose M.S., 1969b. Fruitfulness in grapevines: Effects of light intensity and temperature. Botanical Gazette, 130(3), 166  173. doi:10.1086/336486
 Buttrose M.S., 1970. Fruitfulness in grapevines: Development of leaf primordia in buds in relation to bud fruitfulness. Botanical Gazette, 131, 7883. doi:10.1086/336515
 Buttrose M.S., 1974. Climatic factors and fruitfulness in grapevines. Horticultural Abstracts, 44, 319325.
 Buttrose M.S. and Hale C.R., 1973. Effect of temperature on development of the grapevine inflorescence after bud burst. American Journal of Enology and Viticulture, 24(1), 1416.
 Chappell P.R., 2016. The climate and weather of Marlborough. NIWA Science and Technology Series, 69, 40.
 Coombe, B. G., 1962. The effect of removing leaves, flowers and shoot tips on fruitset in Vitis vinifera L. J Hort Sci, 37, 1  15. doi:10.1080/00221589.1962.11514023
 Coombe B.G., 1995a. Adoption of a system for identifying grapevine growth stages. Australian Journal of Grape and Wine Research, 1(2), 104110. doi:10.1111/j.17550238.1995.tb00086.x
 Coombe B.G., 1995b. Growth stages of the grapevine: adoption of a system for identifying grapevine growth stages. Australian Journal of Grape and Wine Research, 1, 104  110. doi:10.1111/j.17550238.1995.tb00086.x
 Dunn G.M. and Martin S.R., 2007. A functional association in Vitis vinifera L. cv. Cabernet Sauvignon between the extent of primary branching and the number of flowers formed per inflorescence. Australian Journal of Grape and Wine Research, 13(2), 95100. doi:10.1111/j.17550238.2007.tb00239.x
 Ebadi A., May P. and Coombe B.G., 1996. Effect of shortterm temperature and shading on fruitset, seed and berry development in model vines of V. vinifera, cvs Chardonnay and Shiraz. Australian Journal of Grape and Wine Research, 2(1), 18. doi:10.1111/j.17550238.1996.tb00087.x
 Ebadi A., May P., Sedgley M. and Coombe B.G., 1995. Effect of low temperature near flowering time on ovule development and pollen tube growth in the grapevine (Vitis vinifera L.), cvs Chardonnay and Shiraz. Australian Journal of Grape and Wine Research, 1(1), 1118. doi: doi:10.1111/j.17550238.1995.tb00072.x
 Eltom M., Trought M.C., Agnew R., Parker A. and Winefield C.S., 2017. Pre‐budburst temperature influences the inner and outer arm morphology, phenology, flower number, fruitset, TSS accumulation and variability of Vitis vinifera L. Sauvignon Blanc bunches. Australian Journal of Grape and Wine Research, 23(2), 280286. doi:10.1111/ajgw.12260
 Faluomi V., 2017. Mathematical modelling of vine and grapevine development to evaluate the impact of climate change on wine quality. Paper presented at the Fifth International Congress on Mountain and steep slope viticulture Conegliano, Treviso  Veneto, Italy.
 Fraga H., Garcia de Cortazar Atauri I., Malheiro A.C. and Santos J.A., 2016. Modelling climate change impacts on viticultural yield, phenology and stress conditions in Europe. Glob Chang Biol, 22(11), 37743788. doi:10.1111/gcb.13382
 Friend A., 2005. Berry set and development in Vitis vinifera L. (PhD), PhD thesis, Lincoln University, Lincoln, New Zealand.
 Friend A.P., Trought M.C.T. and Creasy G.L., 2009. The influence of seed weight on the development and growth of berries and live green ovaries in Vitis vinifera L. cvs. Pinot noir and CabernetSauvignon. Australian Journal of Grape and Wine Research, 15(2), 166174. doi:10.1111/j.17550238.2009.00050.x
 García de CortázarAtauri I., Daux V., Garnier E., Yiou P., Viovy, N., Seguin B., . . . Chuine I., 2010. Climate reconstructions from grape harvest dates: Methodology and uncertainties. The Holocene, 20(4), 599608. doi:10.1177/0959683609356585
 Hui F.K.C., Müller S. and Welsh A.H., 2019. Testing random effects in linear mixed models: another look at the Ftest (with discussion). Australian and New Zealand Journal of Statistics, 61(1), 6184. doi:10.1111/anzs.12256
 Jackson D.I. and Cherry N.J., 1988. Prediction of a District's GrapeRipening Capacity Using a LatitudeTemperature Index (LTI). American Journal of Enology and Viticulture, 39(1), 1928.
 Jackson R.S., 2008. Wine science: principles and applications: Academic press.
 Jon D.A., Gregory V.J. Andrew T., Andrew H. and Michael C.T.T., 2012. Analysis of viticulture region climate structure and suitability in New Zealand. Oeno One, 46(3). doi:10.20870/oenoone.2012.46.3.1515
 Keller M., 2015. The Science of Grapevines: Anatomy and Physiology: Elsevier Science.
 Kliewer W.M., 1977. Influence of temperature, solar radiation and nitrogen on coloration and composition of Emperor grapes. American Journal of Enology and Viticulture, 28(2), 96103.
 Kliewer W.M. and Dokoozlian N.K., 2005. Leaf area/crop weight ratios of grapevines: influence on fruit composition and wine quality. American Journal of Enology and Viticulture, 56(2), 170181.
 Kuznetsova A., Brockhoff P.B. and Christensen R.H.B., 2017. lmerTest package: tests in linear mixed effects models. Journal of Statistical Software, 82(13). doi:10.18637/jss.v082.i13
 May P., 2000. From bud to berry, with special reference to inflorescence and bunch morphology in Vitis vinifera L. Australian Journal of Grape and Wine Research, 6(2), 8298. doi:10.1111/j.17550238.2000.tb00166.x
 Mullins M.G., Bouquet A. and Williams L.E., 1992. Biology of the grapevine: Cambridge University Press.
 New Zealand Winegrowers., 2018. NZ Winegrowers Annual Report 2018. Auckland New Zealand: New Zealand Winegrowers.
 Ojeda H., Deloire A. and Carbonneau A., 2001. Influence of water deficits on grape berry growth. Vitis, 40(3), 141145.
 Pagay V., Zufferey V. and Lakso A.N., 2015. The influence of water stress on grapevine (Vitis vinifera L.) shoots in a cool , humid climate : growth , gas exchange and hydraulics. Functional Plant Biology, 43, 827837. doi:10.1071/FP16017
 Parker A.K., Hofmann R.W., van Leeuwen C., McLachlan A.R.G. and Trought M.C.T., 2014. Leaf area to fruit mass ratio determines the time of veraison in Sauvignon Blanc and Pinot Noir grapevines. Australian Journal of Grape and Wine Research, 20(3), 422431. doi:10.1111/ajgw.12092
 Perold A.I., 1927. Treatise on viticulture. London: Macmillan and Co.
 Peterson B.G., Carl P., Boudt K., Bennett R., Ulrich J., Zivot E., . . . Balkissoon K., 2018. Package ‘PerformanceAnalytics’. R Team Cooperation.
 Petrie P.R., Trought M.C.T. and Howell G.S., 2000. Fruit composition and ripening of Pinot Noir (Vitis vinifera L.) in relation to leaf area. Australian Journal of Grape and Wine Research, 6, 4651. doi:10.1111/j.17550238.2000.tb00161.x
 R Development Core Team R., 2019. R: A Language and Environment for Statistical Computing (Version 3.6.0. Retrieved from http://www.rproject.org
 Sadras V.O. and Soar C.J., 2009. Shiraz vines maintain yield in response to a 24 °C increase in maximum temperature using an opentop heating system at key phenostages. European Journal of Agronomy, 31(4), 250258. doi:10.1016/j.eja.2009.09.004
 Sakamoto Y., Ishiguro M. and Kitagawa G., 1986. Akaike information criterion statistics. Dordrecht, The Netherlands: D. Reidel, 81.
 Sanchez L.A. and Dokoozlian N.K., 2005. Bud microclimate and fruitfulness in Vitis vinifera L. American Journal of Enology and Viticulture, 56(4), 319329.
 Santos J.A., Malheiro A.C., Karremann M.K. and Pinto J. G., 2011. Statistical modelling of grapevine yield in the Port Wine region under present and future climate conditions. International Journal of Biometeorology, 55(2), 119131. doi:10.1007/s0048401003180
 Staudt G., 1982. Pollen germination and pollen tube growth in vivo with Vitis [grapes] and their dependence on temperature. Vitis.
 Swanepoel J.J. and Archer E., 1988. The ontogeny and development of Vitis vinifera L. cv. Chenin blanc inflorescence in relation to phonological stages. . Vitis, 27, 133141.
 Trought M.C.T., 2000. Hitting an ndimensional target in moving space: The aim of the viticulturist and the role of the grapevine physiologist. Paper presented at the 5th International Symposium on Cool Climate Viticulture and Oenology, Melbourne, Australia.
 Trought M.C.T., 2005. Fruitset  possible implications on wine quality. Paper presented at the Transforming flowers to fruit, Mildura, Australia.
 Vasconcelos M.C., Greven M., Winefield C.S., Trought M.C.T. and Raw V., 2009. The Flowering Process of Vitis vinifera: A Review. American Journal of Enology and Viticulture, 60(4), 411434.
 Victor S., Martin M. and Paul P., 2017. Resilience of grapevine yield in response to warming. Oeno One, 51(4. doi:10.20870/oenoone.2017.51.4.1913
 Watt A., 2010. Environmental Factors Influencing Inflorescence Differentation and Development and Bunch Architecture, of Vitis vinifera L. cvs. Chardonnay, Shiraz and Sauvignon Blanc. University of Melbourne, Melbourne School of Land and Environment,
 Watt A.M., Dunn G.M., May P.B., Crawford S.A. and Barlow E.W.R., 2008. Development of inflorescence primordia in Vitis vinifera L. cv. Chardonnay from hot and cool climates. Australian Journal of Grape and Wine Research, 14(1), 4653. doi:10.1111/j.17550238.2008.00006.x
 Wilson G., 1996. The influence of site environment and the effects of varying light and temperature on inflorescence development and flowering in grapevines, Vitis vinifera L. CabernetSauvignon. (MSc), MSc thesis, Lincoln University.