TROPOMI reveals dry-season increase of solar-induced chlorophyll fluorescence in the Amazon forest.

Photosynthesis of the Amazon rainforest plays an important role in the regional and global carbon cycles, but, despite considerable in situ and space-based observations, it has been intensely debated whether there is a dry-season increase in greenness and photosynthesis of the moist tropical Amazonian forests. Solar-induced chlorophyll fluorescence (SIF), which is emitted by chlorophyll, has a strong positive linear relationship with photosynthesis at the canopy scale. Recent advancements have allowed us to observe SIF globally with Earth observation satellites. Here we show that forest SIF did not decrease in the early dry season and increased substantially in the late dry season and early part of wet season, using SIF data from the Tropospheric Monitoring Instrument (TROPOMI), which has unprecedented spatial resolution and near-daily global coverage. Using in situ CO2 eddy flux data, we also show that cloud cover rarely affects photosynthesis at TROPOMI's midday overpass, a time when the forest canopy is most often light-saturated. The observed dry-season increases of forest SIF are not strongly affected by sun-sensor geometry, which was attributed as creating a pseudo dry-season green-up in the surface reflectance data. Our results provide strong evidence that greenness, SIF, and photosynthesis of the tropical Amazonian forest increase during the dry season.

Divergent trends of open-surface water body area in the contiguous United States from 1984 to 2016.

The contiguous United States (CONUS), especially the West, faces challenges of increasing water stress and uncertain impacts of climate change. The historical information of surface water body distribution, variation, and multidecadal trends documented in remote-sensing images can aid in water-resource planning and management, yet is not well explored. Here, we detected open-surface water bodies in all Landsat 5, 7, and 8 images (∼370,000 images, >200 TB) of the CONUS and generated 30-meter annual water body frequency maps for 1984-2016. We analyzed the interannual variations and trends of year-long water body area, examined the impacts of climatic and anthropogenic drivers on water body area dynamics, and explored the relationships between water body area and land water storage (LWS). Generally, the western half of the United States is prone to water stress, with small water body area and large interannual variability. During 1984-2016, water-poor regions of the Southwest and Northwest had decreasing trends in water body area, while water-rich regions of the Southeast and far north Great Plains had increasing trends. These divergent trends, mainly driven by climate, enlarged water-resource gaps and are likely to continue according to climate projections. Water body area change is a good indicator of LWS dynamics in 58% of the CONUS. Following the 2012 prolonged drought, LWS in California and the southern Great Plains had a larger decrease than surface water body area, likely caused by massive groundwater withdrawals. Our findings provide valuable information for surface water-resource planning and management across the CONUS.

Tracking annual changes of coastal tidal flats in China during 1986-2016 through analyses of Landsat images with Google Earth Engine.

Tidal flats (non-vegetated area), along with coastal vegetation area, constitute the coastal wetlands (intertidal zone) between high and low water lines, and play an important role in wildlife, biodiversity and biogeochemical cycles. However, accurate annual maps of coastal tidal flats over the last few decades are unavailable and their spatio-temporal changes in China are unknown. In this study, we analyzed all the available Landsat TM/ETM+/OLI imagery (~ 44,528 images) using the Google Earth Engine (GEE) cloud computing platform and a robust decision tree algorithm to generate annual frequency maps of open surface water body and vegetation to produce annual maps of coastal tidal flats in eastern China from 1986 to 2016 at 30-m spatial resolution. The resulting map of coastal tidal flats in 2016 was evaluated using very high-resolution images available in Google Earth. The total area of coastal tidal flats in China in 2016 was about 731,170 ha, mostly distributed in the provinces around Yellow River Delta and Pearl River Delta. The interannual dynamics of coastal tidal flats area in China over the last three decades can be divided into three periods: a stable period during 1986-1992, an increasing period during 1993-2001 and a decreasing period during 2002-2016. The resulting annual coastal tidal flats maps could be used to support sustainable coastal zone management policies that preserve coastal ecosystem services and biodiversity in China.

