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.

Generally, there is an inconsistency between the total regional crop area that was obtained from remote sensing technology and the official statistical data on crop areas. When performing scale conversion and data aggregation of remote sensing-based crop mapping results from different administrative scales, it is difficult to obtain accurate crop planting area that match crop area statistics well at the corresponding administrative level. This problem affects the application of remote sensing-based crop mapping results. In order to solve the above problem, taking Fucheng County of Hebei Province in the Huanghuaihai Plain of China as the study area, based on the Sentinel-2 normalized difference vegetation index (NDVI) time series data covering the whole winter wheat growth period, the statistical data of the regional winter wheat planting area were regarded as reference for the winter wheat planting area extracted by remote sensing, and a new method for winter wheat mapping that is based on similarity measurement indicators and their threshold optimizations (WWM-SMITO) was proposed with the support of the shuffled complex evolution-University of Arizona (SCE-UA) global optimization algorithm. The accuracy of the regional winter wheat mapping results was verified, and accuracy comparisons with different similarity indicators were carried out. The results showed that the total area accuracy of the winter wheat area extraction by the proposed method reached over 99.99%, which achieved a consistency that was between the regional remote sensing-based winter wheat planting area and the statistical data on the winter wheat planting area. The crop recognition accuracy also reached a high level, which showed that the proposed method was effective and feasible. Moreover, in the accuracy comparison of crop mapping results based on six different similarity indicators, the winter wheat distribution that was extracted by root mean square error (RMSE) had the best recognition accuracy, and the overall accuracy and kappa coefficient were 94.5% and 0.8894, respectively. The overall accuracies of winter wheat that were extracted by similarity indicators, such as Euclidean distance (ED), Manhattan distance (MD), spectral angle mapping (SAM), and spectral correlation coefficient (SCC) were 94.1%, 93.9%, 93.3%, and 92.8%, respectively, and the kappa coefficients were 0.8815, 0.8776, 0.8657, and 0.8558, respectively. The accuracy of the winter wheat results extracted by the similarity indicator of dynamic time warping (DTW) was relatively low. The results of this paper could provide guidance and serve as a reference for the selection of similarity indicators in crop distribution extraction and for obtaining large-scale, long-term, and high-precision remote sensing-based information on a regional crop spatial distribution that is highly consistent with statistical crop area data.

In-season crop-type maps are required for a variety of agricultural monitoring and decision-making applications. The earlier the crop type maps of the current growing season are obtained, the more beneficial it is for agricultural decision-making and management. With the availability of a large amount of high spatiotemporal resolution remote sensing data, different data sources are expected to increase the frequency of data acquisition, which can provide more information in the early season. To explore the potential of integrating different data sources, a Dual-1DCNN algorithm was built based on the CNN model in this study. Moreover, an incremental training method was used to attain the network on each data acquisition date and obtain the best detection date for each crop type in the early season. A case study for Hengshui City in China was conducted using time series of Sentinel-1A (S1A) and Sentinel-2 (S2) attained in 2019. To verify this method, the classical methods support vector machine (SVM), random forest (RF), and Mono-1DCNN were implemented. The input for SVM and RF was S1A and S2 data, and the input for Mono-1DCNN was S2 data. The results demonstrated the following: (1) Dual-1DCNN achieved an overall accuracy above 85% at the earliest time.; (2) all four types of models achieved high accuracy (F1s were greater than 90%) on summer maize after sowing one month later; (3) for cotton and common yam rhizomes, Dual-1DCNN performed best, with its F1 reaching 85% within 2 months after cotton sowing, 15 days, 20 days, and 45 days ahead of Mono-1DCNN, SVM, and RF, respectively, and its extraction of the common yam rhizome was achieved 1–2 months earlier than other methods within the acceptable accuracy. These results confirmed that Dual-1DCNN offered significant potential in the in-season detection of crop types.

To improve crop model performance for regional crop yield estimates, a new four-dimensional variational algorithm (POD4DVar) merging the Monte Carlo and proper orthogonal decomposition techniques was introduced to develop a data assimilation strategy using the Crop Environment Resource Synthesis (CERES)-Wheat model. Two winter wheat yield estimation procedures were conducted on a field plot and regional scale to test the feasibility and potential of the POD4DVar-based strategy. Winter wheat yield forecasts for the field plots showed a coefficient of determination (R²) of 0.73, a root mean square error (RMSE) of 319 kg/ha, and a relative error (RE) of 3.49%. An acceptable yield at the regional scale was estimated with an R² of 0.997, RMSE of 7346 tons, and RE of 3.81%. The POD4DVar-based strategy was more accurate and efficient than the EnKF-based strategy. In addition to crop yield, other critical crop variables such as the biomass, harvest index, evapotranspiration, and soil organic carbon may also be estimated. The present study thus introduces a promising approach for operationally monitoring regional crop growth and predicting yield. Successful application of this assimilation model at regional scales must focus on uncertainties derived from the crop model, model inputs, data assimilation algorithm, and assimilated observations.

