Daily Mapping of 30 m LAI and NDVI for Grape Yield Prediction in California Vineyards

Wine grape quality and quantity are affected by vine growing conditions during critical phenological stages. Field observations of vine growth stages are too sparse to fully capture the spatial variability of vine conditions. In addition, traditional grape yield prediction methods are time consuming and require large amount grape samples. Remote sensing data provide detailed spatial and temporal information regarding vine development that is useful for vineyard management. In this study, Landsat surface reflectance products from 2013 and 2014 were used to map satellite-based Normalized Difference Vegetation Index (NDVI) and leaf area index (LAI) over two Vitis vinifera L. cv. Pinot Noir vineyards in California, USA. The spatial correlation between grape yield maps and the interpolated daily time series (LAI and NDVI) was quantified. NDVI and LAI were found to have similar performance as a predictor of spatial yield variability, providing peak correlations of 0.8 at specific times during the growing season, and the timing of this peak correlation differed for the two years of study. In addition, correlations with maximum and seasonal-cumulative vegetation indices were also evaluated, and showed slightly lower correlations with the observed yield maps. Finally, the within-season grape yield predictability was examined using a simple strategy in which the relationship between grape yield and vegetation indices were calibrated with limited ground measurements. This strategy has a strong potential to improve the accuracy and efficiency of yield estimation in comparison with traditional approaches used in the wine grape growing industry.


Introduction
Over 99 percent of grapes grown commercially in the United States come from California, which accounts for nearly 90 percent of American wine production valued at approximately $6 billion from 918,000 acres of vineyards in 2015 [1].The ability to accurately and efficiently monitor vine development and estimate grape yields within season will have significant benefit to the wine industry such as efficient harvest organization, regional pricing negotiations, and tank space allocation for vintage and marketing strategies [2].An over-cropped vine has an insufficient exposed leaf area relative to the retained fruit weight, leading to delayed ripening, lower fruit quality and an eventual decrease in vine size, limiting future yield potential and reducing cold hardiness of both buds and wood.An under-cropped vine has an excess of exposed leaf area relative to the weight of the retained fruit crop.It will be more prone to fungal diseases and will also have reduced fruit quality [3].
Viticulturists have developed several systems for estimating yield, but all are based on sampling: (a) the number of grape-bearing vines per acre; (b) the number of grape clusters per vine, and most notably; and (c) cluster weight [2,4,5].The sample size and distribution is determined by row and vine spacing, which means that the spatial variability should be considered.Growers must determine whether, when and how much fruit they should remove at thinning time.Accordingly, daily monitoring of vineyard conditions is required for optimal grape management, but ground-based measurements are not likely to be representative of whole-field conditions and are costly to install and maintain, especially for large scale and distributed production systems.
Remotely sensed data offer strong advantages over other monitoring techniques by providing a timely, synoptic, and up-to-date overview of actual crop growing conditions over large areas at multiple stages during the growing season.It has been found that there are strong relationships between satellite-based vegetation indices such as NDVI and vine development across vineyards [6][7][8].
In addition, remote sensing data have also been used to extract grapevine shape/size [9], to predict grape phenolics and color [10], and to discriminate different varieties within vineyards [11,12].
In addition to NDVI, a direct biophysical parameter for monitoring vine conditions is leaf area index (LAI) which provides specific physical information regarding canopy functioning; therefore, high-quality LAI retrieved from remotely sensed data is a valuable asset [13,14].Many approaches have been developed to estimate LAI using Landsat.These approaches can be classified into two basic categories: (1) using empirical relationships between LAI and vegetation indices [15,16]; and (2) through the inversion of a physical radiative transfer (RT) model [17].The inversion of a RT model could be implemented by seeking the solution with minimum cost iteratively, using a lookup table (LUT) method [13] or neural network (NN) algorithm [18].Although the NDVI approach is simple, it is usually sensitive to soil and atmospheric conditions as well as measurement geometries and spatial resolutions, and thus no unique relationship between LAI and NDVI is universally applicable [19].The conventional RT model inversion with an iterative process is both time consuming and hard to extend to regional and global scales, although it is more universal and may have a better accuracy.In this study we evaluate a reference-based approach to mapping LAI with Landsat developed by Gao et al. [20], which uses Landsat surface reflectance (SR) maps to downscale Moderate Resolution Imaging Spectroradiometer (MODIS) LAI products.This method has been successfully incorporated into the DisALEXI ET model [21] to map energy fluxes at 30 m resolution [22][23][24] and to provide a MODIS-consistent LAI data product from Landsat imagery.
The objective of this paper is to evaluate the relationships between grape yields and Landsat-derived LAI and NDVI time-series at field scale (30 m) for improving grape yield prediction through combining the limited field sampling with remote sensing data.While previous studies have shown the correlation between crop yields and satellite-based vegetation indices (VI) [6,[25][26][27], vineyard yield estimation is more complex due to intensive management of the grape and vine development to maximize production goals, which may vary from year to year based on climate, economics, and market demand.For example, the vineyard manager might remove some clusters early in the season in order to drive quality in the remaining clusters [28].Although it is difficult to establish a universal model to make year-to-year prediction of at-harvest yield using remote-sensing based indices, it is still very important for the grape growers to predict yield before harvest by the sampling method as described above.
This study focuses on the spatial variability of grape yields for 2013 and 2014 achieved in two adjacent Pinot Noir vineyards located in the Central Valley of California, USA.The spatial correlations between yield, and NDVI and LAI were evaluated and analyzed.Based on the findings, a simple grape yield prediction strategy, combing satellite based VI and limited within-field samples, is proposed for improving grape yield estimation before harvest.

Landsat-Resolution LAI Estimation
Landsat data (30 m) provide a unique opportunity to characterize crop conditions and water stress at field or sub-field scales.In addition to spatial resolution requirements, high temporal resolution is also needed to capture the dynamics of crop growth through time, especially at critical growth stages [29].While a single Landsat satellite has a relatively long revisit cycle of 16 days, at present we have access to both Landsat 7 and 8 data, which together can provide acquisitions at one location as often as every 8 days, depending on cloud cover.
A regression tree method is used for retrieving Landsat LAI using homogeneous and high quality MODIS LAI as a reference [20].The MODIS LAI product (500 m) has been validated and analyzed with independent field measurements [30,31].In this approach, Landsat surface reflectances are aggregated to match coarse-resolution MODIS data, extracted from MODIS LAI products over the Landsat scene of interest.Only high-quality MODIS LAI samples retrieved from the main algorithm (based on RT model) are selected based on the product quality flags.The selected MODIS LAI samples are further screened to exclude pixels with high sub-pixel variability in Landsat reflectance, thus ensuring that only homogeneous MODIS pixels are used in the training process.As LAI has a nonlinear relation to spectral reflectance and NDVI, a multi-variant regression tree approach (cubist by RuleQuest) is used in the processing.LAI and surface reflectance samples at MODIS pixel resolution are used to train the MODIS resolution model.The derived regression trees are finally applied to the Landsat surface reflectances at their native 30-m resolution to produce Landsat-scale LAI maps.

TIMESAT
The TIMESAT program package was designed primarily for smoothing and filtering time-series of satellite data [32,33].TIMESAT has been used successfully to analyze time-series VI data from the Advanced Very High Resolution Radiometer (AVHRR) [33,34] and from MODIS [35][36][37].In this study, we used an adaptive Savitzky-Golay (SG) filtering function supplied by TIMESAT to interpolate Landsat-retrieved NDVI and LAI maps from clear days to daily timesteps.

Yield Correlation Analyses
Gridded yield data collected within-field (as described in Section 3.4) were aggregated to 30 m resolution consistent with the Landsat data grids.The spatial correlation between yield VI was quantified using the Pearson correlation coefficient (R).To identify an optimal window during the growing season where an index is the best correlated to the yield at harvest, the daily time-series indices were accumulated and summarized over different intervals (from 1-day to 120-day) prior to computation of correlations.A 1-day interval gives a daily comparison, while a 120-day interval corresponds to almost the entire growing season, which starts from late March/early April and ends at harvest at the end of August/early September.The nominal growing season was defined as Day of Year (DOY) 91 (1 April) to 273 (30 September).

Grape Yield Prediction Model
Using the satellite index from the time interval identified by the methods in Section 2.3 as providing the optimal correlation with yields, a simple operational yield prediction approach is proposed and evaluated.First, the location of the maximum and the minimum satellite index pixel values from the optimal correlation date are located.At these locations within the field, a sampling box is defined where the vegetation index is measured in situ on the optimal correlation date.Then, grape weight and cluster areal density at the same sampling box can be measured in situ before harvest.Finally, a linear function is fitted based on the index and grape weights measured from the field (only two samples), and used to estimate production field-wide based on the satellite index image on the optimal correlation date.Compared to traditional methods which have to choose enough samples to represent the large area, the grape yield prediction model proposed here only needs two samples, which will greatly improve the grape yield prediction efficiency and reduce sampling costs.In this study, the maximum and minimum index and corresponding yield (grape weight and density) were extracted from the images representing the samples measured from fields, and then used to test the model.

Study Area
The Central Valley of California has abundant sunshine, large day-and-night temperature differences, and an evaporative demand ranging from 889 to 1270 mm of water during the growing season, making it a primary wine grape producing area in the USA.Our study area extended over two Vitis vinifera L. cv.Pinot Noir vineyards located near Lodi in central California as Figure 1 shown [22,38].The northern field (Site 1) has an area of about 35 ha and is a more mature vineyard (9-10 years old in 2016), while the southern field (Site 2) is smaller in size, at about 21 ha and 6-7 years old.The height of the vines ranges between 2 and 2.5 m, with row spacing approximately 3.35 m and an average vine spacing along the row of 1.52 m.Both fields have an east-west row orientation.Variations in soil texture, vine age and health result in a fairly wide range in vine biomass across the two fields.The vines typically budburst in mid-late March and grow through August before the grapes are harvested in early September.When winter rains and soil moisture are adequate, a grass cover crop flourishes early in the growing season in the inter-row until becoming senescent starting in early May, which is typically the beginning of the dry season.The cover crop helps to remove excess soil moisture from winter rains before vine budbreak.Drip irrigation typically commences sometime early in June and is continued until harvest if necessary.Decisions to begin irrigation in a field for a season are typically triggered by different factors, including visual assessment of canopy water stress, spot measurements of leaf water potential with gas pressure chambers, soil moisture, and upcoming weather conditions.
values from the optimal correlation date are located.At these locations within the field, a sampling box is defined where the vegetation index is measured in situ on the optimal correlation date.Then, grape weight and cluster areal density at the same sampling box can be measured in situ before harvest.Finally, a linear function is fitted based on the index and grape weights measured from the field (only two samples), and used to estimate production field-wide based on the satellite index image on the optimal correlation date.Compared to traditional methods which have to choose enough samples to represent the large area, the grape yield prediction model proposed here only needs two samples, which will greatly improve the grape yield prediction efficiency and reduce sampling costs.In this study, the maximum and minimum index and corresponding yield (grape weight and density) were extracted from the images representing the samples measured from fields, and then used to test the model.

Study Area
The Central Valley of California has abundant sunshine, large day-and-night temperature differences, and an evaporative demand ranging from 889 to 1270 mm of water during the growing season, making it a primary wine grape producing area in the USA.Our study area extended over two Vitis vinifera L. cv.Pinot Noir vineyards located near Lodi in central California as Figure 1 shown [22,38].The northern field (Site 1) has an area of about 35 ha and is a more mature vineyard (9-10 years old in 2016), while the southern field (Site 2) is smaller in size, at about 21 ha and 6-7 years old.The height of the vines ranges between 2 and 2.5 m, with row spacing approximately 3.35 m and an average vine spacing along the row of 1.52 m.Both fields have an east-west row orientation.Variations in soil texture, vine age and health result in a fairly wide range in vine biomass across the two fields.The vines typically budburst in mid-late March and grow through August before the grapes are harvested in early September.When winter rains and soil moisture are adequate, a grass cover crop flourishes early in the growing season in the inter-row until becoming senescent starting in early May, which is typically the beginning of the dry season.The cover crop helps to remove excess soil moisture from winter rains before vine budbreak.Drip irrigation typically commences sometime early in June and is continued until harvest if necessary.Decisions to begin irrigation in a field for a season are typically triggered by different factors, including visual assessment of canopy water stress, spot measurements of leaf water potential with gas pressure chambers, soil moisture, and upcoming weather conditions.

Ground Measurement Data
Two selected Pinot Noir vineyards were part of the Grape Remote sensing Atmospheric Profiling and Evapotranspiration eXperiment (GRAPEX) [22,38], which aims to combine in-situ and remotely-sensed data to investigate the effects of canopy structure and row orientation of energy and moisture exchange processes within and above the vine canopy.A critical component of GRAPEX has been to work closely and collaboratively with scientist at Ernest & Julio Gallo Winery to ensure that both the experimental data and models can address critical operational needs, including improving grape quality, yield estimation and irrigation strategies.GRAPEX measurements include ground, airborne and satellite remote sensing, surface energy balance, turbulence and mean profile measurements of wind temperature and water vapor focusing on above, below and between vine canopy, surface and sub-surface soil moisture, along with ground-based biophysical measurements of LAI and leaf-level conductance, transpiration and photosynthesis.Several intensive observation periods (IOPs) were conducted in 2013-2014 to collect biophysical measurements of the vine canopy during different phenological stages.The IOPs were scheduled to be coincident with Landsat 8 overpass times to provide validation of satellite-derived variables.The date and the number of measurements of the IOPs were listed in Table 1.The LAI measurements were made using a Li-Cor LAI-2000 (or 2200) instrument along several transects, collecting four below-canopy samples and one above-canopy reading to obtain an estimate of the effective LAI of the vine/inter-row system.The four below-canopy samples were collected: (1) under the vine and directly in-line with the vine row; (2) one-quarter of the distance from the row and the adjacent row of vines; (3) half way between rows in the center of the inter-row; and (4) three-quarters of the distance from the row to the adjacent row.Measurements were made with and without including the inter-row cover crop, which had a very minor impact on the LAI measurements for the field IOPs in June through August.Vines grow in the early spring and rapidly expand in leaf area mid to late May, and the observed LAI increased throughout the season until harvest in early September.

Landsat and MODIS Data Products
Landsat 7 (L7) and Landsat 8 (L8) atmospherically corrected surface reflectance products over the study area for the period from 2013-2014 were directly downloaded from the USGS EarthExplorer website [39].A total of 166 scenes were processed (see Figure 2).The study area lies in the overlap between path 43/row 33 and path 44/row 33, and this fortuitous location-along with low cloud cover during growing season-makes it possible to collect sufficient clear Landsat scenes to support interpolation to daily maps.Although Landsat 7 has striped gaps in all but the center of each scene missing, due to the SLC (Scan Line Corrector) failure that occurred since May of 2003, these gaps can be filled by interpolating Landsat 7 and 8 data using the TIMESAT program.
MODIS 4-day composite LAI products (MCD15A3, Collection 5; Myneni, 2014) for 2013-2014 over the study area were used.In contrast to the Terra-MODIS (MOD15A2), Aqua-MODIS (MYD15A2), and the combined Terra and Aqua MODIS (MCD15A2) versions of the 8-day LAI products, the 4-day composited MODIS LAI product (MCD15A3) has increased temporal frequency and helps to monitor phenology and associated rapid changes (that occur within less than a week), especially during transition periods (green-up and senescence).One MODIS tile (h08v05) was selected to cover both Landsat scenes.The LAI-reflectance samples from the entire Landsat scene over the entire year were used to train and build the regression trees.The regression trees were then applied to Landsat resolution reflectance to generate Landsat-scale LAI as Section 2.1 described.
Remote Sens. 2017, 9, 317 6 of 18 and helps to monitor phenology and associated rapid changes (that occur within less than a week), especially during transition periods (green-up and senescence).One MODIS tile (h08v05) was selected to cover both Landsat scenes.The LAI-reflectance samples from the entire Landsat scene over the entire year were used to train and build the regression trees.The regression trees were then applied to Landsat resolution reflectance to generate Landsat-scale LAI as Section 2.1 described.

Grape Yields
Yield maps for 2013 and 2014 harvest in both fields were provided by E & J Gallo Winery.Advanced Technology Viticulture (ATV, Joslin, Australia) yield monitor systems were installed on three self-propelled, over-the-row, trunk-shaking mechanical GH9000 AIM harvesters (Agricultural Industrial Manufacturing, Lodi, CA, USA).Yield monitors consist of a combined GPS receiver and antenna, as well as a specialized grape weighing and data recording system consisting of a load cell weight bridge, belt speed sensor, and data logger.Harvest data were transferred into a cleanup script written in R Studio software (RStudio Inc., Boston, MA, USA) to convert mass flow units into tons per hectare, eliminate outliers, and normalize data across harvesters.Yield data more than three standard deviations from the mean were removed to reduce scatter, which usually results in removal of less than 5% of the total yield data.After cleanup, the data were interpolated to 3 × 3 m resolution grids, and then aggregated up to the 30-m Landsat grid resolution.

Landsat NDVI and LAI
Figure 3 shows 30-m resolution maps of NDVI and estimated LAI over the GRAPEX fields on three selected days, with corresponding ground pictures taken from Site 1. Day 93, 173, and 253 separately correspond to the start, middle, and end of the growing season.LAI and NDVI have similar distributions, and the spatial pattern of the vegetation conditions within the fields and surrounding areas are clearly evident, supporting the hypothesis that field-scale management practices can be effectively informed by Landsat-scale resolution imagery.The maps for DOY 93 (early in the growing season) primarily reflect green biomass in the grass cover crop (see photo from Site 1).On this day, LAI in the north field is about 1.5-1.8 and NDVI is about 0.5-0.6,while the south field has LAI values that range from 1.6-2.5 and NDVI values of 0.5-0.7.DOY 173 is around the time of peak vine biomass, with NDVI around 0.6-0.7 and LAI of 1.7-2.2over the vineyards.DOY 253 (10 September) occurred before grape harvest, which is at the end of September, so the vegetation cover on this day is much lower than that on DOY 173.NDVI ranged 0.5-0.6 and LAI ranged 1.3-1.5 over GRAPEX vineyard fields on DOY 253.The lower vine LAI on this day is also reflected in the ground photograph in Figure 3.

Grape Yields
Yield maps for 2013 and 2014 harvest in both fields were provided by E & J Gallo Winery.Advanced Technology Viticulture (ATV, Joslin, Australia) yield monitor systems were installed on three self-propelled, over-the-row, trunk-shaking mechanical GH9000 AIM harvesters (Agricultural Industrial Manufacturing, Lodi, CA, USA).Yield monitors consist of a combined GPS receiver and antenna, as well as a specialized grape weighing and data recording system consisting of a load cell weight bridge, belt speed sensor, and data logger.Harvest data were transferred into a cleanup script written in R Studio software (RStudio Inc., Boston, MA, USA) to convert mass flow units into tons per hectare, eliminate outliers, and normalize data across harvesters.Yield data more than three standard deviations from the mean were removed to reduce scatter, which usually results in removal of less than 5% of the total yield data.After cleanup, the data were interpolated to 3 × 3 m resolution grids, and then aggregated up to the 30-m Landsat grid resolution.

Landsat NDVI and LAI
Figure 3 shows 30-m resolution maps of NDVI and estimated LAI over the GRAPEX fields on three selected days, with corresponding ground pictures taken from Site 1. Day 93, 173, and 253 separately correspond to the start, middle, and end of the growing season.LAI and NDVI have similar distributions, and the spatial pattern of the vegetation conditions within the fields and surrounding areas are clearly evident, supporting the hypothesis that field-scale management practices can be effectively informed by Landsat-scale resolution imagery.The maps for DOY 93 (early in the growing season) primarily reflect green biomass in the grass cover crop (see photo from Site 1).On this day, LAI in the north field is about 1.5-1.8 and NDVI is about 0.5-0.6,while the south field has LAI values that range from 1.6-2.5 and NDVI values of 0.5-0.7.DOY 173 is around the time of peak vine biomass, with NDVI around 0.6-0.7 and LAI of 1.7-2.2over the vineyards.DOY 253 (10 September) occurred before grape harvest, which is at the end of September, so the vegetation cover on this day is much lower than that on DOY 173.NDVI ranged 0.5-0.6 and LAI ranged 1.3-1.5 over GRAPEX vineyard fields on DOY 253.The lower vine LAI on this day is also reflected in the ground photograph in Figure 3.The NDVI and LAI maps in Figure 3 also demonstrate within field spatial variation in vine development over the sampling area.The central part of the northern field shows higher LAI and NDVI, indicating a region of more vigorous vine growth.Lower LAI and NDVI in the upper left of the northern field are associated with soil texture discontinuities.The major soil type in the two fields is Kimball silt loam but there are small areas within the field (such as in the upper left corner of the northern field) with silt loam and gravelly loam soils which have poor to moderate water storage capacity [22].The area of low LAI and NDVI at the top middle of the northern field is associated with a house, parking area, and storage sheds where equipment is kept for maintaining the fields.
NDVI and LAI averaged values (3 × 3 pixels) extracted from the sampling locations at Sites 1 and 2 (indicated in Figure 1) are shown in Figure 4. NDVI and LAI derived from L8 SR data are greater than those from L7, especially in the middle of the season where the difference is up to about 0.1 for NDVI and 0.5-1 for LAI.In addition, this difference is greater at Site 1 than at Site 2. Several studies have pointed out that surface reflectance differences between L7 and L8 result from quite different spectral response functions, resulting in higher values of NDVI from L8 than from L7 [40][41][42], which is consistent with our results.The reasons for these differences are complex, being dependent on the sensor differences in bands reflectance, view zenith angle, etc. [41], and will not be further explored in this study.
The GRAPEX study area is located within the overlap between adjacent WRS paths (43 and 44) and typically has low cloud cover during the grape growing season; therefore, there are sufficient Landsat 7 and 8 data available to provide a seasonable baseline, and then to form a smoothed The NDVI and LAI maps in Figure 3 also demonstrate within field spatial variation in vine development over the sampling area.The central part of the northern field shows higher LAI and NDVI, indicating a region of more vigorous vine growth.Lower LAI and NDVI in the upper left of the northern field are associated with soil texture discontinuities.The major soil type in the two fields is Kimball silt loam but there are small areas within the field (such as in the upper left corner of the northern field) with silt loam and gravelly loam soils which have poor to moderate water storage capacity [22].The area of low LAI and NDVI at the top middle of the northern field is associated with a house, parking area, and storage sheds where equipment is kept for maintaining the fields.
NDVI and LAI averaged values (3 × 3 pixels) extracted from the sampling locations at Sites 1 and 2 (indicated in Figure 1) are shown in Figure 4. NDVI and LAI derived from L8 SR data are greater than those from L7, especially in the middle of the season where the difference is up to about 0.1 for NDVI and 0.5-1 for LAI.In addition, this difference is greater at Site 1 than at Site 2. Several studies have pointed out that surface reflectance differences between L7 and L8 result from quite different spectral response functions, resulting in higher values of NDVI from L8 than from L7 [40][41][42], which is consistent with our results.The reasons for these differences are complex, being dependent on the sensor differences in bands reflectance, view zenith angle, etc. [41], and will not be further explored in this study.
The GRAPEX study area is located within the overlap between adjacent WRS paths (43 and 44) and typically has low cloud cover during the grape growing season; therefore, there are sufficient Landsat 7 and 8 data available to provide a seasonable baseline, and then to form a smoothed growing curve using the Savitzky-Golay (SG) filtering algorithms from TIMESAT.The Landsat 7 ETM+ images had less valid pixels due to the SLC failure especially in 2014.We used all valid pixels from Landsat 7 and 8.The original gaps in the Landsat 7 images were filled through the TIMESAT fitting process.In Figure 4, both NDVI and LAI show similar growth curves.The start of the growing season is around DOY 100 (middle of April), after which the leaf area increases quickly until the peak greenness period around DOY 150 (early of June).There is usually more than one peak in leaf biomass after DOY 150, especially for the south field inboth years.Although we do not have photos at Site 2 from 2013 and 2014 to examine the field conditions, it could be the result of pruning followed by vine re-growth which serves to shade the ripening grapes in the late summer from too much radiation.Overall, the LAI continues to decrease after DOY 150, and the end of the growing season is around DOY 250 (early September).
Remote Sens. 2017, 9, 317 8 of 18 growing curve using the Savitzky-Golay (SG) filtering algorithms from TIMESAT.The Landsat 7 ETM+ images had less valid pixels due to the SLC failure especially in 2014.We used all valid pixels from Landsat 7 and 8.The original gaps in the Landsat 7 images were filled through the TIMESAT fitting process.In Figure 4, both NDVI and LAI show similar growth curves.The start of the growing season is around DOY 100 (middle of April), after which the leaf area increases quickly until the peak greenness period around DOY 150 (early of June).There is usually more than one peak in leaf biomass after DOY 150, especially for the south field inboth years.Although we do not have photos at Site 2 from 2013 and 2014 to examine the field conditions, it could be the result of pruning followed by vine re-growth which serves to shade the ripening grapes in the late summer from too much radiation.Overall, the LAI continues to decrease after DOY 150, and the end of the growing season is around DOY 250 (early September).Observed values of LAI, collected in-field, are also included in Figure 4.In general, the observed LAI is in good agreement with the smoothed growing curves formed by TIMESAT, except the measurement on DOY 269 of 2014.Although there is no valid estimated Landsat LAI on DOY 269 of 2014, Landsat LAI from both DOY 262 and 278 is 1.4-1.5, which is higher than the observed LAI on DOY 269 (0.96 for north site and 0.86 for south site).The leaves during this period had obviously started to senesce compared to peak conditions, meanwhile the grass cover crop is significant which is reflected on the Landsat pixels.However, the observed LAI measured with the Li-Cor LAI-2200 on DOY 269 did not include the grass cover crop resulting in lower observed LAI.Observed values of LAI, collected in-field, are also included in Figure 4.In general, the observed LAI is in good agreement with the smoothed growing curves formed by TIMESAT, except the measurement on DOY 269 of 2014.Although there is no valid estimated Landsat LAI on DOY 269 of 2014, Landsat LAI from both DOY 262 and 278 is 1.4-1.5, which is higher than the observed LAI on DOY 269 (0.96 for north site and 0.86 for south site).The leaves during this period had obviously started to senesce compared to peak conditions, meanwhile the grass cover crop is significant which is reflected on the Landsat pixels.However, the observed LAI measured with the Li-Cor LAI-2200 on DOY 269 did not include the grass cover crop resulting in lower observed LAI.

Spatial Correlation between Yield and Daily NDVI and LAI
The daily NDVI and LAI images were used to compute spatial correlations with yield maps collected for the north and south vineyards in 2013 and 2014.Correlation coefficients are shown in Figure 5 as a function of day of VI retrieval.The R curves for NDVI and LAI have similar trends in all cases except for the south field for 2014, as discussed below.The growing season started later in 2013 than in 2014, so from DOY 91 to DOY 110 of 2013 the correlation of both NDVI and LAI is strongly affected by the grass cover crop resulting in poor correlations with grape yield in this period.In 2013, the NDVI and LAI correlation curve for both vineyards are similar, increasing to 0.8 in the early growing season (in early May as the grass cover crop senesced), then decreasing to 0.6 (north vineyard) and 0.4 (south vineyard) by the end of the growing season.Figure 6 shows the yield maps for 2013 and the evolution of LAI from DOY 110-140, during the period of maximum yield correlation, which occurred around DOY 140 for the north field and DOY 130 for the south field.From DOY 110-140, the LAI spatial pattern became more and more like the yield spatial pattern.This seems very early in the season to have the highest correlation with yield, which may be due to the management of the vineyard later in the season.
In 2014, high winter (especially February-March) temperatures and faster accumulation of growing degree days led to an early start of season.In the north vineyard (Figure 5c), yield-NDVI/LAI correlations increase quickly from DOY 91 (R = 0) to DOY 110 (about R = 0.5), then more slowly to about 0.7 at the end of growing season.The south vineyard (Figure 5d) exhibits a very different correlation behavior, with negative correlations (R = −0.5)at DOY 130, and a peak positive correlation of R = 0.5-0.6 around DOY 180.The negative correlations early in the season relate to patterns in the development and management of the cover crop in the south vineyard.As apparent in Figure 7, the grass LAI was higher in the northern part of that field-opposite to trends in grape production.It is clear that NDVI/LAI-based predictions of yield variability should be limited to the period of vine growth, where vine LAI dominates the surface reflectance signal.

Spatial Correlation between Yield and Daily NDVI and LAI
The daily NDVI and LAI images were used to compute spatial correlations with yield maps collected for the north and south vineyards in 2013 and 2014.Correlation coefficients are shown in Figure 5 as a function of day of VI retrieval.The R curves for NDVI and LAI have similar trends in all cases except for the south field for 2014, as discussed below.The growing season started later in 2013 than in 2014, so from DOY 91 to DOY 110 of 2013 the correlation of both NDVI and LAI is strongly affected by the grass cover crop resulting in poor correlations with grape yield in this period.In 2013, the NDVI and LAI correlation curve for both vineyards are similar, increasing to 0.8 in the early growing season (in early May as the grass cover crop senesced), then decreasing to 0.6 (north vineyard) and 0.4 (south vineyard) by the end of the growing season.Figure 6 shows the yield maps for 2013 and the evolution of LAI from DOY 110-140, during the period of maximum yield correlation, which occurred around DOY 140 for the north field and DOY 130 for the south field.From DOY 110-140, the LAI spatial pattern became more and more like the yield spatial pattern.This seems very early in the season to have the highest correlation with yield, which may be due to the management of the vineyard later in the season.
In 2014, high winter (especially February-March) temperatures and faster accumulation of growing degree days led to an early start of season.In the north vineyard (Figure 5c), yield-NDVI/LAI correlations increase quickly from DOY 91 (R = 0) to DOY 110 (about R = 0.5), then more slowly to about 0.7 at the end of growing season.The south vineyard (Figure 5d) exhibits a very different correlation behavior, with negative correlations (R = −0.5)at DOY 130, and a peak positive correlation of R = 0.5-0.6 around DOY 180.The negative correlations early in the season relate to patterns in the development and management of the cover crop in the south vineyard.As apparent in Figure 7, the grass LAI was higher in the northern part of that field-opposite to trends in grape production.It is clear that NDVI/LAI-based predictions of yield variability should be limited to the period of vine growth, where vine LAI dominates the surface reflectance signal.While LAI and NDVI have similar performance in terms of yield correlations in most cases, Figure 5 shows that NDVI correlations are higher for the southern field in 2014.The spatial distribution for both indices is similar during the period of peak correlation, but with a slightly higher level of noise/speckle in LAI likely due to artifacts in the MODIS downscaling.In general, NDVI may be a smoother and more reliable indicator of yield in vineyards, which are not likely to suffer from significant saturation effects due to the strong row structure.

Optimal Temporal Window for Yield Prediction Using NDVI and LAI Timeseries
To better identify the temporal window during which NDVI and LAI are most predictive of grape yields in this study, 2-D correlations plots are shown in Figure 8, computed between yields and indices averaged over windows ranging from one-day to 120-days (y-axis), with the window end date varying over the growing season (x-axis).Correlation strength is summarized by color, with red tones indicating index average windows with stronger correlations.The positive slope of the contour of maximum correlation indicates that longer averaging windows are accompanied by a   While LAI and NDVI have similar performance in terms of yield correlations in most cases, Figure 5 shows that NDVI correlations are higher for the southern field in 2014.The spatial distribution for both indices is similar during the period of peak correlation, but with a slightly higher level of noise/speckle in LAI likely due to artifacts in the MODIS downscaling.In general, NDVI may be a smoother and more reliable indicator of yield in vineyards, which are not likely to suffer from significant saturation effects due to the strong row structure.

Optimal Temporal Window for Yield Prediction Using NDVI and LAI Timeseries
To better identify the temporal window during which NDVI and LAI are most predictive of grape yields in this study, 2-D correlations plots are shown in Figure 8, computed between yields and indices averaged over windows ranging from one-day to 120-days (y-axis), with the window end date varying over the growing season (x-axis).Correlation strength is summarized by color, with red tones indicating index average windows with stronger correlations.The positive slope of the contour of maximum correlation indicates that longer averaging windows are accompanied by a While LAI and NDVI have similar performance in terms of yield correlations in most cases, Figure 5 shows that NDVI correlations are higher for the southern field in 2014.The spatial distribution for both indices is similar during the period of peak correlation, but with a slightly higher level of noise/speckle in LAI likely due to artifacts in the MODIS downscaling.In general, NDVI may be a smoother and more reliable indicator of yield in vineyards, which are not likely to suffer from significant saturation effects due to the strong row structure.

Optimal Temporal Window for Yield Prediction Using NDVI and LAI Timeseries
To better identify the temporal window during which NDVI and LAI are most predictive of grape yields in this study, 2-D correlations plots are shown in Figure 8, computed between yields and indices averaged over windows ranging from one-day to 120-days (y-axis), with the window end date varying over the growing season (x-axis).Correlation strength is summarized by color, with red tones indicating index average windows with stronger correlations.The positive slope of the contour of maximum correlation indicates that longer averaging windows are accompanied by a delay in the peak correlation signal.In general, there will be a tradeoff between early notification and confidence in yield forecasts [43].In addition, we found that the highest correlation with yields for both NDVI and LAI for both fields and both years came from 1-day interval.The best correlation information statistics are summarized in Table 2. Sibley et al. [25] also found that when predicting maize yield, only using the vegetation index from a single date performed better than the method using two to four dates per season.However, Esquerdo et al. [44] found that a full-season averaging window was optimal for estimating soybean yields in Brazil using NDVI.Mkhabela et al. [45] found the best time for making an accurate maize yield forecast was from the late January through late March depending on the agro-ecological region.Anderson et al. [43] analyzed the relationship between remotely sensed indices and crop yield in Brazil and also found that in some regions, some degree of additional time-averaging (moving upwards in the plots) helps to improve correlations with yields.The results from all these studies indicate that the time interval related to the best correlation depends not only on the crop type, but also on the climate/weather conditions and management practices.Regardless of this dispute, in this study, we find there is a strong relationship between grape yield and remotely sensed vegetation indices.The correlation between NDVI and yield ranged from 0.66 to 0.83, and for LAI the correlation with yield ranged from 0.53 to 0.82.The best correlation time is different for the north and south vineyard from two years due to different management and weather conditions early in the season.Irrigation start date and amount, as well as frequency and distribution throughout the season will have a significant impact on canopy and crop development.However, since the physical and chemical properties of the soil vary throughout the vineyard, the effect of any uniformly applied input will result in different spatial variability patterns.Canopy management practices, such as leafing and hedging, will also influence canopy density.Nonetheless, a high and stable correlation (greater than 0.6) over both years and fields was found over a common area between 25 June and 9 August.delay in the peak correlation signal.In general, there will be a tradeoff between early notification and confidence in yield forecasts [43].In addition, we found that the highest correlation with yields for both NDVI and LAI for both fields and both years came from 1-day interval.The best correlation information statistics are summarized in Table 2. Sibley et al. [25] also found that when predicting maize yield, only using the vegetation index from a single date performed better than the method using two to four dates per season.However, Esquerdo et al. [44] found that a full-season averaging window was optimal for estimating soybean yields in Brazil using NDVI.Mkhabela et al. [45] found the best time for making an accurate maize yield forecast was from the late January through late March depending on the agro-ecological region.Anderson et al. [43] analyzed the relationship between remotely sensed indices and crop yield in Brazil and also found that in some regions, some degree of additional time-averaging (moving upwards in the plots) helps to improve correlations with yields.The results from all these studies indicate that the time interval related to the best correlation depends not only on the crop type, but also on the climate/weather conditions and management practices.Regardless of this dispute, in this study, we find there is a strong relationship between grape yield and remotely sensed vegetation indices.The correlation between NDVI and yield ranged from 0.66 to 0.83, and for LAI the correlation with yield ranged from 0.53 to 0.82.The best correlation time is different for the north and south vineyard from two years due to different management and weather conditions early in the season.Irrigation start date and amount, as well as frequency and distribution throughout the season will have a significant impact on canopy and crop development.However, since the physical and chemical properties of the soil vary throughout the vineyard, the effect of any uniformly applied input will result in different spatial variability patterns.Canopy management practices, such as leafing and hedging, will also influence canopy density.Nonetheless, a high and stable correlation (greater than 0.6) over both years and fields was found over a common area between 25 June and 9 August.Figure 9 shows the scatter plots between NDVI and the yield on the best correlation date corresponding to Table 2.The yield range, especially for the south vineyard changed a lot between the two years, from 10-25 thousand kg/ha in 2013 to 30-50 in 2014.In addition, the slope of the fitting equation apparently increased in 2014 compared to that in 2013 for the both north and south vineyard, indicating that the yield variation became more sensitive to NDVI variation.In other words, the grape cluster density increased in 2014, which likely resulted from different management of the vineyard.However, the slope of the relationships for the north and south block is very close in the same year (66.8 and 60.6 in 2013, and 82 and 89.4 in 2014), indicating that the grape cluster density is similar in each year on both fields.
Remote Sens. 2017, 9, 317 12 of 18 Figure 9 shows the scatter plots between NDVI and the yield on the best correlation date corresponding to Table 2.The yield range, especially for the south vineyard changed a lot between the two years, from 10-25 thousand kg/ha in 2013 to 30-50 in 2014.In addition, the slope of the fitting equation apparently increased in 2014 compared to that in 2013 for the both north and south vineyard, indicating that the yield variation became more sensitive to NDVI variation.In other words, the grape cluster density increased in 2014, which likely resulted from different management of the vineyard.However, the slope of the relationships for the north and south block is very close in the same year (66.8 and 60.6 in 2013, and 82 and 89.4 in 2014), indicating that the grape cluster density is similar in each year on both fields.2. The x-axis represents NDVI and the y-axis represents the yield (tons/acre).Fitting function labeled on each window.

A Simple Calibrated Method for Estimating Field-Scale Yield Variations
The varying slopes in NDVI vs. yield apparent in Figure 9 demonstrates the difficulties involved in using NDVI data alone to predict yield in wine grape vineyards.Management activities and variability in weather conditions interfere with natural relationships between leaf area and fruit production, and these management decisions are influenced by weather, disease, and also by economic and marketing factors that cannot be remotely-sensed.Still, there is the potential to combine remote sensing with standard field sampling that is already integral to the operational yield projection process which may significantly improve on current procedures and efficiently enhance the spatial detail of vineyard yield projections.
As a simple initial experiment, yields were sampled from the full maps at the points of minimum and maximum NDVI determined within a given field for a given year, as determined on the date of optimal correlation.This is analogous to targeted field sampling, using NDVI maps at critical growth stages to identify sampling points.Projected yield maps were then developed using a linear yield-NDVI function, scaled between these limits.Table 3 lists the prediction model coefficients and error statistics, showing low relative errors from 6.9-11.5%.The predicted and observed yield maps show similar spatial patterns (Figure 10), reflecting primarily variations in soil type and vine health.2. The x-axis represents NDVI and the y-axis represents the yield (tons/acre).Fitting function labeled on each window.

A Simple Calibrated Method for Estimating Field-Scale Yield Variations
The varying slopes in NDVI vs. yield apparent in Figure 9 demonstrates the difficulties involved in using NDVI data alone to predict yield in wine grape vineyards.Management activities and variability in weather conditions interfere with natural relationships between leaf area and fruit production, and these management decisions are influenced by weather, disease, and also by economic and marketing factors that cannot be remotely-sensed.Still, there is the potential to combine remote sensing with standard field sampling that is already integral to the operational yield projection process which may significantly improve on current procedures and efficiently enhance the spatial detail of vineyard yield projections.
As a simple initial experiment, yields were sampled from the full maps at the points of minimum and maximum NDVI determined within a given field for a given year, as determined on the date of optimal correlation.This is analogous to targeted field sampling, using NDVI maps at critical growth stages to identify sampling points.Projected yield maps were then developed using a linear yield-NDVI function, scaled between these limits.Table 3 lists the prediction model coefficients and error statistics, showing low relative errors from 6.9-11.5%.The predicted and observed yield maps show similar spatial patterns (Figure 10), reflecting primarily variations in soil type and vine health.In practice, however, the date of optimal correlation will not be known a priori and may vary from year to year.A more generalized summary of the NDVI data must be used.Therefore, the correlations between yield and maximum and cumulative indices during the growing season (over DOY 110 to 250) and during optimal common period in Section 4.4 (over DOY 151 to 221) were also evaluated and compared as shown in Table 4. Neither index combination is as good as the correlation from the best date.For the maximum index, the correlation of the south field in 2014 is negative, because the maximum index is mostly from DOY 140-150 and it is negatively correlated during this period for the south vineyard in 2014 (Figure 5).Thus, the two types of cumulative index (DOY 110-250 and DOY 151-221) are more stable than the maximum index.According to this comparison, another strategy was attempted to predict the grape yield based on the cumulative NDVI.The yield prediction relative errors range from 10 to 18%, which is, as expected, a little larger than that in Table 3, but it is still better than the industry performance of around 30% [46].In practice, however, the date of optimal correlation will not be known a priori and may vary from year to year.A more generalized summary of the NDVI data must be used.Therefore, the correlations between yield and maximum and cumulative indices during the growing season (over DOY 110 to 250) and during optimal common period in Section 4.4 (over DOY 151 to 221) were also evaluated and compared as shown in Table 4. Neither index combination is as good as the correlation from the best date.For the maximum index, the correlation of the south field in 2014 is negative, because the maximum index is mostly from DOY 140-150 and it is negatively correlated during this period for the south vineyard in 2014 (Figure 5).Thus, the two types of cumulative index (DOY 110-250 and DOY 151-221) are more stable than the maximum index.According to this comparison, another strategy was attempted to predict the grape yield based on the cumulative NDVI.The yield prediction relative errors range from 10 to 18%, which is, as expected, a little larger than that in Table 3, but it is still better than the industry performance of around 30% [46].

Importance of High Spatiotemporal Resolution Remote Sensing Data
At the Landsat scale (30 m), phenology and water use can be differentiated by crop type and land management practice [47], making this an optimal scale for assessments of crop yield, water productivity, and moisture conditions [48][49][50].In this study, Landsat enabled us to investigate and analyze variations in vegetation conditions at sub-field scales.Meanwhile, the crop conditions changed quickly in critical growing stages even during a short time period as shown in Figure 5, which required a continuous timeseries of satellite data.To highlight this point, two pictures taken from the north vineyard during the rapid development stage (9 and 12 May 2013) are shown in Figure 11.Even though these are only three days apart, the leaf density of the vines on 12 May is obviously larger than that on 9 May (readers can compare the right row on the picture).

Importance of High Spatiotemporal Resolution Remote Sensing Data
At the Landsat scale (30 m), phenology and water use can be differentiated by crop type and land management practice [47], making this an optimal scale for assessments of crop yield, water productivity, and moisture conditions [48][49][50].In this study, Landsat enabled us to investigate and analyze variations in vegetation conditions at sub-field scales.Meanwhile, the crop conditions changed quickly in critical growing stages even during a short time period as shown in Figure 5, which required a continuous timeseries of satellite data.To highlight this point, two pictures taken from the north vineyard during the rapid development stage (9 and 12 May 2013) are shown in Figure 11.Even though these are only three days apart, the leaf density of the vines on 12 May is obviously larger than that on 9 May (readers can compare the right row on the picture).There are enough clear Landsat overpass images available in this study area to generate the timeseries NDVI/LAI.In other regions with higher frequency of cloud cover, data fusion methodologies, combining the high spatial resolution of Landsat with the daily temporal frequency of MODIS or similar moderate resolution systems, facilitates estimation of daily vegetation indices or evapotranspiration at sub-field scales [22,23,51].In addition, medium resolution data from Sentinel-2A or other sensors can also be incorporated for monitoring crop conditions at field scale using a data normalization approach [52].
Unlike cereals, vineyards are not continuous crop, which means the ground information from the inter-row other than the vine canopy would be captured by Landsat images as mixed pixels.The effect of grass cover crop could be seen from NDVI/LAI maps on DOY 93 in Figure 3.However, based on the investigation of ground pictures and communication with growers, the grass cover crop had little impact on NDVI/LAI estimation during June and August which is the critical growth stage.As described in Section 3.1, the grass cover crop flourishes in very early season and become senescent in early May.From the correlation curve of 2013 in Figure 5 we can see that this is the period (after DOY 120) when the correlation between NDVI/LAI and yield started to increase, indicating that the grass cover crop had less and less impact.The research here presented that 30 m resolution Landsat images have the ability to evaluate variation in vine conditions.Certainly, very high resolution images with meter-level resolution will help us to better understand the vineyard conditions, especially the inter-row ground information.GRAPEX was designed to collect very high resolution (2 m) images vineyards from an airborne platform at critical growth stages.Combining continuous Landsat 30 m data and 2 m airborne information would be useful for vineyard management.There are enough clear Landsat overpass images available in this study area to generate the timeseries NDVI/LAI.In other regions with higher frequency of cloud cover, data fusion methodologies, combining the high spatial resolution of Landsat with the daily temporal frequency of MODIS or similar moderate resolution systems, facilitates estimation of daily vegetation indices or evapotranspiration at sub-field scales [22,23,51].In addition, medium resolution data from Sentinel-2A or other sensors can also be incorporated for monitoring crop conditions at field scale using a data normalization approach [52].
Unlike cereals, vineyards are not continuous crop, which means the ground information from the inter-row other than the vine canopy would be captured by Landsat images as mixed pixels.The effect of grass cover crop could be seen from NDVI/LAI maps on DOY 93 in Figure 3.However, based on the investigation of ground pictures and communication with growers, the grass cover crop had little impact on NDVI/LAI estimation during June and August which is the critical growth stage.As described in Section 3.1, the grass cover crop flourishes in very early season and become senescent in early May.From the correlation curve of 2013 in Figure 5 we can see that this is the period (after DOY 120) when the correlation between NDVI/LAI and yield started to increase, indicating that the grass cover crop had less and less impact.The research here presented that 30 m resolution Landsat images have the ability to evaluate variation in vine conditions.Certainly, very high resolution images with meter-level resolution will help us to better understand the vineyard conditions, especially the inter-row ground information.GRAPEX was designed to collect very high resolution (2 m) images vineyards from an airborne platform at critical growth stages.Combining continuous Landsat 30 m data and 2 m airborne information would be useful for vineyard management.

Grape Yield Prediction
It has been found that there are significant correlations between NDVI and vine canopy vigor [7], and it is possible to develop a relationship between NDVI and grape yield [53].However, wine grape yield prediction is always a challenge not only because of climate factors, but also the trade-off between yield and quality [3], which means maximizing yield is not the only purpose of wine grape growers, but also market demand and quality have to be considered.In this case, there are many vineyard management techniques in the growing season that can affect yields.In the study of Cunha et al. [6], the relationship between NDVI and wine grape yield was significantly negative, explained by excessive vegetation growth of grapevines with dense canopies can causing shading and low bud fruitfulness [54].However, the relationship is positively correlated in this study, indicating that the relationship between remotely sensed vegetation indices and wine grape yield has different performance under different environmental conditions and management strategies.Meanwhile, Cunha et al. [6] found that NDVI during the middle of April of the previous season (initiation of inflorescence) showed the best correlation with wine grape yield and this relationship is stable and consistent over the vineyard in Portugal.However, in our study, the best relationship only exists in 2013 for both the north and south field, but not for 2014.This may result from management effects in 2014, which is also confirmed by the unusual spatial pattern change over the south vineyard in 2014 (Figure 7), which had much greater yield, even greater than the north vineyard.
Accordingly, it is very difficult to predict the best correlation day and then to find a universal grape yield prediction model based on remotely sensed data.Therefore, in this study, we evaluated the relationship between yield and different NDVI/LAI combinations and found that although cumulative NDVI/LAI is not as good as the NDVI/LAI from the optimal date, it is better than the maximum NDVI/LAI and more stable across both years and both vineyards.Thus cumulative NDVI/LAI is possible to use for yield prediction in practice.Moreover, the finding of a common period of strong correlation across both vineyards and both years should be analyzed using more data from 2015 and 2016, which may be a better solution for yield prediction.The simple grape yield prediction strategy proposed in this study is based on the inner-annual NDVI/LAI combination before harvest, which demonstrates potential for practical application in the future.
Grape yield and quality is also associated with water stress [55].Monitoring vineyard water use and improving water management are important purpose of GRAPEX.Preliminary results about spatio-temporal continuous 30 m resolution evapotranspiration (ET) maps of 2013 generated by a multi-satellite fusion program has been evaluated over GRAPEX vineyards, and the ET maps shows general spatial correspondence with yield maps [22].Therefore, in the next study, further insightsand improved grape yield prediction strategies could be expected by combining high spatio-temporal resolution NDVI/LAI and ET maps.

Conclusions
High temporal-spatial resolution remote sensing vegetation index data or imagery is critical to monitor field scale crop conditions.In this study, Landsat reflectance products from 2013 and 2014 were acquired over two vineyards in California, USA.Daily 30 m resolution NDVI and LAI was interpolated using adaptive Savitzky-Golay filtering function supplied by TIMESAT.The results showed that spatial variability of yield can be described by remote sensing data.There is apparently a correlation in the temporal variation curve for both the north and south vineyards and the best correlation (0.66-0.83) shows a specific time which is inconsistent for the two years evaluated in this study.The relationship between yield and maximum and cumulative index (NDVI and LAI) was also investigated.The correlation with yield for the cumulative index (0.64-0.74, except south vineyard in 2014) is better than the maximum index (0.6-0.7, except the south vineyard in 2014), and both are not as good as the index from the best correlation date.These remote sensing metrics could be used for grape yield prediction calibration using a few ground measurements with relative error of less than 20%.This can improve yield estimation and reduce time from traditional grape yield estimation approaches.To determine the optimal dates for grape yield predicting for the region, more vineyards in the region need to be included and analyzed using the same approach.In addition, inter-annual variability of yield is more complex and influenced by management practices in addition to other factors including weather.

Figure 1 .
Figure 1.Study area.The upper field confined by the dashed black line in the middle image is the north vineyard (Site 1) and the lower field confined by the dashed blue line is the south vineyard (Site 2).The North vineyard is planted with 9-10-year-old Pinot Noir vines, and the South vineyard is planted with 6-7-year-old Pinot Noir vines.The red circles correspond to the LAI measurement locations.

Figure 1 .
Figure 1.Study area.The upper field confined by the dashed black line in the middle image is the north vineyard (Site 1) and the lower field confined by the dashed blue line is the south vineyard (Site 2).The North vineyard is planted with 9-10-year-old Pinot Noir vines, and the South vineyard is planted with 6-7-year-old Pinot Noir vines.The red circles correspond to the LAI measurement locations.

Figure 2 .
Figure 2. Landsat 7 and 8 scenes collected over the study area from 2013-2014.Each dot represents a Landsat acquisition date.

Figure 2 .
Figure 2. Landsat 7 and 8 scenes collected over the study area from 2013-2014.Each dot represents a Landsat acquisition date.

Figure 3 .
Figure 3. Landsat NDVI and LAI maps, and ground pictures on the three selected dates in 2014.(North and south field of GRAPEX vineyard are confined by dashed red line.Camera site is located in the north field and labeled as red dot).

Figure 3 .
Figure 3. Landsat NDVI and LAI maps, and ground pictures on the three selected dates in 2014.(North and south field of GRAPEX vineyard are confined by dashed red line.Camera site is located in the north field and labeled as red dot).

Figure 4 .
Figure 4. Landsat NDVI, LAI and ground measured LAI for both sites in 2013 and 2014.

Figure 4 .
Figure 4. Landsat NDVI, LAI and ground measured LAI for both sites in 2013 and 2014.

Figure 5 .
Figure 5. Spatial correlation coefficients (R) between yield and the smoothed daily NDVI/LAI (solid lines in Figure 4) in 2013 (a,b) and 2014 (c,d) using all pixels in the north (a,c) and south (b,d) vineyard.

Figure 5 .
Figure 5. Spatial correlation coefficients (R) between yield and the smoothed daily NDVI/LAI (solid lines in Figure 4) in 2013 (a,b) and 2014 (c,d) using all pixels in the north (a,c) and south (b,d) vineyard.

Figure 8 .
Figure 8. Correlation of yield with NDVI and LAI plotted as a function of index averaging interval (y-axis) and end date (x-axis).The peak correlation area is identified by a contour.

Figure 8 .
Figure 8. Correlation of yield with NDVI and LAI plotted as a function of index averaging interval (y-axis) and end date (x-axis).The peak correlation area is identified by a contour.

Figure 9 .
Figure 9. Scatter plots between NDVI and yield on the best correlation date identified in Table2.The x-axis represents NDVI and the y-axis represents the yield (tons/acre).Fitting function labeled on each window.

Figure 9 .
Figure 9. Scatter plots between NDVI and yield on the best correlation date identified in Table2.The x-axis represents NDVI and the y-axis represents the yield (tons/acre).Fitting function labeled on each window.

Figure 10 .
Figure 10.Grape yield maps for north and south field of 2013 and 2014.

Figure 10 .
Figure 10.Grape yield maps for north and south field of 2013 and 2014.

Figure 11 .
Figure 11.Vine growing conditions on 9 and 12 May 2013 over the north vineyard.

Figure 11 .
Figure 11.Vine growing conditions on 9 and 12 May 2013 over the north vineyard.

Table 2 .
Date of optimal NDVI/LAI-yield correlation for the north and south vineyards.

Table 2 .
Date of optimal NDVI/LAI-yield correlation for the north and south vineyards.

Table 3 .
The prediction model coefficients (yield = a × NDVI + b) and error statistics.

Table 4 .
Correlation between yield and maximum and cumulative indexes (NDVI and LAI) during the entire growing season (DOY 110 to 250) and the optimal period (DOY 156 to 221).

Table 3 .
The prediction model coefficients (yield = a × NDVI + b) and error statistics.

Table 4 .
Correlation between yield and maximum and cumulative indexes (NDVI and LAI) during the entire growing season (DOY 110 to 250) and the optimal period (DOY 156 to 221).