Variation in NDVI values with change in spatial resolution for semi-arid savanna vegetation: a case study in northwestern South Africa

Natural vegetation and crop-greening patterns in semi-arid savannas are commonly monitored using normalized difference vegetation index (NDVI) values from low spatial resolution sensors such as the Advanced Very High Resolution Radiometer (AVHRR) (1 km, 4 km) and Moderate Resolution Imaging Spectroradiometer (MODIS) (250 m, 500 m). However, because semi-arid savannas characteristically have scattered tree cover, the NDVI values at low spatial resolution suffer from the effect of aggregation of near-infrared and red energy from adjacent vegetated and non-vegetated cover types. This effect is seldom taken into consideration or quantified in NDVI analyses of the vegetation of semi-arid lands. This study examined the effect of pixel size on NDVI values of land-cover features for a semi-arid area, using the 1000 m, 250 m and 10 m pixel sizes. A rainy season Système Pour l'Observation de la Terre 5 (SPOT 5) High Resolution Geometric (HRG) image at 10 m spatial resolution was utilized. Following radiometric and geometric preprocessing, the 10 m pixel size of the image was aggregated to 250 m and 1000 m to simulate imagery at these pixel sizes, and then NDVI images at the spatial resolution scales of 10 m (NDVI_{10 m}), 250 m (NDVI_{250 m}), and 1000 m (NDVI_{1000 m}) derived from the respective images. The simulation of the NDVI_{250 m} image was validated against a concurrent 16 day MODIS NDVI composite (MOD13Q1) image, and the accuracy derived from the validation was generalized to the NDVI_{1000 m} image. With change from low to high spatial resolution, extreme magnitude NDVI values shifted towards the centre (mode) of the resulting approximately Gaussian NDVI distributions. There was a statistically significant difference in NDVI values at the three pixel sizes. Low spatial magnitude vegetation sites (woodland, cropland) had reductions of up to 28% in NDVI value between the NDVI_{10 m} and NDVI_{1000 m} scales. The results indicate that vegetation monitoring using low spatial resolution imagery in semi-arid savannas may only be indicative and needs to be supplemented by higher spatial resolution imagery.     

Analysing NDVI for the African continent using the geostationary meteosat second generation SEVIRI sensor

This study presents first results on Normalized Difference Vegetation Index (NDVI), from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) sensor onboard the geostationary satellite Meteosat Second Generation (MSG) covering the African continent. With a temporal resolution of 15 min MSG offers complementary information for NDVI monitoring compared to vegetation monitoring based on polar orbiting satellites. The improved temporal resolution has potential implications for accurate NDVI assessment of the African continent; e.g. the increased amount of available scenes are expected to help overcome problems related to cloud cover which makes the MSG data particularly well suited for early warning systems. Time series of 2004 MSG NDVI was compared to MODIS (Moderate Resolution Imaging Spectroradiometer) Terra and Aqua NDVI for the Dahra site in the Senegalese Sahel, West Africa. It was found that NDVI was available for 82 days with multiple cloud free acquisitions per day during the growing season as compared to 47 days with information from either MODIS Terra or Aqua for that particular site. Differences in MSG SEVIRI and MODIS BRDF on a seasonal scale were found to influence the time series of NDVI for the test site; MSG NDVI being higher than MODIS in July-August and lower in October-November. Preliminary composite analysis suggests that the period of compositing to produce continent scale cloud free products can be reduced to ∼5 days using MSG NDVI as compared to polar orbiting data. With the availability of diurnal reflectance information the significance of differences between the red and near-infrared wavelengths due to anisotropy become evident, causing diurnal variations in observed NDVI. Diurnal MSG NDVI was compared to in situ measured MSG NDVI at the test site in Senegal and the same “bowl-shaped” diurnal curve was found for a medium dense cover of annual grasses. The range in observed NDVI and time of diurnal minimum was different due to different viewing geometry. Daily minimum of in situ measured NDVI was around solar noon whereas minimum MSG NDVI occurs one hour prior to noon due to the test site location 12° west of the satellite sensor. Diurnal variation in observed NDVI was studied for a number of pixels characterized by different sensor view zenith angles and vegetation types. This analysis illustrated the diurnal NDVI dependency of illumination conditions, view angle and vegetation intensity and pinpoints the importance of proper BRDF modeling to produce daily values of MSG NDVI normalized for acquisition time, which will be the subject of a forthcoming paper.     

