Assimilating a synthetic Kalman filter leaf area index series into the WOFOST model to improve regional winter wheat yield estimation

Abstract The scale mismatch between remote sensing observations and state variables simulated by crop growth models decreases the reliability of crop yield estimates. To overcome this problem, we implemented a two-step data-assimilation approach: first, we generated a time series of 30-m-resolution leaf area index (LAI) by combining Moderate Resolution Imaging Spectroradiometer (MODIS) data and three Landsat TM images with a Kalman filter algorithm (the synthetic KF LAI series); second, the time series were assimilated into the WOFOST crop growth model to generate an ensemble Kalman filter LAI time series (the EnKF-assimilated LAI series). The synthetic EnKF LAI series then drove the WOFOST model to simulate winter wheat yields at 1-km resolution for pixels with wheat fractions of at least 50%. The county-level aggregated yield estimates were compared with official statistical yields. The synthetic KF LAI time series produced a more realistic characterization of LAI phenological dynamics. Assimilation of the synthetic KF LAI series produced more accurate estimates of regional winter wheat yield ( R 2  = 0.43; root-mean-square error (RMSE) = 439 kg ha −1 ) than three other approaches: WOFOST without assimilation (determination coefficient R 2  = 0.14; RMSE = 647 kg ha −1 ), assimilation of Landsat TM LAI ( R 2  = 0.37; RMSE = 472 kg ha −1 ), and assimilation of S-G filtered MODIS LAI ( R 2  = 0.49; RMSE = 1355 kg ha −1 ). Thus, assimilating the synthetic KF LAI series into the WOFOST model with the EnKF strategy provides a reliable and promising method for improving regional estimates of winter wheat yield.

[1]  R. E. Kalman,et al.  A New Approach to Linear Filtering and Prediction Problems , 2002 .

[2]  C. A. van Diepen,et al.  Crop model data assimilation with the Ensemble Kalman filter for improving regional crop yield forecasts , 2007 .

[3]  Jun Qin,et al.  Data Assimilation Methods for Land Surface Variable Estimation , 2008 .

[4]  Pierre Defourny,et al.  Potential performances of remotely sensed LAI assimilation in WOFOST model based on an OSS experiment , 2011 .

[5]  Yanxia Zhao,et al.  Assimilating remote sensing information with crop model using Ensemble Kalman Filter for improving LAI monitoring and yield estimation , 2013 .

[6]  W. Verhoef,et al.  Simulated impact of sensor field of view and distance on field measurements of bidirectional reflectance factors for row crops , 2015 .

[7]  Wei Su,et al.  Estimating regional winter wheat yield by assimilation of time series of HJ-1 CCD NDVI into WOFOST-ACRM model with Ensemble Kalman Filter , 2013, Math. Comput. Model..

[8]  Jianxi Huang,et al.  Jointly Assimilating MODIS LAI and ET Products Into the SWAP Model for Winter Wheat Yield Estimation , 2015, IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing.

[9]  L. Dente,et al.  Assimilation of leaf area index derived from ASAR and MERIS data into CERES - wheat model to map wheat yield , 2008 .

[10]  H. Rauch.  Solutions to the linear smoothing problem , 1963 .

[11]  Caiyan Lin,et al.  An Ensemble Kalman Filter for severe dust storm data assimilation over China , 2008 .

[12]  W. G. M. Bastiaanssen,et al.  Assimilation of satellite data into agrohydrological models to improve crop yield forecasts , 2009 .

[13]  Philip Lewis,et al.  Assimilating canopy reflectance data into an ecosystem model with an Ensemble Kalman Filter , 2008 .

[14]  Jianxi Huang,et al.  Assimilation of MODIS-LAI into the WOFOST model for forecasting regional winter wheat yield , 2013, Math. Comput. Model..

[15]  D. Roy,et al.  Multi-temporal MODIS-Landsat data fusion for relative radiometric normalization, gap filling, and prediction of Landsat data , 2008 .

[16]  Gregory Duveiller,et al.  Estimating regional winter wheat yield with WOFOST through the assimilation of green area index retrieved from MODIS observations , 2012 .

[17]  Tomás Martínez-Marín,et al.  Crop Phenology Estimation Using a Multitemporal Model and a Kalman Filtering Strategy , 2014, IEEE Geoscience and Remote Sensing Letters.

[18]  Xiaolin Zhu,et al.  An enhanced spatial and temporal adaptive reflectance fusion model for complex heterogeneous regions , 2010 .

