Synergistic integration of optical and microwave satellite data for crop yield estimation

Developing accurate models of crop stress, phenology and productivity is of paramount importance, given the increasing need of food. Earth observation (EO) remote sensing data provides a unique source of information to monitor crops in a temporally resolved and spatially explicit way. In this study, we propose the combination of multisensor (optical and microwave) remote sensing data for crop yield estimation and forecasting using two novel approaches. We first propose the lag between Enhanced Vegetation Index (EVI) derived from MODIS and Vegetation Optical Depth (VOD) derived from SMAP as a new joint metric combining the information from the two satellite sensors in a unique feature or descriptor. Our second approach avoids summarizing statistics and uses machine learning to combine full time series of EVI and VOD. This study considers two statistical methods, a regularized linear regression and its nonlinear extension called kernel ridge regression to directly estimate the county-level surveyed total production, as well as individual yields of the major crops grown in the region: corn, soybean and wheat. The study area includes the US Corn Belt, and we use agricultural survey data from the National Agricultural Statistics Service (USDA-NASS) for year 2015 for quantitative assessment. Results show that (1) the proposed EVI-VOD lag metric correlates well with crop yield and outperforms common single-sensor metrics for crop yield estimation; (2) the statistical (machine learning) models working directly with the time series largely improve results compared to previously reported estimations; (3) the combined exploitation of information from the optical and microwave data leads to improved predictions over the use of single sensor approaches with coefficient of determination R≥20.76; (4) when models are used for within-season forecasting with limited time information, crop yield prediction is feasible up to four months before harvest (models reach a plateau in accuracy); and (5) the robustness of the approach is confirmed in a multi-year setting, reaching similar performances than when using single-year data. In conclusion, results confirm the value of using both EVI and VOD at the same time, and the advantage of using automatic machine learning models for crop yield/production estimation.

[1]  Arnaud Mialon,et al.  Comparison between SMOS Vegetation Optical Depth products and MODIS vegetation indices over crop zones of the USA , 2014 .

[2]  Yi Y. Liu,et al.  Global long‐term passive microwave satellite‐based retrievals of vegetation optical depth , 2011 .

[3]  G. Donald,et al.  Estimating spatio-temporal patterns of agricultural productivity in fragmented landscapes using AVHRR NDVI time series , 2003 .

[4]  Gary Schnitkey Concentration of Corn and Soybean Production in the U.S. , 2013 .

[5]  Juan Frausto Solís,et al.  Predictive ability of machine learning methods for massive crop yield prediction , 2014 .

[6]  S. Carpenter,et al.  Solutions for a cultivated planet , 2011, Nature.

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

[8]  Jarad Niemi,et al.  Monitoring Crop Growth in the Us Corn Belt with SMOS Level 2 Tau , 2018, 2018 IEEE 15th Specialist Meeting on Microwave Radiometry and Remote Sensing of the Environment (MicroRad).

[9]  C. Hamel,et al.  Crop Yield , 1999, Springer Berlin Heidelberg.

[10]  Vern C. Vanderbilt,et al.  Spectral-agronomic relationships of corn, soybean and wheat canopies , 1981 .

[11]  Gustau Camps-Valls,et al.  Machine Learning Regression Approaches for Colored Dissolved Organic Matter (CDOM) Retrieval with S2-MSI and S3-OLCI Simulated Data , 2018, Remote. Sens..

[12]  Nate G. McDowell,et al.  Interacting Effects of Leaf Water Potential and Biomass on Vegetation Optical Depth , 2017 .

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

[14]  Ranga B. Myneni,et al.  Monitoring crop yield in USA using a satellite-based climate-variability Impact Index , 2010, 2010 IEEE International Geoscience and Remote Sensing Symposium.

[15]  José Luis Rojo-Álvarez,et al.  Digital Signal Processing with Kernel Methods , 2018 .

[16]  Matthew O. Jones,et al.  Satellite passive microwave remote sensing for monitoring global land surface phenology , 2011 .

[17]  Jesslyn F. Brown,et al.  Measuring phenological variability from satellite imagery , 1994 .

[18]  Terry L. Kastens,et al.  Image masking for crop yield forecasting using AVHRR NDVI time series imagery , 2005 .

[19]  J. Clevers,et al.  Combined use of optical and microwave remote sensing data for crop growth monitoring , 1996 .

[20]  A. J. Stern,et al.  Crop Yield Assessment from Remote Sensing , 2003 .

[21]  Gustau Camps-Valls,et al.  Nonlinear Distribution Regression for Remote Sensing Applications , 2019, IEEE Transactions on Geoscience and Remote Sensing.

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

[23]  International Journal of Applied Earth Observation and Geoinformation , 2017 .

[24]  S. Liang,et al.  Estimating crop yield from multi-temporal satellite data using multivariate regression and neural network techniques , 2007 .

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

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

[27]  Martha C. Anderson,et al.  The shared and unique values of optical, fluorescence, thermal and microwave satellite data for estimating large-scale crop yields. , 2016 .

[28]  Le Yu,et al.  Rapid corn and soybean mapping in US Corn Belt and neighboring areas , 2016, Scientific Reports.

[29]  Fengmei Yao,et al.  Improved maize cultivated area estimation over a large scale combining MODIS–EVI time series data and crop phenological information , 2014 .

[30]  Claire Marais-Sicre,et al.  Estimation of corn yield using multi-temporal optical and radar satellite data and artificial neural networks , 2017, Int. J. Appl. Earth Obs. Geoinformation.

[31]  A. Strahler,et al.  Monitoring vegetation phenology using MODIS , 2003 .

