Next Article in Journal
Marker-Less 3d Object Recognition and 6d Pose Estimation for Homogeneous Textureless Objects: An RGB-D Approach
Next Article in Special Issue
Assessment of Smoke Contamination in Grapevine Berries and Taint in Wines Due to Bushfires Using a Low-Cost E-Nose and an Artificial Intelligence Approach
Previous Article in Journal
A Miniature Intermittent-Flow Respirometry System with a 3D-Printed, Palm-Sized Zebrafish Treadmill for Measuring Rest and Activity Metabolic Rates
Previous Article in Special Issue
Artificial Intelligence Applied to a Robotic Dairy Farm to Model Milk Productivity and Quality based on Cow Data and Daily Environmental Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Classification of Smoke Contaminated Cabernet Sauvignon Berries and Leaves Based on Chemical Fingerprinting and Machine Learning Algorithms

1
Digital Agriculture, Food, and Wine Group, Faculty of Veterinary and Agricultural Sciences, The University of Melbourne, Parkville, VIC 3010, Australia
2
School of Agriculture, Food and Wine, The University of Adelaide, Waite Campus, PMB 1, Glen Osmond, SA 5064, Australia
3
The Australian Research Council Training Centre for Innovative Wine Production, PMB 1, Glen Osmond, SA 5064, Australia
4
Department of Wine, Food and Molecular Biosciences, Faculty of Agriculture and Life Sciences, Lincoln University, Lincoln 7647, Canterbury, New Zealand
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(18), 5099; https://doi.org/10.3390/s20185099
Submission received: 22 August 2020 / Revised: 5 September 2020 / Accepted: 5 September 2020 / Published: 7 September 2020

Abstract

:
Wildfires are an increasing problem worldwide, with their number and intensity predicted to rise due to climate change. When fires occur close to vineyards, this can result in grapevine smoke contamination and, subsequently, the development of smoke taint in wine. Currently, there are no in-field detection systems that growers can use to assess whether their grapevines have been contaminated by smoke. This study evaluated the use of near-infrared (NIR) spectroscopy as a chemical fingerprinting tool, coupled with machine learning, to create a rapid, non-destructive in-field detection system for assessing grapevine smoke contamination. Two artificial neural network models were developed using grapevine leaf spectra (Model 1) and grape spectra (Model 2) as inputs, and smoke treatments as targets. Both models displayed high overall accuracies in classifying the spectral readings according to the smoking treatments (Model 1: 98.00%; Model 2: 97.40%). Ultraviolet to visible spectroscopy was also used to assess the physiological performance and senescence of leaves, and the degree of ripening and anthocyanin content of grapes. The results showed that chemical fingerprinting and machine learning might offer a rapid, in-field detection system for grapevine smoke contamination that will enable growers to make timely decisions following a bushfire event, e.g., avoiding harvest of heavily contaminated grapes for winemaking or assisting with a sample collection of grapes for chemical analysis of smoke taint markers.

1. Introduction

The incidence and intensity of wildfires are increasing worldwide, mainly due to the effects of climate change [1,2,3,4,5]. Bushfires that occur near wine regions can result in grapevine smoke exposure, which can alter the chemical composition of grape berries. Wine produced from these smoke-affected grapes may exhibit unpalatable smoky aromas and flavors, such as “burnt wood”, “ashy”, and “burnt rubber” [6,7,8,9]. These undesirable characters have been attributed to smoke-derived volatile phenols (VPs), including guaiacol, 4-methylguaiacol, cresols, and syringol [7,10,11]. It is thought that these VPs accumulate primarily in the skin of grape berries following smoke exposure and, to a lesser extent, in the pulp and seeds [12,13,14,15]. Grapevine smoke exposure, and the resulting smoke taint in wine, have caused significant financial losses for grape growers and winemakers due to discarded grapes and unsaleable wine. For example, the 2009 Black Saturday bushfires in Victoria, Australia, were estimated to have caused AUD 300 million in lost revenue [16,17,18,19]. More recently, the Australian Grape and Wine Incorporated (AWGI) estimated an AUD 40 million loss from the 2019/2020 summer bushfires [20]. Vineyard smoke exposure, therefore, remains a significant issue for the wine industry, particularly given the increasing frequency and severity of bushfires [21].
Grapevine leaves have also been found to accumulate VPs, and a positive correlation has been demonstrated between the levels of smoke compounds detected in leaves and wine when they were included in the primary fermentation [13,22,23]. From a physiological point of view, smoke exposure has also been shown to decrease stomatal conductance in leaves, which may result from the reaction of carbon dioxide (CO2) and carbon monoxide (CO) with water vapor in the substomatal cavity producing carbonic acid (H2CO3) [24,25]. Carbonic acid reduces the pH in the stomata, resulting in partial or complete stomatal closure [25,26]. Damage to leaf surfaces following smoke exposure has also been observed, with the development of necrotic lesions or, in extreme cases, total leaf necrosis [10,22,27]. This may be the result of ozone (O3) present in smoke, which has been linked to chlorophyll destruction and accelerated leaf senescence [28,29].
Some chromatographic techniques such as gas chromatography-mass spectrometry (GC-MS) and high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) have been developed to quantify levels of free and glycosidically bound VPs in grapes and wines [30,31,32,33]. While these techniques are currently used for qualitative and quantitative analysis and may assist growers in determining the level of smoke taint in the final wine, there are numerous shortcomings: sample preparation is time-consuming and destructive, and analyses require expensive reagents, standards, and equipment, as well as trained personnel. Furthermore, following a bushfire event, there may be long delays in the availability of results due to large numbers of samples being submitted to commercial laboratories for analysis [34,35]. Consequently, alternative methods of smoke taint analysis have recently been investigated and may offer non-destructive sample preparation, as well as accurate and rapid results.
The use of spectroscopic techniques has increased in recent years due to their ease of use, rapid results, minimal sample preparation, and non-destructive nature, all of which allow repeated measurements to be taken [34,35,36,37,38,39]. Furthermore, the development of smaller, handheld spectroscopic devices coupled with decreasing costs, has allowed these technologies to be more readily accessible and affordable to growers and farmers, while their portability allows for in-field use, reducing the risk of sample deterioration during transportation [39,40]. Ultraviolet (UV) to visible (Vis) spectroscopy involves the region between 200–780 nm, which can be used to analyze compounds containing organic acids, phenolic compounds, and pigments such as anthocyanins, carotenoids, and chlorophylls [41]. UV-Vis spectroscopy has been used to determine the contribution of chemical compounds towards the composition of extra virgin olive oils to determine the region in the Mediterranean it was produced, to optimize the aging process of Spanish wines, and to assess the impact of heating edible oils and to determine their acid level [42,43,44]. Near-infrared (NIR) spectroscopy between the light spectra regions of 780–2500 nm has been widely used in agricultural and food science applications, with NIR bands corresponding to overtones resulting from the vibrations of O-H, C-H, N-H, and S-H bonds [39,41]. Various spectroscopic techniques, most notably in the NIR region, have been used for numerous applications in viticulture, including the assessment of grape quality and ripeness as well as the authentication of geographical origin [38,45,46,47,48,49,50]. Research has also been conducted on the use of mid-infrared (MIR) spectroscopy (between 2500–25,000 nm) of the electromagnetic spectrum, as well as synchronous two-dimensional MIR correlation spectroscopy (2D-COS) for the classification of smoke tainted wines [34,35]. Both techniques showed potential for screening smoke tainted wine, with MIR spectroscopy achieving 61 and 70% classification rates for control and smoke affected wines, respectively. However, classification rates were affected by the degree of smoke taint, as well as compositional differences arising from the grape variety and oak maturation [34]. While this technology may help to assess wine samples for smoke taint, it does not provide an early, in-field detection system that could help growers identify which grapes may be contaminated before winemaking. At present, there is very little research investigating the in-field use of Vis-NIR spectroscopy for the classification of smoke-affected grapevine leaves and berries. Research by Fuentes and coworkers [19] developed a model using NIR spectroscopy between the region of 700–1100 nm to predict the levels of guaiacol glycoconjugates in berries and wine, and the levels of guaiacol in wine. These models may offer growers a non-destructive in-field detection system for grapevine smoke contamination. However, further research is required to determine the effectiveness of different NIR regions for monitoring smoke contamination.
Several chemometric techniques have been used to analyze spectral data, including partial least squares (PLS) regression, principal component analysis (PCA), and artificial neural networks (ANN), to name a few [41]. Of these techniques, ANNs have increased in popularity as classification, prediction, and clustering tools, particularly since they can better interpret the non-linear patterns of spectral data [51,52,53,54]. Machine learning (ML) modeling based on ANN can be trained from a set of given data known as ‘inputs’ or independent variables and form complex, non-linear relationships with these inputs and the ‘targets’ or dependent variables [54]. For example, preliminary ML models for the classification of smoke tainted grapevines have been developed using infra-red (IR) thermal imagery from canopies, which gave an indication of changes in stomatal conductance for classification of control and smoke-exposed grapevines [25]. In addition to this, another model has been proposed that aims to quantify levels of smoke derived compounds in grapes and wine using NIR spectroscopy measurements as inputs [25]. Furthermore, UV-Vis spectroscopy may offer insights into the degree of physiological performance of leaves as well as fruit ripening and quality through analyzing pigment content, such as chlorophylls, anthocyanins, and carotenoids [55,56,57,58,59].
The objective of this study was to investigate the use of NIR spectroscopy, coupled with ML modeling for the detection of grapevine smoke contamination. Grapevine leaves and berries were analyzed in the vineyard in a smoke trial using a NIR spectrometer, and the absorbance values were used as inputs to train different machine learning algorithms in order to create ANNs with the best classification performances. In addition to this, UV-Vis spectroscopy was used to assess the physiological performance and degree of senescence of leaves, as well as the degree of ripening and anthocyanin content of grapes. This may offer growers a rapid and non-destructive detection system that they can employ themselves to obtain real-time information regarding smoke exposure. This will facilitate timely decision-making around which fruit to sample for chemical analysis and/or to harvest to maintain wine quality.

