Mapping paddy rice agriculture in southern China using multi-temporal MODIS images

Information on the area and spatial distribution of paddy rice fields is needed for trace gas emission estimates, management of water resources, and food security. Paddy rice fields are characterized by an initial period of flooding and transplanting, during which period open canopy (a mixture of surface water and rice crops) exists. The Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard the NASA EOS Terra satellite has visible, near infrared and shortwave infrared bands; and therefore, a number of vegetation indices can be calculated, including Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and Land Surface Water Index (LSWI) that is sensitive to leaf water and soil moisture. In this study, we developed a paddy rice mapping algorithm that uses time series of three vegetation indices (LSWI, EVI, and NDVI) derived from MODIS images to identify that initial period of flooding and transplanting in paddy rice fields, based on the sensitivity of LSWI to the increased surface moisture during the period of flooding and rice transplanting. We ran the algorithm to map paddy rice fields in 13 provinces of southern China, using the 8-day composite MODIS Surface Reflectance products (500-m spatial resolution) in 2002. The resultant MODIS-derived paddy rice map was evaluated, using the National Land Cover Dataset (1:100,000 scale) derived from analysis of Landsat ETM+ images in 1999/2000. There were reasonable agreements in area estimates of paddy rice fields between the MODIS-derived map and the Landsat-based dataset at the provincial and county levels. The results of this study indicated that the MODIS-based paddy rice mapping algorithm could potentially be applied at large spatial scales to monitor paddy rice agriculture on a timely and frequent basis.     

Observation of flooding and rice transplanting of paddy rice fields at the site to landscape scales in China using VEGETATION sensor data

A unique physical feature of paddy rice fields is that rice is grown on flooded soil. During the period of flooding and rice transplanting, there is a large proportion of surface water in a land surface consisting of water, vegetation and soils. The VEGETATION (VGT) sensor has four spectral bands that are equivalent to spectral bands of Landsat TM, and its mid-infrared spectral band is very sensitive to soil moisture and plant canopy water content. In this study we evaluated a VGT-derived normalized difference water index (NDWI VGT =(B3-MIR)/ (B3+MIR)) for describing temporal and spatial dynamics of surface moisture. Twenty-seven 10-day composites (VGT- S10) from 1 March to 30 November 1999 were acquired and analysed for a study area (175 km by 165 km) in eastern Jiangsu Province, China, where a winter wheat and paddy rice double cropping system dominates the landscape. We compared the temporal dynamics and spatial patterns of normalized difference vegetation index (NDVI VGT ) and NDWI VGT . The NDWI VGT temporal dynamics were sensitive enough to capture the substantial increases of surface water due to flooding and rice transplanting at paddy rice fields. A land use thematic map for the timing and location of flooding and rice transplanting was generated for the study area. Our results indicate that NDWI and NDVI temporal anomalies may provide a simple and effective tool for detection of flooding and rice transplanting across the landscape.     

Use of RADARSAT‐1 data and a digital elevation model to assess flood damage and improve rice production in the lower part of the Chi River Basin,Thailand

In Thailand, flooding due to seasonal monsoon conditions frequently destroys a substantial amount of rice production, the most important agricultural activity of the country. Taking the 2001 monsoon flooding that hit the Lower Chi River Basin as an example, we developed a new method for accurately assessing damage to flood‐affected paddies. A RADARSAT‐1 image acquired during peak flooding was combined with a 30‐m digital elevation model (DEM) to develop a ‘flood‐level‐determination’ algorithm for estimating floodwater depth. Based on the elongation capability of the rice varieties, a water depth of 80 cm was used to separate ‘non‐damaged’ from ‘damaged’ paddy areas, indicating that about 60% of the paddy fields in the flooded areas were non‐damaged paddies. To minimize the loss of rice and maximize farmers' incomes, a map of rice varieties appropriate for the damaged paddy areas was produced, combining the flood‐affected paddy map with the flood frequency map. Our results demonstrate the potential of using single‐date RADARSAT‐1 data and a DEM to provide accurate and economic means of assessing flood damage to rice fields that can be used to improve rice production.     

