Assessing the Variability of Corn and Soybean Yields in Central Iowa Using High Spatiotemporal Resolution Multi-Satellite Imagery

The utility of remote sensing data in crop yield modeling has typically been evaluated at the regional or state level using coarse resolution (>250 m) data. The use of medium resolution data (10–100 m) for yield estimation at field scales has been limited due to the low temporal sampling frequency characteristics of these sensors. Temporal sampling at a medium resolution can be significantly improved, however, when multiple remote sensing data sources are used in combination. Furthermore, data fusion approaches have been developed to blend data from different spatial and temporal resolutions. This paper investigates the impacts of improved temporal sampling afforded by multi-source datasets on our ability to explain spatial and temporal variability in crop yields in central Iowa (part of the U.S. Corn Belt). Several metrics derived from vegetation index (VI) time-series were evaluated using Landsat-MODIS fused data from 2001 to 2015 and Landsat-Sentinel2-MODIS fused data from 2016 and 2017. The fused data explained the yield variability better, with a higher coefficient of determination (R2) and a smaller relative mean absolute error than using a single data source alone. In this study area, the best period for the yield prediction for corn and soybean was during the middle of the growing season from day 192 to 236 (early July to late August, 1–3 months before harvest). These findings emphasize the importance of high temporal and spatial resolution remote sensing data in agricultural applications.

[1]  C. Schaaf,et al.  Capturing rapid land surface dynamics with Collection V006 MODIS BRDF/NBAR/Albedo (MCD43) products , 2018 .

[2]  N. Silleos,et al.  The use of multi-temporal NDVI measurements from AVHRR data for crop yield estimation and prediction , 1993 .

[3]  A. Huete,et al.  Development of a two-band enhanced vegetation index without a blue band , 2008 .

[4]  C. Justice,et al.  A generalized regression-based model for forecasting winter wheat yields in Kansas and Ukraine using MODIS data , 2010 .

[5]  James E. McMurtrey,et al.  Relationship of spectral data to grain yield variation , 1980 .

[6]  A. Gitelson.  Wide Dynamic Range Vegetation Index for remote quantification of biophysical characteristics of vegetation. , 2004, Journal of plant physiology.

[7]  C. Justice,et al.  The Harmonized Landsat and Sentinel-2 surface reflectance data set , 2018, Remote Sensing of Environment.

[8]  B. Wardlow,et al.  Large-area crop mapping using time-series MODIS 250 m NDVI data: An assessment for the U.S. Central Great Plains , 2008 .

[9]  Mathew R. Schwaller,et al.  On the blending of the Landsat and MODIS surface reflectance: predicting daily Landsat surface reflectance , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[10]  María Amparo Gilabert,et al.  Use of NOAA-AVHRR NDVI data for environmental monitoring and crop forecasting in the Sahel. Preliminary results , 1992 .

[11]  A. Gitelson,et al.  Remote estimation of crop gross primary production with Landsat data , 2012 .

[12]  M. S. Rasmussen.  Developing simple, operational, consistent NDVI-vegetation models by applying environmental and climatic information. Part II: Crop yield assessment , 1998 .

[13]  Compton J. Tucker,et al.  Satellite remote sensing of total herbaceous biomass production in the Senegalese Sahel - 1980-1984 , 1985 .

[14]  P. C. Doraiswamya,et al.  Crop condition and yield simulations using Landsat and MODIS , 2004 .

[15]  M. S. Rasmussen.  Assessment of millet yields and production in northern Burkina Faso using integrated NDVI from the AVHRR. , 1992 .

[16]  D. Lobell,et al.  Towards fine resolution global maps of crop yields: Testing multiple methods and satellites in three countries , 2017 .

[17]  Chris Funk,et al.  Phenologically-tuned MODIS NDVI-based production anomaly estimates for Zimbabwe , 2009 .

