Search

Accurate and timely crop yield prediction over large spatial regions is critical to national food security and sustainable agricultural development. However, designing a robust model for crop yield prediction over a large spatial region remains challenging due to inadequate surveyed samples and an under-development of deep-learning frameworks. To tackle this issue, we integrated multi-source (remote sensing, weather, and soil properties) data into a dual-stream deep-learning neural network model for winter wheat in China's major planting regions. The model consists of two branches for robust feature learning: one for sequential data (remote sensing and weather series data) and the other for statical data (soil properties). The extracted features by both branches were aggregated through an adaptive fusion model to forecast the final wheat yield. We trained and tested the model by using official county-level statistics of historical winter wheat yields. The model achieved an average R2 of 0.79 and a root-mean-square error of 650.21 kg/ha, superior to the compared methods and outperforming traditional machine-learning methods. The dual-stream deep-learning neural network model provided decent in-season yield prediction, with an error of about 13% compared to official statistics about two months before harvest. By effectively extracting and aggregating features from multi-source datasets, the new approach provides a practical approach to predicting winter wheat yields at the county scale over large spatial regions. [ABSTRACT FROM AUTHOR]

Grassland aboveground biomass is crucial for evaluating grassland desertification, degradation, and grassland and livestock balance. Given the lack of understanding of mechanical processes and limited simulation accuracy for grassland aboveground biomass estimation, especially at the regional scale, this study investigates a new method combining remote sensing data assimilation technology and a grassland process-based model to estimate regional grassland biomass, focusing on improving the simulation accuracy by modeling and revealing the mechanism interpretability of grassland growth processes. Xilinhot City of Inner Mongolia was used as the study area. The ModVege model was selected as the grass dynamic simulation model. A likelihood function was constructed composed of the LAI, grassland aboveground biomass, and daily measurements wherein the accumulated temperature reached ST2 (the temperature sum defining the end of reproductive growth). Then, the Markov chain Monte Carlo (MCMC) methodology was adapted to calibrate the ModVege model by maximizing the likelihood function. The time-series LAI from MOD15A3H was assimilated into the ModVege model, and the model parameters ST2 and BMGV0 (initial biomass and green vegetative tissues, respectively) were optimized at a 500 m pixel scale based on the four-dimensional variational method (4DVar) method. Compared with August 15th, the RMSE and MAPE of aboveground biomass were 242 kg/ha and 10%, respectively, after calibration. Data assimilation improved this accuracy, with the RMSE decreasing to 214 kg/ha. Overall, the aboveground grassland biomass of Xilinhot City shows spatial distribution patterns of high value in the northeast and low value in the central and southeast areas. Generally, the method implemented in this study provides an important reference for the aboveground biomass estimation of regional grassland. [ABSTRACT FROM AUTHOR]

Accurate, timely, and fine-resolution crop phenology is essential for determining the optimal timing of agronomic management practices supporting precision agriculture and food security. Synthetic Aperture Radar (SAR) methods, unaffected by cloud occlusion, have been widely applied in monitoring maize phenology. Nonetheless, their reliance on manual threshold settings, which depend on the user's expertise, limits their applicability. Furthermore, the neglect of SAR's potential for monitoring other phenological periods (e.g., seven-leaves date (V7), jointing date (JD), tassel date (TD), and milky date (MID)) hinders their robustness, particularly for regional-scale applications. To address these issues, this study used an adaptive dynamic threshold to evaluate the ability of the Sentinel-1 cross-polarization ratio (CR) in detecting the three-leaves date (V3), V7, JD, TD, MID, and maturity date (MD) of maize. We analyzed the effect of incidence angle, precipitation, and wind speed on Sentinel-1 features to identify the optimal feature for time series fitting. Then, we employed linear regression to determine the optimal threshold and developed an adaptive dynamic threshold for phenology detection. This approach effectively mitigated the speckle noise of Sentinel-1 and minimized artificial interference caused by customary conventional thresholds. Finally, we mapped phenology across 8.3 million ha in Heilongjiang Province. The results indicated that the approach has a higher ability to detect JD (RMSE = 11.10 d), MID (RMSE = 10.31 d), and MD (RMSE = 9.41 d) than that of V3 (RMSE = 32.07 d), V7 (RMSE = 56.37 d), and TD (RMSE = 43.33 d) in Sentinel-1. Compared with Sentinel-2, the average RMSE of JD, MID, and MD decreased by 4.14%, 35.28%, and 26.48%. Moreover, when compared to different thresholds, the adaptive dynamic threshold can quickly determine the optimal threshold for detecting each phenological stage. CR is least affected by incident angle, precipitation, and wind speed, effectively suppressing noise to reflect phenological development better. This approach supports the rapid and feasible mapping of maize phenology across broad spatial regions with a few samples. [ABSTRACT FROM AUTHOR]

In July 2021, a flooding event, which attracted the attention of the whole country and even the world, broke out in Henan, resulting in dramatic losses across multiple fields (e.g., economic and agricultural). The basin at the junction of Hebi, Xinxiang, and Anyang was the most affected region, as the spread of water from the Wei river submerged surrounding agricultural land (e.g., corn-dominated). To comprehensively evaluate the flooding impacts, we proposed a framework to detect the flooding area and evaluated the degree of loss using satellite time series data. First, we proposed a double-Gaussian model to adaptively determine the threshold for flooding extraction using Synthetic Aperture Radar (SAR) data. Then, we evaluated the disaster levels of flooding with field survey samples and optical satellite images. Finally, given that crops vary in their resilience to flooding, we measured the vegetation index change before and after the flooding event using satellite time series data. We found the proposed double-Gaussian model could accurately extract the flooding area, showing great potential to support in-time flooding evaluation. We also showed that the multispectral satellite images could potentially support the classification of disaster levels (i.e., normal, slight, moderate, and severe), with an overall accuracy of 88%. Although these crops were temporarily affected by this flooding event, most recovered soon, especially for the slightly and moderately affected regions. Overall, the distribution of resilience of these affected crops was basically in line with the results of classified disaster levels. The proposed framework provides a comprehensive aspect to the retrospective study of the flooding process on crops with diverse disaster levels and resilience. It can provide rapid and timely flood damage assessment and support emergency management and disaster verification work. [ABSTRACT FROM AUTHOR]