Mapping of rice-cropping pattern and cultural type using remote-sensing and ancillary data: a case study for South and Southeast Asian countries

A geospatial database on the spatial distribution of rice areas and rice cultural types of major rice-producing countries of South and Southeast Asia has been developed in this study using remote-sensing and ancillary data sets. Multitemporal SPOT VGT normalized difference vegetation index (NDVI) data for the period 2009–2010 were used for the analysis. The classification was performed adopting ISODATA clustering to build a non-agricultural area mask followed by rice area mapping. The derived rice area was stratified by logical modelling of ancillary data sets into five rice cultural types: irrigated wet, upland, flood-prone, drought-prone, and deep-water. The uniqueness of this study is a synergistic approach based solely on single-source, high-temporal remote-sensing data coupled with ancillary data, which demonstrate the application of SPOT VGT NDVI data in building a geospatial database for rice crops over a wide spatial extent. This approach was adopted for cost effectivity as the study extent was vast and thus lacking ground truth information. Comparison of the derived rice area against the reported literature values for validation yielded a good correlation (linear coefficient of determination, R² = 0.95–0.99). The high-temporal resolution NDVI data enabled effective characterization of vegetation phenology. The derived spatial outputs can be used in various studies associated with the assessment of greenhouse gas emissions from paddy fields, change detection, and inputs to crop simulation models, which are significantly related to different rice cultural types.     

Analysis of the applicability of multi-temporal full polarimetric airborne L-band SAR scattering to paddy rice field mapping

ABSTRACT

A Synthetic Aperture Radar (SAR) is an all-weather imaging system that is often used for mapping paddy rice fields and estimating the area. Fully polarimetric SAR is used to detect the microwave scattering property. In this study, a simple threshold analysis of fully polarimetric L-band SAR data was conducted to distinguish paddy rice fields from soybean and other fields. We analysed a set of ten airborne SAR L-band 2 (Pi-SAR-L2) images obtained during the paddy rice growing season (in June, August, and September) from 2012 to 2014 using polarimetric decomposition. Vector data for agricultural land use areas were overlaid on the analysed images and the mean value for each agricultural parcel computed. By quantitatively comparing our data with a reference dataset generated from optical sensor images, effective polarimetric parameters and the ideal observation season were revealed. Double bounce scattering and surface scattering component ratios, derived using a four-component decomposition algorithm, were key to extracting paddy rice fields when the plant stems are vertical with respect to the ground. The alpha angle was also an effective factor for extracting rice fields from an agricultural area. The data obtained during August show maximum agreement with the reference dataset of estimated paddy rice field areas.     

Improved classification of conservation tillage adoption using high temporal and synthetic satellite imagery

Conservation tillage management has been advocated for carbon sequestration and soil quality preservation purposes. Past satellite image analyses have had difficulty in differentiating between no-till (NT) and minimal tillage (MT) conservation classes due to similarities in surface residues, and may have been restricted by the availability of cloud-free satellite imagery. This study hypothesized that the inclusion of high temporal data into the classification process would increase conservation tillage accuracy due to the added likelihood of capturing spectral changes in MT fields following a tillage disturbance. Classification accuracies were evaluated for Random Forest models based on 250-m and 500-m MODIS, 30-m Landsat, and 30-m synthetic reflectance values. Synthetic (30-m) data derived from the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) were evaluated because high frequency Landsat image sets are often unavailable within a cropping season due to cloud issues. Classification results from a five-date Landsat model were substantially better than those reported by previous classification tillage studies, with 94% total and ≥ 88% class producer's accuracies. Landsat-derived models based on individual image scenes (May through August) yielded poor MT classifications, but a monthly increase in accuracy illustrated the importance of temporal sampling for capturing regional tillage disturbance signatures. MODIS-based model accuracies (90% total; ≥ 82% class) were lower than in the five-date Landsat model, but were higher than previous image-based and survey-based tillage classification results. Almost all the STARFM prediction-based models had classification accuracies higher than, or comparable to, the MODIS-based results (> 90% total; ≥ 84% class) but the resulting model accuracies were dependent on the MODIS/Landsat base pairs used to generate the STARFM predictions. Also evident within the STARFM prediction-based models was the ability for high frequency data series to compensate for degraded synthetic spectral values when classifying field-based tillage. The decision to use MODIS or STARFM-based data within conservation tillage analysis is likely situation dependent. A MODIS-based approach requires little data processing and could be more efficient for large-area mapping; however a STARFM-based analysis might be more appropriate in mixed-pixel situations that could potentially compromise classification accuracy.     