Temporal filtering of successive MODIS data in monitoring a locust outbreak

The emergence of high temporal resolution satellite data such as MODIS enables timely monitoring of locust outbreaks from space. This monitoring is hampered by the effect of random atmospheric variations on satellite imagery, which may be suppressed through temporal filtering. This paper aims to evaluate the utility of temporally filtering successive MODIS data in monitoring an outbreak in East China. Of the eight vegetation indices examined, the commonly used NDVI was the most indicative of varying vegetation conditions caused by locust infestation inside the study area. The averaging of three successive days of satellite data improves the R ² value of NDVI regression models by 0.227 over single‐day data. It also outperforms the data averaged from two successive days (a broader window size was not attempted due to the short span of the study period). Temporally, NDVI changed at varying rates daily during the outbreak. Early in the outbreak it increased at a reduced pace until 7.5 days. Afterwards it started to decrease at an accelerated rate. If temporally filtered with a proper window size, successive MODIS data allow the outbreak to be monitored accurately (R ² = 0.696).     

Evaluation of earth observation based long term vegetation trends — Intercomparing NDVI time series trend analysis consistency of Sahel from AVHRR GIMMS, Terra MODIS and SPOT VGT data

AVHRR (Advanced Very High Resolution Radiometer) GIMMS (Global Inventory Modelling and Mapping Studies) NDVI (Normalized Difference vegetation Index) data is available from 1981 to present time. The global coverage 8 km resolution 15-day composite data set has been used for numerous local to global scale vegetation time series studies during recent years. Several aspects however potentially introduce noise in the NDVI data set due to the AVHRR sensor design and data processing. More recent NDVI data sets from both Terra MODIS and SPOT VGT data are considered an improvement over AVHRR and these products in theory provide a possibility to evaluate the accuracy of GIMMS NDVI time series trend analysis for the overlapping period of available data. In this study the accuracy of the GIMMS NDVI time series trend analysis is evaluated by comparison with the 1 km resolution Terra MODIS (MOD13A2) 16-day composite NDVI data, the SPOT Vegetation (VGT) 10-day composite (S10) NDVI data and in situ measurements of a test site in Dahra, Senegal. Linear least squares regression trend analysis on eight years of GIMMS annual average NDVI (2000-2007) has been compared to Terra MODIS (1 km and 8 km resampled) and SPOT VGT NDVI data 1 km (2000-2007). The three data products do not exhibit identical patterns of NDVI trends. SPOT VGT NDVI data are characterised by higher positive regression slopes over the 8-year period as compared to Terra MODIS and AVHRR GIMMS NDVI data, possibly caused by a change in channels 1 and 2 spectral response functions from SPOT VGT1 to SPOT VGT2 in 2003. Trend analysis of AVHRR GIMMS NDVI exhibits a regression slope range in better agreement with Terra MODIS NDVI for semi-arid areas. However, GIMMS NDVI shows a tendency towards higher positive regression slope values than Terra MODIS in more humid areas. Validation of the different NDVI data products against continuous in situ NDVI measurements for the period 2002-2007 in the semi-arid Senegal revealed a good agreement between in situ measurements and all satellite based NDVI products. Using Terra MODIS NDVI as a reference, it is concluded that AVHRR GIMMS coarse resolution NDVI data set is well-suited for long term vegetation studies of the Sahel-Sudanian areas receiving < 1000 mm rainfall, whereas interpretation of GIMMS NDVI trends in more humid areas of the Sudanian-Guinean zones should be done with certain reservations.     

Assessment of MODIS sun-sensor geometry variations effect on observed NDVI using MSG SEVIRI geostationary data

