Crop yield estimation at a regional scale over a long period of time is of great significance to food security. In past decades, the integration of remote sensing observations and crop growth models has been recognized as a promising approach for crop growth monitoring and yield estimation. Optical remote sensing data are susceptible to cloud and rain, while synthetic aperture radar (SAR) can penetrate through clouds and has all-weather capabilities. This allows for more reliable and consistent crop monitoring and yield estimation in terms of radar sensor data. The aim of this study is to improve the accuracy for winter wheat yield estimation by assimilating time series soil moisture images, which are retrieved by a water cloud model using SAR and optical data as input, into the crop model. In this study, SAR images were acquired by C-band SAR sensors boarded on Sentinel-1 satellites and optical images were obtained from a Sentinel-2 multi-spectral instrument (MSI) for Hengshui city of Hebei province in China. Remote sensing data and ground data were all collected during the main growing season of winter wheat. Both the normalized difference vegetation index (NDVI), derived from Sentinel-2, and backscattering coefficients and polarimetric indicators, computed from Sentinel-1, were used in the water cloud model to derive time series soil moisture (SM) images. To improve the prediction of crop yields at the field scale, we incorporated remotely sensed soil moisture into the World Food Studies (WOFOST) model using the Ensemble Kalman Filter (EnKF) algorithm. In general, the trend of soil moisture inversion was consistent with the ground measurements, with the coefficient of determination (R²) equal to 0.45, 0.53, and 0.49, respectively, and RMSE was 9.16%, 7.43%, and 8.53%, respectively, for three observation dates. The winter wheat yield estimation results showed that the assimilation of remotely sensed soil moisture improved the correlation of observed and simulated yields (R² = 0.35; RMSE =934 kg/ha) compared to the situation without data assimilation (R² = 0.21; RMSE = 1330 kg/ha). Consequently, the results of this study demonstrated the potential and usefulness of assimilating SM retrieved from both Sentinel-1 C-band SAR and Sentinel-2 MSI optical remote sensing data into WOFOST model for winter wheat yield estimation and could also provide a reference for crop yield estimation with data assimilation for other crop types.

Massive desert locust swarms have been threatening and devouring natural vegetation and agricultural crops in East Africa and West Asia since 2019, and the event developed into a rare and globally concerning locust upsurge in early 2020. The breeding, maturation, concentration and migration of locusts rely on appropriate environmental factors, mainly precipitation, temperature, vegetation coverage and land-surface soil moisture. Remotely sensed images and long-term meteorological observations across the desert locust invasion area were analyzed to explore the complex drivers, vegetation losses and growing trends during the locust upsurge in this study. The results revealed that (1) the intense precipitation events in the Arabian Peninsula during 2018 provided suitable soil moisture and lush vegetation, thus promoting locust breeding, multiplication and gregarization; (2) the regions affected by the heavy rainfall in 2019 shifted from the Arabian Peninsula to West Asia and Northeast Africa, thus driving the vast locust swarms migrating into those regions and causing enormous vegetation loss; (3) the soil moisture and NDVI anomalies corresponded well with the locust swarm movements; and (4) there was a low chance the eastwardly migrating locust swarms would fly into the Indochina Peninsula and Southwest China.

The crop model data assimilation approach has been acknowledged as an effective tool for monitoring crop growth and estimating yield. However, the choice of assimilated variables and the mismatch in scale between remotely sensed observations and crop model-simulated state variables have various effects on the performance of yield estimation. This study aims to examine the accuracy of crop yield estimation through the joint assimilation of leaf area index (LAI) and soil moisture (SM) and to examine the scale effect between remotely sensed data and crop model simulations. To address these issues, we proposed an improved crop data-model assimilation (CDMA) framework, which integrates LAI and SM, as retrieved from remotely sensed data, into the World Food Studies (WOFOST) model using the ensemble Kalman filter (EnKF) approach for winter wheat yield estimation. The results showed that the yield estimation at a 10 m grid size outperformed that at a 500 m grid size, using the same assimilation strategy. Additionally, the winter wheat yield estimation accuracy was higher when using the bivariate data assimilation method (R² = 0.46, RMSE = 756 kg/ha) compared to the univariate method. In conclusion, our study highlights the advantages of joint assimilating LAI and SM for crop yield estimation and emphasizes the importance of finer spatial resolution in remotely sensed observations for crop yield estimation using the CDMA framework. The proposed approach would help to develop a high-accuracy crop yield monitoring system using optical and SAR retrieved parameters.

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.