Tracking seasonal drought effects on ecosystem light use efficiency with satellite-based PRI in a Mediterranean forest

Gross primary productivity (GPP) changes occur at different time-scales and due to various mechanisms such as variations in leaf area, chlorophyll content, rubisco activity, and stomatal conductance. Diagnostic estimates of primary productivity are obviously error prone when these changes are not accounted for. Additional complications arise when factors inuencing a biome-specific maximum light use efficiency (LUE) must be estimated over a large area. In these cases a direct estimation of ecosystem LUE could reduce uncertainty of GPP estimates. Here, we analyse whether a MODIS-based photochemical reectance index (PRI) is a useful proxy for the light use efficiency of a Mediterranean Quercus ilex forest. As the originally proposed reference band for PRI is not available on MODIS, we tested the reference bands 1 (620-670 nm), 4 (545-565 nm), 12 (546-556 nm), 13 (662-672 nm), and 14 (673-683 nm) using different atmospheric correction algorithms. We repeated the analysis with different temporal resolutions of LUE (half-hourly to daily). The strongest correlation between LUE and PRI was found when considering only a narrow range of viewing angles at a time (especially 0-10° and 30-40°). We found that the MODIS-based PRI was able to track ecosystem LUE even during severe summer time water limitation. For this Mediterranean-type ecosystem we could show that a GPP estimation based on PRI is a huge improvement compared to the MODIS GPP algorithm. In this study, MODIS spectral band 1 turned out to be the most suitable reference band for PRI, followed by the narrow red bands 13 and 14. As to date no universally applicable reference band was identified in MODIS-based PRI studies, we advocate thorough testing for the optimal band combination in future studies.     

A polygon-based clustering and analysis framework for mining spatial datasets

Polygons provide natural representations for many types of geospatial objects, such as countries, buildings, and pollution hotspots. Thus, polygon-based data mining techniques are particularly useful for mining geospatial datasets. In this paper, we propose a polygon-based clustering and analysis framework for mining multiple geospatial datasets that have inherently hidden relations. In this framework, polygons are first generated from multiple geospatial point datasets by using a density-based contouring algorithm called DCONTOUR. Next, a density-based clustering algorithm called Poly-SNN with novel dissimilarity functions is employed to cluster polygons to create meta-clusters of polygons. Finally, post-processing analysis techniques are proposed to extract interesting patterns and user-guided summarized knowledge from meta-clusters. These techniques employ plug-in reward functions that capture a domain expert’s notion of interestingness to guide the extraction of knowledge from meta-clusters. The effectiveness of our framework is tested in a real-world case study involving ozone pollution events in Texas. The experimental results show that our framework can reveal interesting relationships between different ozone hotspots represented by polygons; it can also identify interesting hidden relations between ozone hotspots and several meteorological variables, such as outdoor temperature, solar radiation, and wind speed.     

Annually modelling built-settlements between remotely-sensed observations using relative changes in subnational populations and lights at night