2. Materials and Methods

2.1. Vineyard Site and Experimental Design for the Smoke Trial

The smoke trial was conducted in late January-early February during the 2018/2019 growing season, at the University of Adelaide’s Waite Campus in Urrbrae, South Australia (34°58′ S, 138°38′ E). The trial, described previously by Szeto and colleagues [60], involved the application of smoke and/or in-canopy misting to Cabernet Sauvignon grapevines and comprised five different treatments: a control (C), i.e., neither misting nor smoke exposure; (ii) a control with misting (CM), i.e., in-canopy misting but no smoke exposure; (iii) a high-density smoke treatment (HS); (iv) a high-density smoke treatment with misting (HSM); and (v) a low-density smoke treatment without misting (LS). Treatments were applied to Cabernet Sauvignon grapevines planted in 1998 at 2.0 and 3.3 m vine and row spacings, and trained to a bilateral cordon, vertical shoot positioned trellis system (VSP), hand-pruned to a two-node spur system, with under vine drip irrigation (twice weekly, from fruit set to pre-harvest). Smoke treatments were applied (approximately seven days post-véraison, the period grapes are thought to be most susceptible to smoke contamination [10]) using a purpose-built smoke tent (Figure 1a,b) and experimental conditions reported previously [4,61]: low and high-density smoke treatments were achieved by burning different fuel loads (i.e., ~1.5 and 5 kg of barley straw, respectively). In-canopy misting was evaluated as a method for mitigating the uptake of smoke-derived volatile phenols by grapes and involved the continuous application of fine water droplets (65 µm) to the grapevine bunch zone using a purpose-built sprinkler system (delivering water at 11 L/h), as previously described [62]. Each treatment was applied to six vines from three adjacent panels, except the HS treatment, which comprised only five vines, with treatments separated by at least one buffer vine. LS, HS, and HSM treatments comprised duplicate applications of smoke to 1.5 panels/three vines at a time (except for one HS treatment). The in-canopy sprinkler system was turned on 5 min before the first HSM treatment was applied and off 15 min after the second HSM treatment was completed, such that CM and HSM grapevines were misted for approximately 2.5 h in total. The second and fifth vine from each treatment (the middle vines from smoke treatments) were then selected for physiological and NIR measurements.

2.2. Physiological Measurements

The rate of photosynthesis (A), stomatal conductance (gs), and transpiration (E) were determined using a portable infrared gas analyzer equipped with a broad leaf chamber (LCpro-SD, ADC Bioscientific Ltd., Hoddesdon, UK). Measurements were taken on three leaves of each side of the canopy per vine (n = 12 leaves per treatment) with a photosynthetic photon flux density of 1000 µmol m−2 s−1 supplied by a high efficiency, low heat output, mixed red-blue light-emitting diode (LED) array unit. Water vapor and CO2 concentration in the chamber were set to ambient. Measurements were taken one day (24 h) after smoke treatments were applied, on clear, sunny days.

2.3. Determination of Volatile Phenols and Their Glycoconjugates in Grape Juice/Homogenate

The concentration of volatile phenols and their glycoconjugates were determined (in grape juice and homogenate, respectively) using analytical methods described previously [30,32,33,60]. Volatile phenols were measured by stable isotope dilution analysis (SIDA) [3,30,33], using an Agilent 6890 gas chromatograph coupled to a 5973-mass spectrometer (Agilent Technologies, Forest Hill, Vic., Australia). Isotopically labeled standards, i.e., d4-guaiacol and d3-syringol, were prepared in-house using methods outlined previously [3,30,33]. The limit of quantitation for volatile phenols was 1–2 µg/L. Volatile phenol glycoconjugates were also measured by SIDA [30,32], using an Agilent 1200 high-performance liquid chromatograph (HPLC) equipped with a 1290 binary pump, coupled to an AB SCIEX Triple QuadTM 4500 tandem mass spectrometer, with a Turbo VTM ion source (Framingham, MA, USA). The preparation of the isotopically labeled internal standard, i.e., d3-syringol gentiobioside, has been reported previously [30,32]. The limit of quantitation for volatile phenol glycosides was 1 µg/kg.

2.4. Near-Infrared Data Collection

Grapevine leaf and berry spectra were collected one day after smoke exposure, using a microPHAZIRTM RX Analyzer (Thermo Fisher Scientific, Waltham, MA, USA), which had a spectral range of 1596 to 2396 nm at intervals of 7–9 nm. Prior to undertaking the measurements and after every 10–15 readings, the device was calibrated using a white background calibration standard (included with the device). The white background was placed on top of the leaf while measuring to avoid signal noise inclusion due to variation in light or environmental changes. Leaves and berries were also analyzed using the Lighting Passport ProTM handheld spectrometer (Asensetek Incorporation, Xindian District, New Taipei City, Taiwan), which has a spectral range of 380–780 nm at intervals of 1 nm. Measurements were taken at approximately 3 cm from the leaves and berries. All measurements were conducted at ambient temperature between 9:00 a.m. and 6:00 p.m.
For the leaf spectral measurements, nine sunlit and nine shaded, mature, fully expanded leaves were selected (i.e., 18 leaves per vine, 36 leaves per treatment). Leaves were free of any visible signs of disease or blemishes. Each leaf was measured in three areas, in triplicate, using the microPHAZIRTM RX Analyzer, while three measurements per leaf were taken with the Lighting Passport ProTM handheld spectrometer. For the berry spectra, two bunches were selected per vine, and nine berries (three from the top, middle, and bottom of each bunch) were measured, in triplicates using the microPHAZIRTM RX Analyzer (n = 540). On the other hand, twelve berries per treatment were analyzed using the Lighting Passport ProTM (n = 180) while still attached to the bunch.

2.5. Calculating Spectral Indices

Spectral indices for the analysis of pigment content were calculated for both leaves and berries. Leaf spectra taken using the Lighting Passport ProTM were used to calculate the normalized difference vegetation index (NDVI), normalized anthocyanin index (NAI), plant senescence reflectance index (PSRI), and carotenoid reflectance index (CRI) [56,57,59,63,64,65]. Berry spectra were used to calculate the NAI and PSRI. The calculations and wavelengths used for determining these indices are given in Table 1.

2.6. Statistical Analysis

Physiological measurements, spectral indices, volatile phenols, and their glycoconjugates were analyzed by one-way analysis of variance (ANOVA) using Minitab® version 18.1 (Minitab Inc., State College, PA, USA). Mean comparisons were performed using the Fisher least significant difference (LSD) method as a post-hoc test at α = 0.05. Near-infrared data were analyzed using The Unscrambler X version 10.3 software (CAMO Software, Oslo, Norway). Absorbance values for all wavelengths were plotted for both the microPHAZIRTM RX Analyzer and Lighting Passport ProTM leaf and berry readings. Principal component analysis (PCA) was also performed using The Unscrambler X program. All microPHAZIRTM RX Analyzer measurements were pre-processed using the second derivative transformation, Savitzky–Golay derivation, and smoothing using The Unscrambler X version 10.3 software prior to the plotting of graphs and statistical analysis.

2.7. Artificial Neural Network Modeling

