Integrated multiblock data analysis for improved understanding of grape maturity and vineyard site contributions to wine composition and sensory domains

A large varietal collection including over 1700 varieties was maintained in Conegliano, ITA, since the 1950s. Phenological data on a subset of 400 grape varieties including wine grapes, table grapes, and raisins were acquired at bud break, flowering, veraison, and ripening since 1964. Despite the efforts in maintaining and acquiring data over such an extensive collection, the data set has varying degrees of missing cases depending on the variety and the year. This is ubiquitous in phenology datasets with significant size and length. In this work, we evaluated four state-of-the-art methods to estimate missing values in this phenological series: k-Nearest Neighbour (kNN), Multivariate Imputation by Chained Equations (mice), MissForest, and Bidirectional Recurrent Imputation for Time Series (BRITS). For each phenological stage, we evaluated the performance of the methods in two ways. 1) On the full dataset, we randomly hold-out 10% of the true values for use as a test set and repeated the process 1000 times (Monte Carlo cross-validation). 2) On a reduced and almost complete subset of varieties, we varied the percentage of missing values from 10% to 70% by random deletion. In all cases, we evaluated the performance on the original values using normalized root mean squared error. For the full dataset we also obtained performance statistics by variety and by year. MissForest provided average errors of 17% (3 days) at budbreak, 14% (4 days) at flowering, 14.5% (7 days) at veraison, and 17% (3 days) at maturity. We completed the imputations of the Conegliano dataset, one of the world’s most extensive and varied phenological time series and a steppingstone for future climate change studies in grapes. The dataset is now ready for further analysis, and a rigorous evaluation of imputation errors is included.

In the study of terroir, a separate analysis of its many component factors can be of great help in accurately identifying a vineyard’s natural elements that impact wine quality and typicity. This research used a dedicated pluri-disciplinary approach to investigate the ecological characteristics, including geology and geographical features, of 14 vineyards that produce Gewürztraminer and Sauvignon Blanc cultivars in the alpine Alto Adige DOC wine region. Both the geopedological method using Vineyards Geological Identity (VGI) and the new Solar Radiaton Identity (SRI) topoclimatic classification method were used to provide analytical measurements and qualitative/quantitative characterisations. In addition, wide-ranging targeted and untargeted oenological and chemical analyses were carried out on grapes, musts and wines to correlate the soils’ geomineral and physical conditions with the biochemical properties of their fruits and wines. The research identified strong correlations between vineyard geo-identity and wine biofingerprint, confirming a mineral traceability of strontium rubidium ratio and some minerals distinctive to the local geology, such as K, Ca, Ag, Ba and Mn.  The study also discovered that particular geomineral and physical soil conditions of the studied vineyards are related to the different amount of amino acids, primary varietal aromas and polyphenols found in grapes, musts and wines. The research confirmed that winemaking technologies support oenological quality, although in some cases, human practices can overpower certain characteristic elements in wine, erasing the typical imprint left by the vineyards’ natural terroir, which becomes less traceable. Terroir abiotic ecological factors and vineyard identity can be classified in detail using the new VGI and SRI analysis methods to discover interrelationships between geo-pedological and topoclimatic conditions that impact wine quality. These methods are also helpful in identifying which ecological elements are exclusive to a particular vineyard or wine sub-region.

In wine growing regions around the world, climate change has the potential to affect vine transpiration and overall vineyard water use due to related changes in atmospheric demand and soil water deficits. Grapevines control their transpiration in response to a changing environment by regulating conductance of water through the soil-plant-atmosphere continuum. Most vineyard water use models currently estimate vine transpiration by applying generic crop coefficients to estimates of reference evapotranspiration, but this does not account for changes in vine conductance associated with water stress, nor differences thought to exist between varieties. The response of bulk stomatal conductance to daily weather variability and seasonal drought stress was studied on Cabernet-Sauvignon, Merlot, Tempranillo, Ugni blanc, and Semillon vines in a non-irrigated vineyard in Bordeaux France. Whole vine sap flow, temperature and humidity in the vine canopy, and net radiation absorbed by the vine canopy were measured on 15-minute intervals from early July through mid-September 2020, together with periodic measurement of leaf area, canopy porosity, and predawn leaf water potential. From this data, bulk stomatal conductance was calculated on 15-minute intervals, and multiple regression analysis was performed to identify key variables and their relative effect on conductance. Attention was focused on addressing multicollinearity and time-dependency in the explanatory variables and developing regression models that were readily interpretable. Variability of vapor pressure deficit over the day, and predawn water potential over the season explained much of the variability in conductance, with relative differences in response coefficients observed across the five varieties. By characterizing this conductance response, the dynamics of vine transpiration can be better parameterized in vineyard water use modeling of current and future climate scenarios.

Most European vineyards are managed under rainfed conditions, where seasonal water deficit has become increasingly important. The flowering-veraison phenophase represents an important period for vine response to water stress, which is seldomly thoroughly evaluated. Therefore, we aim to quantify the flowering-veraison water stress levels using Crop Water Stress Indicator (CWSI) over 1986–2015 for important European wine regions, and to assess the respective potential Yield Lose Rate (YLR). Additionally, we also investigate whether an advanced flowering-veraison phase may help alleviating the water stress with improved yield. A process-based grapevine model STICS is employed, which has been extensively calibrated for flowering and veraison stages using observed data at 38 locations with 10 different grapevine varieties. Subsequently, the model is being implemented at the regional level, considering site-specific calibration results and gridded climate and soil datasets. The findings suggest wine regions with stronger flowering-veraison CWSI tend to have higher potential YLR. However, contrasting patterns are found between wine regions in France-Germany-Luxembourg and Italy-Portugal-Spain. The former tends to have slight-to-moderate drought conditions (CWSI<0.5) and a negligible-to-moderate YLR (<30%), whereas the latter possesses severe-to-extreme CWSI (>0.5) and substantial YLR (>40%). Wine regions prone to a high drought risk (CWSI>0.75) are also identified, which are concentrated in southern Mediterranean Europe. An advanced flowering-veraison phase may have benefited from cooler temperatures and a higher fraction of spring precipitation in wine regions of Italy-Portugal-Spain, resulting in alleviated CWSI and moderate reductions of YLR. For those of France-Germany-Luxembourg, this can have reduced flowering-veraison precipitation, but prevalent alleviations of YLR are also found, possibly because of shifted phase towards a cooler growing season with reduced evaporative demands. Overall, such a retrospective analysis might provide new insights towards better management of seasonal water deficit for conventionally vulnerable Mediterranean wine regions, but also for relatively cooler and wetter Central European regions.

Current climate change is increasing inter- and intra-annual variability in atmospheric conditions leading to grapevine phenological shifts as well altered grape ripening and composition at ripeness. This study aims to (i) detect weather anomalies within a long-term time series, (ii) model grape ripening revealing altered traits in time to target specific ripeness thresholds for four Vitis vinifera cultivars, and (iii) establish empirical relationships between ripening and weather anomalies with forecasting purposes. The Day of the Year (DOY) to reach specific grape ripeness targets was determined from time series of sugar concentrations, total acidity and pH collected from a private company in the period 2009-2021 in North-Eastern Italy. Non-linear models for the DOY to reach the specified ripeness thresholds were assessed for model efficiency (EF) and error of prediction (RMSE) in four grapevine cultivars (Merlot, Cabernet Sauvignon, Glera and Garganega). For each vintage and cultivar, advances or delays in DOY to target specified ripeness thresholds were assessed with respect to the average ripening dynamics. Long-term meteorological series monitored at ground weather station by means of hourly air temperature and rainfall data were analyzed. Climate statistics were obtained and for each time period (month, bimester, quarter and year) weather anomalies were identified. A linear regression analysis was performed to assess a possible correlation that may exist between ripening and weather anomalies. For each cultivar, ripeness advances or delays expressed in number of days to target the specific ripening threshold were assessed in relation to registered weather anomalies and the specific reference time period in the vintage. Precipitation of the warmest month and spring quarter are key to understanding the effect of climate change on sugar ripeness. Minimum temperatures of May-June bimester and maximum temperatures of spring quarter best correlate with altered total acidity evolution and pH increment during the ripening process, respectively.