Mapping urban features/human built-settlement extents at the annual time step has a wide variety of applications in demography, public health, sustainable development, and many other fields. Recently, while more multitemporal urban features/human built-settlement datasets have become available, issues still exist in remotely-sensed imagery due to spatial and temporal coverage, adverse atmospheric conditions, and expenses involved in producing such datasets. Remotely-sensed annual time-series of urban/built-settlement extents therefore do not yet exist and cover more than specific local areas or city-based regions. Moreover, while a few high-resolution global datasets of urban/built-settlement extents exist for key years, the observed date often deviates many years from the assigned one. These challenges make it difficult to increase temporal coverage while maintaining high fidelity in the spatial resolution. Here we describe an interpolative and flexible modelling framework for producing annual built-settlement extents. We use a combined technique of random forest and spatio-temporal dasymetric modelling with open source subnational data to produce annual 100 m × 100 m resolution binary built-settlement datasets in four test countries located in varying environmental and developmental contexts for test periods of five-year gaps. We find that in the majority of years, across all study areas, the model correctly identified between 85 and 99% of pixels that transition to built-settlement. Additionally, with few exceptions, the model substantially out performed a model that gave every pixel equal chance of transitioning to built-settlement in each year. This modelling framework shows strong promise for filling gaps in cross-sectional urban features/built-settlement datasets derived from remotely-sensed imagery, provides a base upon which to create urban future/built-settlement extent projections, and enables further exploration of the relationships between urban/built-settlement area and population dynamics.     

A hybrid classification scheme for mining multisource geospatial data

Supervised learning methods such as Maximum Likelihood (ML) are often used in land cover (thematic) classification of remote sensing imagery. ML classifier relies exclusively on spectral characteristics of thematic classes whose statistical distributions (class conditional probability densities) are often overlapping. The spectral response distributions of thematic classes are dependent on many factors including elevation, soil types, and ecological zones. A second problem with statistical classifiers is the requirement of the large number of accurate training samples (10 to 30 × |dimensions|), which are often costly and time consuming to acquire over large geographic regions. With the increasing availability of geospatial databases, it is possible to exploit the knowledge derived from these ancillary datasets to improve classification accuracies even when the class distributions are highly overlapping. Likewise newer semi-supervised techniques can be adopted to improve the parameter estimates of the statistical model by utilizing a large number of easily available unlabeled training samples. Unfortunately, there is no convenient multivariate statistical model that can be employed for multisource geospatial databases. In this paper we present a hybrid semi-supervised learning algorithm that effectively exploits freely available unlabeled training samples from multispectral remote sensing images and also incorporates ancillary geospatial databases. We have conducted several experiments on Landsat satellite image datasets, and our new hybrid approach shows over 24% to 36% improvement in overall classification accuracy over conventional classification schemes.     

Mapping paddy rice fields by applying machine learning algorithms to multi-temporal Sentinel-1A and Landsat data

Sentinel-1A synthetic aperture radar (SAR) data present an opportunity for acquiring crop information without restrictions caused by weather and illumination conditions, at a spatial resolution appropriate for individual rice fields and a temporal resolution sufficient to capture the growth profiles of different crop species. This study investigated the use of multi-temporal Sentinel-1A SAR data and Landsat-derived normalized difference vegetation index (NDVI) data to map the spatial distribution of paddy rice fields across parts of the Sanjiang plain, in northeast China. The satellite sensor data were acquired throughout the rice crop-growing season (May–October). A co-registered set of 10 dual polarization (VH/VV) SAR and NDVI images depicting crop phenological development were used as inputs to Support Vector Machine (SVM) and Random Forest (RF) machine learning classification algorithms in order to map paddy rice fields. The results showed a significant increase in overall classification when the NDVI time-series data were integrated with the various combinations of multi-temporal polarization channels (i.e. VH, VV, and VH/VV). The highest classification accuracies overall (95.2%) and for paddy rice (96.7%) were generated using the RF algorithm applied to combined multi-temporal VH polarization and NDVI data. The SVM classifier was most effective when applied to the dual polarization (i.e. VH and VV) SAR data alone and this generated overall and paddy rice classification accuracies of 91.6% and 82.5%, respectively. The results demonstrate the practicality of implementing RF or SVM machine learning algorithms to produce 10 m spatial resolution maps of paddy rice fields with limited ground data using a combination of multi-temporal SAR and NDVI data, where available, or SAR data alone. The methodological framework developed in this study is apposite for large-scale implementation across China and other major rice-growing regions of the world.     