In recent years, remote sensing (RS) research on crop growth status monitoring has gradually turned from static spectrum information retrieval in large-scale to meso-scale or micro-scale, timely multi-source data cooperative analysis; this change has presented higher requirements for RS data acquisition and analysis efficiency. How to implement rapid and stable massive RS data extraction and analysis becomes a serious problem. This paper reports on a Raster Dataset Clean & Reconstitution Multi-Grid (RDCRMG) architecture for remote sensing monitoring of vegetation dryness in which different types of raster datasets have been partitioned, organized and systematically applied. First, raster images have been subdivided into several independent blocks and distributed for storage in different data nodes by using the multi-grid as a consistent partition unit. Second, the “no metadata model” ideology has been referenced so that targets raster data can be speedily extracted by directly calculating the data storage path without retrieving metadata records; third, grids that cover the query range can be easily assessed. This assessment allows the query task to be easily split into several sub-tasks and executed in parallel by grouping these grids. Our RDCRMG-based change detection of the spectral reflectance information test and the data extraction efficiency comparative test shows that the RDCRMG is reliable for vegetation dryness monitoring with a slight reflectance information distortion and consistent percentage histograms. Furthermore, the RDCGMG-based data extraction in parallel circumstances has the advantages of high efficiency and excellent stability compared to that of the RDCGMG-based data extraction in serial circumstances and traditional data extraction. At last, an RDCRMG-based vegetation dryness monitoring platform (VDMP) has been constructed to apply RS data inversion in vegetation dryness monitoring. Through actual applications, the RDCRMG architecture is proven to be appropriate for timely vegetation dryness RS automatic monitoring with better performance, more reliability and higher extensibility. Our future works will focus on integrating more kinds of continuously updated RS data into the RDCRMG-based VDMP and integrating more multi-source datasets based collaborative analysis models for agricultural monitoring.

Land surface temperature (LST) has a critical impact on the energy balance of land surface processes and ecosystem stability. Meanwhile, LST is controlled by multiple factors at the surface, resulting in heterogeneity of its spatial distribution. To understand the drivers of LST spatial heterogeneity and their contributions, the effects of air temperature, normalized difference vegetation index (NDVI), soil moisture, net surface radiation, precipitation, aerosol optical depth (AOD), evapotranspiration, water vapor, digital elevation model (DEM), climate type, and land cover type on LST spatial heterogeneity was analyzed in this study with GeoDetector. The results showed that the explanatory ability of air temperature to impact the spatial heterogeneity of LST was the largest in each year with a mean value of 0.74, followed by water vapor with a mean value of 0.7, and the driving effect of the factors on LST showed an increasing trend year by year. However, the land cover type did not have an effect on the spatial heterogeneity of LST for the univariate analysis in this study. In addition, the interaction analysis indicated that the spatial distribution of LST was jointly driven by all the driving factors. Among them, air temperature had the strongest interaction with other factors, with the strength of the effect in the range of 0.73–0.8. In terms of the highly sensitive area of LST for each driver, AOD has the largest driving area, accounting for 15.8% of the total area, followed by WV, TA, and ET at about 11%, and the remaining variables are less than 10%. During the study period, the area of the highly sensitive region of LST for each factor showed an overall decreasing trend, indicating that the influence of the driving factors on LST will be stronger and more concentrated. Generally, this study provides meaningful understanding of the spatial heterogeneity of LST since 2003 and provides a scientific reference for coping with climate change, analyzing surface environmental patterns, and protecting ecological environment.