Area and spatial distribution information of paddy rice are important for understanding of food security, water use, greenhouse gas emission, and disease transmission. Due to climatic warming and increasing food demand, paddy rice has been expanding rapidly in high latitude areas in the last decade, particularly in northeastern (NE) Asia. Current knowledge about paddy rice fields in these cold regions is limited. The phenology- and pixel-based paddy rice mapping (PPPM) algorithm, which identifies the flooding signals in the rice transplanting phase, has been effectively applied in tropical areas, but has not been tested at large scale of cold regions yet. Despite the effects from more snow/ice, paddy rice mapping in high latitude areas is assumed to be more encouraging due to less clouds, lower cropping intensity, and more observations from Landsat sidelaps. Moreover, the enhanced temporal and geographic coverage from Landsat 8 provides an opportunity to acquire phenology information and map paddy rice. This study evaluated the potential of Landsat 8 images on annual paddy rice mapping in NE Asia which was dominated by single cropping system, including Japan, North Korea, South Korea, and NE China. The cloud computing approach was used to process all the available Landsat 8 imagery in 2014 (143 path/rows, ~3290 scenes) with the Google Earth Engine (GEE) platform. The results indicated that the Landsat 8, GEE, and improved PPPM algorithm can effectively support the yearly mapping of paddy rice in NE Asia. The resultant paddy rice map has a high accuracy with the producer (user) accuracy of 73% (92%), based on the validation using very high resolution images and intensive field photos. Geographic characteristics of paddy rice distribution were analyzed from aspects of country, elevation, latitude, and climate. The resultant 30-m paddy rice map is expected to provide unprecedented details about the area, spatial distribution, and landscape pattern of paddy rice fields in NE Asia, which will contribute to food security assessment, water resource management, estimation of greenhouse gas emissions, and disease control.

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...

No abstract

Agricultural land use change substantially affects climate, water, ecosystems, biodiversity, and human welfare. In recent decades, due to increasing population and food demand and the backdrop of global warming, croplands have been expanding into higher latitude regions. One such hotspot is paddy rice expansion in northeast China. However, there are no maps available for documenting the spatial and temporal patterns of continuous paddy rice expansion. In this study, we developed an automated, Landsat-based paddy rice mapping (Landsat-RICE) system that uses time series Landsat images and a phenology-based algorithm based on the unique spectral characteristics of paddy rice during the flooding/transplanting phase. As a pilot study, we analyzed all the available Landsat images from 1986 to 2010 (498 scenes) in one tile (path/row 113/27) of northeast China, which tracked paddy rice expansion in epochs with five-year increments (1986-1990, 1991-1995, 1996-2000, 2001-2005, and 2006-2010). Several maps of land cover types (barren land and built-up land; evergreen, deciduous and sparse vegetation types; and water-related land cover types such as permanent water body, mixed pixels of water and vegetation, spring flooded wetlands and summer flooded land) were generated as masks. Air temperature was used to define phenology timing and crop calendar, which were then used to select Landsat images in the phenology-based algorithms for paddy rice and masks. The resultant maps of paddy rice in the five epochs were evaluated using validation samples from multiple sources, and the overall accuracies and Kappa coefficients ranged from 84 to 95% and 0.6-0.9, respectively. The paddy rice area in the study area substantially increased from 1986 to 2010, particularly after the 1990s. This study demonstrates the potential of the Landsat-RICE system and time series Landsat images for tracking agricultural land use changes at 30-m resolution in the temperate zone with single crop cultivation.

Knowledge of the area and spatial distribution of paddy rice is important for assessment of food security, management of water resources, and estimation of greenhouse gas (methane) emissions. Paddy rice agriculture has expanded rapidly in northeastern China in the last decade, but there are no updated maps of paddy rice fields in the region. Existing algorithms for identifying paddy rice fields are based on the unique physical features of paddy rice during the flooding and transplanting phases and use vegetation indices that are sensitive to the dynamics of the canopy and surface water content. However, the flooding phenomena in high latitude area could also be from spring snowmelt flooding. We used land surface temperature (LST) data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor to determine the temporal window of flooding and rice transplantation over a year to improve the existing phenology-based approach. Other land cover types (e.g., evergreen vegetation, permanent water bodies, and sparse vegetation) with potential influences on paddy rice identification were removed (masked out) due to their different temporal profiles. The accuracy assessment using high-resolution images showed that the resultant MODIS-derived paddy rice map of northeastern China in 2010 had a high accuracy (producer and user accuracies of 92% and 96%, respectively). The MODIS-based map also had a comparable accuracy to the 2010 Landsat-based National Land Cover Dataset (NLCD) of China in terms of both area and spatial pattern. This study demonstrated that our improved algorithm by using both thermal and optical MODIS data, provides a robust, simple and automated approach to identify and map paddy rice fields in temperate and cold temperate zones, the northern frontier of rice planting.