Mapping Crop Cycles in China Using MODIS-EVI Time Series

As the Earth’s population continues to grow and demand for food increases, the need for improved and timely information related to the properties and dynamics of global agricultural systems is becoming increasingly important. Global land cover maps derived from satellite data provide indispensable information regarding the geographic distribution and areal extent of global croplands. However, land use information, such as cropping intensity (defined here as the number of cropping cycles per year), is not routinely available over large areas because mapping this information from remote sensing is challenging. In this study, we present a simple but efficient algorithm for automated mapping of cropping intensity based on data from NASA’s (NASA: The National Aeronautics and Space Administration) MODerate Resolution Imaging Spectroradiometer (MODIS). The proposed algorithm first applies an adaptive Savitzky-Golay filter to smooth Enhanced Vegetation Index (EVI) time series derived from MODIS surface reflectance data. It then uses an iterative moving-window methodology to identify cropping cycles from the smoothed EVI time series. Comparison of results from our algorithm with national survey data at both the provincial and prefectural level in China show that the algorithm provides estimates of gross sown area that agree well with inventory data. Accuracy assessment comparing visually interpreted time series with algorithm results for a random sample of agricultural areas in China indicates an overall accuracy of 91.0% for three classes defined based on the number of cycles observed in EVI time series. The algorithm therefore appears to provide a straightforward and efficient method for mapping cropping intensity from MODIS time series data

Spatio-temporal Analysis of Agriculture in the Vietnamese Mekong Delta using MODIS Imagery