A vegetation index based technique for spatial sharpening of thermal imagery

High spatial resolution (∼ 100 m) thermal infrared band imagery has utility in a variety of applications in environmental monitoring. However, currently such data have limited availability and only at low temporal resolution, while coarser resolution thermal data (∼ 1000 m) are routinely available, but not as useful for identifying environmental features for many landscapes. An algorithm for sharpening thermal imagery (TsHARP) to higher resolutions typically associated with the shorter wavebands (visible and near-infrared) used to compute vegetation indices is examined over an extensive corn/soybean production area in central Iowa during a period of rapid crop growth. This algorithm is based on the assumption that a unique relationship between radiometric surface temperature (T_R) relationship and vegetation index (VI) exists at multiple resolutions. Four different methods for defining a VI − T_R basis function for sharpening were examined, and an optimal form involving a transformation to fractional vegetation cover was identified. The accuracy of the high-resolution temperature retrieval was evaluated using aircraft and Landsat thermal imagery, aggregated to simulate native and target resolutions associated with Landsat, MODIS, and GOES short- and longwave datasets. Applying TsHARP to simulated MODIS thermal maps at 1-km resolution and sharpening down to ∼ 250 m (MODIS VI resolution) yielded root-mean-square errors (RMSE) of 0.67-1.35 °C compared to the ‘observed’ temperature fields, directly aggregated to 250 m. Sharpening simulated Landsat thermal maps (60 and 120 m) to Landsat VI resolution (30 m) yielded errors of 1.8-2.4 °C, while sharpening simulated GOES thermal maps from 5 km to 1 km and 250 m yielded RMSEs of 0.98 and 1.97, respectively. These results demonstrate the potential for improving the spatial resolution of thermal-band satellite imagery over this type of rainfed agricultural region. By combining GOES thermal data with shortwave VI data from polar orbiters, thermal imagery with 250-m spatial resolution and 15-min temporal resolution can be generated with reasonable accuracy. Further research is required to examine the performance of TsHARP over regions with different climatic and land-use characteristics at local and regional scales.     

基于SAR与Landsat TM的小区域稻田的识别研究——以南京市江宁区为例

以南京市江宁区为研究区域,根据区域特征、作物物候期和水稻的生长特点,采用分层分类的方法提取稻田分布信息。通过比较多时相SAR数据、TM和多时相SAR融合与TM和单时相SAR融合数据识别水稻的精度和提取的水稻种植面积,分析了不同数据对区域多云雨,不同种植方式、面积小且分布破碎的水稻稻田的识别程度,并根据野外实地走访调查分析了主要影响因素。结果表明:多时相SAR数据、TM和多时相SAR数据的水稻识别精度都高于72%,高于TM和单时相SAR融合数据的结果;前两者提取的水稻种植面积和稻田分布接近,主要影响因素是地物分布、不同种植方式水稻物候期和水稻稻田面积小且分布破碎。     

不同水分管理下施用稻草对土壤有机酸和养分有效性的影响