Mapping coastal wetlands of China using time series Landsat images in 2018 and Google Earth Engine.

Coastal wetlands, composed of coastal vegetation and non-vegetated tidal flats, play critical roles in biodiversity conservation, food production, and the global economy. Coastal wetlands in China are changing quickly due to land reclamation from sea, aquaculture, industrialization, and urbanization. However, accurate and updated maps of coastal wetlands (including vegetation and tidal flats) in China are unavailable, and the detailed spatial distribution of coastal wetlands are unknown. Here, we developed a new pixel- and phenology-based algorithm to identify and map coastal wetlands in China for 2018 using time series Landsat imagery (2,798 ETM+/OLI images) and the Google Earth Engine (GEE). The resultant map had a very high overall accuracy (98%). There were 7,474.6 km2 of coastal wetlands in China in 2018, which included 5,379.8 km2 of tidal flats, 1,856.4 km2 of deciduous wetlands, and 238.3 km2 of evergreen wetlands. Jiangsu Province had the largest area of coastal wetlands in China, followed by Shandong, Fujian, and Zhejiang Provinces. Our study demonstrates the high potential of time series Landsat images, pixel- and phenology-based algorithm, and GEE for mapping coastal wetlands at large scales. The resultant coastal wetland maps at 30-m spatial resolution serve as the most current dataset for sustainable management, ecological assessments, and conservation of coastal wetlands in China.

Enhanced gross primary production and evapotranspiration in juniper-encroached grasslands.

Woody plant encroachment (WPE) into grasslands has been occurring globally and may be accelerated by climate change in the future. This land cover change is expected to alter the carbon and water cycles, but it remains uncertain how and to what extent the carbon and water cycles may change with WPE into grasslands under current climate. In this study, we examined the difference of vegetation indices (VIs), evapotranspiration (ET), gross primary production (GPP), and solar-induced chlorophyll fluorescence (SIF) during 2000-2010 between grasslands and juniper-encroached grasslands. We also quantitatively assessed the changes of GPP and ET for grasslands with different proportions of juniper encroachment (JWPE). Our results suggested that JWPE increased the GPP, ET, greenness-related VIs, and SIF of grasslands. Mean annual GPP and ET were, respectively, ~55% and ~45% higher when grasslands were completely converted into juniper forests under contemporary climate during 2000-2010. The enhancement of annual GPP and ET for grasslands with JWPE varied over years ranging from about +20% GPP (~+30% for ET) in the wettest year (2007) to about twice as much GPP (~+55% for ET) in the severe drought year (2006) relative to grasslands without encroachment. Additionally, the differences in GPP and ET showed significant seasonal dynamics. During the peak growing season (May-August), GPP and ET for grasslands with JWPE were ~30% and ~40% higher on average. This analysis provided insights into how and to what degree carbon and water cycles were impacted by JWPE, which is vital to understanding how JWPE and ecological succession will affect the regional and global carbon and water budgets in the future.

Exacerbated grassland degradation and desertification in Central Asia during 2000-2014.

Grassland degradation and desertification is a complex process, including both state conversion (e.g., grasslands to deserts) and gradual within-state change (e.g., greenness dynamics). Existing studies hardly separated the two components and analyzed it as a whole based on time series vegetation index data, which cannot provide a clear and comprehensive picture for grassland degradation and desertification. Here we propose an integrated assessment strategy, by considering both state conversion and within-state change of grasslands, to investigate grassland degradation and desertification process in Central Asia. First, annual maps of grasslands and sparsely vegetated land were generated to track the state conversions between them. The results showed increasing grasslands were converted to sparsely vegetated lands from 2000 to 2014, with the desertification region concentrating in the latitude range of 43-48° N. A frequency analysis of grassland vs. sparsely vegetated land classification in the last 15 yr allowed a recognition of persistent desert zone (PDZ), persistent grassland zone (PGZ), and transitional zone (TZ). The TZ was identified in southern Kazakhstan as one hotspot that was unstable and vulnerable to desertification. Furthermore, the trend analysis of Enhanced Vegetation Index during thermal growing season (EVITGS ) was investigated in individual zones using linear regression and Mann-Kendall approaches. An overall degradation across the area was found; moreover, the second desertification hotspot was identified in northern Kazakhstan with significant decreasing in EVITGS , which was located in PGZ. Finally, attribution analyses of grassland degradation and desertification were conducted by considering precipitation, temperature, and three different drought indices. We found persistent droughts were the main factor for grassland degradation and desertification in Central Asia. Considering both state conversion and gradual within-state change processes, this study provided reference information for identification of desertification hotspots to support further grassland degradation and desertification treatment, and the method could be useful to be extended to other regions.

