Field crop yield prediction is crucial to grain storage, agricultural field management, and national agricultural decision-making. Currently, crop models are widely used for crop yield prediction. However, they are hampered by the uncertainty or similarity of input parameters when extrapolated to field scale. Data assimilation methods that combine crop models and remote sensing are the most effective methods for field yield estimation. In this study, the World Food Studies (WOFOST) model is used to simulate the growing process of spring maize. Common assimilation methods face some difficulties due to the scarce, constant, or similar nature of the input parameters. For example, yield spatial heterogeneity simulation, coexistence of common assimilation methods and the nutrient module, and time cost are relatively important limiting factors. To address the yield simulation problems at field scale, a simple yet effective method with fast algorithms is presented for assimilating the time-series HJ-1 A/B data into the WOFOST model in order to improve the spring maize yield simulation. First, the WOFOST model is calibrated and validated to obtain the precise mean yield. Second, the time-series leaf area index (LAI) is calculated from the HJ data using an empirical regression model. Third, some fast algorithms are developed to complete assimilation. Finally, several experiments are conducted in a large farmland (Hongxing) to evaluate the yield simulation results. In general, the results indicate that the proposed method reliably improves spring maize yield estimation in terms of spatial heterogeneity simulation ability and prediction accuracy without affecting the simulation efficiency.

Leaf area index (LAI) is a crucial crop biophysical parameter that has been widely used in a variety of fields. Five state-of-the-art machine learning regression algorithms (MLRAs), namely, artificial neural network (ANN), support vector regression (SVR), Gaussian process regression (GPR), random forest (RF) and gradient boosting regression tree (GBRT), have been used in the retrieval of cotton LAI with Sentinel-2 spectral bands. The performances of the five machine learning models are compared for better applications of MLRAs in remote sensing, since challenging problems remain in the selection of MLRAs for crop LAI retrieval, as well as the decision as to the optimal number for the training sample size and spectral bands to different MLRAs. A comprehensive evaluation was employed with respect to model accuracy, computational efficiency, sensitivity to training sample size and sensitivity to spectral bands. We conducted the comparison of five MLRAs in an agricultural area of Northwest China over three cotton seasons with the corresponding field campaigns for modeling and validation. Results show that the GBRT model outperforms the other models with respect to model accuracy in average ( $\overline{R^2}$ = 0.854, $\overline{RMSE}$ = 0.674 and $\overline{MAE}$ = 0.456). SVR achieves the best performance in computational efficiency, which means it is fast to train, and to validate that it has great potentials to deliver near-real-time operational products for crop management. As for sensitivity to training sample size, GBRT behaves as the most robust model, and provides the best model accuracy on the average among the variations of training sample size, compared with other models ( $\overline{R^2}$ = 0.884, $\overline{RMSE}$ = 0.615 and $\overline{MAE}$ = 0.452). Spectral bands sensitivity analysis with dCor (distance correlation), combined with the backward elimination approach, indicates that SVR, GPR and RF provide relatively robust performance to the spectral bands, while ANN outperforms the other models in terms of model accuracy on the average among the reduction of spectral bands ( $\overline{R^2}$ = 0.881, $\overline{RMSE}$ = 0.625 and $\overline{MAE}$ = 0.480). A comprehensive evaluation indicates that GBRT is an appealing alternative for cotton LAI retrieval, except for its computational efficiency. Despite the different performance of the ML models, all models exhibited considerable potential for cotton LAI retrieval, which could offer accurate crop parameters information timely and accurately for crop fields management and agricultural production decisions.

Rapidly and accurately estimating yields at field scales are very significant. Each type of currently yield estimation model has been well studied, yet all of them have certain limitations. Based on a coupled Carnegie-Ames-Stanford approach (CASA)-World Food Studies (WOFOST) model and time series Sentinel-2 imagery, we achieved daily crop simulations and crop yield estimations at two adjacent farms in China. The results indicated that the coupled model inherited the high computing speed of light use efficiency (LUE) models and the mechanistic advantages of crop growth models. The ${R}^{2}$ of yield simulation was 0.64 for the CASA model, 0.84 for the coupled model, and 0.86 for the WOFOST model, and the root mean square error (RMSE) values were 948.32, 792.11, and 623.64 kg/ha, respectively. The operating times of the CASA model, CASA-WOFOST coupled model, and WOFOST model over the growing season of wheat were 37 min, 48 min, and 1 day 5 h 7 min, respectively. Compared with the WOFOST model, the coupled model provided a much faster running speed in yield simulations and a similar accuracy; therefore, the proposed model can be applied for assessments at large farms with high-spatial-resolution images to obtain accurate yield simulations. Compared with the CASA model, the coupled model provided higher simulation accuracy in mountainous areas and regions of uneven terrain. It can be concluded that the proposed coupled CASA-WOFOST model can improve the precision, reliability, and stability of crop yield estimation; provide theoretical support for crop yield estimation at field scales; and promote the development of precision agriculture.

Predicting crop maturity dates is important for improving crop harvest planning and grain quality. The prediction of crop maturity dates by assimilating remote sensing information into crop growth model has not been fully explored. In this study, a data assimilation framework incorporating the leaf area index (LAI) product from Moderate Resolution Imaging Spectroradiometer (MODIS) into a World Food Studies (WOFOST) model was proposed to predict the maturity dates of winter wheat in Henan province, China. Minimization of normalized cost function was used to obtain the input parameters of the WOFOST model. The WOFOST model was run with the re-initialized parameter to forecast the maturity dates of winter wheat grid by grid, and THORPEX Interactive Grand Global Ensemble (TIGGE) was used as forecasting period weather input in the future 15 days (d) for the WOFOST model. The results demonstrated a promising regional maturity date prediction with determination coefficient (R²) of 0.94 and the root mean square error (RMSE) of 1.86 d. The outcomes also showed that the optimal forecasting starting time for Henan was 30 April, corresponding to a stage from anthesis to grain filling. Our study indicated great potential of using data assimilation approaches in winter wheat maturity date prediction.

The approach of using multispectral remote sensing (RS) to estimate soil available nutrients (SANs) has been recently developed and shows promising results. This method overcomes the limitations of commonly used methods by building a statistical model that connects RS-based crop growth and nutrient content. However, the stability and accuracy of this model require improvement. In this article, we replaced the statistical model by integrating the World Food Studies (WOFOST) model and time series of remote sensing (T-RS) observations to ensure stability and accuracy. Time series of HJ-1 A/B data was assimilated into the WOFOST model to extrapolate crop growth simulations from a single point to a large area using a specific assimilation method. Because nutrient-limited growth within the growing season is required and the SAN parameters can only be used at the end of the growing season in the original model, the WOFOST model was modified. Notably, the calculation order was changed, and new soil nutrient uptake algorithms were implemented in the model for nutrient-limited growth estimation. Finally, experiments were conducted in the spring maize plots of Hongxing Farm to analyze the effects of nutrient stress on crop growth and the SAN simulation accuracy. The results confirm the differences in crop growth status caused by a lack of soil nutrients. The new approach can take advantage of these differences to provide better SAN estimates. In general, the new approach can overcome the limitations of existing methods and simulate the SAN status with reliable accuracy.