为探索稻草还田短期内土壤有机酸积累和养分供应特征,在温室培养条件下,采用多孔聚酯管溶液采集器采集土壤溶液和"通用佳"离子交换树脂球法,分别研究了连续淹水和干湿交替下,施用稻草30天内对土壤有机酸产生和养分有效性的影响。结果表明,(1)无论连续淹水还是干湿交替,施用稻草后土壤中有机酸积累明显增加,但在干湿交替下的积累要低于连续淹水下;施用稻草后产生的有机酸以乙酸为主,且5cm处的有机酸积累明显低于15cm处。(2)施用稻草降低了土壤氮的有效性,但提高了土壤磷、铁、钾、钙和镁的有效性;与连续淹水相比,干湿交替下施用稻草降低土壤氮有效性的程度更大,而提高土壤磷、铁、钾、钙、镁的程度则较小。     

Growth monitoring and classification of rice fields using multitemporal RADARSAT-2 full-polarimetric data

Spaceborne synthetic aperture radar (SAR) can be used for agricultural monitoring. In this study, three single-polarimetric and four full-polarimetric observation data sets were analysed. A rice paddy field in northern Japan was used as the study site; the data for this site were obtained using RADARSAT-2, which carries a full-polarimetric C-band SAR. Soybean and grass fields were also present within the paddy fields. The temporal change in the backscattering coefficient of the rice paddy fields for the single-polarization data agreed with the temporal change obtained for a rice growth model based on radiative transfer theory. A three-component decomposition approach was applied to the full-polarimetric data. With each rice growth stage, the volume scattering component ratio increased, whereas the surface scattering component ratio generally decreased. The soybean and grass fields showed a smaller double-bounce scattering component than the rice fields for all the acquired data. The results of this study show that multitemporal observation by full-polarimetric SAR has great potential to be utilized for estimating rice-planted areas and monitoring rice growth.     

A geospatial tool for delineating streambanks

Ecologists, civil engineers, and conservation managers currently lack an automated software tool for delineating streambanks, which is an important input for various hydrological studies. Therefore, a computational method for automatically delineating streambanks using aerial imagery and stream centerline datasets was developed and incorporated in a stand-alone, user-friendly desktop-based geospatial tool. This interactive tool, titled Streambank Delineator (StreBanD), was tested on the 161-km Lower White River and 158-km L'Anguille River in Eastern Arkansas using 1-m resolution 3-band aerial imagery and National Hydrography Dataset Plus (NHDPlus) stream centerline. The near-infrared band and normalized difference water index (NDWI) were evaluated for feasibility of streambank delineation. Results showed that the NDWI, a multi-band index approach, consistently provided superior delineations when compared to the single band (near-infrared) approach especially for complex channel morphologies, including braided channels and meandering segments, in the L'Anguille River. The geospatial tool successfully delineated the streambanks for both rivers in less than 10 min with a mean error ranging from 0.3 m to 10.3 m when compared with five manually delineated streambanks. Due to the generic and simple nature of this tool, it should assist scientists and conservationists in rapidly delineating streambanks for their study areas.     

Regional scale mapping of fractional rice cropping change using a phenology-based temporal mixture analysis

Monitoring changes of paddy rice is challenging due to its diverse cropping patterns and spectral variation. To investigate the spatio-temporal changes of rice cropping, we used the 10-day composited Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) time series data with a spatial resolution of 250 m to map the sub-pixel rice spatial distributions in the Hunan Province, the top one region in rice planting area in southern of China. A method of improved phenology-based temporal mixture analysis (PTMA) was presented to identify early, middle, and late rice cropping patterns. The results show that the PTMA is effective to extract rice cropping. The nine rice cropping patterns were classified as early, middle, and late rice cropping, and fractional rice cropping within 250 m pixels was obtained to analyse the internal changes. Both the local planting conditions and different forms of rice cultivation were compared with statistical data. Overall, MODIS-estimated fractional rice agreed well with field samples at the pixel level and statistical data at the county level, which demonstrates the effectiveness of the PTMA method for mapping rice in these hilly regions with small-size paddy rice field. The changes show that single-cropping rice and double-cropping rice have been frequently transferred in space, which could be important information to support agricultural decision-making.     

