A Global Irrigated Area Map (GIAM) has been produced for the end of the last millennium using multiple satellite sensor, secondary, Google Earth and groundtruth data. The data included: (a) Advanced Very High Resolution Radiometer (AVHRR) 3-band and Normalized Difference Vegetation Index (NDVI) 10 km monthly time-series for 1997-1999, (b) Syste`me pour l'Observation de la Terre Vegetation (SPOT VGT) NDVI 1 km monthly time series for 1999, (c) East Anglia University Climate Research Unit (CRU) rainfall 50km monthly time series for 1961-2000, (d) Global 30 Arc-Second Elevation Data Set (GTOPO30) 1 km digital elevation data of the World, (e) Japanese Earth Resources Satellite-1 Synthetic Aperture Radar (JERS-1 SAR) data for the rain forests during two seasons in 1996 and (f) University of Maryland Global Tree Cover 1 km data for 1992-1993. A single mega-file data-cube (MFDC) of the World with 159 layers, akin to hyperspectral data, was composed by re-sampling different data types into a common 1 km resolution. The MFDC was segmented based on elevation, temperature and precipitation zones. Classification was performed on the segments. Quantitative spectral matching techniques (SMTs) used in hyperspectral data analysis were adopted to group class spectra derived from unsupervised classification and match them with ideal or target spectra. A rigorous class identification and labelling process involved the use of: (a) space-time spiral curve (ST-SC) plots, (b) brightness-greenness-wetness (BGW) plots, (c) time series NDVI plots, (d) Google Earth very-high-resolution imagery (VHRI) 'zoom-in views' in over 11 000 locations, (e) groundtruth data broadly sourced from the degree confluence project (3 864 sample locations) and from the GIAM project (1 790 sample locations), (f) high-resolution Landsat-ETM+ Geocover 150m mosaic of the World and (g) secondary data (e.g. national and global land use and land cover data). Mixed classes were resolved based on decision tree algorithms and spatial modelling, and when that did not work, the problem class was used to mask and re-classify the MDFC, and the class identification and labelling protocol repeated. The sub-pixel area (SPA) calculations were performed by multiplying full-pixel areas (FPAs) with irrigated area fractions (IAFs) for every class. A 28 class GIAMwas produced and the area statistics reported as: (a) annualized irrigated areas (AIAs), which consider intensity of irrigation (i.e. sum of irrigated areas from different seasons in a year plus continuous year-round irrigation or gross irrigated areas), and (b) total area available for irrigation (TAAI), which does not consider intensity of irrigation (i.e. irrigated areas at any given point of time plus the areas left fallow but 'equipped for irrigation' at the same point of time or net irrigated areas). The AIA of the World at the end of the last millennium was 467million hectares (Mha), which is sum of the non-overlapping areas of: (a) 252Mha from season one, (b) 174Mha from season two and (c) 41Mha from continuous yea...

A satellite-derived cropland extent map at high spatial resolution (30-m or better) is a must for food and water security analysis. Precise and accurate global cropland extent maps, indicating cropland and non-cropland areas, are starting points to develop higher-level products such as crop watering methods (irrigated or rainfed), cropping intensities (e.g., single, double, or continuous cropping), crop types, cropland fallows, as well as for assessment of cropland productivity (productivity per unit of land), and crop water productivity (productivity per unit of water). Uncertainties associated with the cropland extent map have cascading effects on all higher-level cropland products. However, precise and accurate cropland extent maps at high spatial resolution over large areas (e.g., continents or the globe) are challenging to produce due to the small-holder dominant agricultural systems like those found in most of Africa and Asia. Cloud-based geospatial computing platforms and multi-date, multi-sensor satellite image inventories on Google Earth Engine offer opportunities for mapping croplands with precision and accuracy over large areas that satisfy the requirements of broad range of applications. Such maps are expected to provide highly significant improvements compared to existing products, which tend to be coarser in resolution, and often fail to capture fragmented small-holder farms especially in regions with high dynamic change within and across years. To overcome these limitations, in this research we present an approach for cropland extent mapping at high spatial resolution (30-m or better) using the 10-day, 10 to 20-m, Sentinel-2 data in combination with 16-day, 30-m, Landsat-8 data on Google Earth Engine (GEE). First, nominal 30-m resolution satellite imagery composites were created from 36,924 scenes of Sentinel-2 and Landsat-8 images for the entire African continent in 2015-2016. These composites were generated using a median-mosaic of five bands (blue, green, red, near-infrared, NDVI) during each of the two periods (period 1: January-June 2016 and period 2: July-December 2015) plus a 30-m slope layer derived from the Shuttle Radar Topographic Mission (SRTM) elevation dataset. Second, we selected Cropland/Non-cropland training samples (sample size = 9791) from various sources in GEE to create pixel-based classifications. As supervised classification algorithm, Random Forest (RF) was used as the primary classifier because of its efficiency, and when over-fitting issues of

No abstract

Abstract. Our goal is to map the rice areas of six South Asian countries using moderateresolution imaging spectroradiometer (MODIS) time-series data for the time period 2000 to 2001. South Asia accounts for almost 40% of the world's harvested rice area and is also home to 74% of the population that lives on less than $2.00 a day. The population of the region is growing faster than its ability to produce rice. Thus, accurate and timely assessment of where and how rice is cultivated is important to craft food security and poverty alleviation strategies. We used a time series of eight-day, 500-m spatial resolution composite images from the MODIS sensor to produce rice maps and rice characteristics (e.g., intensity of cropping, cropping calendar) taking data for the years 2000 to 2001 and by adopting a suite of methods that include spectral matching techniques, decision trees, and ideal temporal profile data banks to rapidly identify and classify rice areas over large spatial extents. These methods are used in conjunction with ancillary spatial data sets (e.g., elevation, precipitation), national statistics, and maps, and a large volume of field-plot data. The resulting rice maps and statistics are compared against a subset of independent field-plot points and the best available subnational statistics on rice areas for the main crop growing season (kharif season). A fuzzy classification accuracy assessment for the 2000 to 2001 rice-map product, based on field-plot data, demonstrated accuracies from 67% to 100% for individual rice classes, with an overall accuracy of 80% for all classes. Most of the mixing was within rice classes. The derived physical rice area was highly correlated with the subnational statistics with R 2 values of 97% at the district level and 99% at the state level for 2000 to 2001. These results suggest that the methods, approaches, algorithms, and data sets we used are ideal for rapid, accurate, and large-scale mapping of paddy rice as well as for generating their statistics over large areas. C 2011 Society of Photo-Optical Instrumentation Engineers (SPIE).

No abstract

Diverse irrigated areas were mapped in the Krishna River Basin (258,912 km 2 ), southern India, using an irrigated fraction approach and multiple ancillary data sources. Unsupervised classification of a monthly time series of net difference vegetation index (NDVI) images from the Moderate Resolution Imaging Spectrometer (MODIS) over January-December 2002 generated 40 classes. Nine generalized classes included five irrigated classes with distinct NDVI time signatures: continuous irrigation, double-cropped, irrigated with low biomass, minor irrigation, and groundwater irrigation. Areas irrigated by surface water began greening 45 days after groundwater-irrigated areas, which allowed separation of surface and groundwater irrigation in the classification. The fraction of each class area irrigated was determined using three different methods: ground truth data, a linear regression model calibrated to agricultural census data, and visual interpretation of Landsat TM imagery. Irrigated fractions determined by the three methods varied least for the double-cropped irrigated class (0.62-0.79) and rangeland (0.00-0.02), and most for the minor irrigated class (0.06-0.43). Small irrigated patches (,0.1 km 2 ) accounted for more irrigated area than all major surface water irrigated areas combined. The irrigated fractions of the minor and groundwater-irrigated classes differed widely by method, suggesting that mapping patchy and small irrigated areas remains challenging, but comparison of multiple data sources improves confidence in the classification and highlights areas requiring more intensive fieldwork.

Rice is the most consumed staple food in the world and a key crop for food security. Much of the world's rice is produced and consumed in Asia where cropping intensity is often greater than 100% (more than one crop per year), yet this intensity is not sufficiently represented in many land use products. Agricultural practices and investments vary by season due to the different challenges faced, such as drought, salinity, or flooding, and the different requirements such as varietal choice, water source, inputs, and crop establishment methods. Thus, spatial and temporal information on the seasonal extent of rice is an important input to decision making related to increased agricultural productivity and the sustainable use of limited natural resources. The goal of this study was to demonstrate that hyper temporal moderate-resolution imaging spectroradiometer (MODIS) data can be used to map the spatial distribution of the seasonal rice crop extent and area. The study was conducted in Bangladesh where rice can be cropped once, twice, or three times a year.MODIS normalized difference vegetation index (NDVI) maximum value composite (MVC) data at 500 m resolution along with seasonal field-plot information from year 2010 were used to map rice crop extent and area for three seasons, boro (December/January-April), aus (April/May-June/July), and aman (July/ August-November/December), in Bangladesh. A subset of the field-plot information was used to assess the pixel-level accuracy of the MODIS-derived rice area. Seasonal district-level rice area statistics were used to assess the accuracy of the rice area estimates. When compared to field-plot data, the maps of rice versus non-rice exceeded 90% accuracy in all three seasons and the accuracy of the five rice classes varied from 78% to 90% across the three seasons. On average, the MODIS-derived rice area estimates were 6% higher than the sub-national statistics during boro, 7% higher during aus, and 3% higher during the aman season. The MODIS-derived sub-national areas explained (R 2 values) 96%, 93%, and 96% of the variability at the district level for boro, aus, and aman seasons, respectively.The results demonstrated that the methods we applied for analysing and interpreting moderate spatial and high temporal resolution imagery can accurately capture the seasonal variability in rice crop extent and area. We discuss the robustness of the approach and highlight issues that must be addressed before similar methods are used across other areas of Asia where a mix of rainfed, irrigated, or supplemental irrigation permits single, double, and triple cropping in a single calendar year.

The overarching goal of this study was to establish optimal hyperspectral vegetation indices (HVIs) and hyperspectral narrowbands (HNBs) that best characterize, classify, model, and map the world's main agricultural crops. The primary objectives were: (1) crop biophysical modeling through HNBs and HVIs, (2) accuracy assessment of crop type discrimination using Wilks' Lambda through a discriminant model, and (3) meta-analysis to select optimal HNBs and HVIs for applications related to agriculture. The study was conducted using two Earth Observing One (EO-1) Hyperion scenes and other surface hyperspectral data for the eight leading worldwide crops (wheat, corn, rice, barley, soybeans, pulses, cotton, and alfalfa) that occupy 70% of all cropland areas globally. This study integrated data collected from multiple study areas in various agroecosystems of Africa, the Middle East, Central Asia, and India. Data were collected for the eight crop types in six distinct growth stages. These included (a) ¿eld spectroradiometer measurements (350-2500 nm) sampled at 1-nm discrete bandwidths, and (b) ¿eld biophysical variables (e.g., biomass, leaf area index) acquired to correspond with spectroradiometer measurements. The eight crops were described and classi¿ed using 20 HNBs. The accuracy of classifying these 8 crops using HNBs was around 95%, which was 25% better than the multi-spectral results possible from Landsat-7's Enhanced Thematic Mapper+ or EO-1's Advanced Land Imager. Further, based on this research and meta-analysis involving over 100 papers, the study established 33 optimal HNBs and an equal number of speci¿c two-band normalized difference HVIs to best model and study speci¿c biophysical and biochemical quantities of major agricultural crops of the world. Redundant bands identi¿ed in this study will help overcome the Hughes Phenomenon (or "the curse of high dimensionality")