The quality of Earth observation (EO) based vegetation monitoring has improved during recent years, which can be attributed to the enhanced sensor design of new satellites such as MODIS (Moderate Resolution Imaging Spectroradiometer) on Terra and Aqua. It is however expected that sun-sensor geometry variations will have a more visible impact on the Normalized Difference Vegetation Index (NDVI) from MODIS compared to earlier data sources, since noise related to atmosphere and sensor calibration is substantially reduced in the MODIS data stream. For this reason, the effect of varying MODIS viewing geometry on red, near-infrared (NIR) and NDVI needs to be quantified. Data from the geostationary MSG (Meteosat Second Generation) SEVIRI (Spinning Enhanced Visible and Infrared Imager) sensor is well suited for this purpose due to the fixed position of the sensor, the spectral resolution, including a red and NIR band, and the high temporal resolution (15 min) of data, enabling MSG data to be used as a reference for estimating MODIS surface reflectance and NDVI variations caused by varying sun-sensor geometry. The study was performed on data covering West Africa for periods of lowest possible cloud cover for three consecutive years (2004–2006). An analysis covering the entire range of NDVI revealed day-to-day variations in observed MODIS NDVI of 50–60% for medium dense vegetation (NDVI ≈ 0.5) caused by variations in MODIS view zenith angles (VZAs) between nadir and the high forward-scatter view direction. Statistical analysis on red, NIR and NDVI from MODIS and MSG SEVIRI for three transects (characterized by different vegetation densities) showed that both MODIS red and NIR reflectances are highly dependant on MODIS VZA and relative azimuth angle (RAA), due to the anisotropic behaviour of red and NIR reflectances. The anisotropic reflectance in the red and NIR band was to some degree minimized by the ratioing properties of NDVI. The minimization by the NDVI normalization is very dependent on the vegetation density however, since the degree of anisotropy in red and NIR reflectances depends on the amount of vegetation present. MODIS VZA and RAA effects on NDVI were highest for medium dense vegetation (NDVI ≈ 0.5–0.6). The VZA and RAA effects were less for sparsely vegetated areas (NDVI ≈ 0.3–0.35) and the smallest effect on NDVI was found for dense vegetation (NDVI ≈ 0.7). These results have implications for the end users' interpretation of NDVI, and challenge the expediency of the MODIS NDVI compositing technique, which should be refined to distinguish between forward- and backward-scatter viewing direction by taking RAA into account.     

Transitioning from MODIS to VIIRS: an analysis of inter-consistency of NDVI data sets for agricultural monitoring

The Visible/Infrared Imager/Radiometer Suite (VIIRS) aboard the Suomi National Polar-orbiting Partnership (S-NPP) satellite was launched in 2011, in part to provide continuity with the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument aboard National Aeronautics and Space Administration’s (NASA) Terra and Aqua remote-sensing satellites. The VIIRS will eventually replace Aqua MODIS for both land science and applications and add to the coarse-resolution, long-term data record. It is, therefore, important to provide the user community with an assessment of the consistency of equivalent products from the two sensors. For this study, we do this in the context of example agricultural monitoring applications. Surface reflectance that is routinely delivered within the M{O,Y}D09 and VNP09 series of products provides critical input for generating downstream products. Given the range of applications utilizing the normalized difference vegetation index (NDVI) generated from the M{O,Y}D09 and VNP09 products and the inherent differences between MODIS and VIIRS sensors in calibration, spatial sampling, and spectral bands, the main objective of this study is to quantify uncertainties associated with transitioning from using MODIS to VIIRS-based NDVIs. In particular, we compare NDVIs derived from two sets of Level 3 MYD09 and VNP09 products with various spatial-temporal characteristics, namely 8-day composites at 500 m spatial resolution and daily climate modelling grid images at 0.05° spatial resolution. Spectral adjustment of VIIRS I1 (red) and I2 (near infra-red – NIR) bands to match MODIS/Aqua b1 (red) and b2 (NIR) bands is performed to remove a bias between MODIS and VIIRS-based red, NIR, and NDVI estimates. Overall, red reflectance, NIR reflectance, and NDVI uncertainties were 0.014, 0.029, and 0.056, respectively, for the 500 m product and 0.013, 0.016, and 0.032 for the 0.05° product. The study shows that MODIS and VIIRS NDVI data can be used interchangeably for applications with an uncertainty of less than 0.02–0.05, depending on the scale of spatial aggregation, which is typically the uncertainty of the individual data sets.     