[32]  Stefano Ermon,et al.  Deep Gaussian Process for Crop Yield Prediction Based on Remote Sensing Data , 2017, AAAI.

[33]  Ramesh P. Singh,et al.  Crop yield estimation model for Iowa using remote sensing and surface parameters , 2006 .

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

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

[36]  Michael Rabadi,et al.  Kernel Methods for Machine Learning , 2015 .

[37]  P. Ciais,et al.  Coupling of ecosystem-scale plant water storage and leaf phenology observed by satellite , 2018, Nature Ecology & Evolution.

[38]  A. Huete,et al.  Overview of the radiometric and biophysical performance of the MODIS vegetation indices , 2002 .

[39]  C. Field,et al.  Global scale climate–crop yield relationships and the impacts of recent warming , 2007, Environmental Research Letters.

[40]  Woubet G. Alemu,et al.  Comparing Passive Microwave with Visible-To-Near-Infrared Phenometrics in Croplands of Northern Eurasia , 2017, Remote. Sens..

[41]  Andrew E. Suyker,et al.  A Production Efficiency Model-Based Method for Satellite Estimates of Corn and Soybean Yields in the Midwestern US , 2013, Remote. Sens..

[42]  Ian McCallum,et al.  A comparison of global agricultural monitoring systems and current gaps , 2019, Agricultural Systems.

[43]  Zoltán Barcza,et al.  Statistical modelling of crop yield in Central Europe using climate data and remote sensing vegetation indices , 2018, Agricultural and Forest Meteorology.

[44]  Luciano Vieira Dutra,et al.  Mapping croplands, cropping patterns, and crop types using MODIS time-series data , 2018, Int. J. Appl. Earth Obs. Geoinformation.

[45]  A. Al Bitar,et al.  Modelling the Passive Microwave Signature from Land Surfaces: A Review of Recent Results and Application to the L-Band SMOS SMAP Soil Moisture Retrieval Algorithms , 2017 .

[46]  H. Akaike Fitting autoregressive models for prediction , 1969 .

[47]  Dara Entekhabi,et al.  L-band vegetation optical depth and effective scattering albedo estimation from SMAP. , 2017 .

[48]  T. Schmugge,et al.  Vegetation effects on the microwave emission of soils , 1991 .

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

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

[51]  Feng Zhang,et al.  Using time series of SPOT VGT NDVI for crop yield forecasting , 2003, IGARSS 2003. 2003 IEEE International Geoscience and Remote Sensing Symposium. Proceedings (IEEE Cat. No.03CH37477).

[52]  Dara Entekhabi,et al.  Vegetation optical depth and scattering albedo retrieval using time series of dual-polarized L-band radiometer observations , 2016 .

[53]  Sanja Brdar,et al.  Data Mining Approach for Predictive Modeling of Agricultural Yield Data , 2009 .

[54]  Brian K. Hornbuckle,et al.  Initial Validation of SMOS Vegetation Optical Thickness in Iowa , 2012, IEEE Geoscience and Remote Sensing Letters.

[55]  K. Cassman,et al.  Estimating crop yield potential at regional to national scales , 2013 .

[56]  Richard K. Moore,et al.  Microwave Remote Sensing, Active and Passive , 1982 .

[57]  Lorenzo Bruzzone,et al.  Kernel methods for remote sensing data analysis , 2009 .

[58]  S. Idso,et al.  Remote-Sensing of Crop Yields , 1977, Science.

[59]  Martin Brandt,et al.  Remote sensing of vegetation dynamics in drylands: Evaluating vegetation optical depth (VOD) using AVHRR NDVI and in situ green biomass data over West African Sahel , 2016 .

[60]  Markus Reichstein,et al.  Predicting carbon dioxide and energy fluxes across global FLUXNET sites with regression algorithms , 2016 .

[61]  Thomas Jagdhuber,et al.  Remote sensing of vegetation dynamics in agro-ecosystems using smap vegetation optical depth and optical vegetation indices , 2017, 2017 IEEE International Geoscience and Remote Sensing Symposium (IGARSS).

[62]  Nguyen Hieu Trung,et al.  A comparative analysis of multitemporal MODIS EVI and NDVI data for large-scale rice yield estimation , 2014 .

[63]  Adriano Camps,et al.  L-band vegetation optical depth seasonal metrics for crop yield assessment , 2018, Remote Sensing of Environment.

[64]  Clement Atzberger,et al.  Advances in Remote Sensing of Agriculture: Context Description, Existing Operational Monitoring Systems and Major Information Needs , 2013, Remote. Sens..

[65]  Nello Cristianini,et al.  Kernel Methods for Pattern Analysis , 2003, ICTAI.

[66]  Andrew E. Suyker,et al.  SMOS Optical Thickness Changes in Response to the Growth and Development of Crops, Crop Management, and Weather , 2016 .

[67]  Arnaud Mialon,et al.  Comparison of SMOS and AMSR-E vegetation optical depth to four MODIS-based vegetation indices , 2016 .

[68]  Luis Guanter,et al.  Nonlinear Statistical Retrieval of Atmospheric Profiles From MetOp-IASI and MTG-IRS Infrared Sounding Data , 2012, IEEE Transactions on Geoscience and Remote Sensing.

[69]  Rick Mueller,et al.  Intercomparison of Soil Moisture, Evaporative Stress, and Vegetation Indices for Estimating Corn and Soybean Yields Over the U.S. , 2017, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[70]  Dara Entekhabi,et al.  Moisture pulse-reserve in the soil-plant continuum observed across biomes , 2018, Nature Plants.