Satellite-based models have tremendous potential for monitoring crop production because satellite data can provide temporally and spatially continuous crop growth information at large scale. This study used a satellite-based vegetation production model (i.e., eddy covariance light use efficiency, EC-LUE) to estimate national winter wheat gross primary production, and then combined this model with the harvest index (ratio of aboveground biomass to yield) to convert the estimated winter wheat production to yield. Specifically, considering the spatial differences of the harvest index, we used a cross-validation method to invert the harvest index of winter wheat among counties, municipalities and provinces. Using the field-surveyed and statistical yield data, we evaluated the model performance, and found the model could explain more than 50% of the spatial variations of the yield both in field-surveyed regions and most administrative units. Overall, the mean absolute percentage errors of the yield are less than 20% in most counties, municipalities and provinces, and the mean absolute percentage errors for the production of winter wheat at the national scale is 4.06%. This study demonstrates that a satellite-based model is an alternative method for crop yield estimation on a larger scale. [ABSTRACT FROM AUTHOR]

Timely and accurate regional crop-yield estimates are crucial for guiding agronomic practices and policies to improve food security. In this study, a crop-growth model was integrated with time series of remotely sensed data through deep learning (DL) methods to improve the accuracy of regional wheat-yield estimations in Henan Province, China. Firstly, the time series of moderate-resolution imaging spectroradiometer (MODIS) normalized difference vegetation index (NDVI) were input into the long short-term memory network (LSTM) model to identify the wheat-growing region, which was further used to estimate wheat areas at the municipal and county levels. Then, the leaf area index (LAI) and grain-yield time series simulated by the Crop Environment REsource Synthesis for Wheat (CERES-Wheat) model were used to train and evaluate the LSTM, one-dimensional convolutional neural network (1-D CNN) and random forest (RF) models, respectively. Finally, an exponential model of the relationship between the field-measured LAI and MODIS NDVI was applied to obtain the regional LAI, which was input into the trained LSTM, 1-D CNN and RF models to estimate wheat yields within the wheat-growing region. The results showed that the linear correlations between the estimated wheat areas and the statistical areas were significant at both the municipal and county levels. The LSTM model provided more accurate estimates of wheat yields, with higher R2 values and lower root mean square error (RMSE) and mean relative error (MRE) values than the 1-D CNN and RF models. The LSTM model has an inherent advantage in capturing phenological information contained in the time series of the MODIS-derived LAI, which is important for satellite-based crop-yield estimates. [ABSTRACT FROM AUTHOR]

Crop mapping using satellite imagery is crucial for agriculture applications. However, a fundamental challenge that hinders crop mapping progress is the scarcity of samples. The latest foundation model, Segment Anything Model (SAM), provides an opportunity to address this issue, yet few studies have been conducted in this area. This study investigated the parcel segmentation performance of SAM on commonly used medium-resolution satellite imagery (i.e., Sentinel-2 and Landsat-8) and proposed a novel automated sample generation framework based on SAM. The framework comprises three steps. First, an image optimization automatically selects high-quality images as the inputs for SAM. Then, potential samples are generated based on the masks produced by SAM. Finally, the potential samples are subsequently subjected to a sample cleaning procedure to acquire the most reliable samples. Experiments were conducted in Henan Province, China, and southern Ontario, Canada, using six proven effective classifiers. The effectiveness of our method is demonstrated through the combination of field-survey-collected samples and differently proportioned generated samples. Our results indicated that directly using SAM for parcel segmentation remains challenging, unless the parcels are large, regular in shape, and have distinct color differences from surroundings. Additionally, the proposed approach significantly improved the performance of classifiers and alleviated the sample scarcity problem. Compared to classifiers trained only by field-survey-collected samples, our method resulted in an average improvement of 16% and 78.5% in Henan and Ontario, respectively. The random forest achieved relatively good performance, with weighted-average F1 of 0.97 and 0.996 obtained using Sentinel-2 imagery in the two study areas, respectively. Our study contributes insights into solutions for sample scarcity in crop mapping and highlights the promising application of foundation models like SAM. [ABSTRACT FROM AUTHOR]

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 (R2) 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. [ABSTRACT FROM AUTHOR]