Crop condition and yield simulations using Landsat and MODIS

Monitoring crop condition and yields at regional scales using imagery from operational satellites remains a challenge because of the problem in scaling local yield simulations to the regional scales. NOAA AVHRR satellite imagery has been traditionally used to monitor vegetation changes that are used indirectly to assess crop condition and yields. Additionally, the 1-km spatial resolution of NOAA AVHRR is not adequate for monitoring crops at the field level. Imagery from the new MODIS sensor onboard the NASA Terra satellite offers an excellent opportunity for daily coverage at 250-m resolution, which is adequate to monitor field sizes are larger than 25 ha. A field study was conducted in the predominantly corn and soybean area of Iowa to evaluate the applicability of the 8-day MODIS composite imagery in operational assessment of crop condition and yields. Ground-based canopy reflectance and leaf area index (LAI) measurements were used to calibrate the models. The MODIS data was used in a radiative transfer model to estimate LAI through the season. LAI was integrated into a climate-based crop simulation model to scale from local simulation of crop development and responses to a regional scale. Simulations of corn and soybean yields at a 1.6×1.6-km² grid scale were comparable to county yields reported by the USDA-National Agricultural Statistics Service (NASS). Weekly changes in soil moisture for the top 1-m profile were also simulated as part of the crop model as one of the critical parameters influencing crop condition and yields.     

Multi‐sensor data fusion and comparison of total column ozone

With many remote‐sensing instruments onboard satellites exploring the Earth's atmosphere, most data are processed to gridded daily maps. However, differences in the original spatial, temporal, and spectral resolution—as well as format, structure, and temporal and spatial coverage—make the data merging, or fusion, difficult. NASA Goddard Earth Sciences Data and Information Services Center (GES‐DISC) has archived several data products for various sensors in different formats, structures, and multi‐temporal and spatial scales for ocean, land, and atmosphere. In this investigation using Earth science data sets from multiple sources, an attempt was made to develop an optimal technique to merge the atmospheric products and provide interactive, online analysis tools for the user community. The merged/fused measurements provide a more comprehensive view of the atmosphere and improve coverage and accuracy, compared with a single instrument dataset. This paper describes ways of merging/fusing several NASA Earth Observing Systems (EOS) remote‐sensing datasets available at GES‐DISC. The applicability of various methods was investigated for merging total column ozone to implement these methods into Giovanni, the online interactive analysis tool developed by GES‐DISC. Ozone data fusion of MODerate resolution Imaging Spectrometer (MODIS) Terra and Aqua Level‐3 daily data sets was conducted, and the results were found to provide better coverage. Weighted averaging of Terra and Aqua data sets, with the consequent interpolation through the remaining gaps using Optimal Interpolation (OI), also was conducted and found to produce better results. Ozone Monitoring Instrument (OMI) total column ozone is reliable and provides better results than Atmospheric Infrared Sounder (AIRS) and MODIS. However, the agreement among these instruments is reasonable. The correlation is high (0.88) between OMI and AIRS total column ozone, while the correlation between OMI and MODIS Terra/Aqua fused total column ozone is 0.79.     

Evaluating MODIS soil fractional cover for arid regions,using albedo from high-spatial resolution satellite imagery

Broad-scale high-temporal frequency satellite imagery is increasingly used for environmental monitoring. While the normalized difference vegetation index (NDVI) is the most commonly used index to track changes in vegetation cover, newer spectral mixture approaches aim to quantify sub-pixel fractions of photosynthesizing vegetation, non-photosynthesizing vegetation, and exposed soil. Validation of the unmixing products is essential to enable confident use of the products for management and decision-making. The most frequently used validation method is by field data collection, but this is very time consuming and costly, in particular in remote regions where access is difficult.