Three ANN models were developed for berry and leaf NIR readings, which were used as inputs to classify the different smoke treatments using customized code written in MATLAB® (version R2020a, MathWorks Inc., Natick, MA USA) (Figure 2). This code tested a total of 17 training algorithms in a loop to find the optimum in terms of accuracy and performance. Once the optimum training algorithm was identified, further training was performed to develop the most accurate ANN model. For both models, the Levenberg–Marquardt training algorithm was found to be the best algorithm, resulting in models with the highest accuracy and no signs of overfitting.
Overtones within the 1596–1800 nm range were used as inputs for the microPHAZIRTM leaf model (Model 1). This region was selected to avoid water overtones and any classification resulting from the water status of the vines. The entire spectral range was used for the microPHAZIRTM berry model (Model 2) (1596–2396 nm). The two models were developed using a random data division with 70% (n = 1134 for Model 1 and 378 for Model 2) training, 15% (n = 243 for Model 1 and 81 for Model 2) for validation with a mean squared error (MSE) performance algorithm and 15% (n = 243 for Model 1 and 81 for Model 2) for testing with a default derivative function. Ten hidden neurons were selected for each of the two models after conducting a trimming exercise with three, five, and ten neurons.

3. Results

3.1. Physiological Measurements

Results of gas exchange parameters are shown in Table 2. The transpiration rate was lower for the HS treatment (P < 0.005) with a mean rate of 1.43 mmol m−2 s−1, while no differences were observed in the other treatments. The CM and C treatments both had the highest gs values with an average value of 0.15 mol m−2 s−1 for each, while HS and LS treatments had the lowest average gs at 0.056 mol m−2 s−1 and 0.082 mol m−2 s−1 respectively. Mean rates of A were found to be highest in the C and CM treatments (10.77 µmol m−2 s−1 and 9.66 µmol m−2 s−1, respectively), while the LS and HS treatments had the lowest (7.01 µmol m−2 s−1 and 5.59 µmol m−2 s−1, respectively).

3.2. Levels of Smoke Taint Marker Compounds in Grape Juice/Homogenate

Differences in volatile phenol concentrations between HS and HSM treatments were found for guaiacol, 4-methylsyringol, and syringol (P < 0.05; Table S1). In particular, 4-methylsyringol and syringol had the largest differences in concentrations amongst the smoke treatments, with the HS treatment exhibiting the highest mean values (17 and 126 µg/L, respectively) followed by the HSM treatment (9 and 59 µg/L, respectively) while the CM treatments exhibited the lowest mean values (2 and 8 µg/L), which displayed the lowest mean value. There were no differences between the HS and HSM treatments, nor between the C, CM, and LS treatments for 4-methylguaiacol, phenol, and total cresols; however, HS and HSM grapes had significantly higher volatile phenol concentrations than C, CM, and LS grapes.
Some differences in volatile phenol glycoconjugate levels could be seen amongst the five smoke treatments. Some glycoconjugates displayed differences between the HS and HSM treatments. There was no difference in GuRG levels between the LS, HS, and HSM treatments, with no levels detected in the C and CM treatments. The HS smoke treatment had the highest levels of PhRG, PhGG, CrPG, SyGG, and SyPG, followed by the HSM and LS treatments and then the C and CM treatments. Interestingly the C and HS treatments had the highest level of CrGG followed by the CM and HSM treatment, while the LS treatment had the lowest concentration.

3.3. NIR Absorbance Patterns for Leaves and Berries

Absorbance spectra for the averages of replicates for both raw and transformed leaf absorbance spectra are depicted in Figure 3 and Figure 4. For the microPHAZIRTM RX Analyzer leaf absorbances, clear differences in spectral readings were observed for each smoking treatment. A peak was observed at approximately 1784–1793 nm (Figure 3a), while for the transformed data (Figure 3b), large peaks are present between 1596–1647 nm.
Differences in absorption readings were also found for the microPHAZIRTM RX Analyzer berry absorbance spectra (Figure 4a). Peaks were originally observed at approximately 1785 and 1902 nm, but in the transformed data (Figure 4b), large peaks were observed between approximately 1596–1640 nm and 1820–1940 nm.

3.4. Principal Component Analysis

Figure 5a shows the principal component analysis (PCA) for the microPHAZIRTM RX Analyzer leaf spectra with absorbance values between 1600–1800 nm. The first principal component (PC1) accounted for 62% of the data variability, while principal component two (PC2) accounted for 24%. Hence, 86% of the total variability was explained by these PCs. There was no clear separation of the different smoke treatments when modeled with the microPHAZIRTM leaf spectra. PC1 was represented by wavelengths between 1604–1621 nm and between 1621–1647 nm (loadings shown in Figure 5b). PC2 was represented by wavelengths between 1613–1647 nm, as well as 1604 nm.
Figure 6a shows the PCA for the microPHAZIRTM RX Analyzer berry spectra, where 59% of the data variability was described by PC1, while PC2 accounted for 10% of the data variability; thus, a total of 69% of the total data variability was explained by the first two components of the PCA. As with the microPHAZIRTM RX Analyzer leaf spectra, most of the smoke treatments overlapped quadrants. The CM treatment was grouped primarily in the upper right quadrant, while C and LS treatments were grouped primarily in the lower right. The HS treatment was located primarily in the upper right and left quadrants, while the HSM treatment was grouped in the left upper and lower quadrants. PC1 one was represented by the wavelength region 1604–1622. PC2 was represented by the wavelengths between 1630–1647 nm and 2374–2389 nm (loadings shown in Figure 6b).

3.5. Spectral Indices

Results for the spectral indices are shown in Table 3. In the case of the leaf NDVI and NAI, the HS and C treatments had the lowest mean values (0.72 and 0.64 for the HS treatment and 0.84 and 0.74 for the C) (P < 0.05). There were no differences for the remaining treatments. For the leaf PSRI, the HS treatment had the highest mean value at 0.065, with no differences for the remaining treatments. For the leaf CRI500, the LS and HS treatments had the highest values at 1.45 and 1.20, respectively, and for the CRI700, the LS treatments had the highest mean values at 1.76, with no differences for the remaining treatments.
In the case of the berry NAI, the HS and LS treatments had the highest mean values with 0.88 and 0.87, with both the C and LS treatments having the lowest mean values of 0.80 and 0.75. For the PSRI, both the LS and C treatments had the highest mean values of 0.02, while the HSM had the lowest value at −0.02.

3.6. Artificial Neural Network Models

Table 4 shows the confusion matrices for the two models developed using the spectral readings as inputs and the experimental treatments as targets. Both models displayed high accuracy in classifying the spectral readings according to the treatments, with an overall accuracy of 98% for the microPHAZIRTM leaf model (Model 1) and 97.4% for the microPHAZIRTM berry model (Model 2). Models 1 and 2 presented validation accuracies (94% and 93%, respectively) close to those of the training stage (100% both models). Furthermore, performance values for training (Models 1 and 2: MSE < 0.01) were lower than the other stages and validation (Model 1: MSE = 0.02; Model 2: MSE = 0.03) and testing (Model 1: MSE = 0.02; Model 2: MSE = 0.04) were similar; this indicates that there were no signs of overfitting for both Model 1 and Model 2.
Figure 7 depicts the receiver operating characteristic (ROC) curves for the two ANN models developed. All models showed high true-positive rates (sensitivity) and low false-positive rates (specificity) for classifying the spectral readings according to the experimental treatment, which can also be observed in the last column of each confusion matrix. For Model 2, the HS treatment had the highest sensitivity (100%), followed by the CM and HSM treatments (99.1% each) and LS treatment (96.3%). The C treatment had the lowest sensitivity of 92.6% for this model. For Model 1, the C treatment had the highest sensitivity (99.1%), followed by the LS treatment (98.8%), HS treatment (97.8%), and CM treatment (97.5%), while the HSM had the lowest sensitivity of 96.9%.

4. Discussion

4.1. Physiological Measurements

Leaf gas exchange parameters were measured the day after smoking. The three smoke treatments showed significant reductions in gs, in particular, the high-density smoke without misting (HS) treatment, which showed the lowest average reading for gs (Table 2). Stomatal closure is one of the first responses to smoke exposure undertaken by plants [6,26], and a study by Ristic and colleagues [26] found that the time required for gs to recover following one hour of smoke exposure for Cabernet Sauvignon grapevines was approximately 6–10 days. A previous study by Bell et al. [6] found that gs of potted Cabernet Sauvignon grapevines had returned to 60% of pre-smoke exposure rate following fifteen min exposure to smoke using Tasmanian blue gum (Eucalyptus globulus L.) leaves as fuel, while rates had returned to 80% of pre-smoke values following exposure to smoke derived from Coast Live Oak (Quercus agrifolia Née) leaves. This indicates that in addition to the type of fuel used, the intensity of smoke exposure may also affect the extent of stomatal closure and, hence, reduction in gs. It is, therefore, not surprising that the HS treatment had the lowest gs. However, it is interesting that the low smoke treatment (LS) had lower gs than the high smoke with misting treatment (HSM), which indicates that misting may have reduced the effect of smoke exposure on gs. During a bushfire, the type of fuel burnt will vary depending on the region and the type of plant species native to the area, as well as the amount of smoke exposure due to land topography and wind vectors; therefore, the effect on gs may vary [17,18,23,66]. While misting only partially prevented the uptake of volatile phenols and glycoconjugates in grapes [60], it did appear to have a physiological effect. It is evident that misting reduced the effect of smoke exposure on gs. Smoke contains a complex mixture of gases such as sulfur dioxide (SO2), O3, and nitrogen dioxide (NO2), as well as dust particles that have been shown to inhibit photosynthesis and affect stomatal opening [6,26,29]. Stomata are the primary point of entry for these gases and dust particles [6]; therefore, misting may help prevent the uptake of dust and other particles by trapping them in water that has condensed on the leaf surface, preventing their entrance into the stomata. The present water may also act as a solvent for gases such as SO2 and NO2, thereby incorporating them into a solution that then may drip off the leaf surface. In addition to this, smoke exposure may trigger stomatal closure by producing high vapor pressure deficits [26,29]. The presence of misting may help reduce the leaf-to-air vapor pressure difference produced by smoke exposure, thereby reducing the impact on gs. Misting also appeared to reduce the effect of smoke exposure on transpiration rate (E) as there were no differences between the two control treatments and the LS and HSM treatments. Only the HS treatment had significantly reduced E. Mean rates of photosynthesis (A) followed similar patterns to gs, with the HS treatment having the lowest value, followed by the LS treatment and then the HSM treatment, while the control without misting (C) had the highest rate of A. This indicates that while misting may have reduced A in the control treatments, it may also help reduce the effects of smoke exposure on A.

