Remote Sensing Mapping Flooded Rice Paddies Using Time Series of Modis Imagery in the Krishna River Basin, India

Rice is one of the major crops cultivated predominantly in flooded paddies, thus a large amount of water is consumed during its growing season. Accurate paddy rice maps are therefore important inputs for improved estimates of actual evapotranspiration in the agricultural landscape. The main objective of this study was to obtain flooded paddy rice maps using multi-temporal images of Moderate Resolution Imaging Spectroradiometer (MODIS) in the Krishna River Basin, India. First, ground-based spectral samples collected by a field spectroradiometer, CROPSCAN, were used to demonstrate unique contrasts between the Normalized Difference Vegetation Index (NDVI) and the Land Surface Water Index (LSWI) observed during the transplanting season of rice. The contrast between Enhanced Vegetation Index (EVI) and Land Surface Water Index (LSWI) from MODIS time series data was then used to generate classification decision rules to map flooded rice paddies, for the transplanting seasons of Kharif and Rabi rice crops in the Krishna River Basin. Consistent with ground spectral observations, the relationship of the MODIS EVI vs. LSWI of paddy rice fields showed distinct features from other crops during the transplanting seasons. The MODIS-derived maps were validated against extensive reference data collected from multiple land use field surveys. The accuracy of the paddy rice maps, when determined using field plot data, was approximately 78%. The MODIS-derived rice crop areas were also OPEN ACCESS 8859 compared with the areas reported by Department of Agriculture (DOA), Government of India (Government Statistics). The estimated root mean square difference (RMSD) of rice area estimated using MODIS and those reported by the Department of Agriculture over 10 districts varied between 3.4% and 6.6% during 10 years of our study period. Some of the major factors responsible for this difference include high noise of the MODIS images during the prolonged monsoon seasons (typically June–October) and the coarse spatial resolution (500 m) of the MODIS images compared to the small crop fields in the basin. However, this study demonstrates, based on multi-year analysis, that MODIS images can still provide robust and consistent flooded paddy rice extent and areas over a highly heterogeneous large river basin.


