National-Scale Cropland Mapping Based on Phenological Metrics, Environmental Covariates, and Machine Learning on Google Earth Engine

Many challenges prevail in cropland mapping over large areas, including dealing with massive volumes of datasets and computing capabilities. Accordingly, new opportunities have been opened at a breakneck pace with the launch of new satellites, the continuous improvements in data retrieval technology, and the upsurge of cloud computing solutions such as Google Earth Engine (GEE). Therefore, the present work is an attempt to automate the extraction of multi-year (2016–2020) cropland phenological metrics on GEE and use them as inputs with environmental covariates in a trained machine-learning model to generate high-resolution cropland and crop field-probabilities maps in Morocco. The comparison of our phenological retrievals against the MODIS phenology product shows very close agreement, implying that the suggested approach accurately captures crop phenology dynamics, which allows better cropland classification. The entire country is mapped using a large volume of reference samples collected and labelled with a visual interpretation of high-resolution imagery on Collect-Earth-Online, an online platform for systematically collecting geospatial data. The cropland classification product for the nominal year 2019–2020 showed an overall accuracy of 97.86% with a Kappa of 0.95. When compared to Morocco’s utilized agricultural land (SAU) areas, the cropland probabilities maps demonstrated the ability to accurately estimate sub-national SAU areas with an R-value of 0.9. Furthermore, analyzing cropland dynamics reveals a dramatic decrease in the 2019–2020 season by 2% since the 2018–2019 season and by 5% between 2016 and 2020, which is partly driven by climate conditions, but even more so by the novel coronavirus disease 2019 (COVID-19) that impacted the planting and managing of crops due to government measures taken at the national level, like complete lockdown. Such a result proves how much these methods and associated maps are critical for scientific studies and decision-making related to food security and agriculture.

COHERENCE and BACKSCATTER BASED CROPLAND MAPPING USING MULTI-TEMPORAL SENTINEL-1 with DYNAMIC TIME WARPING

Cropland mapping is an important inventory for food security and decision making operated by governments. Crop mapping is used to identify the croplands and their spatial distribution. For a reliable analysis and forecast for projection, multi-temporal data play a key role. Even current open and frequent optical satellite data such as Sentinel-2 and Landsat support monitoring, they are not always operational due to atmospheric conditions (rain, cloud cover, haze, etc.). On the other hand, Synthetic Aperture Radar (SAR) satellites provide alternative data sets compared to optical satellites since they can acquire images under all weather conditions. In this study, an annual cropland monitoring study is conducted using Sentinel-1 SAR. For the investigation, Tokat Province an agricultural region of Turkey, where the main source of income is agriculture, was selected. There are 4 different vegetation species (wheat, sunflower, sugar beet, corn) in the study area. Sentinel-1 data was used to generate time-series of each class and phenological structures of the crops. In this context, backscatter images of both vertical-vertical (VV) and vertical-horizontal (VH) polarized data, and coherence of both VV and VH were produced from Sentinel-1 data. Time-Weighted Dynamic Time-Warping (TWDTW) classification approach was used over cropland. The produced time-series are classified under different scenarios. The results showed that only coherence has provided higher accuracies about 81% compared to using only backscatter images as 49%.

Mapping Croplands in the Granary of the Tibetan Plateau Using All Available Landsat Imagery, A Phenology-Based Approach, and Google Earth Engine

The Tibetan Plateau (TP), known as “The Roof of World”, has expansive alpine grasslands and is a hotspot for climate change studies. However, cropland expansion and increasing anthropogenic activities have been poorly documented, let alone the effects of agricultural activities on food security and environmental change in the TP. The existing cropland mapping products do not depict the spatiotemporal characteristics of the TP due to low accuracies and inconsistent cropland distribution, which is affected by complicated topography and impedes our understanding of cropland expansion and its associated environmental impacts. One of the biggest challenges of cropland mapping in the TP is the diverse crop phenology across a wide range of elevations. To decrease the classification errors due to elevational differences in crop phenology, we developed two pixel- and phenology-based algorithms to map croplands using Landsat imagery and the Google Earth Engine platform along the Brahmaputra River and its two tributaries (BRTT) in the Tibet Autonomous Region, also known as the granary of TP, in 2015–2019. Our first phenology-based cropland mapping algorithm (PCM1) used different thresholds of land surface water index (LSWI) by considering varied crop phenology along different elevations. The second algorithm (PCM2) further offsets the phenological discrepancy along elevational gradients by considering the length and peak of the growing season. We found that PCM2 had a higher accuracy with fewer images compared with PCM1. The number of images for PCM2 was 279 less than PCM1, and the Matthews correlation coefficient for PCM2 was 0.036 higher than PCM1. We also found that the cropland area in BRTT was estimated to be 1979 ± 52 km2 in the late 2010s. Croplands were mainly distributed in the BRTT basins with elevations of 3800–4000 m asl. Our phenology-based methods were effective for mapping croplands in mountainous areas. The spatially explicit information on cropland area and distribution in the TP aid future research into the effects of cropland expansion on food security and environmental change in the TP.

Synergistic Use of Radar and Optical Satellite Data for Improved Monsoon Cropland Mapping in India

Monsoon crops play a critical role in Indian agriculture, hence, monitoring these crops is vital for supporting economic growth and food security for the country. However, monitoring these crops is challenging due to limited availability of optical satellite data due to cloud cover during crop growth stages, landscape heterogeneity, and small field sizes. In this paper, our objective is to develop a robust methodology for high-resolution (10 m) monsoon cropland mapping appropriate for different agro-ecological regions (AER) in India. We adapted a synergistic approach of combining Sentinel-1 Synthetic Aperture Radar (SAR) data with Normalized Difference Vegetation Index (NDVI) derived from Sentinel-2 optical data using the Google Earth Engine platform. We developed a new technique, Radar Optical cross Masking (ROM), for separating cropland from non-cropland by masking out forest, plantation, and other non-dynamic features. The methodology was tested for five different AERs in India, representing a wide diversity in agriculture, soil, and climatic variations. Our findings indicate that the overall accuracy obtained by using the SAR-only approach is 90%, whereas that of the combined approach is 93%. Our proposed methodology is particularly effective in regions with cropland mixed with tree plantation/mixed forest, typical of smallholder dominated tropical countries. The proposed agriculture mask, ROM, has high potential to support the global agriculture monitoring missions of Geo Global Agriculture Monitoring (GEOGLAM) and Sentinel-2 for Agriculture (S2Agri) project for constructing a dynamic monsoon cropland mask.