4.2. Near-Infrared Spectroscopy Patterns and Principal Component Analysis

From the PCA biplots (Figure 5 and Figure 6) and spectra (Figure 3 and Figure 4) generated in the current study, it is evident that smoke exposure alters the NIR spectral signals of grapevine leaves and berries, and this may prove useful for the detection of grapevine smoke contamination. For the microPHAZIRTM RX Analyzer leaf spectra, high loadings (Figure 5b.) were observed for the wavelength regions between 1604–1621, 1621–1647, and 1613–1647 nm, all of which correspond to C-H stretching of sugars and aromatic compounds [67,68,69,70]. For the microPHAZIRTM RX Analyzer berry spectra, high loadings (Figure 6b.) were observed for the wavelength regions between 1604–1622, 1630–1647, and 2374–2389 nm, which correspond to C-H stretching of sugars, such as glucose, as well as aromatic hydrocarbons, which may be due to the presence of smoke-derived volatile phenols, such as guaiacols, cresols, and syringols, and their glycoconjugates [67,68,71,72].

4.3. Spectral Indices

4.3.1. Leaf

The normalized difference vegetation index (NDVI) gives an indication of plant vigor and fruit ripening resulting from relative changes in chlorophyll content. It is based on the variation between the maximum absorption of red by chlorophyll pigments and the maximum reflectance in the infrared caused by leaf cellular structure [56,57,73,74,75]. Similarly, relative changes in anthocyanin content are expressed as the normalized anthocyanin index (NAI). Both the NDVI and NAI are expressed as a normalized value between −1 (lack of green or redness) to +1 (green or red) [56,57]. Not surprisingly, HS leaves had the lowest NDVI and NAI values. Previous studies investigating the effects of pollution on leaf pigments found a decrease in photosynthetic pigments following exposure to pollutants, including sulfur dioxide (SO2), carbon dioxide (CO2), nitrogen dioxide (NO2), and ozone (O3) [59,76,77]. These studies are often used as comparisons for investigating the effects of smoke exposure on leaves as compounds in air pollution can also found in smoke [6,22]. There were no differences in NDVI and NAI values between the LS, HSM, and control treatments (C and CM), indicating that misting may reduce the effects of smoke exposure on leaf pigments, and low levels of smoke exposure for one hour may also have no effect. Longer periods of smoke exposure (days or weeks, as is often the case with wildfires) may be required to cause a noticeable change in leaf pigments.
The plant senescence reflectance index (PSRI) gives an indication of the stage of leaf senescence and fruit ripening through assessing changes in carotenoid accumulation and their proportion to chlorophyll. Values range from −1 to +1, with higher values indicating increased stress and carotenoid accumulation [55,63,64,65,78]. The PSRI was highest for the HS treatment, indicating heightened stress and leaf senescence. This also corresponds with the high CRI500 value for this smoke treatment, indicating increased carotenoid accumulation.

4.3.2. Berries

Research by Noestheden et al. [5] found that smoke exposure induced changes in phenylpropanoid metabolites in Pinot Noir berries and wine, some of which are associated with the color and mouthfeel of the wine. Berries exposed to HS and LS treatments had the highest mean NAI values, indicating that smoke exposure may increase anthocyanin content, possibly due to an increase in phenolic accumulation as a stress response induced by exposure to ozone present in smoke [5,79,80]. The HSM treatment had a low NAI value, indicating that misting may reduce anthocyanin concentrations through increased irrigation. Castellarin et al. [81] found that early (before véraison) and late (after the onset of ripening) season, water deficits increased anthocyanin accumulation during ripening. The application of in-canopy misting may reduce water stress and, therefore, reduce anthocyanin accumulation.
Interestingly the HSM followed by the HS treatments had the lowest PSRI values. As carotenoid concentrations in grapes generally decrease during véraison, this may have resulted in lower PSRI values. Therefore, the PSRI may not be suitable for assessing the degree of ripening in grape berries.

4.4. ANN Modeling

Both ANN models classified leaf and berry readings as a function of smoke exposure with high accuracy. The microPHAZIRTM leaf model (model 1) had the highest positive classification, with 98% accuracy (Table 4). The NIR region selected for use in Model 1 was between 1600–1800 nm in order to minimize any possible interference due to the absorption spectra of water in the region of approximately 1930 nm [69]. Furthermore, the region between 1680–1690 nm is associated with aromatic C-H stretching [67]; as such, any patterns observed by the ANN would most likely be due to the presence of smoke-derived volatile phenols. Research by Kennison [22] found a positive correlation between levels of smoke-derived compounds found in leaves and levels in wine; this ANN model developed may, therefore, offer a rapid, in-field method for assessing grapevine smoke contamination. It also demonstrates great promise for further research into the use of NIR spectroscopy coupled with unmanned aerial vehicles (UAVs) with Global Positioning System (GPS) trackers, which could fly over vineyards to scan grapevine canopies and provide maps of smoke contaminated regions.
The microPHAZIRTM berry model (model 2) also had a high overall accuracy in classifying grape berries according to smoke treatment (97.4%). For Model 2, the entire wavelength range between 1600–2396 nm was used. This includes the C-H stretching of aromatic compounds at 1680 nm, O-H stretching at 1930 nm associated with glucose, cellulose, and water, and C=O second overtone associated with carboxylic acids and water between 1900–1910 nm [67,69]. As NIR measurements were conducted in-field on whole berries, this offers a non-destructive tool for assessing grapevine smoke contamination. Whole grapes may be used for assessment as smoke compounds have been found to occur primarily in grape skins [3,25]. Furthermore, the Lighting PassportTM smart handheld spectrometer may be of interest to growers due to its affordability compared to other spectrometers. It is also very small and lightweight, making it easy to undertake measurements in-field, and it can be connected to smartphones via Bluetooth, where data can be stored and retrieved for later analysis [82].
The two ANN models more accurately differentiated the spectral readings relative to PCA. This may be because ANNs are better suited to handle complex, non-linear data, and more readily find patterns or relationships between data than other forms of analysis [53,83,84,85]. Research by Janik et al. [53] found that the combination of ANNs with partial least squares (PLS) or PCA overcomes issues of non-linearity as well as increasing the accuracy of regression models in predicting total anthocyanin concentrations in red grape homogenates. This may also explain why Model 2 was able to accurately differentiate the berry spectral readings from C, CM, and LS treatments, despite analysis of variance indicating there were no statistically significant differences.
As smoke exposure altered the chemical fingerprinting of grapevine leaves and berries, the ANN models were able to detect changes in the spectral patterns and then classify the readings as a function of experimental treatments. This may offer grape growers a rapid method of assessing the level of smoke contamination in grape berries and leaves, with a high level of accuracy and precision. This may assist growers in deciding which berry samples to send for further chemical analysis to quantify the levels of smoke compounds in grapes and predict the level of smoke taint in the final wine, or they may decide to avoid harvesting heavily contaminated grapes for winemaking. Furthermore, as this method is non-destructive, repeated measurements are possible. By knowing the level of smoke contamination, growers can make informed decisions.
While the ANN models developed were able to classify Cabernet Sauvignon leaf and berry spectra accurately, further research is required to assess whether these models can be used for other grape varieties, as differences in berry composition and leaf physiology may affect the accuracy of classification [6,34]. Previous research evaluated MIR spectroscopy for the classification of smoke tainted wines found compositional differences due to grape variety prevailed over differences resulting from low levels of smoke exposure [34]. Furthermore, the physiological responses of different grape varieties to smoke were found to vary, both in magnitude and in recovery time [6,26]. Thus, further testing of these models using berry and leaf spectra from different grapevine varieties is required.