New methodologies using MODIS time‒series imagery were developed for revealing spatio‒temporal changes of agricultural environments and land use patterns in the Vietnamese Mekong Delta. The following methodologies were proposed:a Wavelet based Filter for Crop Phenology (WFCP), a Wavelet‒based fi lter for evaluating the spatial distribution of Cropping Systems (WFCS), and a Wavelet‒based fi lter for detecting spatio‒temporal changes in Flood Inundation(WFFI). The WFCP algorithm involves smoothing the temporal profi le of the Enhanced Vegetation Index (EVI) using the wavelet transform approach. As a result of validation using the agricultural statistical data in Japan, it was shown that the WFCP was able to estimate rice growing stages, including transplanting date, heading date and harvesting date from the smoothed EVI data, with 9‒12 days accuracy(RMSE). The WFCS algorithm was developed for detecting rice‒cropping patterns in the Vietnamese Mekong delta based on WFCP. It was revealed that the spatial distribution of rice cropping seasons was characterized by both annual fl ood inundation around the upper region in the rainy season and salinity intrusion around the coastal region in the dry season. The WFFI algorithm was developed for estimating start and end dates of fl ood inundation by using time‒series Land Surface Water Index and EVI. Annual intensity of Mekong fl oods was evaluated from 2000 to 2004, at a regional scale. Applying a series of wavelet‒based methodologies to the MODIS data acquired from 2000 to 2006, it was confi rmed that the cropping season for the winter‒spring rice in the fl ood‒prone area fl uctuated depending on the annual change of fl ood scale. It was also confi rmed that the triple rice‒cropped area in the An Giang province expanded from 2000 to 2005, because the construction of a ring‒dike system and water‒resource infrastructure made it possible to sustain a third rice cropping season during the fl ood season. The proposed methodologies(WFCP, WFCS, WFFI) based on MODIS time‒series imagery made it clear that while the rice cropping in the Vietnamese Mekong Delta was quantitatively(annual fl ooding) and qualitatively(salinity intrusion) affected by water‒resource changes, there were some regions where the cultivation system was changed from double rice cropping to triple rice cropping because of the implementation of measures against fl ooding.日本の食料自給率 (2005年時の供給熱量ベース) は、40% と先進7カ国の中で最も低い。日本は、その食料海外依存度の高さから、世界的な食料価格の変動の影響を最も受け易い国と言える。近年の経済発展に伴う中国の大豆輸入量の増加や世界的なエネルギー政策の転換 (バイオエタノール政策) は、世界の穀物需給バランスを不安定にさせつつあり、世界的な問題となっている。さらに、地球温暖化による農業生産影響、増加し続ける世界人口、鈍化する穀物生産性を考えれば、世界の食料需給バランスが将来にわたって安定し続けると言うことはできないだろう。他方、食料増産・生産性向上を目的とした集約的農業の展開は、発展途上国の農業環境にさらなる負荷を与えるかもしれない。世界の食料生産と密接な関係にある日本は、自国の食料安全保障を議論する前提として、急速に変わり行く世界の農業生産現場やそれを取り巻く農業環境を客観的に理解し、世界の農業環境情報を独自の手法によって収集・整理する必要がある。そこで、筆者は、衛星リモートセンシング技術を活用することによって、地球規模の視点で、時間的・空間的な広がりを持って変わり行く農業生産活動とそれを取り巻く農業環境情報を把握・理解するための時系列衛星データ解析手法の確立を目指すこととした。本研究では、インドシナ半島南端に位置するベトナム・メコンデルタを調査対象領域とした。ベトナムは、タイに次ぐ世界第2位のコメ輸出国であり、その輸出米の9割近くが、ベトナム・メコンデルタで生産されたものである。筆者は、ベトナム・メコンデルタを世界の食料安全保障を考える上で重要な食料生産地帯の一つであると考え、本地域における農業環境及び土地利用パターンの時空間変化を明らかにするためのMODIS データを用いた新たな時系列解析手法の開発を行った。　本研究において提案する時系列解析手法は、次の三つである。1. Wavelet‒based Filter for Crop Phenology (WFCP) ,2. Wavelet‒based Filter for evaluating the spatial distribution of Cropping System (WFCS) , 3. Wavelet‒based Filterfor detecting spatio‒temporal changes in Flood Inundation (WFFI) . WFCP は、時系列植生指数 (EVI) を平滑化するためにウェーブレット変換手法を利用しており、日本の農業統計データを用いた検証結果から、水稲生育ステージ (田植日、出穂日、収獲日) を約9－12日 (RMSE) の精度で推定可能であることが示された。WFCP を基に改良されたWFCS は、水稲作付パターンの年次把握を可能にし、ベトナムメコンデルタにおける水稲作付時期の空間分布が、上流部において毎年雨期に発生する洪水と沿岸部において乾季に発生する塩水遡上によって特徴づけられていることを明らかにした。WFFI は、時系列水指数 (LSWI) と植生指数 (EVI) から、湛水期間、湛水開始日・湛水終息日を広域把握し、メコン川洪水強度の年次変化を地域スケールで評価することを可能にする。そして、ウェーブレット変換を利用した一連の手法を、2000～2006年までのMODIS 時系列画像に適用することによって、メコンデルタ上流部の洪水常襲地帯において、冬春米の作付時期が、年次変化する洪水規模に依存していることを明らかにした。また、An Giang 省において、堤防建設 (輪中) や水利施設の建設によって、洪水期における水稲三期作が可能になった地域が、2000～2005年にかけて拡大していることを明らかにした。本研究で提案したMODIS 時系列画像を利用した時系列解析手法 (WFCP、WFCS、WFFI) によって、ベトナムメコンデルタにおける水稲生産が水資源の量的 (洪水) ・質的 (塩水遡上) 変動影響を受ける一方、洪水対策の実施によって、栽培体系を二期作から三期作に変更している地域があることを明らかにした

Mapping of multitemporal rice (Oryza sativa L.) growth stages using remote sensing with multi-sensor and machine learning : a thesis dissertation presented in partial fulfilment of the requirements for the degree of Doctor of Philosophy in Earth Science at Massey University, Manawatū, New Zealand

