Improving winter wheat yield estimation by assimilation of the leaf area index from Landsat TM and MODIS data into the WOFOST model

Abstract To predict regional-scale winter wheat yield, we developed a crop model and data assimilation framework that assimilated leaf area index (LAI) derived from Landsat TM and MODIS data into the WOFOST crop growth model. We measured LAI during seven phenological phases in two agricultural cities in China’s Hebei Province. To reduce cloud contamination, we applied Savitzky–Golay (S–G) filtering to the MODIS LAI products to obtain a filtered LAI. We then regressed field-measured LAI on Landsat TM vegetation indices to derive multi-temporal TM LAIs. We developed a nonlinear method to adjust LAI by accounting for the scale mismatch between the remotely sensed data and the model’s state variables. The TM LAI and scale-adjusted LAI datasets were assimilated into the WOFOST model to allow evaluation of the yield estimation accuracy. We constructed a four-dimensional variational data assimilation (4DVar) cost function to account for the observations and model errors during key phenological stages. We used the shuffled complex evolution–University of Arizona algorithm to minimize the 4DVar cost function between the remotely sensed and modeled LAI and to optimize two important WOFOST parameters. Finally, we simulated winter wheat yield in a 1-km grid for cells with at least 50% of their area occupied by winter wheat using the optimized WOFOST, and aggregated the results at a regional scale. The scale adjustment substantially improved the accuracy of regional wheat yield predictions (R2 = 0.48; RMSE = 151.92 kg ha−1) compared with the unassimilated results (R2 = 0.23; RMSE = 373.6 kg ha−1) and the TM LAI results (R2 = 0.27; RMSE = 191.6 kg ha−1). Thus, the assimilation performance depends strongly on the LAI retrieval accuracy and the scaling correction. Our research provides a scheme to employ remotely sensed data, ground-measured data, and a crop growth model to improve regional crop yield estimates.

[1]  Elizabeth Pattey,et al.  Using Leaf Area Index, retrieved from optical imagery, in the STICS crop model for predicting yield and biomass of field crops , 2012 .

[2]  Hongliang Fang,et al.  Corn‐yield estimation through assimilation of remotely sensed data into the CSM‐CERES‐Maize model , 2008 .

[3]  V. Singh,et al.  The EPIC model. , 1995 .

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

[5]  Harry Vereecken,et al.  Multivariate and Multiscale Data Assimilation in Terrestrial Systems: A Review , 2012, Sensors.

[6]  W. Verhoef,et al.  A spectral directional reflectance model of row crops , 2010 .

[7]  R. Colombo,et al.  Retrieval of leaf area index in different vegetation types using high resolution satellite data , 2003 .

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

[9]  P. Sellers.  Canopy reflectance, photosynthesis, and transpiration. II. the role of biophysics in the linearity of their interdependence , 1987 .

[10]  J. Wolf,et al.  WOFOST: a simulation model of crop production. , 1989 .

[11]  C. A. van Diepen,et al.  User's guide for the WOFOST 7.1 crop growth simulation model and WOFOST Control Center 1.5. , 1998 .

[12]  Kiyoshi Honda,et al.  Soil hydraulic parameters estimated from satellite information through data assimilation , 2011 .

[13]  F. Baret,et al.  Remotely sensed green area index for winter wheat crop monitoring:10-Year assessment at regional scale over a fragmented landscape , 2012 .

[14]  A. Huete.  A soil-adjusted vegetation index (SAVI) , 1988 .

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

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

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

[18]  W. Bastiaanssen,et al.  A new crop yield forecasting model based on satellite measurements applied across the Indus Basin, Pakistan , 2003 .

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

[20]  G. Hoogenboom,et al.  Integration of MODIS LAI and vegetation index products with the CSM–CERES–Maize model for corn yield estimation , 2011 .

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

[22]  D. Raes,et al.  AquaCrop-The FAO Crop Model to Simulate Yield Response to Water: I. Concepts and Underlying Principles , 2009 .

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

[24]  James W. Jones,et al.  The DSSAT cropping system model , 2003 .

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

[26]  Jin Chen,et al.  A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter , 2004 .

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

[28]  Soroosh Sorooshian,et al.  Optimal use of the SCE-UA global optimization method for calibrating watershed models , 1994 .

[29]  X. Moa,et al.  Prediction of crop yield , water consumption and water use efficiency with a SVAT-crop growth model using remotely sensed data on the North China Plain , 2005 .

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

[31]  Frédéric Baret,et al.  Using Thermal Time and Pixel Purity for Enhancing Biophysical Variable Time Series: An Interproduct Comparison , 2013, IEEE Transactions on Geoscience and Remote Sensing.

[32]  A. Bondeau,et al.  Combining agricultural crop models and satellite observations: from field to regional scales , 1998 .

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

[34]  F. Baret,et al.  Retrieving wheat Green Area Index during the growing season from optical time series measurements based on neural network radiative transfer inversion , 2011 .

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

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

[37]  M. Claverie,et al.  Maize and sunflower biomass estimation in southwest France using high spatial and temporal resolution remote sensing data , 2012 .

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

[39]  S. Liang,et al.  Validation of MODIS and CYCLOPES LAI products using global field measurement data , 2012 .

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

[41]  S. Running,et al.  Relationship of thematic mapper simulator data to leaf area index , 1987 .