[19]  Michael E. Schaepman,et al.  A review on reflective remote sensing and data assimilation techniques for enhanced agroecosystem modeling , 2007, Int. J. Appl. Earth Obs. Geoinformation.

[20]  Qingguo Zhou,et al.  Assimilating remote sensing information into a coupled hydrology-crop growth model to estimate regional maize yield in arid regions , 2014 .

[21]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[22]  G. Evensen,et al.  Analysis Scheme in the Ensemble Kalman Filter , 1998 .

[23]  Jianxi Huang,et al.  Regional Crop Yield Assessment by Combination of a Crop Growth Model and Phenology Information Derived from MODIS , 2011 .

[24]  Jianxi Huang,et al.  Comparison of Two Optimization Algorithms for Estimating Regional Winter Wheat Yield by Integrating MODIS Leaf Area Index and World Food Studies Model , 2013 .

[25]  Jing Wang,et al.  Estimating near future regional corn yields by integrating multi-source observations into a crop growth model , 2013 .

[26]  Pieter Kempeneers,et al.  A Kalman Filter-Based Method to Generate Continuous Time Series of Medium-Resolution NDVI Images , 2014, Remote. Sens..

[27]  F. Baret,et al.  Influence of landscape spatial heterogeneity on the non-linear estimation of leaf area index from moderate spatial resolution remote sensing data , 2006 .

[28]  Ling Lu,et al.  Impacts and uncertainties of upscaling of remote-sensing data validation for a semi-arid woodland , 2008 .

[29]  Greg Welch,et al.  Welch & Bishop , An Introduction to the Kalman Filter 2 1 The Discrete Kalman Filter In 1960 , 1994 .

[30]  R. Colombo,et al.  Inversion of a radiative transfer model with hyperspectral observations for LAI mapping in poplar plantations , 2004 .

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

[32]  Bernhard Geiger,et al.  Use of a Kalman filter for the retrieval of surface BRDF coefficients with a time-evolving model based on the ECOCLIMAP land cover classification , 2008 .

[33]  Jianxi Huang,et al.  Improving Timeliness of Winter Wheat Production Forecast in United States of America, Ukraine and China Using MODIS Data and NCAR Growing Degree Day , 2014 .

[34]  Dehai Zhu,et al.  Improving winter wheat yield estimation by assimilation of the leaf area index from Landsat TM and MODIS data into the WOFOST model , 2015 .

[35]  Chad J. Shuey,et al.  Validating MODIS land surface reflectance and albedo products: methods and preliminary results , 2002 .

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

[37]  F. Baret,et al.  Crop specific green area index retrieval from MODIS data at regional scale by controlling pixel-target adequacy , 2011 .

[38]  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 .

[39]  Pierre Defourny,et al.  Estimating regional wheat yield from the shape of decreasing curves of green area index temporal profiles retrieved from MODIS data , 2012, Int. J. Appl. Earth Obs. Geoinformation.

[40]  Peter S. Maybeck,et al.  Stochastic Models, Estimation And Control , 2012 .

[41]  Chunlin Huang,et al.  Experiments of one-dimensional soil moisture assimilation system based on ensemble Kalman filter , 2008 .

[42]  T. Başar,et al.  A New Approach to Linear Filtering and Prediction Problems , 2001 .

[43]  James Hansen,et al.  Assimilation of remotely sensed soil moisture and vegetation with a crop simulation model for maize yield prediction , 2013 .

[44]  Luis Alonso,et al.  Multitemporal fusion of Landsat/TM and ENVISAT/MERIS for crop monitoring , 2013, Int. J. Appl. Earth Obs. Geoinformation.

[45]  Pierre-Philippe Mathieu,et al.  Data assimilation: From photon counts to Earth System forecasts , 2008 .

[46]  Patrick Matgen,et al.  Enhanced biomass prediction by assimilating satellite data into a crop growth model , 2014, Environ. Model. Softw..

[47]  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.

[48]  M. Guérif,et al.  Assimilating remote sensing data into a crop model to improve predictive performance for spatial applications , 2005 .

[49]  Li Zhang,et al.  Monitoring winter wheat growth in North China by combining a crop model and remote sensing data , 2008, Int. J. Appl. Earth Obs. Geoinformation.

[50]  M. Schaepman,et al.  Downscaling time series of MERIS full resolution data to monitor vegetation seasonal dynamics , 2009 .