Mapping paddy rice planting area in northeastern Asia with Landsat 8 images, phenology-based algorithm and Google Earth Engine.

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.

Mapping paddy rice planting area in rice-wetland coexistent areas through analysis of Landsat 8 OLI and MODIS images.

Accurate and up-to-date information on the spatial distribution of paddy rice fields is necessary for the studies of trace gas emissions, water source management, and food security. The phenology-based paddy rice mapping algorithm, which identifies the unique flooding stage of paddy rice, has been widely used. However, identification and mapping of paddy rice in rice-wetland coexistent areas is still a challenging task. In this study, we found that the flooding/transplanting periods of paddy rice and natural wetlands were different. The natural wetlands flood earlier and have a shorter duration than paddy rice in the Panjin Plain, a temperate region in China. We used this asynchronous flooding stage to extract the paddy rice planting area from the rice-wetland coexistent area. MODIS Land Surface Temperature (LST) data was used to derive the temperature-defined plant growing season. Landsat 8 OLI imagery was used to detect the flooding signal and then paddy rice was extracted using the difference in flooding stages between paddy rice and natural wetlands. The resultant paddy rice map was evaluated with in-situ ground-truth data and Google Earth images. The estimated overall accuracy and Kappa coefficient were 95% and 0.90, respectively. The spatial pattern of OLI-derived paddy rice map agrees well with the paddy rice layer from the National Land Cover Dataset from 2010 (NLCD-2010). The differences between RiceLandsat and RiceNLCD are in the range of ±20% for most 1-km grid cell. The results of this study demonstrate the potential of the phenology-based paddy rice mapping algorithm, via integrating MODIS and Landsat 8 OLI images, to map paddy rice fields in complex landscapes of paddy rice and natural wetland in the temperate region.

Mapping paddy rice planting areas through time series analysis of MODIS land surface temperature and vegetation index data.

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.

Mapping paddy rice planting area in cold temperate climate region through analysis of time series Landsat 8 (OLI), Landsat 7 (ETM+) and MODIS imagery.

Accurate and timely rice paddy field maps with a fine spatial resolution would greatly improve our understanding of the effects of paddy rice agriculture on greenhouse gases emissions, food and water security, and human health. Rice paddy field maps were developed using optical images with high temporal resolution and coarse spatial resolution (e.g., Moderate Resolution Imaging Spectroradiometer (MODIS)) or low temporal resolution and high spatial resolution (e.g., Landsat TM/ETM+). In the past, the accuracy and efficiency for rice paddy field mapping at fine spatial resolutions were limited by the poor data availability and image-based algorithms. In this paper, time series MODIS and Landsat ETM+/OLI images, and the pixel- and phenology-based algorithm are used to map paddy rice planting area. The unique physical features of rice paddy fields during the flooding/open-canopy period are captured with the dynamics of vegetation indices, which are then used to identify rice paddy fields. The algorithm is tested in the Sanjiang Plain (path/row 114/27) in China in 2013. The overall accuracy of the resulted map of paddy rice planting area generated by both Landsat ETM+ and OLI is 97.3%, when evaluated with areas of interest (AOIs) derived from geo-referenced field photos. The paddy rice planting area map also agrees reasonably well with the official statistics at the level of state farms (R2 = 0.94). These results demonstrate that the combination of fine spatial resolution images and the phenology-based algorithm can provide a simple, robust, and automated approach to map the distribution of paddy rice agriculture in a year.