5. Conclusions

Results from this study indicate that smoke exposure alters the NIR spectra of Cabernet Sauvignon grapevine leaves and berries. As a result, accurate classification models can be developed using ANN modeling. Artificial neural networks are better at classifying non-linear or complex data than traditional techniques, such as principal component analysis. Furthermore, the use of UV-Vis spectroscopy may offer insights into the physiological performance of leaves and the quality and degree of ripening of grapes. These techniques may assist grape growers in identifying grapevines that have been contaminated by smoke, thereby informing decision-making to avoid harvesting and processing heavily contaminated grapes and/or the need for mitigation techniques to manage the risk of smoke taint in resulting wine. Further testing of the ANN models developed in the current study is required to assess their accuracy in classifying grapevine leaf and berry spectra from other grape varieties.

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-8220/20/18/5099/s1, Table S1: Concentrations of volatile phenols in grape juice (µg/L) and their glycoconjugates in grape homogenate (µg/kg) one hour after smoke treatments.

Author Contributions

Conceptualization, V.S., and S.F.; data curation, V.S., C.G.V., and S.F.; formal analysis, V.S.; funding acquisition, S.F.; investigation, V.S., and C.G.V.; methodology, V.S., C.G.V., C.S., K.L.W., and S.F.; project administration, K.L.W., and S.F.; resources, K.L.W., R.D.B., and S.F.; software, C.G.V., and S.F.; supervision, D.D.T., A.P., and S.F.; validation, C.G.V., K.L.W., and S.F.; visualization, V.S., C.G.V., D.D.T., and S.F.; writing—original draft, V.S.; writing—review and editing, V.S., C.G.V., C.S., K.L.W., D.D.T., A.P., R.D.B., and S.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