A comparison of sampling designs for estimating deforestation from Landsat imagery: A case study of the Brazilian Legal Amazon

Three sampling designs — simple random, stratified random, and systematic sampling — are compared on the basis of precision of estimated loss of intact humid tropical forest area in the Brazilian Legal Amazon from 2000 to 2005. MODIS-derived deforestation is used to partition the study area into strata to intensify sampling within forest clearing hotspots. The precision of the estimator of deforestation area for each design is calculated from a population of wall-to-wall PRODES deforestation data available for the study area. Both systematic and stratified sampling yield smaller standard errors than simple random sampling, and the stratified design has smaller standard errors than the systematic design at each sample size evaluated. The results of this case study demonstrate the utility of a stratified design based on MODIS-derived deforestation data to improve precision of the estimated loss of intact forest area as estimated from sampling Landsat imagery.     

An interactive framework for spatial joins: a statistical approach to data analysis in GIS

Many Geographic Information Systems (GIS) handle a large volume of geospatial data. Spatial joins over two or more geospatial datasets are very common operations in GIS for data analysis and decision support. However, evaluating spatial joins can be very time intensive due to the size of datasets. In this paper, we propose an interactive framework that provides faster approximate answers of spatial joins. The proposed framework utilizes two statistical methods: probabilistic join and sampling based join. The probabilistic join method provides speedup of two orders of magnitude with no correctness guarantee, while the sampling based method provides an order of magnitude improvement over the full indexing tree joins of datasets and also provides running confidence intervals. The framework allows users to trade-off speed versus bounded accuracy, hence it provides truly interactive data exploration. The two methods are evaluated empirically with real and synthetic datasets.     

Mapping evapotranspiration using MODIS and MM5 Four-Dimensional Data Assimilation

Evapotranspiration (ET), the sum of evaporation from soil and transpiration from vegetation, is of vital importance in the hydrologic cycle and must be taken into consideration in assessments of the water resources of any region. The MODerate resolution Imaging Spectroradiometer (MODIS) sensor offers a promising opportunity for estimating daily ET with a 1 km spatial resolution, but is hampered by frequent cloud contamination or data gaps from other factors. In this study, 1) a stand-alone ET model was applied and tested during clear or partial cloudy sky conditions using MODIS-based inputs of land surface and atmospheric data and 2) meteorological simulations by using Four-Dimensional Data Assimilation (FDDA) system between MODIS and the 5th Generation Meso-scale Meteorological Model (MM5) was used in cloudy conditions to facilitate continuous daily ET estimates. The MODIS ET algorithm modified from Mu et al. (2007) is based on the Penman-Monteith equation and was applied to predict ET at flux measurement sites. This algorithm considers both the effects of surface energy partitioning processes and environmental constraints on ET. We devised gap-filling approaches for MODIS aerosol and albedo data that were identified as bottlenecks to determine retrieval rates of insolation and ET. MODIS-derived input variables (i.e., meteorological variables and radiation components) for estimating ET showed a good agreement with flux tower observations at each site. The retrieval rate of MODIS ET doubled at four flux measurement sites after gap-filling with negligible compensation was undertaken for accuracy. In spite of the high accuracy of MODIS-derived input variables, MODIS ET showed meaningful errors at the four flux measurement sites. These errors were mainly associated with errors in the estimated canopy conductance. During clear sky conditions, MODIS was used to calculate ET, while the MODIS-MM5 FDDA system provided input variables for the calculation of ET under cloudy sky conditions. The performance of the MODIS-MM5 FDDA system was evaluated by comparing ET based on MODIS, which showed a good agreement with the MODIS ET for various land cover types. Our results indicate that MODIS can be applied to monitor the land surface energy budget and ET with reasonable accuracy and that MODIS-MM5 FDDA has the potential to provide reasonable input data of ET estimation under cloudy conditions.