It is well known that timely crop growth monitoring and accurate crop yield estimation at a fine scale is of vital importance for agricultural monitoring and crop management. Crop growth models have been widely used for crop growth process description and yield prediction. In particular, the accurate simulation of important state variables, such as leaf area index (LAI) and root zone soil moisture (SM), is of great importance for yield estimation. Data assimilation is a useful tool that combines a crop model and external observations (often derived from remote sensing data) to improve the simulated crop state variables and consequently model outputs like crop total biomass, water use and grain yield. In spite of its effectiveness, applying data assimilation for monitoring crop growth at the regional scale in China remains challenging, due to the lack of high spatiotemporal resolution satellite data that can match the small field sizes which are typical for agriculture in China. With the accessibility of freely available images acquired by Sentinel satellites, it becomes possible to acquire data at high spatiotemporal resolution (10–30 m, 5–6 days), which offers attractive opportunities to characterize crop growth. In this study, we assimilated remotely sensed LAI and SM into the Word Food Studies (WOFOST) model to estimate winter wheat yield using an ensemble Kalman filter (EnKF) algorithm. The LAI was calculated from Sentinel-2 using a lookup table method, and the SM was calculated from Sentinel-1 and Sentinel-2 based on a change detection approach. Through validation with field data, the inverse error was 10% and 35% for LAI and SM, respectively. The open-loop wheat yield estimation, independent assimilations of LAI and SM, and a joint assimilation of LAI + SM were tested and validated using field measurement observation in the city of Hengshui, China, during the 2016–2017 winter wheat growing season. The results indicated that the accuracy of wheat yield simulated by WOFOST was significantly improved after joint assimilation at the field scale. Compared to the open-loop estimation, the yield root mean square error (RMSE) with field observations was decreased by 69 kg/ha for the LAI assimilation, 39 kg/ha for the SM assimilation and 167 kg/ha for the joint LAI + SM assimilation. Yield coefficients of determination (R²) of 0.41, 0.65, 0.50, and 0.76 and mean relative errors (MRE) of 4.87%, 4.32%, 4.45% and 3.17% were obtained for open-loop, LAI assimilation alone, SM assimilation alone and joint LAI + SM assimilation, respectively. The results suggest that LAI was the first-choice variable for crop data assimilation over SM, and when both LAI and SM satellite data are available, the joint data assimilation has a better performance because LAI and SM have interacting effects. Hence, joint assimilation of LAI and SM from Sentinel-1 and Sentinel-2 at a 20 m resolution into the WOFOST provides a robust method to improve crop yield estimations. However, there is still bias between the key soil moisture in the root zone and the Sentinel-1 C band retrieved SM, especially when the vegetation cover is high. By active and passive microwave data fusion, it may be possible to offer a higher accuracy SM for crop yield prediction.

Leaf area index (LAI) and canopy water content (CWC) are important variables for monitoring crop growth and drought, which can be estimated from remotely sensed data. The goal of this study was to evaluate the suitability of the Sentinel-2 multispectral instrument (S2 MSI) data for winter wheat LAI and CWC estimation with three different inversion approaches in the main farming region in North China. During the winter wheat key growth stages in 2017, 22 fields, each with five independent samples, the total number of sample plot is 110, were designed for experimental measurements. In this study, the LAI and CWC were retrieved separately using empirical models through different spectral indices, neural network (NN) algorithms, and lookup table (LUT) methods based on the PROSAIL model. The accuracies of the estimated LAI and CWC were assessed through in situ measurements. The results show that the LUT inversion approach was more suitable for LAI and CWC estimation than the spectral index-based empirical model or the NN algorithm. With the LUT approach, LAI was obtained with a root mean square error (RMSE) of 0.43 m ² ·m ^-2 and a relative RMSE (RRMSE) of 11% using seven S2 MSI bands, and CWC was obtained with an RMSE of 0.41 kg·m ^-2 , and an RRMSE of 32% using five S2 MSI bands. In all the three methods, S2 MSI was sensitive to LAI variation and able to reach higher accuracies when red edge bands were used. However, CWC inversion was still a challenge using S2 MSI data.