Abstract
As farmland systems vary over space and time (season and year), accurate and updated maps of paddy rice are needed for studies of food security and environmental problems. We selected a wheat-rice double-cropped area from fragmented landscapes along the ruralâurban complex (Jiangsu Province, China) and explored the potential utility of integrating time series optical images (Landsat-8, MODIS) and radar images (PALSAR) in mapping paddy rice planting areas. We first identified several main types of non-cropland land cover and then identified paddy rice fields by selecting pixels that were inundated only during paddy rice flooding periods. These key temporal windows were determined based on MODIS Land Surface Temperature and vegetation indices. The resultant paddy rice map was evaluated using regions of interest (ROIs) drawn from multiple high-resolution images, Google Earth and in-situ cropland photos. The estimated overall accuracy and Kappa coefficient were 89.8% and 0.79, respectively. In comparison with the National Land Cover Data (China) from 2010, the resultant map better detected changes in the paddy rice fields and revealed more details about their distribution. These results demonstrate the efficacy of using images from multiple sources to generate paddy rice maps for two-crop rotation systems.
Similar content being viewed by others
Introduction
Studies on paddy rice fields aim to provide direct or indirect information for researches on food security, water resource management and environmental sustainability. Paddy rice fields provide one of the main staple foods for more than half of the worldâs population with 11% of cultivated land1. In Asia, the majority of rice agriculture relies on irrigation, accounting for 70% of global fresh water withdrawals2. Determining the area of paddy rice fields is an important component of obtaining more accurate information about agricultural water use to effectively manage fresh water resources3. In addition, as a kind of cultivated wetland, seasonally flooded paddy fields contribute 10â13% of the atmosphereâs anthropogenic methane4. Meanwhile, paddy rice fields are changing at a breakneck pace due to dramatic encroachment by expanding cities5 and the potentially reduced availability of water resources caused by climate change6. Therefore, it is urgently necessary to update and refine information about paddy rice planting areas in order to efficiently and accurately estimate crop production7, manage water resources8,9 and monitor greenhouse gas emissions10.
At the global and regional scales, several early studies of ecosystems and land cover have involved mapping paddy rice fields based on agricultural census data. In the late 1980s and early 1990s, several paddy rice datasets with coarse spatial resolution were produced to analyze global climate and greenhouse gas emission11,12. In the years following, two global cropland datasets representing crop patterns in the early 1990s and in the year 2000 were created at a spatial resolution of 5 arc minutes (~10âkm)13,14. At the regional scale, a rice dataset for Asia was developed at the beginning of the 1980s15. In recent years, several studies on paddy rice planting areas were conducted by combining agricultural census data at the national scale16,17. Although most of the crop datasets were produced with input from multiple sources, these datasets were developed mainly by relying on statistics with coarse spatial resolution that were supplied by administrative units. Given the limitations of this spatial and temporal information, it is a challenge to apply these datasets to finer spatial research and to update them year to year.
Remote sensing is an efficient technique to acquire temporal and spatial cropland information repeatedly and consistently18. Historically, two main kinds of satellites have been used to map paddy rice fields: microwave and optical. Microwave satellites can penetrate through clouds and are thus superior for mapping paddy rice in regions dominated by long-term cloudy and rainy weather19,20,21,22, but available synthetic aperture radar (SAR) imagery is limited19 or expensive23. Commonly used optical sensors are the Multispectral Scanner System/Thematic Mapper/Enhanced Thematic Mapper Plus (MSS/TM/ETM+), the Moderate Resolution Imaging Spectroradiometer (MODIS)24,25,26,27,28,29, NOAAâs Advanced Very High Resolution Radiometer (AVHRR)30,31 and SPOT High Resolution Geometrical/High Resolution Visible Infrared/VEGETATION (HRG/HRVIR/VGT)32,33,34. Owing to continuous archiving and free to acquire, MODIS and Landsat have been the prevalent data sources of mapping paddy rice fields over the last several years.
Optical image classification techniques often use either individual image(s) or individual pixel(s) of time series data as input, here namely image-based methods and pixel-based methods. Image-based methods quantify the relationships (such as similarities or differences in spectra or texture) among all the pixels in an image for classification or object detection. These image-based methods have been applied to paddy rice mapping using single or multi-temporal optical images (e.g., MODIS and Landsat) at regional scales23,35,36,37. For image-based methods, collection of training samples (pixels) from ground reference data in each corresponding year remains a challenge38. Pixel-based methods rely primarily on the time series data for a pixel. These methods track the seasonal dynamics of a type of land cover and provide phenology information about the land surface. Several pixel-based algorithms have been developed to classify cropland using various optical images (e.g., MODIS, Landsat, SPOT and FORMOSAT-2)32,38,39,40. To take into account the phenology of paddy rice, a pixel- and phenology-based algorithm using time series data of vegetation indices has been proposed. This was successfully applied to VGT and MODIS data for southern China and South and Southeast Asia28,29. Due to the high temporal resolution and continuous observation, MODIS data have frequently been combined with pixel- and phenology-based algorithms to track paddy rice phenology in order to map paddy rice areas or detect the intensity of paddy rice fields24,26,27,41,42,43. However, MODIS-based phenology information cannot capture the sub-pixel dynamics of small paddy rice fields in heterogeneous and fragmented agricultural landscapes26,44. This could be improved by using Landsat images with 30-m spatial resolution. Nevertheless, more research is needed to document the combination of Landsat images with pixel- and phenology-based algorithms for paddy rice mapping, especially in the case of Landsat 8 OLI images.
Reduced Landsat data availability caused by cloud cover or other problems may result in the failure of Landsat time series with 16-day intervals to distinguish the crops and trees. This problem is obvious in regions with double or multiple rotation agricultural systems, where cropland tends to be covered by plants year round, resulting in unavoidable confusion with natural evergreen forest. Fortunately, this problem can be solved by the Phased Array Type L-band Synthetic Aperture Radar (PALSAR) onboard the Advanced Land Observing Satellite (ALOS), given its capacity for forest detection45,46. The current study, which involves mapping paddy rice planting areas at a 30-m spatial resolution, has two aims: (1) to develop a pixel- and phenology-based algorithm by integrating time series optical images (Landsat-8, MODIS) and radar images (PALSAR) to map paddy rice planting areas in wheat-rice agricultural systems; and (2) to evaluate the potential utility of Landsat-8 OLI data in identifying fragmented paddy rice fields in complex agricultural landscapes. The case study area is located at the Yangzi-Huaihe Plain, China, which is characterized by a two-crop rotation (wheat and rice) agricultural system and intermixture of rural and urban landscapes (Figure S1).
Results
Spectral signatures of flooded pixels and other land cover types
Different characteristics of vegetation indices are the basis for distinguishing flooded paddies from other land cover types. Figure S2 shows one example of the maps of Enhanced Vegetation Index (EVI), Normalized Difference Vegetation Index (NDVI), Land Surface Water Index (LSWI), LSWI-EVI and LSWI-NDVI on Julian day 191 (July 10, 2013). At this time, paddy rice fields are in the midst of the flooding periods (including flooding/transplanting period and reviving period of paddy rice calendar), covered by a mixture of water and plants. Some uplands (e.g. corn fields) are going through the seeding or three leaves periods, covered by soil or by a mixture of soil and plants. The scatterplot graphs (Fig. 1(a,b)) show that EVI and NDVI mainly cluster between 0 and 0.3 for both the paddies and the uplands. The tremendous discrepancy between the two is that paddies cluster with LSWI-EVIââ¥â0 or LSWI-NDVIââ¥â0 and uplands cluster with LSWI-EVIâ<â0 or LSWI-NDVIâ<â0. Water bodies have LSWI-EVIââ¥â0 or LSWI-NDVIââ¥â0, but their EVI and NDVI are close to or less than zero. Other vegetation like forests or shrubs have higher EVI or NDVI (larger than 0.4) and lower LSWI-EVI or LSWI-NDVI (less than 0). Figure 1(c) shows that more paddy fields are detected by using LSWI-EVIââ¥â0 than LSWI-NDVIââ¥â0. In this image, about 10% of the pixels are flooded.
(a,b) show the gathering of four main objects (paddy rice fields, uplands, other vegetation and water) in two-dimensional scatter plots: EVI and LSWI-EVI and NDVI and LSWI-NDVI on Julian day 191(July 10, 2013). (c) Frequency histograms of LSWI-EVI and LSWI-NDVI. It shows that LSWI-EVIââ¥â0 detects more paddy rice fields than LSWI-NDVIââ¥â0. (a,b) created in ENVI 5.0, (c) created in Sigmaplot 12.0.
Spatio-temporal dynamics of paddy rice fields
Two Landsat-8 images (June 24 & July 10, 2013) were acquired during the flooding periods. Figure 2(a,b) show the original images with bands combination: R: SWIR1, G: NIR, B: Green. Figure 2(c,d,e,f) represent relevant paddy rice maps identified by the criteria LSWI-NDVIââ¥â0 and by LSWI-EVIââ¥â0, respectively. Figure 2(g,h) are the combined paddy rice maps identified by the criteria LSWI-NDVIââ¥â0 or LSWI-EVIââ¥â0. They also show cloud/cloud shadow, water body, built-up/barren land and forest masks. Because of less cloud/cloud shadow cover, the image quality on July 10 is better than the one on June 24 and more paddy rice fields are detected from it.
(a,b) Landsat-8 images after atmospheric correction on June 24, 2013 and July 10, 2013 (R,G,Bâ=âSWIR, NIR, Green); (c,d) Flooding pixels identified by the criteria LSWI-NDVIââ¥â0; (e,f) Flooding pixels identified by the criteria LSWI â EVIââ¥â0; (g,h) Paddy rice mappings identified by the criteria LSWI-NDVIââ¥â0 or LSWI - EVIââ¥â0. Clouds/cloud shadows, water, built-up/barren land, forest are all shown. Map created in ArcMap 10.1. (a,b) were downloaded from Earth Resources Observation and Science (EROS) Center, USGS ( http://earthexplorer.usgs.gov/).
Spatial distribution of paddy rice fields
Figure 3 shows the spatial distribution of paddy rice fields in 2013 at 30âm spatial resolution, which integrates the paddy rice maps from June 24 and July 10. Few paddy rice fields are detected within the Shandong Province. However, in the Jiangsu Province, paddy rice fields are distributed widely and extensively. The paddy rice planting area is estimated to be approximately 2406.0âkm2, accounting for 37.6% of the study area. Additionally, this map reveals that the rice agriculture here is mainly conducted in small or medium size croplands. We used the Fragstats software 4.2 ( http://www.umass.edu/landeco/research/fragstats/downloads/fragstats_downloads.html#diagnostic) to evaluate the fragmentation of the paddy rice fields. In this map, the total number of paddy field patches is 180,135 and the mean patch size of paddy fields is 0.0238âkm2 (~154âmâÃâ154âm).
The paddy rice planting area map at 30âm spatial resolution, identified through the criteria LSWI - NDVIââ¥â0 or LSWI - EVIââ¥â0 from Landsat-8 images on June 24, 2013 and July 10, 2013. The total number of paddy rice plots is 180,135 and the mean paddy rice field size is 0.0238âkm2 (~154âm*154âm). Map created in ArcMap 10.1.
Evaluation of Landsat-derived rice map
In this work, 3,610 paddy rice pixels (199 ROIs) and 3,113 non-paddy rice pixels (85 ROIs), located in good observation regions, were employed to calculate the confusion matrix (Table 1).
The validation showed the paddy rice/non-paddy rice map has a reasonably high accuracy. The paddy rice planting area in this map yielded 83.1% producer accuracy and 97.5% user accuracy. The non-paddy rice area in this map has 97.6% producer accuracy and 83.3% user accuracy. The overall accuracy and Kappa coefficient of this map are 89.8% and 0.79, respectively.
Comparison to other available datasets
To compare this result with NLCD2010 (1 km gridded 2010 National Land Cover Data, China) Paddy data, the 30âm Landsat-8 paddy rice planting area mapping was resampled to have 1âkm spatial resolution (PaddyLandsat-8, Fig. 4(a)). In general, the spatial aggregation of paddy rice fields in both mappings (Fig. 4(a,c)) was consistent except for two areas marked with blue circles (A, B). Furthermore, PaddyLandsat-8 revealed more details about paddy rice field patterns than PaddyNLCD2010 data (Fig. 4(c)).
(a) 1âkm Landsat-8 paddy rice map (PaddyLandsat-8) and different colors presenting the occupation levels of paddy rice fields in given pixels; (b) The cloud/cloud shadow masks for images on June 24 and July 10, 2013; (c) Paddy rice map of 1âkm NLCD2010 dataset (PaddyNLCD2010); two significant discrepancies between PaddyLandsat-8 and PaddyNLCD2010 are marked with blue circles (A, B); (d) Upland map of 1âkm NLCD2010 dataset; (e) WorldView-2 image from June 20, 2012, the location of which is marked by the black rectangle in (a, b, c, d). It shows abundant flooded signals in the croplands, which should be classified as paddy rice fields just as in the results of this study, rather than uplands in NLCD2010. (f) Pixel-level comparison between these two datasets. The solid line was drawn from the regression analysis of all the points except the grey ones gathering on the X-axis. These grey points reveal that PaddyLandsat-8 detected new paddy rice fields that were absent in PaddyNLCD2010. Maps (a, b, c, d, e) created in ArcMap 10.1 and Map (f) created in SigmaPlot 12.0. Image (e) was provided by NASA for use in the NASA projects.
The total area of paddy rice planting in 2013 estimated by PaddyLandsat-8 was 2406.0âkm2. It was far less than the paddy rice area (4,986.96âkm2) in 2010 estimated by PaddyNLCD2010. Just considering the paddy rice fields in good observation regions, the paddy rice area from PaddyLandsat-8 was 1,311.14âkm2, 12.2% lower than that (1,492.85âkm2) from PaddyNLCD2010. Figure 4(b) shows the cloud/cloud shadow masks of images on June 24 and July 23, 2013. For the difference marked by the blue circle A, one obvious explanation is cloud contamination on June 24, 2013 (Fig. 4(b)). Likewise, visual analysis showed a large difference occurred in the blue circle B (Fig. 4(a,c)). PaddyLandsat-8 detected some paddy rice fields in this area, while NLCD2010 data classified these croplands as uplands (Fig. 4(d)). One WordView-2 image (Fig. 4(e)) on June 20, 2012, its location marked by the black rectangle, was used to judge the performance of these two results. It shows there are abundant flooding signals in the croplands that should be classified as paddies rather than uplands.
At the pixel level, the regression analysis was carried out to compare these two datasets for the paddy rice area estimate. R-square was 0.56 without considering the grey points gathering on X-axis (Fig. 4(f)). One possible reason for this difference was that PaddyLandsat-8 detected new paddy rice fields that were absent in PaddyNLCD2010, just as the paddy rice fields within the blue circle B. In addition, owing to PaddyNLCD2010 referring to the paddy rice planting area in 2010 and PaddyLandsat-8 focusing on that in 2013, there were uncertainties due to paddy rice changes in actual agricultural practice or rapid Land Use/Land cover changes.
Discussion
In this study, multi-source remote sensing data including Landsat-8, MODIS and PALSAR were used to identify paddy rice fields from the rotation agricultural system of winter wheat and paddy rice. 30âm Landsat OLI images provide more details about the distribution of paddy rice fields (Fig. 4(a)) and they are also helpful in identifying paddy rice from heterogeneous crops, which are common in this study area (Figure S5). MODIS-based LSTnight images have the ability to track the growing season, which is consistent with the results of climate observations (Figure S3). The dynamics of MODIS-based vegetation indices give it the ability to detect the flooding signals (including flooding/transplanting and reviving periods) of paddy rice fields (Fig. 5(a)) and yield results consistent with the crop calendar (Figure S4). Therefore, MODIS-based LSTnight and vegetation indices (VIs) dynamics can help automatically select Landsat-8 images for paddy rice planting area mapping. PALSAR data have advantages in mapping forest, which can be used to generate forest masks in paddy rice mapping to solve the mixture of croplands and forest, especially in the areas with long growing seasons. Therefore, the combination of multi-source remote sensing data makes it possible to map paddy rice fields automatically in complicated agricultural systems.
The seasonal dynamics of NDVI, EVI and LSWI, extracted from MOD09A1 product in 2012 and Landsat-8 images with good quality observations in 2013, for two-crop rotation agricultural systems. (a) Winter wheat and paddy rice rotation (34.286â°N, 119.642â°E), (b) winter wheat and corn rotation (34.237â°N, 119.232â°E).
The integration of Landsat-8, MODIS and PALSAR data also showed some advantages in the identification of fragmented paddy rice fields. The mean patch size of paddy fields in the study area was 154âmâÃâ154âm in 2013. In terms of relatively small and even fragmented rice fields, the Landsat-8 images with finer spatial resolution (30âm) were better than the MODIS satellite data (500âm) in reducing mixed-pixel problems. In last several decades, the croplands in this study area were becoming fragmented rapidly, because of high population pressure, rapid economic development, urbanization and limited arable lands. In 1983, the mean farmland area was about 14âkm2, reduced to 6âkm2 in 200047. In neighboring districts, the mean patch size of paddy rice fields also showed rapid diminution, declining from ~12âkm2 in 1990 to ~1.4âkm2 by 200648. At the same time, the minimum patch area was also reduced from 2âÃâ10â4âkm2 in 1990 to 5âÃâ10â5âkm2 in 2006. Cropland fragmentation is common in Asia40. In China, the crop land area per household was 5.3âÃâ10â3âkm2 divided into 6.06 plots on average49. The net cultivated area per capita was lower than 6âÃâ10â4âkm2 in Bangladesh50 and the average rice field was 1.1âÃâ10â2âkm2 in Taiwan32. Therefore, the integration of multi-source remote sensing data has the potential to provide more valuable information for updating and refining paddy rice maps in Asia.
We recognized that the paddy rice map identified in this study was affected by several potential factors. The first source of uncertainty was the limited availability of images from the 16-day revisit Landsat-8 satellite that detected the flooded signals of paddy rice fields. In the study area, the flooding periods usually lasted around three weeks, from mid to late June or early July. Therefore, the flooded signals of paddy rice fields could only be observed by one or two images. This shortage of image data could be remedied in the future by using multi-year Landsat-8 data and other optical sensors, such as Landsat TM/ETM+ and Sentinel-2A/B51,52. Secondly, paddy rice planting area mapping based on Landsat-8 images did not avoid the impacts of clouds and their shadows, just as the other optical sensors, MODIS, TM, AVHRR27,30,53. Figure 2(a,b) show that the original images from June 24 and July 10 are covered with about 30% cloud. In this study, Fmask was used to detect clouds and cloud shadows. Fig. 2(g,h) show that Fmask worked well, but it overestimated the clouds. Therefore, these factors would result in the underestimation of paddy rice fields. Although 30âm Landsat-8 images were used in this study, it was still a challenge to remove the influence of the mixed pixels (e.g. vegetation and water). These pixels had the same characteristics as paddy rice fields during the flooding periods. Therefore, they might be identified as paddies if the ratio of vegetation and water satisfied the extraction algorithm: LSWI â EVIââ¥â0 or LSWI âNDVIââ¥â0. The confusion of paddy fields with other land cover types caused by mixed pixels is a common problem occurring in various sensors, including Landsat TM/ETM, SPOT, MODIS32,33. Various agricultural practices would be another source of uncertainty. This research extracted flooded paddy rice fields with Landsat-8 images from June 24 and July 10, 2013. If some farmers deviated from the regular agricultural practice calendar, their paddy rice fields with earlier or later plants would not be detected from these two images, because flooding/transplanting signals would be weak or nonexistent. In addition, the rapid Land Use/Land Cover changes caused by urbanization and industrialization in this area from 2010 to 2013 was one of the reasons for the discrepancy between PaddyLandsat-8 and PaddyNLCD2010.
The results of this study have demonstrated the potential of multi-source remote sensing data (Landsat-8, MODIS and PALSAR) to map paddy rice planting areas in the wheat-rice double cropping system, using a pixel- and phenology-based algorithm. MODIS-based LSTnight and VIs dynamics make it possible to automatically select Landsat-8 images within key time windows. PALSAR data can solve the mixture of croplands and forest. Landsat-8 images provide more details about the distribution of paddy rice fields, which is useful for the extraction of fragmented ones. As Landsat-8, MODIS and PALSAR Forest/Non-forest product are available to the public, there is a potential to develop 30âm paddy rice planting area maps across the two-crop zone using this approach.
Methods
Maps of non-cropland land cover types
It is necessary to map some major non-cropland land covers, including water bodies, built-up and barren lands, forests, permanently flooded regions during the growing season (Fig. 6). In practice, these land cover types potentially affect the seasonal dynamics of vegetation indices and the accuracy of the paddy rice detection algorithm.
The workflow for mapping paddy rice planting area in a doubling agricultural system (Yangzi-Huaihe Plain) using Landsat-8 images from 2013. PALSAR 50âm FNF product was used as forest mask. Vegetation indices algorithms were used to detect non-croplands and flooded croplands. Figure created in Microsoft Visio 2010.
Water bodies have lower NDVI and EVI and higher LSWI values. Similar to the water body extraction algorithm based on MODIS29, pixels in each image meeting the condition NDVIâ<â0.1 and NDVIâ<âLSWI were extracted as water. Persistent water bodies were then composed of the pixels that were identified as water in all the good quality observations throughout the plant growing season.
Built-up and barren lands have high reflectance at visible and near infrared bands and low moisture content. Based on these physical features, a simple algorithm was put forward for built-up/barren lands, that is, LSWIâ<â0. Then, we calculated the frequency of a pixel identified as barren/built-up lands in the Landsat-8 time series (the total number of good quality observations). Permanent built-up/barren land mask was made up by pixels with a frequency of â¥90%.
Forest cover can be mapped from optical images or microwave images. Because of frequent cloud cover and 16-day revisit cycle, there is no sufficient number of Landsat images available in 2013 to distinguish croplands and forests. Previous studies showed that the cloud-free L-band SAR is the most advantageous for forest detection46. In this research project, the PALSAR-50âm Forest/Non-forest (FNF) classification map from 2010, provided by Japan Aerospace Exploration Agency (JAXA), were resampled to generate a 30âm forest mask. These PALSAR 50âm products are free to the public at the official ALOS Kyoto and Carbon Initiative website ( http://www.eorc.jaxa.jp/ALOS/en/palsar_fnf/fnf_index.htm)54. Compared to the ground truth data, these forest/non-forest products yielded 84% total accuracy on average55. This forest mask includes evergreen and deciduous forests.
Unlike seasonally flooded paddy rice fields, some lands remain flooded during the whole growing season. Therefore, it is necessary to distinguish between these two kinds of flooded lands. First, this work identified the flooded pixels for each Landsat-8 image following the rule LSWI â EVIââ¥â0 or LSWI â NDVIââ¥â0. Then, a map of the permanently flooded lands was produced from the pixels flooded in all the good-quality observations throughout the growing season.
According to the seasonal dynamics of MODIS-based vegetation indices (Fig. 5(a)) and Landsat-8 datasets (Table S1), seasonally flooded pixels were further divided into three phases: flooded from mid-April to early June, from mid-June to early July and from late July to early November. From mid-April to early June, flooded signals mainly appeared in natural wetlands, some aquaculture areas and irrigation channels. From mid-June to early July, paddy rice fields begin flooding and transplanting with significant flooded signals. From late July to early November, flooded signals disappear from paddy rice fields and they may appear in the irrigation channels of croplands.
Algorithms for identifying inundation and paddy rice fields
Paddy rice is usually planted in flooded fields. Three periods can be differentiated during paddy rice growth. In the flooding periods, the land surface of paddy rice fields is covered by water with a depth of 2-15âcm and green paddy rice plants. About 50 to 60 days later, most of the cropland surface is covered by the canopies of paddy rice plants. At the end of the growth period prior to harvesting, there is a decrease in the number of leaves and a decrease of leaf and stem moisture content.
The phenological features of main crops in the study area were investigated through MOD09A1 8-day composite vegetation indices time series (Fig. 5), obtained from the MODIS data portal at the Earth Observation and Modeling Facility (EOMF), University of Oklahoma ( http://www.eomf.ou.edu/visualization/gmap/). To get reliable vegetation index time series data, the bad-quality observations were gap-filled through a three-step gap-filling procedure56. Figure 5(a,b) show that the largest difference between paddy rice and other crops is the flooded signal during the growing season.
The phenological features of main crops were also investigated via the dynamics of three vegetation indices based on Landsat-8 time series images. Figure 5(a) shows that LSWI is larger than EVI in early July. Figure 5(a,b) indicate that despite a 16-day revisit, it is possible to detect the difference between paddy rice and corn from the Landsat-8 images. According to the seasonal dynamics of NDVI, EVI and LSWI and the finer spatial resolution of Landsat-8 images, an algorithm was proposed to extract paddy rice fields from 30âm Landsat-8 images: LSWI â NDVIââ¥â0 or LSWI â EVIââ¥â0. According to the dynamics of MODIS vegetation indices, Landsat-8 images acquired within the flooding periods were used to map paddy rice fields. In this study, the flooding periods of paddy rice fields was from mid-June to early July, composed by the flooding/transplanting period and the reviving period of paddy rice calendar.
Accuracy assessment of resultant maps
We used ground truth data (field photos (Figure S5)), Google Earth (GE) and high-resolution images to locate and digitize ROIs. Google Earth displays high-resolution images, which have been used to validate land cover classification in several studies57,58,59. However, GE images were not enough to visually interpret ROIs as it lacked images within key time windows. We also ordered multiple high-resolution images from 2012 and 2013 from the NASA Goddard Space Flight Center, including WorldView-2 (WV2), OrbView5 (OV5) and QuickBird2 (QB2). According to the reference information, we generated a series of random sampling points and interpreted them into ROIs. In total, 15,751 Landsat-8 pixels were acquired, including 7,388 paddy rice pixels (173 ROIs) and 8,363 non-paddy rice pixels (427 ROIs) (Figure S6).
The accuracy of the paddy rice map produced in this study was assessed by using the âGround Truth ROIsâ method in ENVI software. We obtained a confusion matrix between the paddy rice map and the ROI data and producerâs accuracy, userâs accuracy, overall accuracy and the Kappa coefficient. We estimated the accuracy of a paddy rice map using only good quality observations during the flooding/transplanting period. Finally, according to the cloud/cloud shadow masks from June 24 and July 10, we selected 3,610 paddy rice pixels (199 ROIs) and 3,113 non-paddy rice pixels (85 ROIs) located in good observation regions that were then used to validate the final paddy rice map.
Comparison with other available datasets of paddy rice fields
We compared our results with NLCD2010 to analyze their differences and the dynamics of paddy rice fields. The paddy rice map was compared with NLCD2010 at two scales. At the regional level, we analyzed the variations of paddy rice fields in spatial distribution and planting areas. At the pixel level, correlation analysis was used to compare these two datasets in regard to the estimation of paddy rice planting areas.
Additional Information
How to cite this article: Wang, J. et al. Mapping paddy rice planting area in wheat-rice double-cropped areas through integration of Landsat-8 OLI, MODIS and PALSAR images. Sci. Rep. 5, 10088; doi: 10.1038/srep10088 (2015).
References
Khush, G. S. What it will take to feed 5.0 billion rice consumers in 2030. Plant Mol. Biol. 59, 1â6 (2005).
Samad, M. et al. Irrigation management strategies for improving the performance of irrigated agriculture. Outlook agr. 21, 279â286 (1992).
Brisco, B. et al. Compact polarimetry assessment for rice and wetland mapping. Int. J. Remote Sens. 34, 1949â1964 (2013).
Zhang, W., Yu, Y., Huang, Y., Li, T. & Wang, P. Modeling methane emissions from irrigated rice cultivation in China from 1960 to 2050. Glob. Chang. Biol. 17, 3511â3523 (2011).
Woolston, C. Rice. Nature 514, S49âS49 (2014).
Dayton, L. Agribiotechnology Blue-sky rice. Nature 514, S52âS54 (2014).
Doraiswamy, P. C. et al. Application of MODIS derived parameters for regional crop yield assessment. Remote Sens. Environ. 97, 192â202 (2005).
Wu, W. B., Shibasaki, R., Yang, P., Tang, H. J. & Sugimoto, K. Modeling changes in paddy rice sown areas in Asia. Sustain. Sci. 5, 29â38 (2010).
Rosenzweig, C. et al. Water resources for agriculture in a changing climate: international case studies. Glob. Environ. Change 14, 345â360 (2004).
Li, C. S. et al. Modeling impacts of farming management alternatives on CO2, CH4 and N2O emissions: A case study for water management of rice agriculture of China. Glob. Biogeochem. Cycles 19, GB3010 (2005).
Matthews, E., Fung, I. & Lerner, J. Methane emission from rice cultivation: Geographic and seasonal distribution of cultivated areas and emissions. Glob. Biogeochem. Cycles 5, 3â24 (1991).
Aselmann, I. & Crutzen, P. J. Global distribution of natural freshwater wetlands and rice paddies, their net primary productivity, seasonality and possible methane emissions. J. Atmos. Chem. 8, 307â358 (1989).
Leff, B., Ramankutty, N. & Foley, J. A. Geographic distribution of major crops across the world. Glob. Biogeochem. Cycles 18, GB1009 (2004).
Monfreda, C., Ramankutty, N. & Foley, J. A. Farming the planet: 2. Geographic distribution of crop areas, yields, physiological types and net primary production in the year 2000. Glob. Biogeochem. Cycles 22, GB1022 (2008).
Huke, R. E. Rice area by type of culture: South, Southeast and East Asia. Los Baños/Laguna/Philippine: International Rice Research Institute, 1982).
Liu, Z. et al. Change analysis of rice area and production in China during the past three decades. J. GEOGR. SCI. 23, 1005â1018 (2013).
Frolking, S. et al. Combining remote sensing and ground census data to develop new maps of the distribution of rice agriculture in China. Glob. Biogeochem. Cycles 16, 1091 (2002).
Thenkabail, P. S. et al. Assessing Future Risks to Agricultural Productivity, Water Resources and Food Security: How Can Remote Sensing Help? Photogramm. Eng. Remote Sensing 78, 773â782 (2012).
Zhang, Y., Wang, C. Z., Wu, J. P., Qi, J. G. & Salas, W. A. Mapping paddy rice with multitemporal ALOS/PALSAR imagery in southeast China. Int. J. Remote Sens. 30, 6301â6315 (2009).
Shao, Y. et al. Rice monitoring and production estimation using multitemporal RADARSAT. Remote Sens. Environ. 76, 310â325 (2001).
Miyaoka, K. et al. Rice-Planted Area Mapping Using Small Sets of Multi-Temporal SAR Data. Geosci. Remote Sens. Lett., IEEE 10, 1507â1511 (2013).
Bouvet, A. & Thuy, L. T. Use of ENVISAT/ASAR wide-swath data for timely rice fields mapping in the Mekong River Delta. Remote Sens. Environ. 115, 1090â1101 (2011).
Li, P., Feng, Z., Jiang, L., Liu, Y. & Xiao, X. Changes in rice cropping systems in the Poyang Lake Region, China during 2004â2010. J. GEOGR. SCI. 22, 653â668 (2012).
Bridhikitti, A. & Overcamp, T. J. Estimation of Southeast Asian rice paddy areas with different ecosystems from moderate-resolution satellite imagery. Agric. Ecosyst. Environ. 146, 113â120 (2012).
Gumma, M. K., Thenkabail, P. S., Maunahan, A., Islam, S. & Nelson, A. Mapping seasonal rice cropland extent and area in the high cropping intensity environment of Bangladesh using MODIS 500 m data for the year 2010. ISPRS J. Photogramm. Remote Sens. 91, 98â113 (2014).
Son, N.-T., Chen, C.-F., Chen, C.-R., Duc, H.-N. & Chang, L.-Y. A Phenology-Based Classification of Time-Series MODIS Data for Rice Crop Monitoring in Mekong Delta, Vietnam. Remote Sens. 6, 135â156 (2014).
Sun, H. S., Huang, J. F., Huete, A. R., Peng, D. L. & Zhang, F. Mapping paddy rice with multi-date moderate-resolution imaging spectroradiometer (MODIS) data in China. J. ZHEJIANG UNIV-SC. A. 10, 1509â1522 (2009).
Xiao, X. et al. Mapping paddy rice agriculture in South and Southeast Asia using multi-temporal MODIS images. Remote Sens. Environ. 100, 95â113 (2006).
Xiao, X. M. et al. Mapping paddy rice agriculture in southern China using multi-temporal MODIS images. Remote Sens. Environ. 95, 480â492 (2005).
Fang, H. L., Wu, B. F., Liu, H. Y. & Huang, X. Using NOAA AVHRR and landsat TM to estimate rice area year-by-year. Int. J. Remote Sens. 19, 521â525 (1998).
Bachelet, D. Rice paddy inventory in a few provinces of China using AVHRR data. Geocarto Int. 10, 23â38 (1995).
Chen, C. F., Huang, S. W., Son, N. T. & Chang, L. Y. Mapping double-cropped irrigated rice fields in Taiwan using time-series Satellite Pour IâObservation de la Terre data. J. Appl. Remote Sens. 5, 053528 (2011).
Thi, T. H. N., De Bie, C. A. J. M., Ali, A., Smaling, E. M. A. & Chu, T. H. Mapping the irrigated rice cropping patterns of the Mekong delta, Vietnam, through hyper-temporal SPOT NDVI image analysis. Int. J. Remote Sens. 33, 415â434 (2012).
Xiao, X. et al. Observation of flooding and rice transplanting of paddy rice fields at the site to landscape scales in China using VEGETATION sensor data. Int. J. Remote Sens. 23, 3009â3022 (2002).
Akbari, M. et al. Crop and land cover classification in Iran using Landsat 7 imagery. Int. J. Remote Sens. 27, 4117â4135 (2006).
Gumma, M. K., Nelson, A., Thenkabail, P. S. & Singh, A. N. Mapping rice areas of South Asia using MODIS multitemporal data. J. Appl. Remote Sens. 5, 053547 (2011).
Panigrahy, S. & Parihar, J. S. Role of Middle Infrared Bands of Landsat Thematic Mapper in Determining the Classification Accuracy of Rice. Int. J. Remote Sens. 13, 2943â2949 (1992).
Zhong, L., Gong, P. & Biging, G. S. Efficient corn and soybean mapping with temporal extendability: A multi-year experiment using Landsat imagery. Remote Sens. Environ. 140, 1â13 (2014).
Zhong, L., Hawkins, T., Biging, G. & Gong, P. A phenology-based approach to map crop types in the San Joaquin Valley, California. Int. J. Remote Sens. 32, 7777â7804 (2011).
Shiu, Y. S., Lin, M. L., Huang, C. H. & Chu, T. H. Mapping paddy rice agriculture in a highly fragmented area using a geographic information system object-based post classification process. J. Appl. Remote Sens. 6, 063526 (2012).
Sakamoto, T., Van Phung, C., Kotera, A., Nguyen, K. D. & Yokozawa, M. Analysis of rapid expansion of inland aquaculture and triple rice-cropping areas in a coastal area of the Vietnamese Mekong Delta using MODIS time-series imagery. Landscape Urban Plann. 92, 34â46 (2009).
Peng, D. L., Huete, A. R., Huang, J. F., Wang, F. M. & Sun, H. S. Detection and estimation of mixed paddy rice cropping patterns with MODIS data. Int. J. Appl. Earth Obs. Geoinf. 13, 13â23 (2011).
Torbick, N., Salas, W. A., Hagen, S. & Xiao, X. M. Monitoring Rice Agriculture in the Sacramento Valley, USA With Multitemporal PALSAR and MODIS Imagery. IEEE J-STARS 4, 451â457 (2011).
Biradar, C. M. & Xiao, X. Quantifying the area and spatial distribution of double- and triple-cropping croplands in India with multi-temporal MODIS imagery in 2005. Int. J. Remote Sens. 32, 367â386 (2011).
Shimada, M. et al. New global forest/non-forest maps from ALOS PALSAR data (2007â2010). Remote Sens. Environ. 155, 13â31 (2014).
Dong, J. et al. Mapping deciduous rubber plantations through integration of PALSAR and multi-temporal Landsat imagery. Remote Sens. Environ. 134, 392â402 (2013).
Ou, W. Y., Guishan; LI Hengpeng; Yv, Xingxiu. Spatio-temporal Variation and Driving Forces of Landscape Patterns in the Coastal Zone of Yancheng, Jiangsu. SCIENTIA GEOGRAPHICA SINICA 24, 610â615 (2004).
Long, H., Liu, Y., Wu, X. & Dong, G. Spatio-temporal dynamic patterns of farmland and rural settlements in SuâXiâChang region: Implications for building a new countryside in coastal China. Land Use Policy 26, 322â333 (2009).
Tan, S., Heerink, N. & Qu, F. Land fragmentation and its driving forces in China. Land use policy 23, 272â285 (2006).
Rahman, S. & Rahman, M. Impact of land fragmentation and resource ownership on productivity and efficiency: The case of rice producers in Bangladesh. Land Use Policy 26, 95â103 (2009).
Drusch, M. et al. Sentinel-2: ESAâs Optical High-Resolution Mission for GMES Operational Services. Remote Sens. Environ. 120, 25â36 (2012).
Berger, M., Moreno, J., Johannessen, J. A., Levelt, P. F. & Hanssen, R. F. ESAâs sentinel missions in support of Earth system science. Remote Sens. Environ. 120, 84â90 (2012).
Motohka, T., Nasahara, K. N., Miyata, A., Mano, M. & Tsuchida, S. Evaluation of optical satellite remote sensing for rice paddy phenology in monsoon Asia using a continuous in situ dataset. Int. J. Remote Sens. 30, 4343â4357 (2009).
Dong, J., Xiao, X., Sheldon, S., Biradar, C. & Xie, G. Mapping tropical forests and rubber plantations in complex landscapes by integrating PALSAR and MODIS imagery. ISPRS J. Photogramm. Remote Sens. 74, 20â33 (2012).
Shimada, M. et al. Generation of 10m resolution PALSAR and JERS-SAR mosaic and forest/non-forest maps for forest carbon tracking. International Geoscience and Remote Sensing Symposium, Vancouver, BC, Canada, IEEE, 3510â3513 (2011).
Jin, C. et al. Phenology and gross primary production of two dominant savanna woodland ecosystems in Southern Africa. Remote Sens. Environ. 135, 189â201 (2013).
Cohen, W. B., Yang, Z. G. & Kennedy, R. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 2. TimeSync - Tools for calibration and validation. Remote Sens. Environ. 114, 2911â2924 (2010).
Montesano, P. M. et al. MODIS tree cover validation for the circumpolar taiga-tundra transition zone. Remote Sens. Environ. 113, 2130â2141 (2009).
Luedeling, E. & Buerkert, A. Typology of oases in northern Oman based on Landsat and SRTM imagery and geological survey data. Remote Sens. Environ. 112, 1181â1195 (2008).
Acknowledgements
This study was supported in part by research grants from the National Aeronautics and Space Administration Land Cover and Land Use Change program (NNX11AJ35G, NNX14AD78G), US National Science Foundation EPSCoR program (IIA-1301789), National Institutes of Health (NIH) NIAID (1R01AI101028-01A1). The high resolution images were provided by NASA for use in the NASA projects with the terms of the National Geospatial-Intelligence Agencyâs (NGA) Nextview License Agreement, NASAâs NGA Commercial Data Archive web site (cad4nasa.gsfc.nasa.gov). We thank Ms. Sarah Xiao for the English editing.
Author information
Authors and Affiliations
Contributions
X.-M.X. and J.âW. designed the research. J.W., Y.-W. Q. and Y.-T.Z. collected field data. J.W., J.-W.D, G.-L.Z., W.-L.K. and C.âJ. performed data analysis. Y.âZ. contributed to the writing of the manuscript. All authors commented on the manuscript.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Wang, J., Xiao, X., Qin, Y. et al. Mapping paddy rice planting area in wheat-rice double-cropped areas through integration of Landsat-8 OLI, MODIS and PALSAR images. Sci Rep 5, 10088 (2015). https://doi.org/10.1038/srep10088
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep10088
This article is cited by
-
Effects of no-till and rice varieties on nitrogen and phosphorus balance in rice fields
Environmental Sciences Europe (2020)
-
Fingerprint of rice paddies in spatialâtemporal dynamics of atmospheric methane concentration in monsoon Asia
Nature Communications (2020)
-
Special issue: Urban green infrastructure and the ecological functions
Landscape and Ecological Engineering (2019)
-
No-tillage effects on grain yield, N use efficiency, and nutrient runoff losses in paddy fields
Environmental Science and Pollution Research (2016)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.