This research was supported through the Australian Government Research Training Program Scholarship, as well as the Digital Viticulture program funded by the University of Melbourne’s Networked Society Institute, Australia. C.S. was supported by the Australian Research Council Training Centre for Innovative Wine Production (www.arcwinecentre.org.au), which is funded as part of the ARC’s Industrial Transformation Research Program (Project No. ICI70100008), with support from Wine Australia and industry partners. The authors greatly acknowledge the Digital Agriculture, Food, and Wine Group.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Cain, N.; Hancock, F.; Rogers, P.; Downey, M. The effect of grape variety and smoking duration on the accumulation of smoke taint compounds in wine. Wine Vitic. J. 2013, 28, 48–49. [Google Scholar]
  2. CSIRO. Australian Government Bureau of Meteorology. State Clim. 2018, 2018, 5. [Google Scholar]
  3. Dungey, K.A.; Hayasaka, Y.; Wilkinson, K.L. Quantitative analysis of glycoconjugate precursors of guaiacol in smoke-affected grapes using liquid chromatography–tandem mass spectrometry based stable isotope dilution analysis. Food Chem. 2011, 126, 801–806. [Google Scholar] [CrossRef]
  4. Kennison, K.; Gibberd, M.; Pollnitz, A.; Wilkinson, K. Smoke-derived taint in wine: The release of smoke-derived volatile phenols during fermentation of merlot juice following grapevine exposure to smoke. J. Agric. Food Chem. 2008, 56, 7379–7383. [Google Scholar] [CrossRef] [PubMed]
  5. Noestheden, M.; Noyovitz, B.; Riordan-Short, S.; Dennis, E.G.; Zandberg, W.F. Smoke from simulated forest fire alters secondary metabolites in Vitis vinifera L. Berries and wine. Planta 2018, 248, 1537–1550. [Google Scholar] [CrossRef]
  6. Bell, T.; Stephens, S.; Moritz, M. Short-term physiological effects of smoke on grapevine leaves. Int. J. Wildland Fire 2013, 22, 933–946. [Google Scholar] [CrossRef] [Green Version]
  7. De Vries, C.; Mokwena, L.; Buica, A.; McKay, M. Determination of volatile phenol in Cabernet sauvignon wines, made from smoke-affected grapes, by using hs-spme GC-MS. S. Afr. J. Enol. Vitic. 2016, 37, 15–21. [Google Scholar] [CrossRef] [Green Version]
  8. Noestheden, M.; Dennis, E.G.; Romero-Montalvo, E.; DiLabio, G.A.; Zandberg, W.F. Detailed characterization of glycosylated sensory-active volatile phenols in smoke-exposed grapes and wine. Food Chem. 2018, 259, 147–156. [Google Scholar] [CrossRef]
  9. Parker, M.; Osidacz, P.; Baldock, G.A.; Hayasaka, Y.; Black, C.A.; Pardon, K.H.; Jeffery, D.W.; Geue, J.P.; Herderich, M.J.; Francis, I.L. Contribution of several volatile phenols and their glycoconjugates to smoke-related sensory properties of red wine. J. Agric. Food Chem. 2012, 60, 2629–2637. [Google Scholar] [CrossRef]
  10. Kennison, K.; Wilkinson, K.; Pollnitz, A.; Williams, H.; Gibberd, M. Effect of smoke application to field-grown Merlot grapevines at key phenological growth stages on wine sensory and chemical properties. Aust. J. Grape Wine Res. 2011, 17, 5–12. [Google Scholar] [CrossRef]
  11. Ristic, R.; van der Hulst, L.; Capone, D.; Wilkinson, K. Impact of bottle aging on smoke-tainted wines from different grape cultivars. J. Agric. Food Chem. 2017, 65, 4146–4152. [Google Scholar] [CrossRef] [PubMed]
  12. Härtl, K.; Schwab, W. Smoke taint in wine-how smoke-derived volatiles accumulate in grapevines. Wines Vines 2018, 99, 62–64. [Google Scholar]
  13. Hayasaka, Y.; Baldock, G.; Pardon, K.; Jeffery, D.; Herderich, M. Investigation into the formation of guaiacol conjugates in berries and leaves of grapevine Vitis vinifera L. Cv. Cabernet sauvignon using stable isotope tracers combined with hplc-ms and ms/ms analysis. J. Agric. Food Chem. 2010, 58, 2076–2081. [Google Scholar] [CrossRef] [PubMed]
  14. Hoj, P.; Pretorius, I.; Blair, R. Investigations conducted into the nature and amelioration of taints in grapes and wine, caused by smoke resulting from bushfires. Aust. Wine Res. Inst. Annu. Rep. 2003, 37–39. [Google Scholar]
  15. Kelly, D.; Zerihun, A.; Singh, D.; Vitzthum von Eckstaedt, C.; Gibberd, M.; Grice, K.; Downey, M. Exposure of grapes to smoke of vegetation with varying lignin composition and accretion of lignin derived putative smoke taint compounds in wine. Food Chem. 2012, 135, 787–798. [Google Scholar] [CrossRef] [Green Version]
  16. Singh, D.; Chong, H.; Pitt, K.; Cleary, M.; Dokoozlian, N.; Downey, M. Guaiacol and 4-methylguaiacol accumulate in wines made from smoke-affected fruit because of hydrolysis of their conjugates. Aust. J. Grape Wine Res. 2011, 17, S13–S21. [Google Scholar] [CrossRef]
  17. Noestheden, M.; Thiessen, K.; Dennis, E.G.; Tiet, B.; Zandberg, W.F. Quantitating organoleptic volatile phenols in smoke-exposed Vitis vinifera berries. J. Agric. Food Chem. 2017, 65, 8418–8425. [Google Scholar] [CrossRef]
  18. Krstic, M.; Johnson, D.; Herderich, M. Review of smoke taint in wine: Smoke-derived volatile phenols and their glycosidic metabolites in grapes and vines as biomarkers for smoke exposure and their role in the sensory perception of smoke taint. Aust. J. Grape Wine Res. 2015, 21, 537–553. [Google Scholar] [CrossRef]
  19. Kennison, K.; Wilkinson, K.; Williams, H.; Smith, J.; Gibberd, M. Smoke-derived taint in wine: Effect of postharvest smoke exposure of grapes on the chemical composition and sensory characteristics of wine. J. Agric. Food Chem. 2007, 55, 10897–10901. [Google Scholar] [CrossRef]
  20. Claughton, C.; Jeffery, C.; Pritchard, M.; Hough, C.; Wheaton, C. Wine Industry’s ‘Black Summer’ as Cost of Smoke Taint, Burnt Vineyards, and Lost Sales Add up. ABC News, 28 February 2020. [Google Scholar]
  21. Dutta, R.; Das, A.; Aryal, J. Big data integration shows Australian bushfire frequency is increasing significantly. R. Soc. Open Sci. 2016, 3, 150241. [Google Scholar] [CrossRef] [Green Version]
  22. Kennison, K. Bushfire Generated Smoke Taint in Grapes and Wine. Final Report to Grape and Wine Research and Development Corporation; RD 05/02–3; Department of Agriculture and Food Western Australia: Kalgoorlie, Australia, 2009. [Google Scholar]
  23. Simos, C. The implications of smoke taint and management practices. Aust. Vitic. Jan./Feb. 2008, 12, 77–80. [Google Scholar]
  24. Fuentes, S.; Tongson, E. Advances in smoke contamination detection systems for grapevine canopies and berries. Wine Vitic. J. 2017, 32, 36. [Google Scholar]
  25. Fuentes, S.; Tongson, E.J.; De Bei, R.; Gonzalez Viejo, C.; Ristic, R.; Tyerman, S.; Wilkinson, K. Non-invasive tools to detect smoke contamination in grapevine canopies, berries and wine: A remote sensing and machine learning modeling approach. Sensors 2019, 19, 3335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Ristic, R.; Fudge, A.; Pinchbeck, K.; De Bei, R.; Fuentes, S.; Hayasaka, Y.; Tyerman, S.; Wilkinson, K. Impact of grapevine exposure to smoke on vine physiology and the composition and sensory properties of wine. Theor. Exp. Plant Physiol. 2016, 28, 67–83. [Google Scholar] [CrossRef]
  27. Collins, C.; Gao, H.; Wilkinson, K. An observational study into the recovery of grapevines (Vitis vinifera L.) following a bushfire. Am. J. Enol. Vitic. 2014, 65, 285–292. [Google Scholar] [CrossRef]
  28. Mikkelsen, T.N.; Heide-Jørgensen, H.S. Acceleration of leaf senescence infagus sylvatica l. By low levels of tropospheric ozone demonstrated by leaf colour, chlorophyll fluorescence and chloroplast ultrastructure. Trees 1996, 10, 145–156. [Google Scholar] [CrossRef]
  29. Calder, J.; Lifferth, G.; Moritz, M.; St Clair, S. Physiological effects of smoke exposure on deciduous and conifer tree species. Int. J. For. Res. 2010, 2010, 438930. [Google Scholar] [CrossRef] [Green Version]
  30. Hayasaka, Y.; Baldock, G.; Parker, M.; Pardon, K.; Black, C.; Herderich, M.; Jeffery, D. Glycosylation of smoke-derived volatile phenols in grapes as a consequence of grapevine exposure to bushfire smoke. J. Agric. Food Chem. 2010, 58, 10989–10998. [Google Scholar] [CrossRef]
  31. Hayasaka, Y.; Dungey, K.; Baldock, G.; Kennison, K.; Wilkinson, K. Identification of a β-d-glucopyranoside precursor to guaiacol in grape juice following grapevine exposure to smoke. Anal. Chim. Acta 2010, 660, 143–148. [Google Scholar] [CrossRef]
  32. Hayasaka, Y.; Parker, M.; Baldock, G.A.; Pardon, K.H.; Black, C.A.; Jeffery, D.W.; Herderich, M.J. Assessing the impact of smoke exposure in grapes: Development and validation of a HPLC-MS/MS method for the quantitative analysis of smoke-derived phenolic glycosides in grapes and wine. J. Agric. Food Chem. 2013, 61, 25–33. [Google Scholar] [CrossRef]
  33. Pollnitz, A.P.; Pardon, K.H.; Sykes, M.; Sefton, M.A. The effects of sample preparation and gas chromatograph injection techniques on the accuracy of measuring guaiacol, 4-methylguaiacol and other volatile oak compounds in oak extracts by stable isotope dilution analyses. J. Agric. Food Chem. 2004, 52, 3244–3252. [Google Scholar] [CrossRef] [PubMed]
  34. Fudge, A.; Wilkinson, K.; Ristic, R.; Cozzolino, D. Classification of smoke tainted wines using mid-infrared spectroscopy and chemometrics. J. Agric. Food Chem. 2012, 60, 52–59. [Google Scholar] [CrossRef] [PubMed]
  35. Fudge, A.; Wilkinson, K.; Ristic, R.; Cozzolino, D. Synchronous two-dimensional MIR correlation spectroscopy (2d-cos) as a novel method for screening smoke tainted wine. Food Chem. 2013, 139, 115–119. [Google Scholar] [CrossRef] [PubMed]
  36. van der Hulst, L.; Munguia, P.; Culbert, J.A.; Ford, C.M.; Burton, R.A.; Wilkinson, K.L. Accumulation of volatile phenol glycoconjugates in grapes following grapevine exposure to smoke and potential mitigation of smoke taint by foliar application of kaolin. Planta 2019, 249, 941–952. [Google Scholar] [CrossRef] [PubMed]
  37. Barbin, D.F.; Felicio, A.L.D.S.M.; Sun, D.-W.; Nixdorf, S.L.; Hirooka, E.Y. Application of infrared spectral techniques on quality and compositional attributes of coffee: An overview. Food Res. Int. 2014, 61, 23–32. [Google Scholar] [CrossRef] [Green Version]
  38. Urraca, R.; Sanz-Garcia, A.; Tardaguila, J.; Diago, M.P. Estimation of total soluble solids in grape berries using a handheld NIR spectrometer under field conditions. J. Sci. Food Agric. 2016, 96, 3007–3016. [Google Scholar] [CrossRef] [PubMed]
  39. Dos Santos, C.A.T.; Lopo, M.; Ricardo, N.; Lopes, J. A review on the applications of portable near-infrared spectrometers in the agro-food industry. Appl. Spectrosc. 2013, 67, 1215–1233. [Google Scholar] [CrossRef]
  40. Hall, A. Remote sensing applications for viticultural terroir analysis. Elements 2018, 14, 185–190. [Google Scholar] [CrossRef]
  41. Yu, J.; Wang, H.; Zhan, J.; Huang, W. Review of recent UV–Vis and infrared spectroscopy researches on wine detection and discrimination. Appl. Spectrosc. Rev. 2018, 53, 65–86. [Google Scholar] [CrossRef]
  42. Zhang, W.; Li, N.; Feng, Y.; Su, S.; Li, T.; Liang, B. A unique quantitative method of acid value of edible oils and studying the impact of heating on edible oils by UV–Vis spectrometry. Food Chem. 2015, 185, 326–332. [Google Scholar] [CrossRef]
  43. Ferreiro-González, M.; Ruiz-Rodríguez, A.; Barbero, G.F.; Ayuso, J.; Álvarez, J.A.; Palma, M.; Barroso, C.G. FT-IR, Vis spectroscopy, color and multivariate analysis for the control of ageing processes in distinctive spanish wines. Food Chem. 2019, 277, 6–11. [Google Scholar] [CrossRef]
  44. Alves, F.C.G.B.S.; Coqueiro, A.; Março, P.H.; Valderrama, P. Evaluation of olive oils from the mediterranean region by UV–Vis spectroscopy and independent component analysis. Food Chem. 2019, 273, 124–129. [Google Scholar] [CrossRef] [PubMed]
  45. Mandrile, L.; Zeppa, G.; Giovannozzi, A.M.; Rossi, A.M. Controlling protected designation of origin of wine by raman spectroscopy. Food Chem. 2016, 211, 260–267. [Google Scholar] [CrossRef] [Green Version]
  46. Larraín, M.; Guesalaga, A.R.; Agosín, E. A multipurpose portable instrument for determining ripeness in wine grapes using NIR spectroscopy. IEEE Trans. Instrum. Meas. 2008, 57, 294–302. [Google Scholar] [CrossRef]
  47. González-Caballero, V.; Pérez-Marín, D.; López, M.-I.; Sánchez, M.-T. Optimization of NIR spectral data management for quality control of grape bunches during on-vine ripening. Sensors 2011, 11, 6109–6124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Barnaba, F.E.; Bellincontro, A.; Mencarelli, F. Portable NIR-AOTF spectroscopy combined with winery FTIR spectroscopy for an easy, rapid, in-field monitoring of sangiovese grape quality. J. Sci. Food Agric. 2014, 94, 1071–1077. [Google Scholar] [CrossRef]
  49. Fernández-Novales, J.; López, M.-I.; Sánchez, M.-T.; Morales, J.; González-Caballero, V. Shortwave-near infrared spectroscopy for determination of reducing sugar content during grape ripening, winemaking, and aging of white and red wines. Food Res. Int. 2009, 42, 285–291. [Google Scholar] [CrossRef]
  50. Cao, F.; Wu, D.; He, Y. Soluble solids content and ph prediction and varieties discrimination of grapes based on visible–near infrared spectroscopy. Comput. Electron. Agric. 2010, 71, S15–S18. [Google Scholar] [CrossRef]
  51. Dolatabadi, Z.; Elhami Rad, A.H.; Farzaneh, V.; Akhlaghi Feizabad, S.H.; Estiri, S.H.; Bakhshabadi, H. Modeling of the lycopene extraction from tomato pulps. Food Chem. 2016, 190, 968–973. [Google Scholar] [CrossRef]
  52. Gumus, Z.P.; Ertas, H.; Yasar, E.; Gumus, O. Classification of olive oils using chromatography, principal component analysis and artificial neural network modelling. J. Food Meas. Charact. 2018, 12, 1325–1333. [Google Scholar] [CrossRef]
  53. Janik, L.J.; Cozzolino, D.; Dambergs, R.; Cynkar, W.; Gishen, M. The prediction of total anthocyanin concentration in red-grape homogenates using visible-near-infrared spectroscopy and artificial neural networks. Anal. Chim. Acta 2007, 594, 107–118. [Google Scholar] [CrossRef] [PubMed]
  54. Pero, M.; Askari, G.; Skåra, T.; Skipnes, D.; Kiani, H. Change in the color of heat-treated, vacuum-packed broccoli stems and florets during storage: Effects of process conditions and modeling by an artificial neural network. J. Sci. Food Agric. 2018, 98, 4151–4159. [Google Scholar] [CrossRef] [PubMed]
  55. Merzlyak, M.N.; Gitelson, A.A.; Chivkunova, O.B.; Rakitin, V.Y. Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening. Physiol. Plant. 1999, 106, 135–141. [Google Scholar] [CrossRef] [Green Version]
  56. Overbeck, V.; Schmitz, M.; Blanke, M. Non-destructive sensor-based prediction of maturity and optimum harvest date of sweet cherry fruit. Sensors 2017, 17, 277. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Solomakhin, A.A.; Blanke, M.M. Overcoming adverse effects of hailnets on fruit quality and microclimate in an apple orchard. J. Sci. Food Agric. 2007, 87, 2625–2637. [Google Scholar] [CrossRef]
  58. Merzlyak, M.N.; Solovchenko, A.E.; Gitelson, A.A. Reflectance spectral features and non-destructive estimation of chlorophyll, carotenoid and anthocyanin content in apple fruit. Postharvest Biol. Technol. 2003, 27, 197–211. [Google Scholar] [CrossRef]
  59. Sims, D.A.; Gamon, J.A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 2002, 81, 337–354. [Google Scholar] [CrossRef]
  60. Szeto, C.; Ristic, R.; Capone, D.; Puglisi, C.; Pagay, V.; Culbert, J.; Jiang, W.; Herderich, M.; Tuke, J.; Wilkinson, K. Uptake and glycosylation of smoke-derived volatile phenols by Cabernet sauvignon grapes and their subsequent fate during winemaking. Molecules 2020, 25, 3720. [Google Scholar] [CrossRef]
  61. Ristic, R.; Osidacz, P.; Pinchbeck, K.; Hayasaka, Y.; Fudge, A.; Wilkinson, K. The effect of winemaking techniques on the intensity of smoke taint in wine. Aust. J. Grape Wine Res. 2011, 17, S29–S40. [Google Scholar] [CrossRef]
  62. Caravia, L.; Pagay, V.; Collins, C.; Tyerman, S.D. Application of sprinkler cooling within the bunch zone during ripening of Cabernet sauvignon berries to reduce the impact of high temperature. Aust. J. Grape Wine Res. 2017, 23, 48–57. [Google Scholar] [CrossRef]
  63. Mate, A.R.; Deshmukh, R.R. Analysis of effects of air pollution on chlorophyll, water, carotenoid and anthocyanin content of tree leaves using spectral indices. Int. J. Eng. Sci 2016, 6, 5465–5474. [Google Scholar]
  64. Gitelson, A.A.; Zur, Y.; Chivkunova, O.B.; Merzlyak, M.N. Assessing carotenoid content in plant leaves with reflectance spectroscopy. Photochem. Photobiol. 2002, 75, 272–281. [Google Scholar] [CrossRef]
  65. García-Estévez, I.; Quijada-Morín, N.; Rivas-Gonzalo, J.C.; Martínez-Fernández, J.; Sánchez, N.; Herrero-Jiménez, C.M.; Escribano-Bailón, M.T. Relationship between hyperspectral indices, agronomic parameters and phenolic composition of Vitis vinifera cv Tempranillo grapes. J. Sci. Food Agric. 2017, 97, 4066–4074. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. De Vries, C.; Buica, A.; McKay, J.B.M. The impact of smoke from vegetation fires on sensory characteristics of Cabernet sauvignon wines made from affected grapes. S. Afr. J. Enol. Vitic. 2016, 37, 22–30. [Google Scholar] [CrossRef] [Green Version]
  67. Boido, E.; Fariña, L.; Carrau, F.; Dellacassa, E.; Cozzolino, D. Characterization of glycosylated aroma compounds in tannat grapes and feasibility of the near infrared spectroscopy application for their prediction. Food Anal. Methods 2013, 6, 100–111. [Google Scholar] [CrossRef]
  68. Cynkar, W.; Cozzolino, D.; Dambergs, R.; Janik, L.; Gishen, M. Effect of variety, vintage and winery on the prediction by visible and near infrared spectroscopy of the concentration of glycosylated compounds (g-g) in white grape juice. Aust. J. Grape Wine Res. 2007, 13, 101–105. [Google Scholar] [CrossRef]
  69. Burns, D.; Ciurczak, E. Handbook of Near-Infrared Analysis; CRC Press: Boca Raton, FL, USA, 2007. [Google Scholar]
  70. Osborne, B.G.; Fearn, T.; Hindle, P.H. Practical NIR Spectroscopy with Applications in Food and Beverage Analysis; Longman Scientific & Technical: Harlow, UK, 1993; Volume 2. [Google Scholar]
  71. Cozzolino, D.; Cynkar, W.; Dambergs, R.; Janik, L.; Gishen, M. Effect of both homogenisation and storage on the spectra of red grapes and on the measurement of total anthocyanins, total soluble solids and pH by visual near infrared spectroscopy. J. Near Infrared Spectrosc. 2005, 13, 213–223. [Google Scholar] [CrossRef]
  72. Martelo-Vidal, M.J.; Vazquez, M. Evaluation of ultraviolet, visible, and near infrared spectroscopy for the analysis of wine compounds. Czech J. Food Sci. 2014, 32, 37–47. [Google Scholar] [CrossRef] [Green Version]
  73. Haboudane, D.; Miller, J.R.; Pattey, E.; Zarco-Tejada, P.J.; Strachan, I.B. Hyperspectral vegetation indices and novel algorithms for predicting green lai of crop canopies: Modeling and validation in the context of precision agriculture. Remote Sens. Environ. 2004, 90, 337–352. [Google Scholar] [CrossRef]
  74. Hall, A.; Lamb, D.; Holzapfel, B.; Louis, J. Optical remote sensing applications in viticulture—A review. Aust. J. Grape Wine Res. 2002, 8, 36–47. [Google Scholar] [CrossRef]
  75. Lamb, D.W.; Weedon, M.M.; Bramley, R.G.V. Using remote sensing to predict grape phenolics and colour at harvest in a Cabernet sauvignon vineyard: Timing observations against vine phenology and optimising image resolution. Aust. J. Grape Wine Res. 2004, 10, 46–54. [Google Scholar] [CrossRef] [Green Version]
  76. Iqbal, M.; Jura-Morawiec, J.; WŁoch, W.; Mahmooduzzafar. Foliar characteristics, cambial activity and wood formation in Azadirachta indica A. Juss. As affected by coal–smoke pollution. Flora Morphol. Distrib. Funct. Ecol. Plants 2010, 205, 61–71. [Google Scholar] [CrossRef]
  77. Nighat, F.; Iqbal, M. Stomatal conductance, photosynthetic rate, and pigment content in Ruellia tuberosa leaves as affected by coal-smoke pollution. Biol. Plant. 2000, 43, 263–267. [Google Scholar] [CrossRef]
  78. Ren, S.; Chen, X.; An, S. Assessing plant senescence reflectance index-retrieved vegetation phenology and its spatiotemporal response to climate change in the inner mongolian grassland. Int. J. Biometeorol. 2017, 61, 601–612. [Google Scholar] [CrossRef] [PubMed]
  79. Sandermann, H.; Ernst, D.; Heller, W.; Langebartels, C. Ozone: An abiotic elicitor of plant defence reactions. Trends Plant Sci. 1998, 3, 47–50. [Google Scholar] [CrossRef]
  80. Bellincontro, A.; Catelli, C.; Cotarella, R.; Mencarelli, F. Postharvest ozone fumigation of Petit Verdot grapes to prevent the use of sulfites and to increase anthocyanin in wine. Aust. J. Grape Wine Res. 2017, 23, 200–206. [Google Scholar] [CrossRef]
  81. Castellarin, S.D.; Matthews, M.A.; Di Gaspero, G.; Gambetta, G.A. Water deficits accelerate ripening and induce changes in gene expression regulating flavonoid biosynthesis in grape berries. Planta 2007, 227, 101–112. [Google Scholar] [CrossRef]
  82. Allied Scientific Pro. Lighting Passport the World’s First Smart Handheld Spectrometer. 2017. Available online: https://lightingpassport.alliedscientificpro.com (accessed on 10 August 2020).
  83. Chandraratne, M.; Kulasiri, D.; Samarasinghe, S. Classification of lamb carcass using machine vision: Comparison of statistical and neural network analyses. J. Food Eng. 2007, 82, 26–34. [Google Scholar] [CrossRef]
  84. Al-Alawi, S.M.; Abdul-Wahab, S.A.; Bakheit, C.S. Combining principal component regression and artificial neural networks for more accurate predictions of ground-level ozone. Environ. Model. Softw. 2008, 23, 396–403. [Google Scholar] [CrossRef]
  85. Diamantopoulou, M.J.; Milios, E. Modelling total volume of dominant pine trees in reforestations via multivariate analysis and artificial neural network models. Biosyst. Eng. 2010, 105, 306–315. [Google Scholar] [CrossRef]