Introduction
Evapotranspiration (ET) from irrigated crop areas comprises the biggest anthropogenic consumption of fresh water in the world.Irrigation is known to consume about 70-80 percent of all water used by humans [1].It is therefore important to estimate current irrigation water use accurately, not only for water resources assessment but also for the prediction of future water demand.Among all the irrigated crops, flooded paddy rice has drawn particular interest due to the large population living on rice consumption and its unique cultivation method, especially in Asia.Rice feeds more than half of the entire population of the world [1,2] and the cultivation of paddy rice on flooded soils requires large quantities of fresh water, which has significant implications to water security and terrestrial water balance.Moreover, waterlogged soils are one of the largest sources of methane gas emissions [3][4][5].Mapping rice paddies is thus a critical component for developing environmentally sustainable water management [6,7].
Currently available land-cover and land-use maps based on satellite sensors, such as global land cover maps based on Landsat and Moderate Resolution Imaging Spectroradiometer (MODIS) produced by the Global Land Cover Facility (GLCF, http://www.landcover.org),do not provide detailed sub-classes within the class of agricultural fields.However, the unique flooded farming method of rice can be utilized to discriminate the paddy rice fields from other agricultural fields.Particularly, during the transplanting season of rice, a significant fraction of the rice paddies consists of ponded surface water, resulting in strong water reflectance.Even though detecting ponded surfaces using only visible (red) and near infrared (NIR) bands that are commonly used to classify vegetation types is not straightforward, it can be aided by using the Land Surface Water Index (LSWI) [8].The LSWI is calculated using spectral signals in shortwave infrared (SWIR) and near infrared (NIR) ranges to detect water at the soil surface [9] and the vegetation water content [10,11].
Previous works [8,[12][13][14] have compared NDVI or the Enhanced Vegetation Index (EVI) with LSWI during the transplant season to map paddy rice fields in China, South Asia and Southeast Asia.Xiao et al. [8] analyzed the multi-temporal SPOT-4 satellite imagery (Vegetation 10-day composite images, 1-km resolution) collected over Jiangsu, China to map flooded paddy rice fields, where pixels with LSWI greater than NDVI were classified as paddy rice.Xiao et al. [12,13] used multi-temporal MODIS images (8-day composites, MOD09A1, 500-m resolution) to extend the approach to larger regions in South China, South Asia and Southeast Asia.They used a relaxed set of criteria, LSWI + 0.05 > EVI or LSWI + 0.05 > NDVI, to identify flooded paddy rice fields.The mapped paddy rice fields were overall in good agreement with census statistics and regional land cover maps derived from the Landsat ETM+, although some discrepancies in land classes existed.The errors were attributed to cloud contamination of optical images, topographic effects, and low spatio-temporal resolution of the images.Sun et al. [14] provided more detailed multi-year LSWI vs. EVI relationships in paddy rice fields depending on the planting season of the crop.Recently, Gumma et al. [15] mapped rice area in South Asia using spectral matching techniques, which were fully based on MODIS NDVI time-series data.
Previous studies reported that the coarse spatial and temporal resolution of satellite observations used is an important source of error in identifying flooded rice fields.The ground coverage of rice increases steeply within a few weeks from transplanting and NDVI quickly saturates to its near-maximum values [16].Moreover, due to the typically small size of paddies in South and Southeast Asian countries, most 500-m MODIS pixels contain mixed flooded rice, non-flooded crops and other vegetation.The MODIS SWIR bands required to calculate LSWI are available for 500-m spatial resolution.Consequently, the sub-pixel-scale heterogeneity can seriously degrade the performance of the NDVI (or EVI)/LSWI-based paddy-rice-mapping.In order to examine NDVI-vs.-LSWIbehavior free from the effect of the sub-pixel-scale heterogeneity, it is therefore necessary to analyze field spectral samples collected from various crop covers and growth stages.Ground-based spectral samples can play a critical role [17,18] in developing and calibrating an algorithm to map rice paddies.Note that the previous works reviewed (e.g., [12][13][14]) developed algorithms directly on the assumed relationship between satellite-based NDVI (and EVI) and LSWI without presenting their actual relationship over homogeneous fields.
In this study, spectral samples of flooded paddy rice at various growth stages, along with those of non-paddy crops, were collected from the Musi sub-basin of the Krishna River Basin in India (Figure 1) using an MSR 16 CROPSCAN spectroradiometer to demonstrate behavioral features of NDVI and LSWI.After confirming a distinctive NDVI vs. LSWI relationship for paddy rice and non-paddy crops, flooded rice paddies are mapped using the time series of three indices-NDVI, LSWI and EVI-derived from the 8-day composite of MODIS data (MOD09A1).Consequently, decision rules for Kharif rice and Rabi rice were developed from a subset of the observed data using LSWI and EVI, the rules were then applied on MODIS data over the whole study area, and the spatial distribution of paddy rice maps for the Kharif and Rabi seasons finally produced over the Basin.

Study Area
The Krishna River Basin (Figure 1) is the fourth largest river basin in India in terms of drainage area, located in the southern part of India between longitudes 73°15′E and 81°15′E, and latitudes 13°05′N and 19°20′N.The total geographical area is 265,752 km 2 , covering the three states of Karnataka, Maharashtra and Andhra Pradesh.The River Krishna originates in the Western Ghats Mountains, flows east across the Deccan Plateau, and discharges into the Bay of Bengal.It has three main tributaries that drain from the northwest, west, and southwest (Figure 1).The climate is dominantly semi-arid, with some dry, sub-humid areas in the eastern delta and humid areas in the Western Ghats.Annual precipitation is 800 mm on average.Most of the rainfall occurs during the Indian monsoon from June to October [19,20].Crops are grown in three seasons, during the monsoon (June to mid-December), called Kharif, in the post-monsoon dry season (mid-December to March) called Rabi, and dry summer period (April and May) [19,20].Compared to Kharif and Rabi, the short dry summer period, called Zaid, comprises a very small fraction of the total rice production thus will not be considered for the land cover classification of this work.Major irrigated crops include single/double cropping rice, single cropping sugarcane, chili, cotton, pulses, fodder grass, and some light irrigated corn, sorghum and sunflower.Rainfed crops include grains (sorghum, millet), pulses (red and green gram, chickpea), and oilseeds (sunflower, groundnut).In the large irrigation command areas, such as Almatti, Tungabhadra, Ujjani-Bhima, Jurala, Nagarjuna Sagar, and Prakasam (Krishna Delta), crop fields are irrigated during the Kharif and Rabi seasons.Major canal irrigation schemes occur along each of the major reservoirs in the upper basin, and along the main drainage in the lower basin and in the delta (Figure 1), which are predominant rice growing areas.Minor irrigated systems include small tanks, small riparian lift schemes, and ground water irrigation.Ground water sources include dug wells, shallow tube wells, and deep tube wells.We conducted intensive land use surveys at different crop growth stages in the Nagajuna Sagar irrigation command (Figure 1), which is located in the eastern part of the basin.Ground-based spectral samples were collected from the Musi sub-basin located in the northeastern part of the basin.

Field Spectral Samples
In order to demonstrate and confirm the unique NDVI vs. LSWI relationship of flooded paddy rice, which were utilized to map rice paddies in previous works, a spectral library of paddy rice at different growth stages were taken from ground-based spectral samples using a portable multispectral radiometer (model MSR16) manufactured by CROPSCAN, Inc.The MSR16 used in this work collects 12 bands, six bands from sky and six bands from ground targets, providing on-site calibrated surface reflectance at six wavelengths ranging from 530 nm to 1640 nm.It is assumed that the irradiance flux density incident on the sky-facing radiometers is identical to the flux density incident to the target surface.Specifications of the spectral data are summarized in Table 1.An intensive ground sampling was conducted as a part of the Australia-India Land Surface Parameterisation Experiment for Remote Sensing (AILSPEX-RS) during February 2011 (AILSPEX-11) and September 2012 (AILSPEX-12) in the Musi sub-basin of the Krishna River Basin in India [17].The first experiment period captured the transplanting season of Rabi (winter, dry cropping season) in the region and the second experiment captured the peak growing season of Kharif.A comprehensive database was constructed for the field spectral sample survey.Information collected for each site includes: 1. Location information (GPS position, location name, date of collection).2. Land use (crop type).3. Vegetation height.4. 12-band reflectance measurements by CROPSCAN.5. Land surface temperature measured by the thermal infrared scanner.6. Soil temperature at the depths of 1 cm, 5 cm and 10 cm. 7. Soil moisture content in the top 5 cm by the theta probe soil moisture sensor (model ML2 by Delta-T Devices Ltd.).8. Digital photographs.
The data collected from this experiment was used to identify unique features of rice paddies and to justify the rice paddy identification criteria from MODIS imagery.LSWI and NDVI were produced over the paddy rice fields and surrounding non-flooded crop fields, and analyzed with varying crop height.

MODIS Surface Reflectance Data
The MODIS 8-day composite product from Terra (MOD09A1) was used in this work.From the seven bands of MOD09A1 product at 500 m resolution (Table 2), four bands (blue, red, NIR and SWIR) captured the seasonal variations in vegetation vigor, soil and vegetation moisture and surface water that characterize the key stages of rice cultivation [13,21].The reflectance data have undergone several pre-processing steps, including algorithms for atmospheric correction.The 8-day composite MODIS product with 500-m spatial resolution (MOD09A1) was downloaded from http://modis-land.gsfc.nasa.govfor the period of June 2000 to May 2010.A total of 46 composites were available each year.For each 8-day composite, the following indices were calculated using surface reflectance values from the blue, red, near infrared (NIR) and shortwave infrared (SWIR) bands (see Table 2) NDVI, EVI [21], and LSWI [12,13].

Field Survey of Land Use
In order to obtain characteristic EVI vs. LSWI for flooded paddy rice and non-paddy crops, and also to validate the paddy rice map of this study, we have used information collected from 91 sampling sites in the northeastern part of the Krishna River Basin were used (see Figure 1, Land use survey locations) during Kharif and Rabi in 2006-2007.Ground-based satellite validation can be complicated by the subpixel-scale surface heterogeneity.In order to reduce the impact of the heterogeneity, field survey sites were chosen over relatively homogeneous cropping fields.The sampling sites include various crop fields: irrigated rice, cotton, chili, sugarcane, turmeric, grams, etc.All of the sites are located within a major surface water irrigated area called the Nagarjuna Sagar irrigation command.A total of eight field surveys were conducted at different crop growth stages.A comprehensive database constructed from the field survey includes: 1. Location information (GPS position, location name, date of collection).2. Land use/land cover type (class name).3. Fraction of individual land cover types such as cropped canopy and no canopy area (water, fallow lands and weeds) within each crop field.4. Crop types, cropping pattern and cropping calendar (i.e., Kharif, Rabi and summer seasons).5. Agricultural intensification, sources of water and presence of irrigated, rain-fed and supplemental irrigation.6. Characteristics of crops such as plant height, soil texture, and density of plants per square meter; 7. Digital photographs.Some of the photos showing growth stages of two paddy rice fields (Sites 25 and 79) in the Nagarjuna Sagar irrigation command area are shown in Figure 2. The first photo in Figure 2a,b shows the rice paddies during the transplanting season, where a large fraction of ponded field was exposed when compared to other growth stages.The fraction of ponded surface is gradually decreasing with the increase in vegetation cover.

Agricultural Statistical Data
District level agricultural census data for the current study area (in 2000-2010) were obtained from the Directorate of Economics and Statistics (DES), Department of Agriculture (DOA), Government of India (http://apy.dacnet.nic.in).The database contains district level statistics on cropping area, production and productivity of each crop.The district-wise statistics of paddy rice area in each year, along with the field survey data, were used to validate the paddy rice maps produced from the satellite imagery in this work.

Growth Stages of Rice
The cropping period of rice varies between 4 and 6 months in the Krishna River Basin, depending on the rice variety.In general, paddy fields are ploughed and flooded before rice transplanting.Rice grows rapidly after transplantation with the temporal dynamics of paddy rice fields characterized by three main periods: flooding/transplanting period, cropping period (vegetative growth, reproductive and ripening stages) and the fallow period after harvest.During the rice transplant period, the land surface is a mixture of surface water and green plants with water depth ranging from 5 to 20 cm.About 50 to 60 days after transplanting, the rice canopy grows and covers most of the cultivated surface area.From the start of end of the growth period, leaf and stem moisture content start decreasing.The NDVI and EVI of paddy fields decrease during the transplanting period and then increase for the rest of the growing season, whereas LSWI steeply increases during the transplanting season.According to time-series data of the spectral reflectance of paddy fields [22], the NDVI and EVI reaches a maximum around the heading date as the rice crop changes its growth phase from vegetative growth to reproductive growth on reaching the heading date, and the leaves begin to wither.NDVI and EVI then decrease abruptly because of harvesting.

Identification of Flooded Paddy Rice
Rice fields are filled with water at the beginning of the transplanting season and remain ponded until approximately 1-2 weeks before harvest.However, due to the large fraction of water surface during the transplanting season, reflectance in the short-wavelength infrared (SWIR) range, which is sensitive to water, can be used to distinguish the flooded rice paddies from other crops.Thus, in order to detect the transplanting date, which follows flooding in paddy rice fields, using the time profile of EVI or NDVI and LSWI, existence of the period when a temporary inversion of NDVI/EVI and where LSWI becomes higher than EVI or NDVI values or approaches closer, can be used.
Due to the locally and regionally varying irrigation schedules in the basin, the transplanting occurs over a period of nearly 60 days from July through September for Kharif and from December through February for Rabi rice.Consequently, at least six to eight MODIS 8-day composite images are needed for each season to discriminate rice paddies, which are based on our field surveys, ground based spectral data collected from field experiments (see Section 4.1 for more details) and agriculture statistical information collected from different parts of the river basin.
Identification of changes in the mixture of surface water and vegetation canopies in paddy rice fields over time requires spectral bands or vegetation indices that are sensitive to both water and vegetation.NDVI is closely correlated with the Leaf Area Index (LAI) of paddy rice fields [11] but is sensitive to aerosols.However, EVI directly adjusts the reflectance in the red band as a function of the reflectance in the blue band, accounting for residual atmospheric contamination and variable soil and canopy background reflectance [21].Consequently, EVI is known to be much less sensitive to aerosols than NDVI and so in this work, EVI is chosen instead of NDVI for the vegetation change detection.The SWIR band is sensitive to leaf water content and soil moisture and is used to develop improved water sensitive indices such as LSWI, which combines NIR with SWIR bands and has the capability of retrieving canopy water content [23,24].The strong light absorption by liquid water in the SWIR (Band 6 of MODIS) range makes LSWI sensitive to the total amount of liquid water in vegetation and its soil background [9].The Band 5 of MODIS (1230-1250 nm) is also sensitive to soil water content, but it is not used to detect LSWI because of weaker sensitivity than Band 6.Although the 8-day composite surface reflectance products of MODIS have been pre-processed to reduce the impacts of clouds, shadows, and aerosols, there remain residual noises due to atmospheric effect and heavy clouds in some regions.Algorithms used for removing noise in time-series methods include moving median or average and Fourier-based algorithms [25,26].In this work, contaminated pixels were identified by abrupt changes in NDVI with time and were substituted by one of the values in the temporally adjacent images based on validity or by the average of the adjacent images when both of them are valid.If the pixels in both the adjacent images are contaminated, then those pixels are discarded.
Due to the prolonged heavy monsoon rains in the study area, temporary flooding can occur during the Kharif season.Xiao et al. [12,13] used the conditions LSWI + 0.05 ≥ EVI or LSWI + 0.05 ≥ NDVI for identifying rice paddy fields on the basis of transplanting period, irrespective of the acquired date of the images.An auxiliary condition was also used: the EVI value reaches half of the maximum EVI value within five 8-day composites following the flooding and transplanting period [12,13].Moreover, the non-rice area was further removed using the digital elevation model (DEM) data and masks of cloud, permanent water, and permanent/evergreen vegetation.This method is expected to avoid the need of identifying various transplanting schedules over a region.However, it can overestimate the rice fields due to the other irrigated crops or temporary flooding due to rainfall.Moreover, Kharif rice and Rabi rice cannot be separated if they are planted in the same fields in a hydrological year.
To resolve these problems, MODIS data were integrated with prior information on the rice growth calendar obtained from field observations in different zones of the basin.The data were obtained from the Department of Agriculture (DOA), the Command Area Development Authority (CADA) and the Reservoir Monitoring Authorities.The planting calendar is relatively similar to the rice growth calendar with small inter-annual variations, thus all potential dates of the flooding and transplanting period in different zones are selected based on the calendars of 2000-2010.As the calendars of the plants in other croplands might be different from the rice calendar, so the noise in the images beyond the transplanting period (e.g., strong water signal) could be eliminated.After harvesting the Kharif rice, the rice fields are kept fallow for approximately 2-4 weeks before the Rabi season starts, or they are cultivated with shortterm dry crops.Based on the cropping patterns of the study region, the available MODIS 8-day composite images in a hydrological year were partitioned into two seasons to generate Kharif and Rabi rice maps, consequently to produce separate categories for single crop (Kharif or Rabi) and double cropped rice (Kharif and Rabi) fields.
The reflectance over these areas is the mixture of paddy rice fields (composed of water, soil, seedlings, etc.) and the background (roads, weeds, other crops, etc.) in the rice-transplanting period.The characteristics of high water and low coverage of vegetation can be detected by comparing LSWI and EVI, since they feature relatively high LSWI and low EVI in this period.Areas of high EVI in this period are regarded as non-rice vegetation fields (trees, shrubs, grass, etc.).In order to obtain characteristic EVI vs. LSWI for flooded paddy rice and non-paddy crops, we examined time series of 8-day EVI and LSWI over the 20 rice and 20 non-rice sites in different regions of our study area including double crop rice (Figure 3) and single crop rice (Figure 4) grown areas.Average EVI and LSWI at each test site during the flooding and transplanting period (for Kharif and Rabi) were calculated.The decision rules for mapping rice are LSWI > 0.12, EVI < 0.27, and LSWI > (EVI − 0.05) for Kharif rice, and LSWI > 0.10, EVI < 0.29, and (LSWI + 0.12) > EVI for Rabi rice (with increased the relaxation from 0.05 to 0.12, see Section 4.2).An additional condition was applied to the initial rice-paddy pixels to remove permanent water bodies: the pixels with EVI ≤ 0.35 in the 6th~11th 8-day MODIS EVI composites following the transplanting period were flagged as permanent water bodies.These conditions were used to generate rice paddy maps for both Kharif and Rabi crops for the entire study region.An overview of the overall approach for mapping rice paddies for Kharif and Rabi seasons is shown in Figure 5 in the form of a flow diagram.

Spectral Characteristics from Ground Reflectance
The main aim of collecting ground reflectance data was to compare the distinct features of NDVI vs. LSWI between the paddy rice and other non-flooded crop fields.Even though a steep increase of LSWI paired with a decrease of NDVI during the transplant season was inferred from their general behavior and used to map rice paddies in previous works, it has not been shown explicitly using ground samples over the systematically varied homogeneous fields.The spectral responses from two categories of fields are summarized in Figure 6.The results show that LSWI of non-rice crops increases linearly with NDVI in general.The strong linear relationship between the indices is due to the dominant impact of vegetation water content and the non-flooded nature on the observed LSWI in the non-rice fields.In contrast with the other non-flooded crops, LSWI of the paddy rice fields showed decreasing trend with NDVI.Even though LSWI is scattered widely around the linear fit, LSWI decreased linearly with NDVI until the NDVI of rice reached approximately 0.75.The two linear models intersected at the point where the NDVI value is 0.85.It also appears that LSWI of paddy rice joins the increasing trend with NDVI beyond the threshold NDVI value.At full canopy cover the influence of ponded water on reflectivity becomes insignificant resulting in both LSWI and NDVI showing high values.The decreasing trend of LSWI is This contrast between LSWI and NDVI during initial growth stages in ponded paddy fields provides contrasting supplementary information, leading to an improved ability in separation of flooded rice from other crops.For paddy rice fields, the decreasing LSWI with increasing NDVI is highly correlated with the height of the rice (Figure 6).The results from the analysis in Figure 6 also shows that examining the trends between NDVI and LSWI are more informative for discriminating paddy rice fields than simply comparing the magnitude of the normalized indices, as was adopted by the previous works.The LSWI values are greater than NDVI only up until the height of the rice reaches approximately 20 cm (Figure 6).However, the spectral samples from the paddy rice exhibit clear distinction from the values of other fields up to NDVI = 0.6, which is equivalent to 35 cm of crop height.The observed difference in the crop height ranges can be used to discriminate paddy rice fields, but only when the rice crop is less than 35 cm tall.With the given limits on temporal scales of satellite images when using the 8-day composite imagery of MODIS for mapping, a larger time window for mapping is gained by explicit LSWI-NDVI analysis when it is assumed that the rice grows linearly with time just after transplanting. Figure 6 shows that plots of paddy rice can be better separated from the other fields while a number of LSWI and NDVI values are mixed when the height of rice is greater than 10 cm.

Decision Rules for Mapping Kharif and Rabi Rice Paddy Using MODIS Indices
The results from the analysis in Figures 7 and 8 show the thresholds of MODIS EVI and MODIS LSWI for the paddy rice fields during transplantation period for Kharif and Rabi seasons.Based on these thresholds, we obtained the final set of decision rules for mapping Kharif and Rabi rice using MODIS 500 m datasets.The decision rules obtained for mapping rice are LSWI > 0.12, EVI < 0.27, and LSWI > (EVI − 0.05) (Figure 7) for Kharif rice, and LSWI > 0.10, EVI < 0.29, and (LSWI + 0.12) > EVI for Rabi rice (with increased the relaxation from 0.05 to 0.12, Figure 8).These conditions were then used to generate rice paddy maps for both Kharif and Rabi seasons using MODIS 500 m datasets.Also an additional condition was applied to the initial rice-paddy pixels to remove permanent water bodies: the pixels with EVI ≤ 0.35 in the 6th~11th 8-day MODIS EVI composites following the transplanting period, were flagged as permanent water bodies.

Spatial Distribution of MODIS-Derived Rice Paddies
The maps of Kharif and Rabi rice paddies were derived from MODIS data using the algorithm explained above.The classification results of 2006-2007 Kharif and Rabi rice fields were shown in Figure 9.The map shows four classes.Blue represents only Kharif-rice, red represents only Rabi-rice and green represents double crop rice (i.e., both Kharif and Rabi rice).The maps were validated using two types of data: field survey data and national statistics of rice cropping areas.

Vegetation Phenology of Various Rice Pixels from Different Zones/Regions in the Basin
Rice crop phenology was studied using NDVI time-series plots from observation points selected from land use survey data (Figures 10 and 11).These NDVI time series profiles provided information on (a) Cropping intensities (e.g., double or single crop); (b) Crop calendar (i.e., when a crop begins and when it is harvested); and (c) Crop health and vigor (indicated by magnitude of NDVI).
Each rice class has a distinctly different phenology depicted by the NDVI magnitude and/or seasonality (Figures 10 and 11).The NDVI time-series also allows the separation of short period rice and long period rice based on factors such as when a crop calendar begins, and the magnitude of NDVI.

Accuracy Assessment from Field Observations
A qualitative assessment was performed to check whether a known rice area was classified as rice or non-rice.This process was done using 91 independent observation sites (see Figure 1) where observations were conducted every 16 days during the Kharif season in 2006-2007.The error matrix (see Table 3) was constructed based on the theoretical description given by Jensen [27].The columns of the error matrix contain the field-plot data points and the rows represent the results of the classified rice maps [28].The error matrix is a multidimensional table in which the cells contain changes from one class to another [28,29].Each field point is categorized as either rice or non-rice.Among the 91 ground sampling points, 44 samples were collected from rice fields in different zones with different crop calendar.The classification results show that 29 of 44 rice-sampling points (66% accuracy) were identified as rice (Table 3) and the remaining sample were identified as non-rice.There are five locations, which are identified as rice but the corresponding ground observation shows that they are from other crops such as irrigated chili, cotton, turmeric and sugarcane crops.These crops are also irrigated by surface water, which can cause high water signals ultimately making it fall into the rice category.The producer accuracy obtained for rice crop is 66% and for non-rice crop is 89%.The user accuracy for rice crop and non-rice crop are 85% and 74%, respectively.The overall accuracy for this classification is 78%.These errors are mainly the result of the MODIS 500 m × 500 m pixel size (22 hectares per pixel) and land parcels are usually relatively smaller than one pixel resulting in pixel heterogeneity of fields in the area.Also, the classification is binary, i.e., either rice or non-rice.Hence if we have a mixed class category such as rice mixed with other crops then accuracy would be expected to improve.

Rice Area Fractions (RAFs)
Each rice pixel contains several different land cover classes in the study site; however, all the rice samples from the field survey, collected from different zones and used to calibrate the thresholds of our classification, have more than 80% rice area fraction.So, all the classified rice pixels are assumed to have relatively high rice fractions.Although, if the rice area fraction is not included, the resulting classified area will have overestimated rice cover.Hence, in order to reduce the potential bias caused by the rice fraction, RAF of individual pixels or its average value needs to be considered when calculating actual rice areas.
The MOD09A1 pixels have a size of approximately 22 hectares, which is larger than many individual agricultural fields of the study area.Thus, many pixels contain more than one land cover classes.Any rice classification derived from such data will provide only the full pixel area (FPA), whereas the actual rice area per pixel can be obtained only by computing sub-pixel area (SPA) [19,[30][31][32], defined as: where SPA is the sub-pixel area of rice pixel, FPA the full-pixel area of the MODIS pixel and RAF the rice area fraction.The rice area fraction, RAF = 0.8472, is derived from the 44 field survey sites distributed over the Nagarjuna Sagar canal command area (Figure 1), by measuring the areas of rice within each sample plot at each site.It is assumed that the areal composition of rice paddies over the field survey areas represents the typical sub-pixel-scale rice paddy fraction of the rice-dominant MODIS pixels for the entire study region.RAF is applied uniformly to the classified pixels to estimate actual rice area for MODIS derived rice paddy pixels for each season.The SPA of each location is computed by multiplying the FPA of that class with RAF of the class.Later, the SPAs of all samples are summed to obtain the rice area fraction from all the zones.

Accuracy Assessment from Agriculture Statistics
Since spatially distributed maps of Kharif or Rabi rice paddies for ground truth are not available, the agricultural statistics data from the Department of Agriculture (http://apy.dacnet.nic.in) are used to evaluate the MODIS derived rice paddy maps.To test the areal accuracy of the MODIS derived rice cropped areas, they were compared with district level agriculture statistics during the period 2000-2001 to 2009-2010.Ten districts (10 spatial units) were chosen, which are completely (not partly covered districts) covered within the extent of the study area with significant rice cropped area.District level rice cropped area were compared with MODIS derived crop area and is shown in the scatter plots (Figures 12 and 13).It is obvious that there is some discrepancy between MODIS derived rice cropped area and the agriculture statistics.Therefore, further analysis of the results by statistical hypothesis significance t-test was undertaken to show whether the discrepancies can be acceptable statistically.The two groups were tested by paired-sample t-test at the level of 0.05.The results are in Table 4 showing that the MODIS derived results and the district level agriculture statistics have significant correlations except for the years 2002-2003 and 2003-2004.MODIS based estimates of rice areas during these two years are lower than the agriculture census data, which might be attributed to the over reporting of crop areas in the agriculture census data.These two years were relatively dry years when compared to others [33,34].The total differences observed between district statistics and MODIS derived rice cropped areas from the selected 10 districts ranged from 0.2% to 24% during 2000-2001 to 2009-2010 (10 consecutive hydrological years).
The estimated RMSD of rice cropped area of the 10-districts ranged from 47,000 ha to 74,500 ha during 2000-2010, which is equivalent to 3.4% to 6.6% (Table 5) of the rice cropped area according to the district wise rice statistics.Discrepancy between the MODIS-based rice area estimates and agricultural census estimates in some administrative units can be attributed to: (1) the methods used in aggregation of rice cropped area per year from the small farm holdings and inconsistent approach and estimation methods from state to state; and (2) limitations of the 500-m resolution MODIS-based algorithm in identifying small patches of agricultural field sizes.For example, relevant to the potential issues with the agricultural census estimates (#1 above), errors in the census data of Raichur might have contributed to the particularly large discrepancy between the MODIS-derived and the census estimates (Figure 13).When the results of Raichur are excluded from evaluation, R 2 of 0.745 increases to 0.841.The poorest performance of 2009-2010 results, R 2 = 0.565, is also associated with the limited accuracy in Raichur.When the Raichur results are removed from the summary evaluation, R 2 in 2009-2010 improves significantly from 0.565 to 0.803.Due to the large-scale irrigated rice paddies of Raichur (large double-cropping fields near the center of the Krishna River Basin in Figure 9), it is unlikely that MODIS-based estimates are more susceptible to errors in the region.Failure to account for double cropping fields potentially could have led to the large underestimation of the rice cropping areas in Raichur.

Summary and Conclusions
This study demonstrates unique contrasts between the NDVI and the LSWI observed during the transplant season of rice based on ground-based spectral samples surveyed by a field spectroradiometer.Also, this study identified a new set of decision rules to map Kharif rice and Rabi rice based on land use survey.Finally, mapped flooded paddy rice for Kharif and Rabi seasons by analyzing the time series of MODIS (MOD09A1) data in the Krishna River Basin, India.Accuracy was determined by correlating MODIS-derived classification using field plot data and MODIS-derived rice areas with subnational statics obtained from the Ministry of Agriculture, Government of India.Overall a 78% of accuracy was obtained for the rice classification using field plot data and 3.4% to 6.6% of error (RMSD) observed in the MODIS derived rice area with sub national statistics over ten years.
There exist some discrepancies in MODIS derived results.We have tried to analyze the factors that might cause errors in the results.We found that clouds have a significant impact on results because all cloud contaminated pixels in the flooding and transplantation period were eliminated in the rice identifying algorithms due to their unreliability, so cloud contamination is a major source that can cause error of omission.The cloud contamination problem impacts classification, especially during Kharif season as the flooding and transplantation period of Kharif rice occurs along with monsoons.Therefore alternative datasets are needed in the regions with frequent and heavy clouds during the flooding and transplantation period.The MODIS derived results in the study have some errors due to the factors like mixed pixels caused by the coarse resolution (relative to the land parcel sizes).However, they are still useful for obtaining spatial distribution maps over a large scale and they might provide valuable information for studies such as estimation of improved evapotranspiration by accounting the open water evaporation from flooded irrigation.

Figure 1 .
Figure 1.Ground survey locations of paddy rice fields used for rice algorithm where land use survey and the CROPSCAN-based spectral sampling were conducted in the Krishna River Basin.

Figure 2 .
Figure 2. (a) Photos of the paddy rice site #25 (transplantation date of 12 August) at different growth stages in 2006-2007 Kharif season.(b) Photos of the paddy rice site #79 (transplantation date of 3 September) at different growth stages in 2006-2007 Kharif season.

Figure 3 .
Figure 3.Time series of LSWI and EVI of observed rice fields in the double crop rice zones.The vertical lines represent transplantation season, during this period we can observe inversion of indices for flooded paddy rice pixels.

Figure 4 .Figure 5 .
Figure 4. Time series of LSWI and EVI of observed rice fields in the single crop rice zones.The vertical lines represent transplantation season, during this period we can observe inversion of indices for flooded paddy rice pixels.

Final
12 , EVI<0.27 and LSWI +0.05>EVI during transplantation and Average EVI>0.35 in the 6 th to 11 th 8-day composites after transplantation Rabi rice: if LSWI>0.10 , EVI<0.29 and LSWI +0.12>EVI during transplantation and Average EVI>0.35 in the 6 th to 11 th 8-day composites after transplantationmainly due to that the fact that the spectral signal coming directly from the ponded surface in the back ground decreases with the growth of rice.During initial growth stages (canopy cover < 100%), the flooded paddy fields absorb both NIR and red heavily, resulting in low NDVI values whereas in LSWI the SWIR band absorption is relatively low (compared to NIR or red), resulting in high LSWI values.

Figure 7 .
Figure 7. Characteristics of paddy rice fields reflected by EVI and LSWI in the transplanting period according to ground observation sites during 2006-2007 Kharif season.The red line is the threshold for potential rice in the transplantation season based on our analysis.

Figure 8 .
Figure 8. Characteristics of paddy rice fields reflected by EVI and LSWI in the transplanting period according to ground observation sites during 2006-2007 Rabi season.The red line is the threshold for potential rice in the transplantation season of Rabi, whereas the green line was threshold used for Kharif.

Figure 9 .
Figure 9. Spatial distribution of Kharif and Rabi paddy rice areas derived using MODIS 500-m datasets in the Krishna River Basin during 2006-2007.Boundaries of the districts chosen for statistical evaluation (used in Section 4.7) are overlaid.

Figure 10 .
Figure 10.NDVI time series of observed rice fields, which are classified as rice pixels in the double crop rice zone.

Figure 11 . 22 -
Figure 11.NDVI time series of observed rice fields, which are classified as rice pixels in the single crop rice zone.

Figure 12 .
Figure 12.The district level comparison of rice cropped areas between agriculture census dataset and MODIS-500 m derived dataset (year wise 2000-2001 to 2009-2010).Numbers of samples are 10 districts per year and number of years are 10.Areas shown on X-and Y-axes are in hectares.

Figure 13 .
Figure 13.Correlations between rice areas derived using MODIS-500 m dataset and district wise agriculture statistics (each series represents one district from 2000-2001 to 2009-2010).Numbers of samples are 10 districts per year and number of years are 10.Areas shown on X-and Y-axes are in hectares.

Table 2 .
Specifics of the MODIS-Terra bands (MOD09A1) used in this study1.
MODIS Bands

2 Band Width (nm) Centre Wavelength (nm) Sub Division Potential Application 3
Note:1Out of the 36 MODIS bands, the first seven bands are specially processed for Land studies; 2 MODIS bands are rearranged to follow the electromagnetic spectrum (e.g., blue band 3 followed by red band 1);3Source: http://modis-land.gsfc.nasa.gov.

Table 3 .
Accuracy assessment using field surveyed data.

Table 4 .
Significance analyses between the MODIS-derived results and agriculture statistics on district scale by paired-sample t-test at the level of 0.05.

Table 5 .
Correlation between districts wise MODIS derived rice cropped areas and agriculture statistics for the selected 10 spatial units.