Large spatial variation and stagnation of cropland gross primary production increases the challenges of sustainable grain production and food security in China.

Sustainable crop grain production and food security is a grand societal challenge. Substantial investments in China's agriculture have been made in the past decades, but our knowledge on cropland gross primary production in China remains limited. Here we analyzed gross primary production (GPP), solar-induced chlorophyll fluorescence (SIF), terrestrial water storage, crop grain production, and agricultural investment and policy during 2000-2018. We found that based on croplands in 2000, approximately 52 × 106 ha (~37%) had continuous increasing trends in GPP during 2000-2018, which were mainly located in northern China. GPP for 63% of croplands was stagnant, declined, or had no significant change. At the national scale, annual cropland GPP increased during 2000-2008 but became stagnant in 2009-2018, which was inconsistent with the interannual trend in the crop grain production data for 2009-2018. The spatial mismatch between crop production and water availability became worse. The major grain exporting provinces, mostly located in water-stressed regions, experienced increased water resource constraints, which posed a challenge for sustainable grain production. The stagnant cropland GPP and increasing water resource constraints highlight the urgent need for sustainable management for crop production and food security in China.

Fingerprint of rice paddies in spatial-temporal dynamics of atmospheric methane concentration in monsoon Asia.

Agriculture (e.g., rice paddies) has been considered one of the main emission sources responsible for the sudden rise of atmospheric methane concentration (XCH4) since 2007, but remains debated. Here we use satellite-based rice paddy and XCH4 data to investigate the spatial-temporal relationships between rice paddy area, rice plant growth, and XCH4 in monsoon Asia, which accounts for ~87% of the global rice area. We find strong spatial consistencies between rice paddy area and XCH4 and seasonal consistencies between rice plant growth and XCH4. Our results also show a decreasing trend in rice paddy area in monsoon Asia since 2007, which suggests that the change in rice paddy area could not be one of the major drivers for the renewed XCH4 growth, thus other sources and sinks should be further investigated. Our findings highlight the importance of satellite-based paddy rice datasets in understanding the spatial-temporal dynamics of XCH4 in monsoon Asia.

Gainers and losers of surface and terrestrial water resources in China during 1989-2016.

Data and knowledge of the spatial-temporal dynamics of surface water area (SWA) and terrestrial water storage (TWS) in China are critical for sustainable management of water resources but remain very limited. Here we report annual maps of surface water bodies in China during 1989-2016 at 30m spatial resolution. We find that SWA decreases in water-poor northern China but increases in water-rich southern China during 1989-2016. Our results also reveal the spatial-temporal divergence and consistency between TWS and SWA during 2002-2016. In North China, extensive and continued losses of TWS, together with small to moderate changes of SWA, indicate long-term water stress in the region. Approximately 569 million people live in those areas with deceasing SWA or TWS trends in 2015. Our data set and the findings from this study could be used to support the government and the public to address increasing challenges of water resources and security in China.

Area extraction and spatiotemporal characteristics of winter wheat-summer maize in Shandong Province using NDVI time series.