This study developed and demonstrates an alternative method for quantifying land-cover fractions using high-spatial resolution satellite imagery. The research aimed to evaluate the bare soil fraction in a sub-pixel product, MODIS Fract-G, for the natural arid landscapes of the far west of South Australia. Twenty-two sample regions, of 3400 sampling points each, were investigated across several arid land types in the study area. Albedo thresholds were carefully determined in Advanced Land Observing Satellite Panchromatic Remote-sensing Instrument Stereo Mapping (ALOS PRISM) images (2.5 m spatial resolution), which separated predominantly bare soil from predominantly vegetated or covered soil, and created classified images. Correlation analysis was carried out between MODIS Fract-G bare soil fractional cover and ALOS PRISM bare soil proportions for the same areas. Results showed much lower correlations than expected, though limited agreement was found in some specific areas. It is posited that the Moderate Resolution Imaging Spectroradiometer (MODIS) fractional cover product, which is based on unmixing using the NDVI and a cellulose absorption index (CAI) proxy, may be generally unable to separate soil from vegetation in situations where both indices are low. In addition, separation is hampered by the lack of ‘pure pixels’ in this heterogeneous landscape. This suggests that the MODIS fractional cover product, at least in its present form, is unsuited to monitor sparsely vegetated arid landscapes.     

Monitoring of forage conditions with MODIS imagery in the Xilingol steppe,Inner Mongolia

A study was conducted to determine the potential suitability of Terra/MODIS imagery for monitoring short‐term phenological changes in forage conditions in a semi‐arid region. The study sites included four meadow steppes and six typical steppes in the Xilingol steppe in central Inner Mongolia, China. The live biomass, dead standing biomass, total biomass, crude protein (CP) concentration and standing CP were estimated from early April to late October using the Enhanced Vegetation Index (EVI) values from Terra imagery (500?m?pixels). Applying regression models, the EVI accounted for 80% of the variation in live biomass, 42% of the dead biomass, 77% of the total biomass, 11% of the CP concentration and 74% of the standing CP. MODIS/EVI is superior to AVHRR/NDVI when estimating forage quantity. Applying these results, the seasonal changes in live biomass and the standing CP could be described in the selected four sites with different degrees of grazing intensity. Generally, the increase in grazing intensity tended to decrease live biomass and standing CP. It was suggested that the EVI obtained from Terra imagery was an available predictor of the forage condition as measured by live biomass and standing CP. The MODIS/EVI values could provide information on the suitable timing of cutting for hay‐making and nutritive value to range managers.     

Sensitivity of vegetation phenology detection to the temporal resolution of satellite data

Vegetation phenology derived from satellite data has increasingly received attention for applications in environmental monitoring and modelling. The accuracy of phenological estimates, however, is unknown at the regional and global level because field validation data are insufficient. To assess the accuracy of satellite‐derived phenology, this study investigates the sensitivity of phenology detection to both the temporal resolution of sampling and the number of consecutive missing values (usually representing cloud cover) in the time series of satellite data. To do this, time series of daily vegetation index data for various ecosystems are modelled and simulated using data from Moderate‐Resolution Imaging Spectroradiometer (MODIS) data. The annual temporal data are then fitted using piecewise logistic functions, which are employed to calculate curvature change rate for detecting phenological transition dates. The results show that vegetation phenology can be estimated with a high precision from time series with temporal resolutions of 6–16 days even if daily data contains some uncertainties. If the temporal resolution is no coarser than 16 days for time series sampled using an average composite, the absolute errors are less than 3 days. On the other hand, the phase shift of temporal sampling is shown to have limited impacts on phenology detection. However, the accuracy of phenology detection may be reduced greatly if missing values in the time series of 16‐day MODIS data occur around the onsets of phenological transition dates. Even so, the probability that the absolute error in phenological estimates is greater than 5 days is less than 4% when only one period is missing in the time series of 16‐day data during vegetation growing seasons; this probability increases to 20% if there are two consecutive missing values.     

Evaluation of the onset of green-up in temperate deciduous broadleaf forests derived from Moderate Resolution Imaging Spectroradiometer (MODIS) data

Vegetation phenology is the chronology of periodic phases of development. It constitutes an efficient bio-indicator of impacts of climate changes and a key parameter for understanding and modelling vegetation-climate interactions and their implications on carbon cycling. Numerous studies were devoted to the remote sensing of vegetation phenology. Most of these were carried out using data acquired by AVHRR instrument onboard NOAA meteorological satellites. Since 1999, multispectral images were acquired over the whole earth surface every one to two days by MODIS instrument onboard Terra and Aqua platforms. In comparison with AVHRR, MODIS constitutes a significant technical improvement in terms of spatial resolution, spectral resolution, geolocation accuracy, atmospheric corrections scheme and cloud screening and sensor calibration. In this study, 250 m daily MODIS data were used to derive precise vegetation phenological dates over deciduous forest stands. Phenological markers derived from MODIS time-series and provided by MODIS Global Land Cover Dynamics product (MOD12Q2) were compared to field measurements carried out over the main deciduous forest stands across France and over five years. We show that the inflexion point of the asymmetric double-sigmoid function fitted to NDVI temporal profile is a good marker of the onset of green-up in deciduous stands. At plot level, the prediction uncertainty is 8.5 days and the bias is 3.5 days. MODIS Global Land Cover Dynamics MOD12Q2 provides estimates of onset of green-up dates which deviate substantially from in situ observations and do not perform better than the null model. RMSE values are 20.5 days (bias -17 days) using the onset of greenness increase and 36.5 days (bias 34.5 days) using the onset of greenness maximum. An improvement of prediction quality is obtained if we consider the average of MOD12Q2 onset of greenness increase and maximum as marker of onset of green-up date. RMSE decreases to 16.5 days and bias to 7.5 days.     

Application of MODIS derived parameters for regional crop yield assessment

NOAA AVHRR has been used extensively for monitoring vegetation condition and changes across the United States. Integration of crop growth models with MODIS imagery at 250 m resolution from the Terra Satellite potentially offers an opportunity for operational assessment of the crop condition and yield at both field and regional scales. The primary objective of this research was to evaluate the quality of the MODIS 250 m resolution data for retrieval of crop biophysical parameters that could be integrated in crop yield simulation models. A secondary objective was evaluating the potential use of MODIS 250 m resolution data for crop classification. A field study (24 fields) was conducted during the 2000 crop season in McLean County, Illinois, in the U.S. Midwest to evaluate the applicability of the MODIS 8-day, 250 m resolution composite imagery (version 4) for operational assessment of crop condition and yields. Ground-based canopy and leaf reflectance and leaf area index (LAI) measurements were used to calibrate a radiative transfer model to create a look up table (LUT) that was used to simulate LAI. The seasonal trend of MODIS derived LAI was used to find crop model parameters by adjusting the LAI simulated from the climate-based crop yield model. Other intermediate products such as crop phenological events were adjusted from the LAI seasonal profile. Corn (Zea mays L.) and soybean (Glycine max (L.) Merr.) yield simulations were conducted on a 1.6 × 1.6 km² spatial resolution grid and the results integrated to the county level. The results were within 10% of county yields reported by the USDA National Agricultural Statistics Service (NASS).     

Feasibility of using Geostationary Ocean Colour Imager (GOCI) data for land applications after atmospheric correction and bidirectional reflectance distribution function modelling

Korea's Geostationary Ocean Colour Imager (GOCI) has very high temporal resolution as well as wide spatial coverage. There is thus great interest in testing its applicability for monitoring land areas in addition to ocean areas. GOCI has eight spectral bands, from blue to near-infrared. These bands can be sensitive to vegetation change, but their wavelength ranges are slightly different from those of the extensively studied Moderate Resolution Imaging Spectroradiometer (MODIS). This study examines whether GOCI data can be applied for land monitoring and how GOCI data should be processed so as to reflect the spectral characteristics of land surfaces as detected by polar-orbit satellite sensors. Several image processing steps were performed for the GOCI data, including atmospheric correction and semi-empirical bidirectional reflectance distribution function modelling, before the results were compared with the MODIS land-surface product. Among the four GOCI normalized difference vegetation index (NDVI) products tested in this study, the GOCI NDVI with viewing-angle-adjusted reflectance showed the best agreement with MODIS NDVI calculated from normalized reflectance, with the lowest root mean square error of 0.126. Additionally, its temporal trends over forest and mixed vegetation areas were similar to those of MODIS NDVI during the study period from September to December.     

Land surface response to precipitation events using MODIS and NEXRAD data

We have developed a wavelet‐based information theoretic approach to examine the interaction between precipitation (PPT) forcing events and the land surface response. Combining Next Generation Weather Radar (NEXRAD) PPT with Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation (NDVI) and surface temperature (T _s) data over the Missouri Basin in the north‐central USA, we are able to address the spatial and temporal fluctuations surrounding the hydrometeorology of grassland ecosystems. Information theory metrics of entropy and mutual information content are combined with a wavelet multi‐resolution analysis to examine to what extent the observed PPT signal directly determines the spatial distribution of the land surface temperature and vegetation and how this relationship varies with spatial scale. Results indicate that (1) there is a reduction in the temporal variance of the wavelet coefficients as the signal is transferred from the PPT into the surface temperature and finally the vegetation signal, (2) there are significant correlations as a function of spatial resolution between PPT–NDVI and PPT–T _s signals which generally increase with spatial resolution, while there is little correlation between the NDVI and T _s signals as a function of resolution, and (3) the scale‐wise entropy and the mutual information content of the signals increase for all fields as the spatial resolution increases. This provides a methodology for determining the relative impact of regional climatology and local land–atmosphere interactions as a function of spatial scale.     

2001~2010年松木希错流域植被动态变化遥感研究

遥感在区域植被变化研究中具有十分重要的作用,能为大面积监测植被状况的演化过程提供技术支持。NDVI在高植被覆盖地区存在过饱和现象,对稀疏地区的植被变化尤其敏感。以古里雅冰帽南部的松木希错流域植被相对稀疏区域为研究区,基于MODIS NDVI数据和逐月气象观测数据,以及RS和GIS平台,对该区域2001~2010年主要植被变化趋势进行了初步研究,并对植被变化与气候驱动因子的关系进行了分析和探讨。结果表明:① 2001~2010年间该区域的植被活动有加强趋势;② NDVI表明研究区植被生长季较短(5~9月),NDVI浮动区间为0.11～0.13,低于全国水平(0.3～0.35),也低于全球稀疏灌丛的平均水平(0.2～0.4);③NDVI与年均气温整体上呈正相关,而与年降水量相关性不强。表明近年来持续升温是影响该区域植被活动加强的最主要原因。     

Evaluation of MODIS and NOAA AVHRR vegetation indices with in situ measurements in a semi‐arid environment

Much effort has been made in recent years to improve the spectral and spatial resolution of satellite sensors to develop improved vegetation indices reflecting surface conditions. In this study satellite vegetation indices from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Advanced Very High Resolution Radiometer (AVHRR) are evaluated against two years of in situ measurements of vegetation indices in Senegal. The in situ measurements are obtained using four masts equipped with self‐registrating multispectral radiometers designed for the same wavelengths as the satellite sensor channels. In situ measurements of the MODIS Normalized Difference Vegetation Index (NDVI) and AVHRR NDVI are equally sensitive to vegetation; however, the MODIS NDVI is consistently higher than the AVHRR NDVI. The MODIS Enhanced Vegetation Index (EVI) proved more sensitive to dense vegetation than both AVHRR NDVI and MODIS NDVI. EVI and NDVI based on the MODIS 16‐day constrained view angle maximum value composite (CV‐MVC) product captured the seasonal dynamics of the field observations satisfactorily but a standard 16‐day MVC product estimated from the daily MODIS surface reflectance data without view angle constraints yielded higher correlations between the satellite indices and field measurements (R ² values ranging from 0.74 to 0.98). The standard MVC regressions furthermore approach a 1?:?1 line with in situ measured values compared to the CV‐MVC regressions. The 16‐day MVC AVHRR data did not satisfactorily reflect the variation in the in situ data. Seasonal variation in the in situ measurements is captured reasonably with R ² values of 0.75 in 2001 and 0.64 in 2002, but the dynamic range of the AVHRR satellite data is very low—about a third to a half of the values from in situ measurements. Consequently the in situ vegetation indices were emulated much better by the MODIS indices than by the AVHRR NDVI.     

Assessing vegetation condition in the presence of biomass burning smoke by applying the Aerosol‐free Vegetation Index (AFRI) on MODIS images

Vegetation indices (VIs) such as the Normalized Difference Vegetation Index (NDVI) are widely used for assessing vegetation cover and condition. One of the NDVI's significant disadvantages is its sensitivity to aerosols in the atmosphere, hence several atmospherically resistant VIs were formulated using the difference in the radiance between the blue and the red spectral bands. The state‐of‐the‐art atmospherically resistant VI, which is a standard Moderate Resolution Imaging Spectroradiometer (MODIS) product, together with the NDVI, is the Enhanced Vegetation Index (EVI). A different approach introduced the Aerosol‐free Vegetation Index (AFRI) that is based on the correlation between the shortwave infrared (SWIR) and the visible red bands. The AFRI main advantage is in penetrating an opaque atmosphere influenced by biomass burning smoke, without the need for explicit correction for the aerosol effect. The objective of this research was to compare the performance of these three VIs under smoke conditions. The AFRI was applied to the 2.1 µm SWIR channel of the MODIS sensor onboard the Earth Observing System (EOS) Terra and Aqua satellites in order to assess its functionality on these imaging platforms. The AFRI performance was compared with those of NDVI and EVI. All VIs were calculated on images with and without present smoke, using the surface‐reflectance MODIS product, for three case studies of fires in Arizona, California, and Zambia. The MODIS Fire Product was embedded on the images in order to identify the exact location of the active fires. Although good correlations were observed between all VIs in the absence of smoke (in the Arizona case R ² = 0.86, 0.77, 0.88 for the NDVI–EVI, AFRI–EVI, and AFRI–NDVI, respectively) under smoke conditions a high correlation was maintained between the NDVI and the EVI, while low correlations were found for the AFRI–EVI and AFRI–NDVI (0.21 and 0.16, for the Arizona case, respectively). A time series of MODIS images recorded over Zambia during the summer of 2000 was tested and showed high NDVI fluctuations during the study period due to oscillations in aerosol optical thickness values despite application of aerosol corrections on the images. In contrast, the AFRI showed smoother variations and managed to better assess the vegetation condition. It is concluded that, beneath the biomass burning smoke, the AFRI is more effective than the EVI in observing the vegetation conditions.     

Obtaining crop-specific time profiles of NDVI: the use of unmixing approaches for serving the continuity between SPOT-VGT and PROBA-V time series

The study examined the potential of two unmixing approaches for deriving crop-specific normalized difference vegetation index (NDVI) profiles so that upon availability of Project for On-Board Autonomy – Vegetation (PROBA-V) imagery in winter 2013, this new data set can be combined with existing Satellite Pour l’Observation de la Terre – VEGETATION (SPOT-VGT) data despite the differences in spatial resolution (300 m of PROBA-V versus 1 km of SPOT-VGT). To study the problem, two data sets were analysed: (1) a set of 10 temporal NDVI images, with 300 and 1000 m spatial resolution, from the state of São Paulo (Brazil) synthesized from 30 m Landsat Thematic Mapper (TM) images, and (2) a corresponding set of 10 observed Moderate Resolution Imaging Spectroradiometer (MODIS) images (250 m spatial resolution). To mimic the influence of noise on the retrieval accuracy, different sensor/atmospheric noise levels were applied to the first data set. For the unmixing analysis, a high-resolution land-cover (LC) map was used. The LC map was derived beforehand using a different set of Landsat TM images. The map distinguishes nine classes, with four different sugarcane stages, two agricultural sub-classes, plus forest, pasture, and urban/water. Unmixing aiming at the retrieval of crop-specific NDVI profiles was done at administrative level. For the synthesized data set it was demonstrated that the ‘true’ NDVI temporal profiles of different land-cover classes (from 30 m TM data) can generally be retrieved with high accuracy. The two simulated sensors (PROBA-V and SPOT-VGT) and the two unmixing algorithms gave similar results. Analysing the MODIS data set, we also found a good correspondence between the modelled NDVI profiles (both approaches) and the (true) Landsat temporal endmembers.     

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.