[18]  Yang Yang,et al.  Remote Sensing of Irrigated Agriculture: Opportunities and Challenges , 2010, Remote. Sens..

[19]  C. Daughtry,et al.  Spectral estimates of absorbed radiation and phytomass production in corn and soybean canopies , 1992 .

[20]  Qingyuan Zhang,et al.  Monitoring interannual variation in global crop yield using long-term AVHRR and MODIS observations. , 2016, ISPRS journal of photogrammetry and remote sensing : official publication of the International Society for Photogrammetry and Remote Sensing.

[21]  Pierre-Majorique Léger,et al.  The early explanatory power of NDVI in crop yield modelling , 2008 .

[22]  Bernardo Friedrich Theodor Rudorff,et al.  Wheat yield estimation at the farm level using TM Landsat and agrometeorological data , 1991 .

[23]  M. S. Rasmussen.  Operational yield forecast using AVHRR NDVI data: reduction of environmental and inter-annual variability , 1997 .

[24]  J. I. Ortiz-Monasterio,et al.  Combining Field Surveys, Remote Sensing, and Regression Trees to Understand Yield Variations in an Irrigated Wheat Landscape , 2005, Agronomy Journal.

[25]  Xiaoxiong Xiong,et al.  Impact of Spatial Sampling on Continuity of MODIS–VIIRS Land Surface Reflectance Products: A Simulation Approach , 2017, IEEE Transactions on Geoscience and Remote Sensing.

[26]  Mutlu Ozdogan,et al.  A new methodology to map irrigated areas using multi-temporal MODIS and ancillary data: An application example in the continental US , 2008 .

[27]  Feng Gao,et al.  Relationships between the evaporative stress index and winter wheat and spring barley yield anomalies in the Czech Republic , 2016 .

[28]  Donghui Xie,et al.  Mapping Paddy Rice Area and Yields Over Thai Binh Province in Viet Nam From MODIS, Landsat, and ALOS-2/PALSAR-2 , 2018, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[29]  Marvin E. Bauer,et al.  Spectral estimates of solar radiation intercepted by corn canopies. , 1983 .

[30]  Martha C. Anderson,et al.  Toward mapping crop progress at field scales through fusion of Landsat and MODIS imagery , 2017 .

[31]  Anatoly A. Gitelson,et al.  MODIS-based corn grain yield estimation model incorporating crop phenology information , 2013 .

[32]  R. D. Jackson,et al.  Multidate spectral reflectance as predictors of yield in water stressed wheat and barley , 1981 .

[33]  C. Justice,et al.  Atmospheric correction of MODIS data in the visible to middle infrared: first results , 2002 .

[34]  D. Roy,et al.  Conterminous United States crop field size quantification from multi-temporal Landsat data , 2015 .

[35]  Martha C. Anderson,et al.  Landsat-8: Science and Product Vision for Terrestrial Global Change Research , 2014 .

[36]  Jianxi Huang,et al.  Improving the timeliness of winter wheat production forecast in the United States of America, Ukraine and China using MODIS data and NCAR Growing Degree Day information , 2015 .

[37]  Anatoly A. Gitelson,et al.  An evaluation of MODIS 8- and 16-day composite products for monitoring maize green leaf area index , 2012 .

[38]  Xiaolin Zhu,et al.  Spatiotemporal Fusion of Multisource Remote Sensing Data: Literature Survey, Taxonomy, Principles, Applications, and Future Directions , 2018, Remote. Sens..

[39]  E. Vermote,et al.  Combined use of Landsat-8 and Sentinel-2A images for winter crop mapping and winter wheat yield assessment at regional scale. , 2017, AIMS geosciences.

[40]  Matthias Drusch,et al.  Sentinel-2: ESA's Optical High-Resolution Mission for GMES Operational Services , 2012 .

[41]  C. Tucker,et al.  Remote Sensing of Total Dry-Matter Accumulation in Winter Wheat , 1981 .