The use of remote sensing to rapidly and accurately obtain information on the spatiotemporal distribution of large-scale wheat and maize acreage is of great significance for improving the level of food production management and ensuring food security. We constructed a MODIS-NDVI time series dataset, combined linear interpolation and the Harmonic Analysis of Time Series algorithm to smooth the time series data curve, and classified the data with random forest algorithms. The results show that winter wheat-summer maize planting areas were mainly distributed in the western plains, southern region, and north-eastern part of the middle mountainous regions while the eastern hilly regions were less distributed and scattered. The winter wheat-summer maize planting areas in the study area continued to grow from 2004-2016, with the most significant growth in the northern part of the western plains and Yellow River Delta. The spatial planting probability reflected the planting core area and showed an intensive planting pattern. During the study period, the peak value and time for the NDVI of the winter wheat were significantly different and showed an increasing trend, while these parameters for the summer maize were relatively stable with little change. Therefore, we mapped a spatial distribution of the winter wheat and summer maize, using the time series data pre-processing synthesis and phenology curve random forest classification methods. Through precision analysis, we obtained satisfactory results, which provided a straightforward and efficient method to monitor the winter wheat and summer maize.

Satellite-observed pantropical carbon dynamics.

Changes in terrestrial tropical carbon stocks have an important role in the global carbon budget. However, current observational tools do not allow accurate and large-scale monitoring of the spatial distribution and dynamics of carbon stocks1. Here, we used low-frequency L-band passive microwave observations to compute a direct and spatially explicit quantification of annual aboveground carbon (AGC) fluxes and show that the tropical net AGC budget was approximately in balance during 2010 to 2017, the net budget being composed of gross losses of -2.86 PgC yr-1 offset by gross gains of -2.97 PgC yr-1 between continents. Large interannual and spatial fluctuations of tropical AGC were quantified during the wet 2011 La Niña year and throughout the extreme dry and warm 2015-2016 El Niño episode. These interannual fluctuations, controlled predominantly by semiarid biomes, were shown to be closely related to independent global atmospheric CO2 growth-rate anomalies (Pearson's r = 0.86), highlighting the pivotal role of tropical AGC in the global carbon budget.

A global moderate resolution dataset of gross primary production of vegetation for 2000-2016.

Accurate estimation of the gross primary production (GPP) of terrestrial vegetation is vital for understanding the global carbon cycle and predicting future climate change. Multiple GPP products are currently available based on different methods, but their performances vary substantially when validated against GPP estimates from eddy covariance data. This paper provides a new GPP dataset at moderate spatial (500 m) and temporal (8-day) resolutions over the entire globe for 2000-2016. This GPP dataset is based on an improved light use efficiency theory and is driven by satellite data from MODIS and climate data from NCEP Reanalysis II. It also employs a state-of-the-art vegetation index (VI) gap-filling and smoothing algorithm and a separate treatment for C3/C4 photosynthesis pathways. All these improvements aim to solve several critical problems existing in current GPP products. With a satisfactory performance when validated against in situ GPP estimates, this dataset offers an alternative GPP estimate for regional to global carbon cycle studies.

Continued decrease of open surface water body area in Oklahoma during 1984-2015.

Oklahoma contains the largest number of manmade lakes and reservoirs in the United States. Despite the importance of these open surface water bodies to public water supply, agriculture, thermoelectric power, tourism and recreation, it is unclear how these water bodies have responded to climate change and anthropogenic water exploitation in past decades. In this study, we used all available Landsat 5 and 7 images (16,000 scenes) from 1984 through 2015 and a water index- and pixel-based approach to analyze the spatial-temporal variability of open surface water bodies and its relationship with climate and water exploitation. Specifically, the areas and numbers of four water body extents (the maximum, year-long, seasonal, and average extents) were analyzed to capture variations in water body area and number. Statistically significant downward trends were found in the maximum, year-long, and annual average water body areas from 1984 through 2015. Furthermore, these decreases were mainly attributed to the continued shrinking of large water bodies (>1km2). There were also significant decreases in maximum and year-long water body numbers, which suggested that some of the water bodies were vanishing year by year. However, remarkable inter-annual variations of water body area and number were also found. Both water body area and number were positively related to precipitation, and negatively related to temperature. Surface water withdrawals mainly influenced the year-long water bodies. The smaller water bodies have a higher risk of drying under a drier climate, which suggests that small water bodies are more vulnerable under climate-warming senarios.

Temporal consistency between gross primary production and solar-induced chlorophyll fluorescence in the ten most populous megacity areas over years.

The gross primary production (GPP) of vegetation in urban areas plays an important role in the study of urban ecology. It is difficult however, to accurately estimate GPP in urban areas, mostly due to the complexity of impervious land surfaces, buildings, vegetation, and management. Recently, we used the Vegetation Photosynthesis Model (VPM), climate data, and satellite images to estimate the GPP of terrestrial ecosystems including urban areas. Here, we report VPM-based GPP (GPPvpm) estimates for the world's ten most populous megacities during 2000-2014. The seasonal dynamics of GPPvpm during 2007-2014 in the ten megacities track well that of the solar-induced chlorophyll fluorescence (SIF) data from GOME-2 at 0.5° × 0.5° resolution. Annual GPPvpm during 2000-2014 also shows substantial variation among the ten megacities, and year-to-year trends show increases, no change, and decreases. Urban expansion and vegetation collectively impact GPP variations in these megacities. The results of this study demonstrate the potential of a satellite-based vegetation photosynthesis model for diagnostic studies of GPP and the terrestrial carbon cycle in urban areas.

Spatiotemporal patterns of paddy rice croplands in China and India from 2000 to 2015.

Due to rapid population growth and urbanization, paddy rice agriculture is experiencing substantial changes in the spatiotemporal pattern of planting areas in the two most populous countries-China and India-where food security is always the primary concern. However, there is no spatially explicit and continuous rice-planting information in either country. This knowledge gap clearly hinders our ability to understand the effects of spatial paddy rice area dynamics on the environment, such as food and water security, climate change, and zoonotic infectious disease transmission. To resolve this problem, we first generated annual maps of paddy rice planting areas for both countries from 2000 to 2015, which are derived from time series Moderate Resolution Imaging Spectroradiometer (MODIS) data and the phenology- and pixel-based rice mapping platform (RICE-MODIS), and analyzed the spatiotemporal pattern of paddy rice dynamics in the two countries. We found that China experienced a general decrease in paddy rice planting area with a rate of 0.72 million (m) ha/yr from 2000 to 2015, while a significant increase at a rate of 0.27mha/yr for the same time period happened in India. The spatial pattern of paddy rice agriculture in China shifted northeastward significantly, due to simultaneous expansions in paddy rice planting areas in northeastern China and contractions in southern China. India showed an expansion of paddy rice areas across the entire country, particularly in the northwestern region of the Indo-Gangetic Plain located in north India and the central and south plateau of India. In general, there has been a northwesterly shift in the spatial pattern of paddy rice agriculture in India. These changes in the spatiotemporal patterns of paddy rice planting area have raised new concerns on how the shift may affect national food security and environmental issues relevant to water, climate, and biodiversity.

Mapping paddy rice distribution using multi-temporal Landsat imagery in the Sanjiang Plain, northeast China.

Information of paddy rice distribution is essential for food production and methane emission calculation. Phenology-based algorithms have been utilized in the mapping of paddy rice fields by identifying the unique flooding and seedling transplanting phases using multi-temporal moderate resolution (500 m to 1 km) images. In this study, we developed simple algorithms to identify paddy rice at a fine resolution at the regional scale using multi-temporal Landsat imagery. Sixteen Landsat images from 2010-2012 were used to generate the 30 m paddy rice map in the Sanjiang Plain, northeast China-one of the major paddy rice cultivation regions in China. Three vegetation indices, Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), and Land Surface Water Index (LSWI), were used to identify rice fields during the flooding/transplanting and ripening phases. The user and producer accuracies of paddy rice on the resultant Landsat-based paddy rice map were 90% and 94%, respectively. The Landsat-based paddy rice map was an improvement over the paddy rice layer on the National Land Cover Dataset, which was generated through visual interpretation and digitalization on the fine-resolution images. The agricultural census data substantially underreported paddy rice area, raising serious concern about its use for studies on food security.