Figure 2.1 is adapted and re-used under a Creative Commons Attribution 4.0 International (CC BY 4.0) license.Rice (Oryza Sativa) plays a pivotal role in food security for Asian countries, especially in Indonesia. Due to the increasing pressure of environmental changes, such as land use and climate, rice cultivation areas need to be monitored regularly and spatially to ensure sustainable rice production. Moreover, timely information of rice growth stages (RGS) can lead to more efficient of inputs distribution from water, seed, fertilizer, and pesticide. One of the efficient solutions for regularly mapping the rice crop is using Earth observation satellites. Moreover, the increasing availability of open access satellite images such as Landsat-8, Sentinel-1, and Sentinel-2 provides ample opportunities to map continuous and high-resolution rice growth stages with greater accuracy. The majority of the literature has focused on mapping rice area, cropping patterns and relied mainly on the phenology of vegetation. However, the mapping process of RGS was difficult to assess the accuracy, time-consuming, and depended on only one sensor. In this work, we discuss the use of machine learning algorithms (MLA) for mapping paddy RGS with multiple remote sensing data in near-real-time. The study area was Java Island, which is the primary rice producer in Indonesia. This study has investigated: (1) the mapping of RGS using Landsat-8 imagery and different MLAs, and their rigorous performance was evaluated by conducting a multitemporal analysis; (2) the temporal consistency of predicting RGS using Sentinel-2, MOD13Q1, and Sentinel-1 data; (3) evaluating the correlation of local statistics data and paddy RGS using Sentinel-2, PROBA-V, and Sentinel-1 with MLAs. The ground truth datasets were collected from multi-year web camera data (2014-2016) and three months of the field campaign in different regions of Java (2018). The study considered the RGS in the analysis to be vegetative, reproductive, ripening, bare land, and flooding, and MLAs such as support vector machines (SVMs), random forest (RF), and artificial neural network (ANN) were used. The temporal consistency matrix was used to compare the classification maps within three sensor datasets (Landsat-8 OLI, Sentinel-2, and Sentinel-2, MOD13Q1, Sentinel-1) and in four periods (5, 10, 15, 16 days). Moreover, the result of the RGS map was also compared with monthly data from local statistics within each sub-district using cross-correlation analysis. The result from the analysis shows that SVM with a radial base function outperformed the RF and ANN and proved to be a robust method for small-size datasets (< 1,000 points). Compared to Sentinel-2, Landsat-8 OLI gives less accuracy due to the lack of a red-edge band and larger pixel size (30 x 30 m). Integration of Sentinel-2, MOD13Q1, and Sentinel-1 improved the classification performance and increased the temporal availability of cloud-free maps. The integration of PROBA-V and Sentinel-1 improved the classification accuracy from the Landsat-8 result, consistent with the monthly rice planting area statistics at the sub-district level. The western area of Java has the highest accuracy and consistency since the cropping pattern only relied on rice cultivation. In contrast, less accuracy was noticed in the eastern area because of upland rice cultivation due to limited irrigation facilities and mixed cropping. In addition, the cultivation of shallots to the north of Nganjuk Regency interferes with the model predictions because the cultivation of shallots resembles the vegetative phase due to the water banks. One future research idea is the auto-detection of the cropping index in the complex landscape to be able to use it for mapping RGS on a global scale. Detection of the rice area and RGS using Google Earth Engine (GEE) can be an action plan to disseminate the information quickly on a planetary scale. Our results show that the multitemporal Sentinel-1 combined with RF can detect rice areas with high accuracy (>91%). Similarly, accurate RGS maps can be detected by integrating multiple remote sensing (Sentinel-2, Landsat-8 OLI, and MOD13Q1) data with acceptable accuracy (76.4%), with high temporal frequency and lower cloud interference (every 16 days). Overall, this study shows that remote sensing combined with the machine learning methodology can deliver information on RGS in a timely fashion, which is easy to scale up and consistent both in time and space and matches the local statistics. This thesis is also in line with the existing rice monitoring projects such as Crop Monitor, Crop Watch, AMIS, and Sen4Agri to support disseminating information over a large area. To sum up, the proposed workflow and detailed map provide a more accurate method and information in near real-time for stakeholders, such as governmental agencies against the existing mapping method. This method can be introduced to provide accurate information to rice farmers promptly with sufficient inputs such as irrigation, seeds, and fertilisers for ensuring national food security from the shifting planting time due to climate change