The crop harvest index (HI) is of great significance for research on the application of crop variety breeding, crop growth simulation, crop management in precision agriculture and crop yield estimation, among other topics. To obtain spatial information on the crop dynamic HI (D-HI), taking winter wheat as the research object and fully considering the changes in crop biomass and the grain filling process from the flowering period to the maturity period, the dynamic f_G (D-f_G) parameter was estimated as the ratio between the aboveground biomass accumulated in different growth periods, from the flowering stage to the maturity stage, and the aboveground biomass in the corresponding periods. Based on the D-f_G parameter estimation using unmanned aerial vehicle (UAV) hyperspectral remote sensing data, a technical method for obtaining spatial information on the winter wheat D-HI was proposed and the accuracy of the proposed method was verified. A correlation analysis was performed between the normalized difference spectral index (NDSI), which was calculated using pairs of any two bands of the UAV hyperspectral spectrum, and the measured D-f_G. Based on this correlation analysis, the center of gravity of the local maximum region of R² was used to determine the sensitive band center to accurately estimate D-f_G. On this basis, remote sensing estimation of the D-f_G was realized by using the NDSI constructed by the sensitive hyperspectral band centers. Finally, based on the D-f_G remote sensing parameters and the D-HI estimation model, spatial information on the D-HI of winter wheat was accurately obtained. The results revealed five pairs of sensitive hyperspectral band centers (i.e., λ(476 nm, 508 nm), λ(444 nm, 644 nm), λ(608 nm, 788 nm), λ(724 nm, 784 nm) and λ(816 nm, 908 nm)) for D-f_G estimation, and the results of the D-f_G remote sensing estimation showed high precision. The root mean square error (RMSE) was between 0.0436 and 0.0604, the normalized RMSE (NRMSE) was between 10.31% and 14.27% and the mean relative error (MRE) was between 8.28% and 12.55%. In addition, the D-f_G parameter estimation, using the NDSI constructed by the above five sensitive remote sensing band centers, yielded highly accurate spatial D-HI information with an RMSE between 0.0429 and 0.0546, an NRMSE between 9.87% and 12.57% and an MRE between 8.33% and 10.90%. The D-HI estimation results based on the hyperspectral sensitive band centers λ(724 nm, 784 nm) had the highest accuracy, with RMSE, NRMSE and MRE values of 0.0429, 9.87% and 8.33%, respectively. The proposed method of acquiring spatial information on the winter wheat D-HI in this study was shown to be feasible, and it might provide a technical reference toward developing satellite-based indices to monitor large-scale crop HI information.

Timely and accurate estimation of the winter wheat planting area and its spatial distribution is essential for the implementation of crop growth monitoring and yield estimation, and hence for the development of national agricultural production and food security. In remotely sensed winter wheat mapping based on spectral similarity, the reference curve is obtained by averaging multiple standard curves, which limits mapping accuracy. We propose a spectral reconstruction method based on singular value decomposition (SR-SVD) for winter wheat mapping based on the unique growth characteristics of crops. Using Sentinel-2 A/B satellite data, we tested the SR-SVD method in Puyang County, and Shenzhou City, China. Performance was increased, with the optimal overall accuracy and the Kappa of Puyang County and Shenzhou City were 99.52% and 0.99, and 98.26% and 0.97, respectively. We selected the spectral angle mapper (SAM) and Euclidean Distance (ED) as the similarity measures. Compared to spectral similarity methods, the SR-SVD method significantly improves mapping accuracy, as it avoids excessive extraction, can identify more detailed information, and is advantageous in distinguishing non-winter wheat pixels. Three commonly used supervised classification methods, support vector machine (SVM), maximum likelihood (ML), and minimum distance (MD) were used for comparison. Results indicate that SR-SVD has the highest mapping accuracy and greatly reduces the number of misidentified pixels. Therefore, the SR-SVD method can achieve high-precision crop mapping and provide technical support for monitoring regional crop planting structure information.