Within-Field Relationships between Satellite-Derived Vegetation Indices, Grain Yield and Spike Number of Winter Wheat and Triticale

The aims of this study were to: (i) evaluate the relationships between vegetation indices (VIs) derived from Sentinel-2 imagery and grain yield (GY) and the number of spikes per square meter (SN) of winter wheat and triticale; (ii) determine the dates and plant growth stages when the above relationships were the strongest at individual field scale, thus allowing for accurate yield prediction. Observations of GY and SN were performed at harvest on six fields (three locations in two seasons: 2017 and 2018) in three regions of Poland, i.e., northeastern (A—Brożówka), central (B—Zdziechów) and southeastern Poland (C—Kryłów). Vegetation indices (Normalized Difference Vegetation Index (NDVI), Soil-Adjusted Vegetation Index (SAVI), modified SAVI (mSAVI), modified SAVI 2 (mSAVI2), Infrared Percentage Vegetation Index (IPVI), Global Environmental Monitoring Index (GEMI), and Ratio Vegetation Index (RVI)) calculated for sampling points from mid-March until mid-July, covering within-field soil and topographical variability, were included in the analysis. Depending on the location, the highest correlation coefficients (of about 0.6–0.9) for most of VIs with GY and SN were obtained about 4–6 weeks before harvest (from the beginning of shooting to milk maturity). Therefore, satellite-derived VIs are useful for the prediction of within-field cereal GY as well as SN variability. Information on GY, predicted together with the results for soil nutrient availability, is the basis for the formulation of variable fertilize rates in precision agriculture. All examined VIs were similarly correlated with GY and SN via the commonly used NDVI. The increase in NDVI by 0.1 unit was related to an average increase in GY by about 2 t ha−1.