Figure 1. Smoke treatments were applied to grapevines using a purpose-built smoke tent; grapevines were enclosed in the tent and exposed to smoke derived from the combustion of barley straw (a,b).
Figure 1. Smoke treatments were applied to grapevines using a purpose-built smoke tent; grapevines were enclosed in the tent and exposed to smoke derived from the combustion of barley straw (a,b).
Sensors 20 05099 g001
Figure 2. Two-layer feedforward network with ten hidden neurons and sigmoid function for the three classification models: (a). microPHAZIRTM leaf model (Model 1) and (b). microPHAZIRTM berry model (Model 2). Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Figure 2. Two-layer feedforward network with ten hidden neurons and sigmoid function for the three classification models: (a). microPHAZIRTM leaf model (Model 1) and (b). microPHAZIRTM berry model (Model 2). Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Sensors 20 05099 g002
Figure 3. Raw leaf absorbance (a) and second derivative spectra (b) measured with the microPHAZIRTM near-infrared (NIR) analyzer for the different smoke and misting treatments. Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Figure 3. Raw leaf absorbance (a) and second derivative spectra (b) measured with the microPHAZIRTM near-infrared (NIR) analyzer for the different smoke and misting treatments. Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Sensors 20 05099 g003
Figure 4. Raw berry absorbance (a) and second derivative spectra (b) measured with the microPHAZIRTM NIR analyzer for the different smoke and misting treatments. Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Figure 4. Raw berry absorbance (a) and second derivative spectra (b) measured with the microPHAZIRTM NIR analyzer for the different smoke and misting treatments. Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Sensors 20 05099 g004aSensors 20 05099 g004b
Figure 5. Principal component analysis (PCA) for the microPHAZIRTM leaf absorbance values between 1600–1800 nm (a) and loadings (b). Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Figure 5. Principal component analysis (PCA) for the microPHAZIRTM leaf absorbance values between 1600–1800 nm (a) and loadings (b). Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Sensors 20 05099 g005
Figure 6. Principal component analysis (PCA) for the microPHAZIRTM berry absorbance values between 1600–2396 nm (a) and loadings (b). Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Figure 6. Principal component analysis (PCA) for the microPHAZIRTM berry absorbance values between 1600–2396 nm (a) and loadings (b). Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Sensors 20 05099 g006
Figure 7. Receiver operating characteristic (ROC) curves for the two models developed (a) the microPHAZIRTM leaf model, (b) the microPHAZIRTM berry model. Colored lines represent the different smoking treatments. Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Figure 7. Receiver operating characteristic (ROC) curves for the two models developed (a) the microPHAZIRTM leaf model, (b) the microPHAZIRTM berry model. Colored lines represent the different smoking treatments. Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke.
Sensors 20 05099 g007
Table 1. Calculations for the spectral indices investigated in this study.
Table 1. Calculations for the spectral indices investigated in this study.
Index NameIndex AbbreviationEquationReferences
Normalized difference vegetation indexNDVI ( I 780 I 660 ) ( I 780 + I 660 ) [56,57]
Normalized anthocyanin indexNAI ( I 780 I 570 ) ( I 780 + I 570 ) [56,57]
Carotenoid reflectance indexCRI550 1 I 510 1 I 550 [63,64]
Carotenoid reflectance indexCRI700 1 I 510 1 I 700 [65]
Plant senescence reflectance indexPSRI I 680 I 500 I 750 [59]
Table 2. Gas exchange parameters measured for the different smoke treatments.
Table 2. Gas exchange parameters measured for the different smoke treatments.
Smoke TreatmentE (mmol m−2 s−1)gs (mol m−2 s−1)A (µmol m−2 s−1)
MeanSDMeanSDMeanSD
C2.48 a0.700.15 a0.0510.77 a3.46
CM2.31 a0.540.15 a0.059.66 ab2.31
HS1.43 b0.620.06 c0.035.59 d2.8
HSM2.06 a0.440.10 b0.038.15 bc1.97
LS2.18 a0.780.08 bc0.037.01 cd2.42
Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke; SD = standard deviation. Means followed by different letters are significantly different based on Fisher least significant difference (LSD) post hoc test (α = 0.05).
Table 3. Means and standard deviation (SD) of spectral indices calculated for leaves and berries.
Table 3. Means and standard deviation (SD) of spectral indices calculated for leaves and berries.
TreatmentLeaf Berry
NDVINAIPSRICRI500CRI700 NAIPSRI
MeanSDMeanSDMeanSDMeanSDMeanSDMeanSDMeanSD
CM0.85 a0.100.77 a0.110.00 b0.010.70 b0.640.82 b0.79----
C0.84 ab0.0820.74 ab0.110.01 b0.020.67 b0.780.77 b0.870.80 b0.070.02 a0.02
HS 0.72 b0.500.64 b0.490.07 a0.191.20 a0.240.82 b0.620.88 a0.040.00 b0.00
HSM0.87 a0.110.79 a0.110.00 b0.020.48 b0.060.58 b0.450.75 b0.10−0.02 c0.00
LS0.92 a0.040.84 a0.080.00 b0.011.45 a1.081.76 a1.400.87 a0.050.02 a0.01
Abbreviations: C = control without misting; CM = control with misting; HS = high density smoke without misting; HSM = high density smoke with misting; and LS = low density smoke. Means followed by different letters are statistically significant based on Fisher’s least significant difference (LSD) post hoc test (α = 0.05).
Table 4. Statistical results for the artificial neural networks pattern recognition models. Model 1: microPHAZIRTM for leaves, and Model 2: microPHAZIRTM for berries. Performance is based on means squared error (MSE).
Table 4. Statistical results for the artificial neural networks pattern recognition models. Model 1: microPHAZIRTM for leaves, and Model 2: microPHAZIRTM for berries. Performance is based on means squared error (MSE).
StageSamples (n)Accuracy %Error %Performance (MSE)
Model 1
Training113110000.00
Validation 24394.25.80.02
Testing24392.67.40.02
Overall161798.02-
Model 2
Training37810000.00
Validation 8192.67.40.03
Testing8190.19.90.04
Overall54097.42.6-