Timely and accurate forecasting of crop yields is crucial to food security and sustainable development in the agricultural sector. However, winter wheat yield estimation and forecasting on a regional scale still remains challenging. In this study, we established a two-branch deep learning model to predict winter wheat yield in the main producing regions of China at the county level. The first branch of the model was constructed based on the Long Short-Term Memory (LSTM) networks with inputs from meteorological and remote sensing data. Another branch was constructed using Convolution Neural Networks (CNN) to model static soil features. The model was then trained using the detrended statistical yield data during 1982 to 2015 and evaluated by leave-one-year-out-validation. The evaluation results showed a promising performance of the model with the overall R 2 and RMSE of 0.77 and 721 kg/ha, respectively. We further conducted yield prediction and uncertainty analysis based on the two-branch model and obtained the forecast accuracy in one month prior to harvest of 0.75 and 732 kg/ha. Results also showed that while yield detrending could potentially introduce higher uncertainty, it had the advantage of improving the model performance in yield prediction. [ABSTRACT FROM AUTHOR]

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 (R2) 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 (R2 = 0.35; RMSE =934 kg/ha) compared to the situation without data assimilation (R2 = 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. [ABSTRACT FROM AUTHOR]

Leaf angle is a critical structural parameter for retrieving canopy leaf area index (LAI) using the PROSAIL model. However, the traditional method using default leaf angle distribution in the PROSAIL model does not capture the phenological dynamics of canopy growth. This study presents a LAI retrieval method for corn canopies using PROSAIL model with leaf angle distribution functions referred from terrestrial laser scanning points at four phenological stages during the growing season. Specifically, four inferred maximum-probability leaf angles were used in the Campbell ellipsoid leaf angle distribution function of PROSAIL. A Lookup table (LUT) is generated by running the PROSAIL model with inferred leaf angles, and the cost function is minimized to retrieve LAI. The results show that the leaf angle distribution functions are different for the corn plants at different phenological growing stages, and the incorporation of derived specific corn leaf angle distribution functions distribute the improvement of LAI retrieval using the PROSAIL model. This validation is done using in-situ LAI measurements and MODIS LAI in Baoding City, Hebei Province, China, and compared with the LAI retrieved using default leaf angle distribution function at the same time. The root-mean-square error (RMSE) between the retrieved LAI on 4 September 2014, using the modified PROSAIL model and the in-situ measured LAI was 0.31 m2/m2, with a strong and significant correlation (R2 = 0.82, residual range = 0 to 0.6 m2/m2, p < 0.001). Comparatively, the accuracy of LAI retrieved results using default leaf angle distribution is lower, the RMSE of which is 0.56 with R2 = 0.76 and residual range = 0 to 1.0 m2/m2, p < 0.001. This validation reveals that the introduction of inferred leaf angle distributions from TLS data points can improve the LAI retrieval accuracy using the PROSAIL model. Moreover, the comparations of LAI retrieval results on 10 July, 26 July, 19 August and 4 September with default and inferred corn leaf angle distribution functions are all compared with MODIS LAI products in the whole study area. This validation reveals that improvement exists in a wide spatial range and temporal range. All the comparisons demonstrate the potential of the modified PROSAIL model for retrieving corn canopy LAI from Landsat imagery by inferring leaf orientation from terrestrial laser scanning data. [ABSTRACT FROM AUTHOR]

Crop residue cover is vital for reducing soil erosion and improving soil fertility, which is an important way of conserving tillage to protect the black soil in Northeast China. How much the crop residue covers on cropland is of significance for black soil protection. Landsat-8 and Sentinel-2 images were used to estimate corn residue coverage (CRC) in Northeast China in this study. The estimation model of CRC was established for improving CRC estimation accuracy by the optimal combination of spectral indices and textural features, based on soil texture zoning, using the random forest regression method. Our results revealed that (1) the optimization C5 of spectral indices and textural features improves the CRC estimation accuracy after harvesting and before sowing with determination coefficients (R2) of 0.78 and 0.73, respectively; (2) the random forest improves the CRC estimation accuracy after harvesting and before sowing with R2 of 0.81 and 0.77, respectively; (3) considering the spatial heterogeneity of the soil background and the usage of soil texture zoning models increase the accuracy of CRC estimation after harvesting and before sowing with R2 of 0.84 and 0.81, respectively. In general, the CRC estimation accuracy after harvesting was better than that before sowing. The results revealed that the corn residue coverage in most of the study area was 0.3 to 0.6 and was mainly distributed in the Songnen Plain. By the estimated corn residue coverage results, the implementation of conservation tillage practices is identified, which is vital for protecting the black soil in Northeast China. [ABSTRACT FROM AUTHOR]

Crop-type mapping is the foundation of grain security and digital agricultural management. Accuracy, efficiency and large-scale scene consistency are required to perform crop classification from remote sensing images. Many current remote-sensing crop extraction methods based on deep learning cannot account for adaptation effects in large-scale, complex scenes. Therefore, this study proposes a novel adaptive feature-fusion network for crop classification using single-temporal Sentinel-2 images. The selective patch module implemented in the network can adaptively integrate the features of different patch sizes to assess complex scenes better. TabNet was used simultaneously to extract spectral information from the center pixels of the patches. Multitask learning was used to supervise the extraction process to improve the weight of the spectral characteristics while mitigating the negative impact of a small sample size. In the network, superpixel optimization was applied to post-process the classification results to improve the crop edges. By conducting the crop classification of peanut, rice, and corn based on Sentinel-2 images in 2022 in Henan Province, China, the novel method proposed in this paper was more accurate, indicated by an F1 score of 96.53%, than other mainstream methods. This indicates our model's potential for application in crop classification in large scenes. [ABSTRACT FROM AUTHOR]

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 (R2 = 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. [ABSTRACT FROM AUTHOR]

Efficient, fast, and accurate crop lodging monitoring is urgent for farmers, agronomists, insurance loss adjusters, and policymakers. This study aims to explore the potential of Chinese GF-1 PMS high-spatial-resolution images for corn lodging monitoring and to find a robust and efficient way to identify corn lodging accurately and efficiently. Three groups of image features and five machine-learning approaches are used for classifying non-lodged, moderately lodged, and severely lodged areas. Our results reveal that (1) the combination of spectral bands, optimized vegetation indexes, and texture features classify corn lodging with an overall accuracy of 93.81% and a Kappa coefficient of 0.91. (2) The random forest is an efficient, robust, and easy classifier to identify corn lodging with the F1-score of 0.95, 0.92, and 0.95 for non-lodged, moderately lodged, and severely lodged areas, respectively. (3) The GF-1 PMS image has great potential for identifying corn lodging on a regional scale. [ABSTRACT FROM AUTHOR]

Sun-induced chlorophyll fluorescence (SIF) has shown potential in quantifying plant responses to environmental changes by which abiotic drivers are dominated. However, SIF is a mixed signal influenced by factors such as leaf physiology, canopy structure, and sun-sensor geometry. Whether the physiological information contained in SIF can better quantify crop disease stresses dominated by biological drivers, and clearly explain the physiological variability of stressed crops, has not yet been sufficiently explored. On this basis, we took winter wheat naturally infected with stripe rust as the research object and conducted a study on the responses of physiological signals and reflectivity spectrum signals to crop disease stress dominated by biological drivers, based on in situ canopy-scale and leaf-scale data. Physiological signals include SIF, SIFyield (normalized by absorbed photosynthetically active radiation), fluorescence yield (ΦF) retrieved by NIRvP (non-physiological components of canopy SIF) and relative fluorescence yield (ΦF-r) retrieved by near-infrared radiance of vegetation (NIRvR). Reflectance spectrum signals include normalized difference vegetation index (NDVI) and near-infrared reflectance of vegetation (NIRv). At the canopy scale, six signals reached extremely significant correlations (P < 0.001) with disease severity levels (SL) under comprehensive experimental conditions (SL without dividing the experimental samples) and light disease conditions (SL < 20%). The strongest correlation between NDVI and SL (R = 0.69) was observed under the comprehensive experimental conditions, followed by NIRv (R = 0.56), ΦF-r (R = 0.53) and SIF (R = 0.51), and the response of ΦF (R = 0.45) and SIFyield (R = 0.34) to SL was weak. Under lightly diseased conditions, ΦF-r (R = 0.62) showed the strongest response to disease, followed by SIFyield (R = 0.60), SIF (R = 0.56) and NIRv (R = 0.54). The weakest correlation was observed between ΦF and SL (R = 0.51), which also showed a result approximating NDVI (R = 0.52). In the case of a high level of crop disease severity, NDVI showed advantages in disease monitoring. In the early stage of crop diseases, which we pay more attention to, compared with SIF and reflectivity spectrum signals, ΦF-r estimated by the newly proposed 'NIRvR approach' (which uses SIF together with NIRvR (i.e., SIF/ NIRvR) as a substitute for ΦF) showed superior ability to monitor crop physiological stress, and was more sensitive to plant physiological variation. At the leaf scale, the response of SIF to SL was stronger than that of NDVI. These results validate the potential of ΦF-r estimated by the NIRvR approach to monitoring disease stress dominated by biological drivers, thus providing a new research avenue for quantifying crop responses to disease stress. [ABSTRACT FROM AUTHOR]

Lodging is one of the typical abiotic adversities during paddy rice growth. In addition to affecting photosynthesis, it can seriously damage crop growth and development, such as reducing rice quality and hindering automated harvesting. It is, therefore, imperative to accurately and in good time acquire crop-lodging areas for yield prediction, agricultural insurance claims, and disaster-management decisions. However, the accuracy requirements for crop-lodging monitoring remain challenging due to complicated impact factors. Aiming at identifying paddy rice lodging on Shazai Island, Guangdong, China, caused by heavy rainfall and strong wind, a decision-tree model was constructed using multiple-parameter information from Sentinel-1 SAR images and the in situ lodging samples. The model innovatively combined the five backscattering coefficients with five polarization decomposition parameters and quantified the importance of each parameter feature. It was found that the decision-tree method coupled with polarization decomposition can be used to obtain an accurate distribution of paddy rice-lodging areas. The results showed that: (1) Radar parameters can capture the changes in lodged paddy rice. The radar parameters that best distinguish paddy rice lodging are VV, VV+VH, VH/VV, and Span. (2) Span is the parameter with the strongest feature importance, which shows the necessity of adding polarization parameters to the classification model. (3) The dual-polarized Sentinel-1 database classification model can effectively extract the area of lodging paddy rice with an overall accuracy of 84.38%, and a total area precision of 93.18%. These observations can guide the future use of SAR-based information for crop-lodging assessment and post-disaster management. [ABSTRACT FROM AUTHOR]

Accurate and timely regional crop yield information, particularly field-level yield estimation, is essential for commodity traders and producers in planning production, growing, harvesting, and other interconnected marketing activities. In this study, we propose a novel data assimilation framework. Firstly, we construct the likelihood constraints for a process-based crop growth model based on the previous year's statistical yield and the current year's field observations. Then, we infer the posterior sets of model-simulated time-series LAI and the final yield of winter wheat with a Markov chain Monte Carlo (MCMC) method for each meteorological data grid of the European Centre for Medium-Range Weather Forecasts Reanalysis (v5ERA5). Finally, we estimate the winter wheat yield at the spatial resolution of 10 m by combining Sentinel-2 LAI and the WOFOST model in Hengshui, the prefecture-level city of Hebei province of China. The results show that the proposed framework can estimate the winter wheat yield with a coefficient of determination R2 equal to 0.29 and mean absolute percentage error MAPE equal to 7.20% compared to within-field measurements. However, the agricultural stress that crop growth models cannot quantitatively simulate, such as lodging, can greatly reduce the accuracy of yield estimates. [ABSTRACT FROM AUTHOR]

Urban forms are closely related to the urban environment, providing great potential to analyze human socioeconomic activities. However, limited studies have investigated the impacts of three-dimensional (3-D) urban forms on socioeconomic activities across cities. In this paper, we explored the relationship between urban form and socioeconomic activities using 3-D building height data from 38 cities in China. First, we aggregated the building footprint data and calculated three building indicators at the grid scale, based on which the spatial patterns of building height and road density were analyzed. Then, we examined the capacities of two-dimensional (2D)/3D urban forms in characterizing socioeconomic activities using satellite-derived nighttime light (NTL) data. Finally, we analyzed the relationship between road density distributions and building heights across 38 cities in China. Our results suggest that the building height information can improve the correlation between urban form and NTL. Different patterns of road distribution were revealed according to the distribution of road density change from the building hotspots, showing the capacity of 3-D building height data in helping characterize socioeconomic activities. Our study indicates that the 3-D building height information is of great potential to support a variety of studies in urban domains, such as population distribution and carbon emissions, with significantly improved capacities. [ABSTRACT FROM AUTHOR]

Grasslands, as an important part of terrestrial ecosystems, are facing serious threats of land degradation. Therefore, the remote monitoring of grasslands is an important tool to control degradation and protect grasslands. However, the existing methods are often disturbed by clouds and fog, which makes it difficult to achieve all-weather and all-time grassland remote sensing monitoring. Synthetic aperture radar (SAR) data can penetrate clouds, which is helpful for solving this problem. In this study, we verified the advantages of the fusion of multi-spectral (MS) and SAR data for improving classification accuracy, especially for cloud-covered areas. We also proposed an adaptive feature fusion method (the SK-like method) based on an attention mechanism, and tested two types of patch construction strategies, single-size and multi-size patches. Experiments have shown that the proposed SK-like method with single-size patches obtains the best results, with 93.12% accuracy and a 0.91 average f1-score, which is a 1.02% accuracy improvement and a 0.01 average f1-score improvement compared with the commonly used feature concatenation method. Our results show that the all-weather, all-time remote sensing monitoring of grassland is possible through the fusion of MS and SAR data with suitable feature fusion methods, which will effectively enhance the regulatory capability of grassland resources. [ABSTRACT FROM AUTHOR]

Garlic and winter wheat are major economic and grain crops in China, and their boundaries have increased substantially in recent decades. Updated and accurate garlic and winter wheat maps are critical for assessing their impacts on society and the environment. Remote sensing imagery can be used to monitor spatial and temporal changes in croplands such as winter wheat and maize. However, to our knowledge, few studies are focusing on garlic area mapping. Here, we proposed a method for coupling active and passive satellite imagery for the identification of both garlic and winter wheat in Northern China. First, we used passive satellite imagery (Sentinel-2 and Landsat-8 images) to extract winter crops (garlic and winter wheat) with high accuracy. Second, we applied active satellite imagery (Sentinel-1 images) to distinguish garlic from winter wheat. Third, we generated a map of the garlic and winter wheat by coupling the above two classification results. For the evaluation of classification, the overall accuracy was 95.97%, with a kappa coefficient of 0.94 by eighteen validation quadrats (3 km by 3 km). The user's and producer's accuracies of garlic are 95.83% and 95.85%, respectively; and for the winter wheat, these two accuracies are 97.20% and 97.45%, respectively. This study provides a practical exploration of targeted crop identification in mixed planting areas using multisource remote sensing data. [ABSTRACT FROM AUTHOR]

In the Loess Plateau (LP) of China, the vegetation degradation and soil erosion problems have been shown to be curbed after the implementation of the Grain for Green program. In this study, the LP is divided into the northwestern semi-arid area and the southeastern semi-humid area using the 400 mm isohyet. The spatial–temporal evolution of the vegetation NDVI during 2000–2015 are analyzed, and the driving forces (including factors of climate, environment, and human activities) of the evolution are quantitatively identified using the geographical detector model (GDM). The results showed that the annual mean NDVI in the entire LP was 0.529, and it decreased from the semi-humid area (0.619) to the semi-arid area (0.346). The mean value of the coefficient of variation of the NDVI was 0.1406, and it increased from the semi-humid area (0.1165) to the semi-arid area (0.1926). The annual NDVI growth rate in the entire LP was 0.0079, with the NDVI growing faster in the semi-humid area (0.0093) than in the semi-arid area (0.0049). The largest increments of the NDVI were from grassland, farmland, and woodland. The GDM results revealed that changes in the spatial distribution of the NDVI could be primarily explained by the climatic and environmental factors in the semi-arid area, such as precipitation, soil type, and vegetation type, while the changes were mainly explained by the anthropogenic factors in the semi-humid area, such as the GDP density, land-use type, and population density. The interactive analysis showed that interactions between factors strengthened the impacts on the vegetation change compared with an individual factor. Furthermore, the ranges/types of factors suitable for vegetation growth were determined. The conclusions of this study have important implications for the formulation and implementation of ecological conservation and restoration strategies in different regions of the LP. [ABSTRACT FROM AUTHOR]

Black soil is one of the most productive soils with high organic matter content. Crop residue covering is important for protecting black soil from alleviating soil erosion and increasing soil organic carbon. Mapping crop residue covered areas accurately using remote sensing images can monitor the protection of black soil in regional areas. Considering the inhomogeneity and randomness, resulting from human management difference, the high spatial resolution Chinese GF-1 B/D image and developed MSCU-net+C deep learning method are used to mapping corn residue covered area (CRCA) in this study. The developed MSCU-net+C is joined by a multiscale convolution group (MSCG), the global loss function, and Convolutional Block Attention Module (CBAM) based on U-net and the full connected conditional random field (FCCRF). The effectiveness of the proposed MSCU-net+C is validated by the ablation experiment and comparison experiment for mapping CRCA in Lishu County, Jilin Province, China. The accuracy assessment results show that the developed MSCU-net+C improve the CRCA classification accuracy from IOUAVG = 0.8604 and KappaAVG = 0.8864 to IOUAVG = 0.9081 and KappaAVG = 0.9258 compared with U-net. Our developed and other deep semantic segmentation networks (MU-net, GU-net, MSCU-net, SegNet, and Dlv3+) improve the classification accuracy of IOUAVG/KappaAVG with 0.0091/0.0058, 0.0133/0.0091, 0.044/0.0345, 0.0104/0.0069, and 0.0107/0.0072 compared with U-net, respectively. The classification accuracies of IOUAVG/KappaAVG of traditional machine learning methods, including support vector machine (SVM) and neural network (NN), are 0.576/0.5526 and 0.6417/0.6482, respectively. These results reveal that the developed MSCU-net+C can be used to map CRCA for monitoring black soil protection. [ABSTRACT FROM AUTHOR]

Light pollution, a phenomenon in which artificial nighttime light (NTL) changes the form of brightness and darkness in natural areas such as protected areas (PAs), has become a global concern due to its threat to global biodiversity. With ongoing global urbanization and climate change, the light pollution status in global PAs deserves attention for mitigation and adaptation. In this study, we developed a framework to evaluate the light pollution status in global PAs, using the global NTL time series data. First, we classified global PAs (30,624) into three pollution categories: non-polluted (5974), continuously polluted (8141), and discontinuously polluted (16,509), according to the time of occurrence of lit pixels in/around PAs from 1992 to 2018. Then, we explored the NTL intensity (e.g., digital numbers) and its trend in those polluted PAs and identified those hotspots of PAs at the global scale with consideration of global urbanization. Our study shows that global light pollution is mainly distributed within the range of 30°N and 60°N, including Europe, north America, and East Asia. Although the temporal trend of NTL intensity in global PAs is increasing, Japan and the United States of America (USA) have opposite trends due to the implementation of well-planned ecological conservation policies and declining population growth. For most polluted PAs, the lit pixels are close to their boundaries (i.e., less than 10 km), and the NTL in/around these lit areas has become stronger over the past decades. The identified hotspots of PAs (e.g., Europe, the USA, and East Asia) help support decisions on global biodiversity conservation, particularly with global urbanization and climate change. [ABSTRACT FROM AUTHOR]

Global biophysical products at decametric resolution derived from Sentinel-2 imagery have emerged as a promising dataset for fine-scale ecosystem modeling and agricultural monitoring. Evaluating uncertainties of different Sentinel-2 biophysical products over various regions and vegetation types is pivotal in the application of land surface models. In this study, we quantified the performance of Sentinel-2-derived Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), and Fractional Vegetation Cover (FVC) estimates using global ground observations with consistent measurement criteria. Our results show that the accuracy of vegetation and non-vegetated classification based on Sentinel-2 surface reflectance products is greater than 95%, which indicates the vegetation identification is favorable for the practical application of biophysical estimates, as several LAI, FAPAR, and FVC retrievals were derived for non-vegetated pixels. The rate of best retrievals is similar between LAI and FAPAR estimates, both accounting for 87% of all vegetation pixels, while it is almost 100% for FVC estimates. Additionally, the Sentinel-2 FAPAR and FVC estimates agree well with ground-measurements-derived (GMD) reference maps, whereas a large discrepancy is observed for Sentinel-2 LAI estimates by comparing with both GMD effective LAI (LAIe) and actual LAI (LAI) reference maps. Furthermore, the uncertainties of Sentinel-2 LAI, FAPAR and FVC estimates are 1.09 m2/m2, 1.14 m2/m2, 0.13 and 0.17 through comparisons to ground LAIe, LAI, FAPAR, and FVC measurements, respectively. Given the temporal difference between Sentinel-2 observations and ground measurements, Sentinel-2 LAI estimates are more consistent with LAIe than LAI values. The robustness of evaluation results can be further improved as long as more multi-temporal ground measurements across different regions are obtained. Overall, this study provides fundamental information about the performance of Sentinel-2 LAI, FAPAR, and FVC estimates, which imbues our confidence in the broad applications of these decametric products. [ABSTRACT FROM AUTHOR]

Understanding spatio-temporal changes in winter wheat (Triticum aestivum L) phenology and its response to temperature will be vital for adapting to climate change in the coming years. For this purpose, the heading date (HD), maturity date (MD), and length of the reproductive growth period (LRGP) were detected from the remotely sensed leaf area index (LAI) data by a threshold-based method during the harvest year 2003 to 2018 across the North China Plain. The results show that there was high spatial heterogeneity of winter wheat phenology in pixel scale across the whole area, which could not be detected in previous site-based studies. The results also verified that climate warming could explain part of the change in the HD. However, for the LRGP, the potential impact of non-climate effects should be further investigated. This study presents the spatio-temporal changes both in winter wheat phenology and corresponding mean temperature and then analyzes their relationships in pixel scale. Additionally, this study further discusses the potential impact of non-climate effects on the LRGP. [ABSTRACT FROM AUTHOR]

Nowadays, GF-1 (GF is the acronym for GaoFen which means high-resolution in Chinese) remote sensing images are widely utilized in agriculture because of their high spatio-temporal resolution and free availability. However, due to the transferrable rationale of optical satellites, the GF-1 remote sensing images are inevitably impacted by clouds, which leads to a lack of ground object's information of crop areas and adds noises to research datasets. Therefore, it is crucial to efficiently detect the cloud pixel of GF-1 imagery of crop areas with powerful performance both in time consumption and accuracy when it comes to large-scale agricultural processing and application. To solve the above problems, this paper proposed a cloud detection approach based on hybrid multispectral features (HMF) with dynamic thresholds. This approach combined three spectral features, namely the Normalized Difference Vegetation Index (NDVI), WHITENESS and the Haze-Optimized Transformation (HOT), to detect the cloud pixels, which can take advantage of the hybrid Multispectral Features. Meanwhile, in order to meet the variety of the threshold values in different seasons, a dynamic threshold adjustment method was adopted, which builds a relationship between the features and a solar altitude angle to acquire a group of specific thresholds for an image. With the test of GF-1 remote sensing datasets and comparative trials with Random Forest (RF), the results show that the method proposed in this paper not only has high accuracy, but also has advantages in terms of time consumption. The average accuracy of cloud detection can reach 90.8% and time consumption for each GF-1 imagery can reach to 5 min, which has been reduced by 83.27% compared with RF method. Therefore, the approach presented in this work could serve as a reference for those who are interested in the cloud detection of remote sensing images. [ABSTRACT FROM AUTHOR]

Accurate crop planting area information is of significance for understanding regional food security and agricultural development planning. While increasing numbers of medium resolution satellite imagery and improved classification algorithms have been used for crop mapping, limited efforts have been made in feature selection, despite its vital impacts on crop classification. Furthermore, different crop types have their unique spectral and phenology characteristics; however, the different features of individual crop types have not been well understood and considered in previous studies of crop mapping. Here, we examined an optimized strategy to integrate specific features of individual crop types for mapping an improved crop type layer in the Sanjiang Plain, a new food bowl in China, by using all Sentinel-2 time series images in 2018. First, an automatic spectro-temporal feature selection (ASTFS) method was used to obtain optimal features for individual crops (rice, corn, and soybean), including sorting all features by the global separability indices for each crop and removing redundant features by accuracy changes when adding new features. Second, the ASTFS-based optimized feature sets for individual crops were used to produce three crop probability maps with the Random Forest classifier. Third, the probability maps were then composited into the final crop layer by considering the probability of each crop at every pixel. The resultant crop layer showed an improved accuracy (overall accuracy = 93.94%, Kappa coefficient = 0.92) than the other classifications without such a feature optimizing process. Our results indicate the potential of the ASTFS method for improving regional crop mapping. [ABSTRACT FROM AUTHOR]

Continuous and accurate estimates of crop canopy leaf area index (LAI) and chlorophyll content are of great importance for crop growth monitoring. These estimates can be useful for precision agricultural management and agricultural planning. Our objectives were to investigate the joint retrieval of corn canopy LAI and chlorophyll content using filtered reflectances from Sentinel-2 and MODIS data acquired during the corn growing season, which, being generally hot and rainy, results in few cloud-free Sentinel-2 images. In addition, the retrieved time series of LAI and chlorophyll content results were used to monitor the corn growth behavior in the study area. Our results showed that: (1) the joint retrieval of LAI and chlorophyll content using the proposed joint probability distribution method improved the estimation accuracy of both corn canopy LAI and chlorophyll content. Corn canopy LAI and chlorophyll content were retrieved jointly and accurately using the PROSAIL model with fused Kalman filtered (KF) reflectance images. The relation between retrieved and field measured LAI and chlorophyll content of four corn-growing stages had a coefficient of determination (R2) of about 0.6, and root mean square errors (RMSEs) ranges of mainly 0.1–0.2 and 0.0–0.3, respectively. (2) Kalman filtering is a good way to produce continuous high-resolution reflectance images by synthesizing Sentinel-2 and MODIS reflectances. The correlation between fused KF and Sentinel-2 reflectances had an R2 value of 0.98 and RMSE of 0.0133, and the correlation between KF and field-measured reflectances had an R2 value of 0.8598 and RMSE of 0.0404. (3) The derived continuous KF reflectances captured the crop behavior well. Our analysis showed that the LAI increased from day of year (DOY) 181 (trefoil stage) to DOY 236 (filling stage), and then increased continuously until harvest, while the chlorophyll content first also increased from DOY 181 to DOY 236, and then remained stable until harvest. These results revealed that the jointly retrieved continuous LAI and chlorophyll content could be used to monitor corn growth conditions. [ABSTRACT FROM AUTHOR]

Phenotyping provides important support for corn breeding. Unfortunately, the rapid detection of phenotypes has been the major limiting factor in estimating and predicting the outcomes of breeding programs. This study was focused on the potential of phenotyping to support corn breeding using unmanned aerial vehicle (UAV) images, aiming at mining and deepening UAV techniques for comparing phenotypes and screening new corn varieties. Two geometric traits (plant height, canopy leaf area index (LAI)) and one lodging resistance trait (lodging area) were estimated in this study. It was found that stereoscopic and photogrammetric methods were promising ways to calculate a digital surface model (DSM) for estimating corn plant height from UAV images, with R2 = 0.7833 (p < 0.001) and a root mean square error (RMSE) = 0.1677. In addition to a height estimation, the height variation was analyzed for depicting and validating the corn canopy uniformity stability for different varieties. For the lodging area estimation, the normalized DSM (nDSM) method was more promising than the gray-level co-occurrence matrix (GLCM) textural features method. The estimation error using the nDSM ranged from 0.8% to 5.3%, and the estimation error using the GLCM ranged from 10.0% to 16.2%. Associations between the height estimation and lodging area estimation were done to find the corn varieties with optimal plant heights and lodging resistance. For the LAI estimation, the physical radiative transfer PROSAIL model offered both an accurate and robust estimation performance both at the middle (R2 = 0.7490, RMSE = 0.3443) and later growing stages (R2 = 0.7450, RMSE = 0.3154). What was more exciting was that the estimated sequential time series LAIs revealed a corn variety with poor resistance to lodging in a study area of Baogaofeng Farm. Overall, UAVs appear to provide a promising method to support phenotyping for crop breeding, and the phenotyping of corn breeding in this study validated this application. [ABSTRACT FROM AUTHOR]

Large-scale crop mapping provides important information in agricultural applications. However, it is a challenging task due to the inconsistent availability of remote sensing data caused by the irregular time series and limited coverage of the images, together with the low spatial resolution of the classification results. In this study, we proposed a new efficient method based on grids to address the inconsistent availability of the high-medium resolution images for large-scale crop classification. First, we proposed a method to block the remote sensing data into grids to solve the problem of temporal inconsistency. Then, a parallel computing technique was introduced to improve the calculation efficiency on the grid scale. Experiments were designed to evaluate the applicability of this method for different high-medium spatial resolution remote sensing images and different machine learning algorithms and to compare the results with the widely used nonparallel method. The computational experiments showed that the proposed method was successful at identifying large-scale crop distribution using common high-medium resolution remote sensing images (GF-1 WFV images and Sentinel-2) and common machine learning classifiers (the random forest algorithm and support vector machine). Finally, we mapped the croplands in Heilongjiang Province in 2015, 2016, 2017, which used a random forest classifier with the time series GF-1 WFV images spectral features, the enhanced vegetation index (EVI) and normalized difference water index (NDWI). Ultimately, the accuracy was assessed using a confusion matrix. The results showed that the classification accuracy reached 88%, 82%, and 85% in 2015, 2016, and 2017, respectively. In addition, with the help of parallel computing, the calculation speed was significantly improved by at least seven-fold. This indicates that using the grid framework to block the data for classification is feasible for crop mapping in large areas and has great application potential in the future. [ABSTRACT FROM AUTHOR]

Continuous monitoring of crop growth status using time-series remote sensing image is essential for crop management and yield prediction. The growing season of summer corn in the North China Plain with the period of rain and hot, which makes the acquisition of cloud-free satellite imagery very difficult. Therefore, we focused on developing image datasets with both a high temporal resolution and medium spatial resolution by harmonizing the time-series of MOD09GA Normalized Difference Vegetation Index (NDVI) images and 30-m-resolution GF-1 WFV images using the improved Kalman filter model. The harmonized images, GF-1 images, and Landsat 8 images were then combined and used to monitor the summer corn growth from 5th June to 6th October, 2014, in three counties of Hebei Province, China, in conjunction with meteorological data and MODIS Evapotranspiration Data Set. The prediction residuals ( Δ P R K ) in NDVI between the GF-1 observations and the harmonized images was in the range of −0.2 to 0.2 with Gauss distribution. Moreover, the obtained phenological curves manifested distinctive growth features for summer corn at field scales. Changes in NDVI over time were more effectively evaluated and represented corn growth trends, when considered in conjunction with meteorological data and MODIS Evapotranspiration Data Set. We observed that the NDVI of summer corn showed a process of first decreasing and then rising in the early growing stage and discuss how the temperature and moisture of the environment changed with the growth stage. The study demonstrated that the synthesized dataset constructed using this methodology was highly accurate, with high temporal resolution and medium spatial resolution and it was possible to harmonize multi-source remote sensing imagery by the improved Kalman filter for long-term field monitoring. [ABSTRACT FROM AUTHOR]

More than 90% of the sugar production in China comes from sugarcane, which is widely grown in South China. Optical image time series have proven to be efficient for sugarcane mapping. There are, however, two limitations associated with previous research: one is that the critical observations during the sugarcane growing season are limited due to frequent cloudy weather in South China; the other is that the classification method requires imagery time series covering the entire growing season, which reduces the time efficiency. The Sentinel-1A (S1A) synthetic aperture radar (SAR) data featuring relatively high spatial-temporal resolution provides an ideal data source for all-weather observations. In this study, we attempted to develop a method for the early season mapping of sugarcane. First, we proposed a framework consisting of two procedures: initial sugarcane mapping using the S1A SAR imagery time series, followed by non-vegetation removal using Sentinel-2 optical imagery. Second, we tested the framework using an incremental classification strategy based on S1A imagery covering the entire 2017–2018 sugarcane season. The study area was in Suixi and Leizhou counties of Zhanjiang city, China. Results indicated that an acceptable accuracy, in terms of Kappa coefficient, can be achieved to a level above 0.902 using time series three months before sugarcane harvest. In general, sugarcane mapping utilizing the combination of VH + VV as well as VH polarization alone outperformed mapping using VV alone. Although the XGBoost classifier with VH + VV polarization achieved a maximum accuracy that was slightly lower than the random forest (RF) classifier, the XGBoost shows promising performance in that it was more robust to overfitting with noisy VV time series and the computation speed was 7.7 times faster than RF classifier. The total sugarcane areas in Suixi and Leizhou for the 2017–2018 harvest year estimated by this study were approximately 598.95 km2 and 497.65 km2, respectively. The relative accuracy of the total sugarcane mapping area was approximately 86.3%. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]