Predicting wheat grain yield and spatial variability at field scale using a simple regression or a crop model in conjunction with Landsat images

• Published in: Computers and electronics in agriculture Vol. 159; pp. 75 - 83
• Main Authors: Gaso, Deborah V., Berger, Andrés G., Ciganda, Verónica S.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.04.2019; Elsevier BV
• Subjects: Accuracy; Agricultural practices; Agricultural production; Crop growth model; Crop yield; Harvesters; Identification methods; Image acquisition; Image enhancement; Land acquisition; Landsat; Landsat satellites; Leaf area index; Model accuracy; Optimization; Predictions; Regression models; Remote sensing; Satellite imagery; Temporal resolution; Thematic Mappers (LANDSAT); Vegetation; Wheat; Yield; HI; Nl (gN m−2); fr; NDVI; Crop growth model; SRM; ᵨNIR; ᵨRED; TB (g m−2 d−1); ei (MJ m−2 d−1); LUE; CMM; CI; Landsat; HI_ini; LAI (m−2 m−2); ᵨGREEN; RMSE; Rg (MJ m−2 d−1); VI; ec (MJ m−2 d−1); RRMSE; Yield; Wheat; Leaf area index
• ISSN: 0168-1699; 1872-7107
• DOI: 10.1016/j.compag.2019.02.026

•Model estimated yield was not more accurate than a simple regression estimation.•Assimilating remote sensing data into a crop model correctly estimated yield spatial variability.•Remote sensing-model data assimilation could be useful to recreate yield maps from archive data.•Adjusted model parameters during assimilation could provide useful derived metrics of the site. Early prediction of crop yields has been a challenge frequently resolved through the combination of remote sensing data and crop models. The aim of this study was to evaluate two different methods based on remote sensing data for predicting winter wheat (Triticum aestivum L.) yield at field scale. We compared the accuracy of: (i) a simple regression method between different vegetation indices at anthesis and grain yield, and (ii) a crop model method based on optimization of two parameters (specific leaf nitrogen and initial aboveground-biomass) using time series of vegetation indices. Vegetation indices were derived from Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Landsat-8 Operational Land Imager (OLI) images acquired for two growing seasons (2013, 2014) across 22 fields in south western Uruguay with an average size of 128 ha. At all sites, leaf area index (LAI) was measured during a field campaign, and grain yield was measured with yield monitors on harvesters. The simple regression method (SRM) achieved higher accuracy than the model-based method (CMM) for the estimation of yield at field scale (RMSE = 966 kg ha−1 and RMSE = 1532 kg ha−1, respectively). When deviations between observed and estimated yields were evaluated at pixel (30 × 30 m) level, the model-based method was better at detecting existing spatial variability in grain yield and at identifying areas of different yield potential. Even though both methods have limited utility to estimate yield at field scale with very high accuracy due to large RMSE, the methodologies are suitable to predict harvest volumes at large agricultural areas or at country level, and to construct synthetic yield maps reflecting within field variability. Higher temporal resolution of images would improve accuracy in estimating yield and spatial variability at field scale.