Share and Cite

MDPI and ACS Style

Summerson, V.; Gonzalez Viejo, C.; Szeto, C.; Wilkinson, K.L.; Torrico, D.D.; Pang, A.; De Bei, R.; Fuentes, S. Classification of Smoke Contaminated Cabernet Sauvignon Berries and Leaves Based on Chemical Fingerprinting and Machine Learning Algorithms. Sensors 2020, 20, 5099. https://doi.org/10.3390/s20185099

AMA Style

Summerson V, Gonzalez Viejo C, Szeto C, Wilkinson KL, Torrico DD, Pang A, De Bei R, Fuentes S. Classification of Smoke Contaminated Cabernet Sauvignon Berries and Leaves Based on Chemical Fingerprinting and Machine Learning Algorithms. Sensors. 2020; 20(18):5099. https://doi.org/10.3390/s20185099

Chicago/Turabian Style

Summerson, Vasiliki, Claudia Gonzalez Viejo, Colleen Szeto, Kerry L. Wilkinson, Damir D. Torrico, Alexis Pang, Roberta De Bei, and Sigfredo Fuentes. 2020. "Classification of Smoke Contaminated Cabernet Sauvignon Berries and Leaves Based on Chemical Fingerprinting and Machine Learning Algorithms" Sensors 20, no. 18: 5099. https://doi.org/10.3390/s20185099

APA Style

Summerson, V., Gonzalez Viejo, C., Szeto, C., Wilkinson, K. L., Torrico, D. D., Pang, A., De Bei, R., & Fuentes, S. (2020). Classification of Smoke Contaminated Cabernet Sauvignon Berries and Leaves Based on Chemical Fingerprinting and Machine Learning Algorithms. Sensors, 20(18), 5099. https://doi.org/10.3390/s20185099

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop