IVAS 2022 banner
IVES 9 IVES Conference Series 9 IVAS 9 IVAS 2022 9 Use of mathematical modelling and multivariate statistical process control during alcoholic fermentation of red wine

Use of mathematical modelling and multivariate statistical process control during alcoholic fermentation of red wine

Abstract

Cyberphysical systems can be seen in the wine industry in the form of precision oenology. Currently, limitations exist with established infrared chemometric models and first principle mathematical models in that they require a high degree of sample preparation, making it inappropriate for use in-line, or that few oenological parameters are considered. To our knowledge, a system which incorporates a more comprehensive mathematical model as well as in-line spectroscopic monitoring for the purpose of precision oenology has not yet been presented.

The use of first principle mathematical modelling was employed to predict the trends of alcoholic fermentation and oenological parameters in a four-phase model based on initial conditions. The components of interest were sugars, alcohol, biomass, nitrogen, carbon dioxide, phenolic parameters, and pH. The phases considered included the lees, the cap, the must, and an intermediate liquid phase present in the cap. For each phase, a system of ordinary differential equations was developed to describe the change of each of the components listed. Parameters such as mass transfer coefficients and partition coefficients need to be determined via regression during the model development stage. To obtain the necessary data, fermentations using three different cultivars (Shiraz, Merlot, and Cabernet Sauvignon) were conducted using three different temperatures (20oC, 25oC, and 28oC). Samples were taken once per day and chemical analysis took place for each of the components. A functional mathematical model capable of generating accurate forecasts for different oenological components using the chemical composition of grapes was attempted. Additionally, the model should describe the change in parameters due to cap mixing and increasing ethanol concentration. The model includes the boundary conditions which can be used to determine if a fermentation is deviating from desired progression.

To complete this process control system, it is still necessary to utilize partial least squares (PLS) calibration models for real time monitoring. Additionally, outlier identification, caused by abnormal spectra, was performed using statistical analysis allowing samples to be re-analysed. The use of machine learning techniques and the development of local and incremental models was explored to assess a live updating of the PLS models. The expected outcome of this study is a combined system using dynamic modelling to predict the fermentation and extraction trends and the monitoring with real time predictions generated by PLS models

DOI:

Publication date: June 23, 2022

Issue: IVAS 2022

Type: Article

Authors

Lambrecht Kiera Nareece¹, Du Toit Prof. W.J.¹, Louw Prof. T.M.²and Aleixandre Tudo Dr. J.L.¹,³

¹Stellenbosch University, South African Grape and Wine Research Institute, Department of Viticulture and Oenology
²Stellenbosch University, Department of Process Engineering
³Universitat Politecnica de Valencia, Instituto de Ingenieria de Alimentos para el Desarrollo (IIAD), Departamento de Tecnología de Alimentos

Contact the author

Keywords

In-line monitoring, process control, dynamic modelling, chemometrics, live modelling

Tags

IVAS 2022 | IVES Conference Series

Citation

Related articles…

Ecophysiological performance of Vitis rootstocks under water stress

The use of rootstocks tolerant to soil water deficit is an interesting strategy to cope with limited water availability. Currently, several nurseries are breeding new genotypes, but the physiological basis of its responses under water stress are largely unknown. To this end, an ecophysiological assessment of the conventional 110-Richter (110R) and SO4, and the new M1 and M4 rootstocks was carried out in potted ungrafted plants. During one season, these Vitis genotypes were grown under greenhouse conditions and subjected to two water regimes, well-watered and water deficit. Water potentials of plants under water deficit down to < -1.4 MPa, and net photosynthesis (AN) <5 μmol m-2 s-1 did not cause leaf oxidative stress damage compared to well-watered conditions in any of the genotypes. The antioxidant capacity was sufficient to neutralize the mild oxidative stress suffered. Under both treatments, gravimetric differences in daily water use were observed among genotypes, leading to differences in the biomass of root, shoot and leaf. Under well-watered conditions, SO4 and 110R were the most vigorous and M1 and M4 the least. However, under water stress, SO4 exhibited the greatest reduction in biomass while M4 showed the lowest. Remarkably, under these conditions, SO4 reached the least negative stem water potential (Ψstem), while M1 reduced stomatal conductance (gs) and AN the most. In addition, SO4 and M1 genotypes also showed the highest and lowest hydraulic conductance values, respectively. Our results suggest that there are differences in water use regulation among genotypes, not only attributed to differences in stomatal regulation or intrinsic water use efficiency at the leaf level. Therefore, because no differences in canopy-to-root ratio were achieved, it is hypothesized that xylem vessel anatomical differences may be driving the reported differences among rootstocks performance. Results demonstrate that each Vitis rootstock differs in its ecophysiological responses under water stress.

Using δ13C and hydroscapes as a tool for discriminating cultivar specific drought response

Measurement of carbon isotope discrimination in berry juice sugars at maturity (δ13C) provides an integrated assessment of water use efficiency (WUE) during the period of berry ripening, and when collected over multiple seasons can be used as an indication of drought stress response. Berry juice δ13C measurements were carried out on 48 different varieties planted in a common garden experiment in Bordeaux, France from 2014 through 2021 and were paired with midday and predawn leaf water potential measurements on the same vines in a subset of six varieties. The aim was to discriminate a large panel of varieties based on their stomatal behaviour and potentially identify hydraulic traits characterizing drought tolerance by comparing δ13C and hydroscapes (the visualisation of plant stomatal behaviour as a response to predawn water potential). Cluster analysis found that δ13C values are likely affected by the differing phenology of each variety, resulting in berry ripening of different varieties taking place under different stress conditions within the same year. We accounted for these phenological differences and found that cluster analysis based on specific δ13C metrics created a classification of varieties that corresponds well to our current empirical understanding of their relative drought tolerances. In addition, we analysed the water potential regulation of the subset of six varieties (using the hydroscape approach) and found that it was well correlated with some δ13C metrics. Surprisingly, a variety’s water potential regulation (specifically its minimum critical leaf water potential under water deficit) was strongly correlated to δ13C values under well-watered conditions, suggesting that base WUE may have a stronger impact on drought tolerance than WUE under water deficit. These results give strong insights on the innate WUE of a very large panel of varieties and suggest that studies of drought tolerance should include traits expressed under non-limiting conditions.

Assessing the relationship between cordon strangulation, dieback, and fungal trunk disease symptom expression

Grapevine trunk diseases including Eutypa dieback are a major factor in the decline of vineyards and may lead to loss of productivity, reduced income, and premature reworking or replanting. Several studies have yielded results indicating that vines may be more likely to express symptoms of vascular disease if their health is already compromised by stress. In Australia and many other wine-growing regions it is a common practice for canes to be wrapped tightly around the cordon wire during the establishment of permanent cordon arms. It is likely that this practice may have a negative effect on health and longevity, as older cordons that have been trained in this manner often display signs of decay and dieback, with the wire often visibly embedded within the wood of the cordon. It is possible that adopting a training method which avoids constriction of the vasculature of the cordon may help to limit the onset of vascular disease symptom expression. A survey was conducted during the spring of two consecutive growing seasons on vineyards in South Australia displaying symptoms of Eutypa lata infection when symptomless shoots were 50–100 cm long. Vines were assessed as follows: (i) the proportion of cordon exhibiting dieback was rated using a 0–100% scale; (ii) the proportion of canopy exhibiting foliar symptoms of Eutypa dieback was rated using a 0–100% scale; (iii) the severity of strangulation was rated using a 0–4 point scale. Images were also taken of each vine for the purpose of measuring plant area index (PAI) using the VitiCanopy App. The goal of the survey was to determine if and to what extent any correlation exists between severity of strangulation and cordon dieback, in addition to Eutypa dieback foliar symptom expression.

A better understanding of the climate effect on anthocyanin accumulation in grapes using a machine learning approach

The current climate changes are directly threatening the balance of the vineyard at harvest time. The maturation period of the grapes is shifted to the middle of the summer, at a time when radiation and air temperature are at their maximum. In this context, the implementation of corrective practices becomes problematic. Unfortunately, our knowledge of the climate effect on the quality of different grape varieties remains very incomplete to guide these choices. During the Innovine project, original experiments were carried out on Syrah to study the combined effects of normal or high air temperature and varying degrees of exposure of the berries to the sun. Berries subjected to these different conditions were sampled and analyzed throughout the maturation period. Several quality characteristics were determined, including anthocyanin content. The objective of the experiments was to investigate which climatic determinants were most important for anthocyanin accumulation in the berries. Temperature and irradiance data, observed over time with a very thin discretization step, are called functional data in statistics. We developed the procedure SpiceFP (Sparse and Structured Procedure to Identify Combined Effects of Functional Predictors) to explain the variations of a scalar response variable (a grape berry quality variable for example) by two or three functional predictors (as temperature and irradiance) in a context of joint influence of these predictors. Particular attention was paid to the interpretability of the results. Analysis of the data using SpiceFP identified a negative impact of morning combinations of low irradiance (lower than about 100 μmol m−2 s−1 or 45 μmol m−2 s−1 depending on the advanced-delayed state of the berries) and high temperature (higher than 25oC). A slight difference associated with overnight temperature occurred between these effects identified in the morning.

Comparison of imputation methods in long and varied phenological series. Application to the Conegliano dataset, including observations from 1964 over 400 grape varieties

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.