PhenoRice:A method for automatic extraction of spatio-temporal information on rice crops using satellite data time series

Agricultural monitoring systems require spatio-temporal information on widely cultivated staple crops like rice. More emphasis has been made on area estimation and crop detection than on the temporal aspects of crop cultivation, but seasonal and temporal information such as i) crop duration, ii) date of crop establishment and iii) cropping intensity are as important as area for understanding crop production. Rice cropping systems are diverse because genetic, environmental and management factors (G × E × M combinations) influence the spatio-temporal patterns of cultivation. We present a rule based algorithm called PhenoRice for automatic extraction of temporal information on the rice crop using moderate resolution hypertemporal optical imagery from MODIS. Performance of PhenoRice against spatially and temporally explicit reference information was tested in three diverse sites: rice-fallow (Italy), rice-other crop (India) and rice-rice (Philippines) systems. Regional product accuracy assessments showed that PhenoRice made a conservative, spatially representative and robust detection of rice cultivation in all sites (r2 between 0.75 and 0.92) and crop establishment dates were in close agreement with the reference data (r2 = 0.98, Mean Error = 4.07 days, Mean Absolute Error = 9.95 days, p < 0.01). Variability in algorithm performance in different conditions in each site (irrigated vs rainfed, direct seeding vs transplanting, fragmented vs clustered rice landscapes and the impact of cloud contamination) was analysed and discussed. Analysis of the maps revealed that cropping intensity and season length per site matched well with local information on agro-practices and cultivated varieties. The results show that PhenoRice is robust for deriving essential temporal descriptions of rice systems in both temperate and tropical regions at a level of spatial and temporal detail that is suitable for regional crop monitoring on a seasonal basis

Monitoring Rice Cropping Pattern and Fallows in Central and Western Part of India

India has the largest area under rice cultivation and holds the second position all over the world as it is one of the principal food crops. Rice-fallow croplands areas are those areas where rice is grown during the Kharif growing season (June- October) followed by fallow during Rabi season (November-February). These croplands are not suitable to grow in Rabi season rice due to their high water needs, but are suitable for short season (≤ 3months). According to national statistics there is an increase in the rice areas in Central and Western states of India. The goal of this project is to monitor the rice-fallow cropland areas & mapping the expansion of rice areas. This study is conducted in Central and Western states of India where different rice eco-systems exist. Time series Moderate Resolution Imaging Spectroradiometer (MODIS) 16days Normalized Difference Vegetation Index (NDVI) at 250m spatial resolution and season wise intensive ground survey data was used. We have applied hierarchical classification and Spectral Matching Techniques (SMT) to map rice areas and the fallows there after (rabi-fallows), in Central and Western states of India. And change detection was carried during 2000-2015 and 2010-2015. The resultant rice maps are compared with available national and sub-national level statistics

Assessment of MODIS spectral indices for determining rice paddy agricultural practices and hydroperiod

The aims of this study were to assess the dynamics of the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI(1) and NDWI(2)) and Shortwave Angle Slope Index (SASI) in relation to rice agricultural practices and hydroperiod, and (2) to assess the capability for these indices to detect phenometrics in rice under different flooding regimes

Spectral Characteristics and Mapping of Rice Fields using Multi-Temporal Landsat and MODIS Data: A Case of District Narowal

Availability of remote sensed data provides powerful access to the spatial and temporal information of the earth surface Real-time earth observation data acquired during a cropping season can assist in assessing crop growth and development performance As remote sensed data is generally available at large scale rather than at field-plot level use of this information would help to improve crop management at broad-scale Utilizing the Landsat TM ETM ISODATA clustering algorithm and MODIS Terra the normalized difference vegetation index NDVI and enhanced vegetation index EVI datasets allowed the capturing of relevant rice cropping differences In this study we tried to analyze the MODIS Terra EVI NDVI February 2000 to February 2013 datasets for rice fractional yield estimation in Narowal Punjab province of Pakistan For large scale applications time integrated series of EVI NDVI 250-m spatial resolution offer a practical approach to measure crop production as they relate to the overall plant vigor and photosynthetic activity during the growing season The required data preparation for the integration of MODIS data into GIS is described with a focus on the projection from the MODIS Sinusoidal to the national coordinate systems However its low spatial resolution has been an impediment to researchers pursuing more accurate classification results and will support environmental planning to develop sustainable land-use practices These results have important implications for parameterization of land surface process models using biophysical variables estimated from remotely sensed data and assist for forthcoming rice fractional yield assessmen

APPLICATIONS OF MODERATE-RESOLUTION REMOTE SENSING TECHNOLOGIES FOR SURFACE AIR POLLUTION MONITORING IN SOUTHEAST ASIA

Retrievals from Earth observation satellites are widely used for many applications, including analyzing dynamic lands and measuring atmospheric components. This research aims to evaluate appropriateness of using satellite retrievals to facilitate understanding characteristics of Southeast Asian (SEA) surface air pollution, attributed to regional biomass burnings and urban activities. The studies in this dissertation focused on using satellite retrievals for 1) mapping potential SEA air pollution sources; which are forests, rice paddies, and urban areas, 2) understanding dynamic optical characteristics of SEA biomass-burning aerosols, and 3) inferring surface ozone level. Data used in this study were from three NASA\u27s Earth Observing System (EOS) satellites, which are Terra, Aqua, and Aura. These retrievals have spatial resolution ranging from hundred meters to ten kilometers. Algorithms used for the SEA land cover classification were developed using time-series analyses of surface reflectance in multiple wavelength bands from Moderate Resolution Imaging Spectroradiometer (MODIS) on the Terra satellite. Comparing the results to national statistical databases, good agreement was obtained for spatial estimation of forest areas after correction with plantation areas. For estimation of rice paddies areas, the agreement depended on the rice ecosystems. It was good for rainfed rice and poor for deepwater rice. Models for irrigated and upland rice areas showed overall high coefficients of determination, suggesting that they effectively simulated the spatial distribution of those rice paddies; but were prone to underestimate and overestimate, respectively. Estimated SEA regional rice area was 42×106 ha, which agrees with previous published values. Analysis of the satellite retrieval could identify large urban areas. However, the satellite-derived urban areas also incorrectly included large sandy beaches. Optical properties of SEA background aerosols were investigated through the multivariate analyses of long-term ground-based aerosol measurements acquired from Aerosol Robotic Network (AERONET). The results in this study showed that from mid-September to December, the aerosol had both fine size and high light scattering efficiency. It was assumed to be largely urban/industrial aerosols, possibly coming from eastern China. From January to April, the aerosol had fine size and had single scattering albedo (SSA at 440 nm) of approximately 0.9. It was assumed to be smoke from local biomass burning. From October to January, when seasonal winds are strongest, more SEA urban aerosol was observed. This aerosol had coarser size and had SSA of ~0.9 or less. The appropriateness of using Ozone Monitoring Instrument (OMI) aerosol retrieval to facilitate understanding SEA biomass-burning aerosol properties was evaluated through three lines of evidence. These are 1) comparisons between the results obtained from multivariate analyses of the OMI aerosol retrieval and those obtained from the ground-measured AERONET data, 2) from Atmospheric Infrared Sounder (AIRS) total column CO product, and 3) from MODIS active fire detections. The results showed that the OMI retrieval used for large-scale SEA biomass-burning aerosol characterization was consistent with these alternative measures only when 1 \u3c OMI aerosol optical depth (442 nm) \u3c 3. The OMI aerosol retrieval was then used for the study on dynamic characteristics of biomass burning aerosol. This study considered the aerosols from two forest-fire episodes, 2007 SEA continent and 2008 Indonesian fires. Dependence of the aerosol optical properties on four variables was investigated. These variables were 1) wind speed/direction, 2) relative humidity (RH), 3) land use/cover as a surrogate of fuel type estimated from time-series analysis of MODIS surface reflectance, and 4) age of aerosol estimated from spatial-temporal analysis of MODIS active fire and the wind characteristics. Results from Pearson Chi-square test for independence showed that the dependence between aerosol group memberships with different optical properties and the limiting variables was significant for most cases, except for Indonesian aerosol age factor. These results agree with prior knowledge on regional burning conditions (types of fuel and relative humidity) and aerosol chemical/physical properties (chemical composition related to aerosol optical properties and hygroscopicity). Using EOS-Aura tropospheric column ozone (TCO) to infer surface ozone level was evaluated through analyses of linear relationships between TCO estimated from OMI and Microwave Limb Sounder (MLS) retrievals and coincident TCO from balloon-based ozonesonde measurements. This evaluation was for different tropospheric ozone profile shapes and for different geographical regions (for low, mid, and high latitudes and for Pacific and Atlantic regions). Results indicate that inference on ozone level derived from the satellite-based TCO requires corresponding information about tropospheric ozone profile shape. The use of satellite-based TCO was more appropriate for polluted low-latitude locations where upper troposphere ozone is rare and surface enhanced ozone is high

Paddy field classification with MODIS-terra multi-temporal image transformation using phenological approach in Java Island

This paper presents the paddy field classification model using the approach based on periodic plant life cycle events and how these elevations in climate as well as habitat factors, such as elevation. The data used are MODIS-Terra two tiles of H28v09 and H29v09 of 2016, consist of 46 series of 8-daily data, with 500 meter resolution in Java region. The paddy field classification method based on the phenological model is done by Maximum Likelihood on the transformed annual multi-temporal image of the reflectance data, index data, and the combination of reflectance and index data. The results of the study showed that, with the reference of the Paddy Field Map from the Ministry of Agriculture (MoA), the overall accuracies of the paddy field classification results using the combination of reflectance and index data provide the highest (85.4%) among the reflectance data (83.5%) and index data (81.7%). The accuracy levels were varied; these depend on the slope and the types of paddy fields. Paddy fields on the slopes of 0-2% could be well identified by MODIS-Terra data, whereas it was difficult to identify the paddy fields on the slope &gt;2%. Rain-fed lowland paddy field type has a lower user accuracy than irrigated paddy fields. This study also performed correlation (r2) between the analysis results and the statistical data based on district and provincial boundaries were &gt;0.85 and &gt;0.99 respectively. These correlations were much higher than the previous study results, which reached 0.49-0.65 (hilly-flat areas of county-level), and 0.80-0.88 (hilly-flat areas of provincial level) for China, and reached 0.44 for Indonesia

Automated cropping intensity extraction from isolines of wavelet spectra

Timely and accurate monitoring of cropping intensity (CI) is essential to help us understand changes in food production. This paper aims to develop an automatic Cropping Intensity extraction method based on the Isolines of Wavelet Spectra (CIIWS) with consideration of intra- class variability. The CIIWS method involves the following procedures: (1) characterizing vegetation dynamics from time–frequency dimensions through a continuous wavelet transform performed on vegetation index temporal profiles; (2) deriving three main features, the skeleton width, maximum number of strong brightness centers and the intersection of their scale intervals, through computing a series of wavelet isolines from the wavelet spectra; and (3) developing an automatic cropping intensity classifier based on these three features. The proposed CIIWS method improves the understanding in the spectral–temporal properties of vegetation dynamic processes. To test its efficiency, the CIIWS method is applied to China’s Henan province using 250 m 8 days composite Moderate Resolution Imaging Spectroradiometer (MODIS) Enhanced Vegetation Index (EVI) time series datasets. An overall accuracy of 88.9% is achieved when compared with in-situ observation data. The mapping result is also evaluated with 30 m Chinese Environmental Disaster Reduction Satellite (HJ-1)-derived data and an overall accuracy of 86.7% is obtained. At county level, the MODIS-derived sown areas and agricultural statistical data are well correlated (r2 = 0.85). The merit and uniqueness of the CIIWS method is the ability to cope with the complex intra-class variability through continuous wavelet transform and efficient feature extraction based on wavelet isolines. As an objective and meaningful algorithm, it guarantees easy applications and greatly contributes to satellite observations of vegetation dynamics and food security efforts