Sentinel-2 is a wide-swath and fine spatial resolution satellite imaging mission designed for data continuity and enhancement of the Landsat and other missions. The Sentinel-2 data are freely available at the global scale, and have similar wavelengths and the same geographic coordinate system as the Landsat data, which provides an excellent opportunity to fuse these two types of satellite sensor data together. In this paper, a new approach is presented for the fusion of Landsat 8 Operational Land Imager and Sentinel-2 Multispectral Imager data to coordinate their spatial resolutions for continuous global monitoring. The 30 m spatial resolution Landsat 8 bands are downscaled to 10 m using available 10 m Sentinel-2 bands. To account for the land-cover/land-use (LCLU) changes that may have occurred between the Landsat 8 and Sentinel-2 images, the Landsat 8 panchromatic (PAN) band was also incorporated in the fusion process. The experimental results showed that the proposed approach is effective for fusing Landsat 8 with Sentinel-2 data, and the use of the PAN band can decrease the errors introduced by LCLU changes. By fusion of Landsat 8 and Sentinel-2 data, more frequent observations can be produced for continuous monitoring (this is particularly valuable for areas that can be covered easily by clouds, thereby, contaminating some Landsat or Sentinel-2 observations), and the observations are at a consistent fine spatial resolution of 10 m. The products have great potential for timely monitoring of rapid changes.

Sentinel-1A synthetic aperture radar (SAR) data present an opportunity for acquiring crop information without restrictions caused by weather and illumination conditions, at a spatial resolution appropriate for individual rice fields and a temporal resolution sufficient to capture the growth profiles of different crop species. This study investigated the use of multi-temporal Sentinel-1A SAR data and Landsat-derived normalized difference vegetation index (NDVI) data to map the spatial distribution of paddy rice fields across parts of the Sanjiang plain, in northeast China. The satellite sensor data were acquired throughout the rice crop-growing season (May–October). A co-registered set of 10 dual polarization (VH/VV) SAR and NDVI images depicting crop phenological development were used as inputs to Support Vector Machine (SVM) and Random Forest (RF) machine learning classification algorithms in order to map paddy rice fields. The results showed a significant increase in overall classification when the NDVI time-series data were integrated with the various combinations of multi-temporal polarization channels (i.e. VH, VV, and VH/VV). The highest classification accuracies overall (95.2%) and for paddy rice (96.7%) were generated using the RF algorithm applied to combined multi-temporal VH polarization and NDVI data. The SVM classifier was most effective when applied to the dual polarization (i.e. VH and VV) SAR data alone and this generated overall and paddy rice classification accuracies of 91.6% and 82.5%, respectively. The results demonstrate the practicality of implementing RF or SVM machine learning algorithms to produce 10 m spatial resolution maps of paddy rice fields with limited ground data using a combination of multi-temporal SAR and NDVI data, where available, or SAR data alone. The methodological framework developed in this study is apposite for large-scale implementation across China and other major rice-growing regions of the world.

Spatio-temporal fusion of MODIS and Landsat data aims to produce new data that have simultaneously the Landsat spatial resolution and MODIS temporal resolution. It is an ill-posed problem involving large uncertainty, especially for reproduction of abrupt changes and heterogeneous landscapes. In this paper, we proposed to incorporate the freely available 250 m MODIS images into spatio-temporal fusion to increase prediction accuracy. The 250 m MODIS bands 1 and 2 are fused with 500 m MODIS bands 3-7 using the advanced area-to-point regression kriging approach. Based on a standard spatio-temporal fusion approach, the interim 250 m fused MODIS data are then downscaled to 30 m with the aid of the available 30 m Landsat data on temporally close days. The 250 m data can provide more information for the abrupt changes and heterogeneous landscapes than the original 500 m MODIS data, thus increasing the accuracy of spatio-temporal fusion predictions. The effectiveness of the proposed scheme was demonstrated using two datasets.

Rice is a globally significant staple food crop. Therefore, it is crucial to have adequate tools for monitoring changes in the extent of rice paddy cultivation. Such a system would require a sustainable and operational workflow that employs open-source medium to high spatial and temporal resolution satellite imagery and efficient classification techniques. This study used similar phenological data from Sentinel-2 (S2) optical and Sentinel-1 (S1) Synthetic Aperture Radar (SAR) satellite imagery to identify paddy rice distribution with deep learning (DL) techniques. Using Google Earth Engine (GEE) and U-Net Convolutional Neural Networks (CNN) segmentation, a workflow that accurately delineates smallholder paddy rice fields using multi-temporal S1 SAR and S2 optical imagery was investigated. The study′s accuracy assessment results showed that the optimal dataset for paddy rice mapping was a fusion of S2 multispectral bands (visible and near infra-red (VNIR), red edge (RE) and short-wave infrared (SWIR)), and S1-SAR dual polarization bands (VH and VV) captured within the crop growing season (i.e., vegetative, reproductive, and ripening). Compared to the random forest (RF) classification, the DL model (i.e., ResU-Net) had an overall accuracy of 94% (three percent higher than the RF prediction). The ResU-Net paddy rice prediction had an F1-Score of 0.92 compared to 0.84 for the RF classification generated using 500 trees in the model. Using the optimal U-Net classified paddy rice maps for the dates analyzed (i.e., 2016–2020), a change detection analysis over two epochs (2016 to 2018 and 2018 to 2020) provided a better understanding of the spatial–temporal dynamics of paddy rice agriculture in the study area. The results indicated that 377,895 and 8551 hectares of paddy rice fields were converted to other land-use over the first (2016–2018) and second (2018–2020) epochs. These statistics provided valuable insight into the paddy rice field distribution changes across the selected districts analyzed. The proposed DL framework has the potential to be upscaled and transferred to other regions. The results indicated that the approach could accurately identify paddy rice fields locally, improve decision making, and support food security in the region.

Satellite data holds considerable potential as a source of information on rice crop growth which can be used to inform agronomy. However, given the typical field sizes in many rice-growing countries such as China, data from coarse spatial resolution satellite systems such as the Moderate Resolution Imaging Spectroradiometer (MODIS) are inadequate for resolving crop growth variability at the field scale. Nevertheless, systems such as MODIS do provide images with sufficient frequency to be able to capture the detail of rice crop growth trajectories throughout a growing season. In order to generate high spatial and temporal resolution data suitable for mapping rice crop phenology, this study fused MODIS data with lower frequency, higher spatial resolution Landsat data. An overall workflow was developed which began with image preprocessing, calculation of multi-temporal normalized difference vegetation index (NDVI) images, and spatiotemporal fusion of data from the two sensors. The Spatial and Temporal Adaptive Reflectance Fusion Model was used to effectively downscale the MODIS data to deliver a time-series of 30 m spatial resolution NDVI data at 8-day intervals throughout the rice-growing season. Zonal statistical analysis was used to extract NDVI time-series for individual fields and signal filtering was applied to the time-series to generate rice phenology curves. The downscaled MODIS NDVI products were able to characterize the development of paddy rice at fine spatial and temporal resolutions, across wide spatial extents over multiple growing seasons. These data permitted the extraction of key crop seasonality parameters that quantified inter-annual growth variability for a whole agricultural region and enabled mapping of the variability in crop performance between and within fields. Hence, this approach can provide rice crop growth data that is suitable for informing agronomic policy and practice across a wide range of scales.