Local-scale agricultural drought monitoring with satellite-based multi-sensor time-series

ABSTRACT

Globally, drought constitutes a serious threat to food and water security. The complexity and multivariate nature of drought challenges its assessment, especially at local scales. The study aimed to assess spatiotemporal patterns of crop condition and drought impact at the spatial scale of field management units with a combined use of time-series from optical (Landsat, MODIS, Sentinel-2) and Synthetic Aperture Radar (SAR) (Sentinel 1) data. Several indicators were derived such as Normalized Difference Vegetation Index (NDVI), Normalized Difference Moisture Index (NDMI), Land Surface Temperature (LST), Tasseled cap indices and Sentinel-1 based backscattering intensity and relative surface moisture. We used logistic regression to evaluate the drought-induced variability of remotely sensed parameters estimated for different phases of crop growth. The parameters with the highest prediction rate were further used to estimate thresholds for drought/non-drought classification. The models were evaluated using the area under the receiver operating characteristic curve and validated with in-situ data. The results revealed that not all remotely sensed variables respond in the same manner to drought conditions. Growing season maximum NDVI and NDMI (70–75%) and SAR derived metrics (60%) reflect specifically the impact of agricultural drought. These metrics also depict stress affected areas with a larger spatial extent. LST was a useful indicator of crop condition especially for maize and sunflower with prediction rates of 86% and 71%, respectively. The developed approach can be further used to assess crop condition and to support decision-making in areas which are more susceptible and vulnerable to drought.     

基于MODIS数据的环北京地区土地资源监测研究

本文基于MODIS 16天合成的NDVI时间序列数据及其他辅助数据,首先用PCA方法对NDVI时间序列数据进行信息增强与压缩处理,结合LST数据、DEM数据及降雨温度数据,利用模糊K-均值非监督分类法,进行环北京地区的土地覆盖分类,得到土地资源现状情况。然后利用变化矢量(CVA)分析方法对环北京地区的土地利用及植被覆盖的多年变化状况进行了分析。结果表明,MODIS数据能很好的应用于大范围的土地资源监测中,并能得到较好的结果。     

Crop classification and acreage estimation in North Korea using phenology features

In North Korea, reliable and timely information on crop acreage and spatial distribution is hard to obtain. In this study, we developed a fast and robust method to estimate crop acreage in North Korea using time-series normalized difference vegetation index (NDVI) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) data. We proposed a method to identify crop type based on NDVI phenology features using data collected in other areas with similar agri-environmental conditions to mitigate the shortage of ground truth data. Eventually the classification map (MODIScrop) was assessed using the Food and Agriculture Organization (FAO) statistical data and high-resolution crop classification maps derived from one Landsat scene (LScrop). The Pareto boundary method was used to assess the accuracy and crop distribution of the MODIScrop maps. Results showed that acreage derived from the MODIScrop maps was generally consistent with that reported in the FAO data (a relative error <4.1% for rice and <6.1% for maize, and <9.0% for soybean except for in 2004, 2008, and 2009) and the maps derived from the LScrop (a relative error about 5% in 2013, and 7% in 2008 and 2014). The classification accuracy reached 74.4%, 69.8%, and 73.1% of the areas covered by the Landsat images in 2008, 2013, and 2014, respectively. This indicates that features derived from NDVI profiles were able to characterize major crops, and the approaches developed in this study are feasible for crop mapping and acreage estimation in regions with limited ground truth data.     

Land surface temperature retrieval and its validation using NOAA AVHRR thermal data

The retrieval of land (soil-vegetation complex) surface temperature (LST) was carried out over semi-arid mixed agriculture landscape of Gujarat using thermal bands (channel 4 and 5) and ground emissivity from atmospherically corrected NDVI of NOAA AVHRR LAC images. The atmospheric correction of Visible and NIR band reflectance was done using SMAC model. The LST computed from split-window method and subsequently corrected with fractional vegetation cover were then compared with near synchronous ground observations of soil and air temperatures made during 13–17 January and April, 1997 at five Land Surface Processes Experiment (LASPEX) sites of Anand, Sanand, Derol, Arnej and Khandha covering 100 km x 100 km. The fractional vegetation cover corrected LST at noon hrs. varied from 301.6 – 311.9K in January and from 315.8 – 325.6K in April. The LST_corr were found to lie in the mid way between AT and ST during January. But in April, LST were found to be more close to ST which may be due to relatively poor vegetation growth as indicated by lower NDVI values in April indicating more contribution to LST from exposed soil surface.     

Investigation and validation of algorithms for estimating land surface temperature from Sentinel-3 SLSTR data

Land surface temperature (LST) is an important indicator of global ecological environment and climate change. The Sea and Land Surface Temperature Radiometer (SLSTR) onboard the recently launched Sentinel-3 satellites provides high-quality observations for estimating global LST. The algorithm of the official SLSTR LST product is a split-window algorithm (SWA) that implicitly assumes and utilizes knowledge of land surface emissivity (LSE). The main objective of this study is to investigate alternative SLSTR LST retrieval algorithms with an explicit use of LSE. Seventeen widely accepted SWAs, which explicitly utilize LSE, were selected as candidate algorithms. First, the SWAs were trained using a comprehensive global simulation dataset. Then, using simulation data as well as in-situ LST, the SWAs were evaluated according to their sensitivity and accuracy: eleven algorithms showed good training accuracy and nine of them exhibited low sensitivity to uncertainties in LSE and column water vapor content. Evaluation based on two global simulation datasets and a regional simulation dataset showed that these nine SWAs had similar accuracy with negligible systematic errors and RMSEs lower than 1.0 K. Validation based on in-situ LST obtained for six sites further confirmed the similar accuracies of the SWAs, with the lowest RMSE ranges of 1.57–1.62 K and 0.49−0.61 K for Gobabeb and Lake Constance, respectively. While the best two SWAs usually yielded good accuracy, the official SLSTR LST generally had lower accuracy. The SWAs identified and described in this study may serve as alternative algorithms for retrieving LST products from SLSTR data.     

Methods for improving accuracy and extending results beyond periods covered by traditional ground-truth in remote sensing classification of a complex landscape

Successful development of approaches to quantify impacts of diverse landuse and associated agricultural management practices on ecosystem services is frequently limited by lack of historical and contemporary landuse data. We hypothesized that ground truth data from one year could be used to extrapolate previous or future landuse in a complex landscape where cropping systems do not generally change greatly from year to year because the majority of crops are established perennials or the same annual crops grown on the same fields over multiple years. Prior to testing this hypothesis, it was first necessary to classify 57 major landuses in the Willamette Valley of western Oregon from 2005 to 2011 using normal same year ground-truth, elaborating on previously published work and traditional sources such as Cropland Data Layers (CDL) to more fully include minor crops grown in the region. Available remote sensing data included Landsat, MODIS 16-day composites, and National Aerial Imagery Program (NAIP) imagery, all of which were resampled to a common 30 m resolution. The frequent presence of clouds and Landsat7 scan line gaps forced us to conduct of series of separate classifications in each year, which were then merged by choosing whichever classification used the highest number of cloud- and gap-free bands at any given pixel. Procedures adopted to improve accuracy beyond that achieved by maximum likelihood pixel classification included majority-rule reclassification of pixels within 91,442 Common Land Unit (CLU) polygons, smoothing and aggregation of areas outside the CLU polygons, and majority-rule reclassification over time of forest and urban development areas. Final classifications in all seven years separated annually disturbed agriculture, established perennial crops, forest, and urban development from each other at 90 to 95% overall 4-class validation accuracy. In the most successful use of subsequent year ground-truth data to classify prior year landuse, an overall 57-class accuracy of 75% was achieved despite the omission of 10 entire classes, most of which were annually disturbed or perennial crops grown on very few fields. Synthetic ground-truth data for the 2004 harvest year based on the most common landuse classes over the following 7 years classified 49 of 57 categories at an overall accuracy of 96% in a final version that included CLU polygon majority rule, default smoothing and aggregation, and forcing of urban development and forest from multi-year majority-rule.     

Mapping of Cropping System for the Indo-Gangetic Plain Using Multi-Date SPOT NDVI-VGT Data

The present study has been carried out to delineate the existing cropping systems in the Indo-Gangetic Plains (IGP) using 10 day composite SPOT VEGETATION (VGT) NDVI data acquired over a crop year (June–May). Results showed that it is feasible to identify the major crops like rice, wheat, sugarcane, potato, and cotton in the dominant growing areas with good accuracy. Double cropping pattern is the most prevalent. Rice-wheat, sugarcane based, cotton-wheat, rice-potato, rice-rice, maize/millet-wheat are some of the major rotations followed. Rice-wheat is the dominant rotation accounting for around 40% of the net sown area. Triple crop rotations was less than 5% of the area and observed in some parts of Uttar Pradesh, Bihar and West Bengal. Single crop rotation of rice-fallow is significant only in West Bengal.     

Land surface temperature and its impact factors in Western Sichuan Plateau,China

The understanding influence of multiple factors variations on land surface temperature (LST) remains elusive. LST was retrieved by the atmospheric correction algorithms. Based on the correlation coefficients, stepwise regression analysis was developed to examine how multiple factors variability led to LST variations. The differences in LST between impact factors vary depending on time in a day. The elevation and land use types significantly affect the LST in sunny slope or shadow areas has a significantly quadratic curve correlation or a negative linear correlation with it, the influence of slope and aspect is not very significant. LST for forestland, grassland and bare land in the sunny slope and shadow area was the cubic polynomial related to its elevation. Normalized difference vegetation index (NDVI) and normalized difference moisture index (NDMI) effectively express LST in mountainous. LST and NDMI or NDVI have a significantly negative correlation, NDMI is more effective and more applicable for the expression of LST.     

Crop area mapping in West Africa using landscape stratification of MODIS time series and comparison with existing global land products

In Africa, food security early warning systems use satellite-derived data concerning crop conditions and agricultural production. Such systems can be improved if they are provided with a more reliable estimation of the cultivated area at national scale. This paper evaluates the potential of using time series from the MODerate resolution Imaging Spectroradiometer MOD13Q1 (16-day composite of normalized difference vegetation index at 250 m resolution) to extract cultivated areas in the fragmented rural landscapes of Mali. To this end, we first stratified Southern Mali into 13 rural landscapes based on the spatio-temporal variability of NDVI and textural indices, using an object-oriented classification scheme.The accuracy of the resulting map (MODIS_crop) and how it compares with existing coarse-resolution global land products (GLC2000 Africa, GLOBCOVER, MODIS V05 and ECOCLIMAP-II), was then assessed against six crop/non-crop maps derived from SPOT 2.5 m resolution images used as references. For crop areal coverage, the MODIS_crop cultivated map was successful in assessing the overall cultivated area at five out of the six validation sites (less than 6% of the absolute difference), while in terms of crop spatial distribution, the producer accuracy was between 33.1% and 80.8%. This accuracy was linearly correlated with the mean patch size index calculated on the SPOT crop maps (r² = 0.8). Using the Pareto boundary as an accuracy assessment method at the study sites, we showed that (i) 20-40% of the classification crop error was due to the spatial resolution of the MODIS sensor (250 m), and that (ii) compared to MODIS V05, which otherwise performed better than the other existing products, MODIS_crop generally minimized omission-commission errors. A spatial validation of the different products was carried out using SPOT image classifications as reference. In the corresponding error matrices, the fraction of correctly classified pixels for our product was 70%, compared to 58% for MODIS V05, while it ranged between 40% and 51% for the GLC2000, the ECOCLIMAP-II and the GLOBCOVER.     

Effectiveness of SID as Spectral Similarity Measure to Develop Crop Spectra from Hyperspectral Image

The present study was undertaken with the objective to check effectiveness of spectral information divergence (SID) to develop spectra from image for crop classes based on spectral similarity with field spectra. In multispectral and hyperspectral remote sensing, classification of pixels is obtained by statistical comparison (by means of spectral similarity) of known field or library spectra to unknown image spectra. Though these algorithms are readily used, little emphasis has been placed on use of various spectral similarity measures to develop crop spectra from the image itself. Hence, in this study methodology suggested to develop spectra for crops based on SID. Absorption features are unique and distinct; hence, validation of the developed spectra is carried out using absorption features by comparing it with field spectra and finding average correlation coefficient r?=?0.982 and computed SID equivalent r?=?0.989. Effectiveness of developed spectra for image classification was computed by probability of spectral discrimination (PSD) and resulted in higher probability for the spectra developed based on SID. Image classification was carried out using field spectra and spectra assigned by SID. Overall classification accuracy of the image classified by field spectra is 78.30% and for the image classified by spectra assigned through SID-based approach is 91.82%. Z test shows that image classification carried out using spectra developed by SID is better than classification carried out using field spectra and significantly different. Validation by absorption features, effectiveness by PSD and higher classification accuracy show possibility of new approach for spectra development based on SID spectral similarity measure.     

Detection and estimation of mixed paddy rice cropping patterns with MODIS data

In this paper, we developed a more sophisticated method for detection and estimation of mixed paddy rice agriculture from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. Previous research demonstrated that MODIS data can be used to map paddy rice fields and to distinguish rice from other crops at large, continental scales with combined Enhanced Vegetation Index (EVI) and Land Surface Water Index (LSWI) analysis during the flooding and rice transplanting stage. Our approach improves upon this methodology by incorporating mixed rice cropping patterns that include single-season rice crops, early-season rice, and late-season rice cropping systems. A variable EVI/LSWI threshold function, calibrated to more local rice management practices, was used to recognize rice fields at the flooding stage. We developed our approach with MODIS data in Hunan Province, China, an area with significant flooded paddy rice agriculture and mixed rice cropping patterns. We further mapped the aerial coverage and distribution of early, late, and single paddy rice crops for several years from 2000 to 2007 in order to quantify temporal trends in rice crop coverage, growth and management systems. Our results were validated with finer resolution (2.5 m) Satellite Pour l’Observation de la Terre 5 High Resolution Geometric (SPOT 5 HRG) data, land-use data at the scale of 1/10,000 and with county-level rice area statistical data. The results showed that all three paddy rice crop patterns could be discriminated and their spatial distribution quantified. We show the area of single crop rice to have increased annually and almost doubling in extent from 2000 to 2007, with simultaneous, but unique declines in the extent of early and late paddy rice. These results were significantly positive correlated and consistent with agricultural statistical data at the county level (P < 0.01).     

Remote sensing research issues of the National Land Use Change Program of China

In 1999, the Ministry of Land and Resources (MLR) of China launched the National Land Use Change Program especially to monitor the scale and distribution of urban expansion and the decrease in cultivated land through remote sensing technology. This Program has been carried out annually and continuously for seven years since then and played an important role in the policy-making of MLR about land management and planning. This paper gives an overview about this Program and discusses several research issues. First, the remote sensing data sources and other ancillary data used in this Program are presented. The approaches for image preprocessing, i.e. radiometric normalization, image geometric rectification and image fusion are then introduced with an emphasis on the algorithm development for image registration. Second, land use change detection technique is the most critical and complex aspect of the Program. The methodologies for change detection using either bi-temporal image pair or one existing land use map and one remotely sensed image are detailed. Third, since the data of land use changes derived from remote sensing will be operationally used for local and central government, field validation and accuracy assessment are crucial to ensure the reliability of change detection results. The strategy of field work and the resulting accuracy evaluations is presented. The land use and change information derived from remotely sensed data has wide applications for land management, including land use database updating, verification of land use planning and monitoring of national high-tech parks. Last, suggestions on how to make full use of the images and change detection result, to improve the consistency of land use classification and to develop change detection algorithms for diverse and complex remote sensing data are given.     

Land cover classification and economic assessment of citrus groves using remote sensing

The citrus industry has the second largest impact on Florida's economy, following tourism. Estimation of citrus area coverage and annual forecasts of Florida's citrus production are currently dependent on labor-intensive interpretation of aerial photographs. Remotely sensed data from satellites has been widely applied in agricultural yield estimation and cropland management. Satellite data can potentially be obtained throughout the year, making it especially suitable for the detection of land cover change in agriculture and horticulture, plant health status, soil and moisture conditions, and effects of crop management practices. In this study, we analyzed land cover of citrus crops in Florida using Landsat Enhanced Thematic Mapper Plus (ETM+) imagery from the University of Maryland Global Land Cover Facility (GLCF). We hypothesized that an interdisciplinary approach combining citrus production (economic) data with citrus land cover area per county would yield a correlation between observable spectral reflectance throughout the year, and the fiscal impact of citrus on local economies. While the data from official sources based on aerial photography were positively correlated, there were serious discrepancies between agriculture census data and satellite-derived cropland area using medium-resolution satellite imagery. If these discrepancies can be resolved by using imagery of higher spatial resolution, a stronger correlation would be observed for citrus production based on satellite data. This would allow us to predict the economic impact of citrus from satellite-derived spectral data analysis to determine final crop harvests.     

卫星遥感地表温度降尺度的光谱归一化指数法

针对卫星遥感技术监测地表温度(land surface temperature,LST)存在时空分辨率矛盾这一难题,以TsHARP温度降尺度算法为基础,根据地表覆盖类型的不同,分别选择与LST相关性更好的光谱指数(归一化植被指数,NDVI;归一化建造指数,NDBI;改进的归一化水体指数,MNDWI;增强型裸土指数,EBSI)提出了新的转换模型,并从定性和定量两个角度评价了TsHARP法和新模型的降尺度精度。结果表明:两种模型在提高LST空间分辨率的同时又能较好地保持MODIS LST影像热特征的空间分布格局,消除了原始1km影像中的马赛克效应,两种模型均能够达到较好的降尺度效果;全局尺度分析表明,不管是在降尺度结果的空间变异性还是精度方面,本文提出的模型(RMSE:1.635℃)均要优于TsHARP法(RMSE:2.736℃);TsHARP法在水体、裸地和建筑用地这些低植被覆盖区表现出较差的降尺度结果,尤其对于裸地和建筑用地更为明显(|MBE|3℃),新模型提高了低植被覆盖区地物的降尺度精度;不同季节的降尺度结果表明,两种模型都是夏、秋季的降尺度结果优于春、冬季,新模型的降尺度结果四季均好于TsHARP法,其中春、冬季的降尺度精度提升效果要优于夏、秋季。     

基于多源遥感数据的玛纳斯河流域荒漠化监测

利用遥感手段对玛纳斯河流域绿洲-荒漠带的荒漠化状况进行监测。利用不同时期TM数据,在小尺度上分析该区域2000~2008年荒漠化土地的分布特征;基于MODIS陆地遥感产品(NDVI、Albedo及LST),利用多源信息复合法计算荒漠化程度指数;基于TM荒漠化监测结果对该指数进行分级处理,初步建立并验证基于MODIS数据的荒漠化监测评价模型。结果表明:近10 a来,绿洲外围荒漠化土地变化较小,大部分变化发生在绿洲内部,主要为局部荒漠化被开发或因盐碱化过重导致的撂荒;荒漠化监测评价模型对研究区荒漠化程度定量评价精度可达91.25%。     

Landsat8和MODIS融合构建高时空分辨率数据识别秋粮作物

本文利用Wu等人提出的遥感数据时空融合方法 STDFA(Spatial Temporal Data Fusion Approach)以Landsat 8和MODIS为数据源构建高时间、空间分辨率的遥感影像数据。以此为基础,构建15种30 m分辨率分类数据集,然后利用支持向量机SVM(Support Vector Machine)进行秋粮作物识别,验证不同维度分类数据集进行秋粮作物识别的适用性。实验结果显示,不同分类数据集的秋粮作物分类结果均达到了较高的识别精度。综合各项精度指标分析,Red+Phenology数据组合对秋粮识别效果最好,水稻识别的制图精度和用户精度分别达到91.76%和82.49%,玉米识别的制图精度和用户精度分别达到85.80%和74.97%,水稻和玉米识别的总体精度达到86.90%。     

A Comparative Analysis of Phenological Curves for Major Crops in Kansas

The goal of this research was to conduct an initial investigation into whether a time-series NDVI reference curve library for crops over a growing season for one year could be used to map crops for a different year. Time-series NDVI libraries of curves for 2001 and 2005 were investigated to ascertain whether or not the 2001 dataset could be used to map crops for 2005. The 2005 16-day composite MODIS 250 m NDVI data were used to extract NDVI values from 1,615 field sites representing alfalfa, corn, sorghum, soybeans, and winter wheat. A k-means cluster analysis of NDVI values from the field sites was performed to identify validation sites with time-series NDVI spectral profiles characteristic of the major crop types grown in Kansas. After completing the field site refinement process, there were 1,254 field sites retained for further analysis, referred to as "final" field sites. The methods employed to evaluate whether the MODIS-based NDVI profiles for major crops in Kansas are stable from year-to-year involved both graphical and statistical analyses. First, the time-series NDVI values for 2005 from the final field sites were aggregated by crop type and the crop NDVI profiles were then visually assessed and compared to the profiles of 2001 to ascertain if each crop's unique phenological pattern was consistent between the two years. Second, separability within each crop class in the time-series NDVI data between 2001 and 2005 was investigated numerically using the Jeffries-Matusita (JM) distance statistic. The results seem to suggest that time-series NDVI response curves for crops over a growing period for one year of valid ground reference data may be useful for mapping crops for a different year when minor temporal shifts in the NDVI values (resulting from inter-annual climate variations or changes in agricultural management practices) are taken into account.     

Landsat ETM＋数据的武汉市地表温度反演研究

以Landsat ETM＋为主要遥感数据源,辅助以气象数据资料,利用NDVI计算出植被覆盖率,继而计算出地表比辐射率,并利用单窗算法精确反演得到了武汉市2002年夏季地表温度。研究表明,武汉城市热岛由武昌、汉口、汉阳三大城区热岛构成,且几乎集中在城区;城区大热岛区中还存在多数热岛效应更为严重的＂岛中岛＂;武汉热岛分布与...     

基于Landsat数据地表特征参数与地表温度关系

以成都市为研究区,定量分析了各地表特征参数与地表温度之间的线性关系。通过对地表温度与归一化植被指数(NDVI)、归一化建筑指数(NDBI)、归一化水汽指数(NDMI)进行局部区域逐像元分析和总体区域统计分析,结果表明NDVI,NDMI,NDBI与地表温度间都存在明显的线性关系,可用于说明地表温度的动态变化,在3月份,NDMI与地温的相关性更优于NDVI。对传统城市热现象研究中,NDMI与NDBI能够用来以NDVI作为分析地表温度随季节而变化的互补的度量标准。     

Land cover classification of the North China Plain using MODIS_EVI time series

Moderate Resolution Imaging Spectroradiometer (MODIS) data have played an important role in global environmental and resource research. However, its low spatial resolution has been an impediment to researchers pursuing more accurate classification results. In this research, the high temporal resolution of MODIS was employed to improve the accuracy of land cover classification of the North China Plain using MODIS_EVI time series from 2003. Harmonic Analysis of Time Series (HANTS) was performed on the MODIS_EVI image time series to reduce cloud and other noise effects. The improved MODIS_EVI time series was then classified into 100 clusters by the Iterative Self-Organizing Data Analysis Technique (ISODATA). To distinguish ambiguous land cover classes, a decision tree was built on five phenological features derived from EVI profiles, Land Surface Temperature (LST) and topographic slope. The overall accuracy of the final land cover map was 75.5%, indicating the promise of using MODIS EVI time series and decision trees for broad area land cover classification.