[42]  R. Myneni,et al.  On the relationship between FAPAR and NDVI , 1994 .

[43]  David B. Lobell,et al.  The use of satellite data for crop yield gap analysis , 2013 .

[44]  J. Malingreau.  Global vegetation dynamics - Satellite observations over Asia , 1986 .

[45]  John R. Miller,et al.  Estimating crop stresses, aboveground dry biomass and yield of corn using multi-temporal optical data combined with a radiation use efficiency model , 2010 .

[46]  Feng Gao,et al.  The Evaporative Stress Index as an indicator of agricultural drought in Brazil: An assessment based on crop yield impacts , 2016 .

[47]  D. Lobell,et al.  Satellite-based assessment of yield variation and its determinants in smallholder African systems , 2017, Proceedings of the National Academy of Sciences.

[48]  David M. Johnson.  An assessment of pre- and within-season remotely sensed variables for forecasting corn and soybean yields in the United States , 2014 .

[49]  F. Hall,et al.  Global Crop Forecasting , 1980, Science.

[50]  Roberto Benedetti,et al.  On the use of NDVI profiles as a tool for agricultural statistics: The case study of wheat yield estimate and forecast in Emilia Romagna , 1993 .

[51]  Michele Meroni,et al.  Towards regional grain yield forecasting with 1km-resolution EO biophysical products: Strengths and limitations at pan-European level , 2015 .

[52]  Douglas K. Bolton,et al.  Forecasting crop yield using remotely sensed vegetation indices and crop phenology metrics , 2013 .

[53]  Shusen Wang,et al.  Crop yield forecasting on the Canadian Prairies using MODIS NDVI data , 2011 .

[54]  Alan H. Strahler,et al.  Global land cover mapping from MODIS: algorithms and early results , 2002 .

[55]  N. C. Strugnell,et al.  First operational BRDF, albedo nadir reflectance products from MODIS , 2002 .

[56]  A. Gitelson,et al.  Near real-time prediction of U.S. corn yields based on time-series MODIS data , 2014 .

[57]  G. Timár,et al.  Crop yield estimation by satellite remote sensing , 2004 .

[58]  Michael Dixon,et al.  Google Earth Engine: Planetary-scale geospatial analysis for everyone , 2017 .

[59]  Donghui Xie,et al.  Improving Spatial-Temporal Data Fusion by Choosing Optimal Input Image Pairs , 2018, Remote. Sens..

[60]  Per Jönsson,et al.  TIMESAT - a program for analyzing time-series of satellite sensor data , 2004, Comput. Geosci..

[61]  Felix Rembold,et al.  Analysis of GAC NDVI Data for Cropland Identification and Yield Forecasting in Mediterranean African Countries , 2001 .

[62]  Ray Jackson,et al.  Development of Agrometeorological Crop Model Inputs from Remotely Sensed Information , 1986, IEEE Transactions on Geoscience and Remote Sensing.

[63]  Prasad S. Thenkabail,et al.  Landsat-5 Thematic Mapper models of soybean and corn crop characteristics , 1994 .

[64]  S. Running,et al.  Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data , 2002 .

[65]  Peijuan Wang,et al.  Fusing Landsat and MODIS Data for Vegetation Monitoring , 2015, IEEE Geoscience and Remote Sensing Magazine.

[66]  Martha C. Anderson,et al.  Field-scale mapping of evaporative stress indicators of crop yield: An application over Mead, NE, USA , 2018, Remote Sensing of Environment.

[67]  D. Lobell,et al.  A scalable satellite-based crop yield mapper , 2015 .

[68]  Zhengwei Yang,et al.  Monitoring US agriculture: the US Department of Agriculture, National Agricultural Statistics Service, Cropland Data Layer Program , 2011 .

[69]  Huajun Tang,et al.  Regional yield estimation for winter wheat with MODIS-NDVI data in Shandong, China , 2008, Int. J. Appl. Earth Obs. Geoinformation.