Validating MODIS surface reflectance based on field spectral measurements

The 16-day revisit cycle of Landsat has long limited its use for studying global biophysical processes, which evolve rapidly during the growing season. In cloudy areas of the Earth, the problem is compounded, and researchers are fortunate to get two to three clear images per year. At the same time, the coarse resolution of sensors such as the Advanced Very High Resolution Radiometer and Moderate Resolution Imaging Spectroradiometer (MODIS) limits the sensors' ability to quantify biophysical processes in heterogeneous landscapes. In this paper, the authors present a new spatial and temporal adaptive reflectance fusion model (STARFM) algorithm to blend Landsat and MODIS surface reflectance. Using this approach, high-frequency temporal information from MODIS and high-resolution spatial information from Landsat can be blended for applications that require high resolution in both time and space. The MODIS 500-m surface reflectance and the 16-day repeat cycle Landsat Enhanced Thematic Mapper Plus (ETM+) 30-m surface reflectance are used to produce a synthetic daily surface reflectance product at ETM+ spatial resolution. The authors present results both with simulated (model) data and actual Landsat/MODIS acquisitions. In general, the STARFM accurately predicts surface reflectance at an effective resolution close to that of the ETM+. However, the performance depends on the characteristic patch size of the landscape and degrades somewhat when used on extremely heterogeneous fine-grained landscapes

The Multi-Angle Imaging Spectroradiometer (MISR) is a powerful sensor for leaf area index (LAI) mapping with its simultaneous multi-angle observations. However, the LAI product derived from MISR observations has low temporal resolution, which is unsatisfactory for many applications. This paper presents an algorithm that expands the MISR LAI product to high temporal resolution with the aid of Moderate Resolution Imaging Spectroradiometer (MODIS) data. The algorithm establishes relationships between the MISR LAI and the MODIS red/near-infrared band ratio (simple ratio (SR)) pixel by pixel using coincident data of these two sensors for the past nine years. Using these pixel-based SR-LAI relationships, a new LAI product with the merits of the original MISR product and high temporal resolution is obtained from MODIS surface reflectance. The expanded LAI series was compared with the original MISR and MODIS LAI products, as well as field LAI measurements made at the Baohe and Maoershan forest sites and the Hulunbeier grassland site, to assess the algorithm's performance. The results show that the temporal coverage of the MISR LAI improved from 15.5% to 65.2% in an 8-day composite, and the mean root-mean-square error is 0.74 for the vegetated pixels. This LAI product has similar temporal consistency and seasonal dynamics to the existing MODIS LAI product generated from the main algorithm, but is more robust against the low quality of reflectance inputs. The expanded LAI product differs with field measurements by about 11.5%, with agreement to field observations at all three sites within an accuracy of 0.8 LAI.

This paper explores snow bidirectional reflectance distribution function (BRDF) properties over some snow covered regions using Moderate Resolution Imaging Spectroradiometer (MODIS) and Multi-angle Imaging SpectroRadiometer (MISR) surface reflectance products. In the visible and near infrared (NIR) region, MODIS and MISR surface bidirectional reflectance factors (BRFs) over snow are accumulated to extract snow BRDF properties. Five surface BRDF models are concerned to simulate snow surface reflectance shape. All the models capture the distribution of snow BRFs with limit of accuracy. The simulated BRFs from several models have similar distribution trend and different details. The BRDF properties discussed can be used as background in the snow BRDF retrieval from spaceborne measurements.

The impact of size variations in the ground instantaneous field of view of pixels on MODIS BRDF modelling

The suitability of using Moderate Resolution Imaging Spectroradiometer (MODIS) images for surface soil moisture estimation to investigate the importance of soil moisture in different applications, such as agriculture, hydrology, meteorology and natural disaster management, is evaluated in this study. Soil moisture field measurements and MODIS images of relevant dates have been acquired. Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and Normalized Difference Water Index (NDWI) are calculated from MODIS images. In addition, MODIS Land Surface Temperature (LST) data (MOD11A1) are used in this analysis. Four different soil moisture estimation models, which are based on NDVI–LST, EVI–LST, NDVI–LST–NDWI and EVI–LST–NDWI, are developed and their accuracies are assessed. Statistical analysis shows that replacing EVI with NDVI in the model that is based on LST and NDVI increases the accuracy of soil moisture estimation. Accuracy evaluation of soil moisture estimation using check points shows that the model based on LST, EVI and NDWI values gives a higher accuracy than that based on LST and EVI values. It is concluded that the model based on the three indices is a suitable model to estimate soil moisture through MODIS imagery.

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.

Knowledge of the variation and regulation of the diffuse attenuation coefficient is important for management of lake water quality and ecological environments. Data were collected seven times (378 samples) from 2006 to 2011 in Lake Taihu to improve the existing semianalytical model of diffuse attenuation coefficient at 490 nm (Kd(490)) by measuring in situ Kd(490) and remote sensing reflectance (Rrs). This semianalytical model was applied to determine the long‐term variation of Kd(490) and its driving factors from Moderate Resolution Imaging Spectroradiometer (MODIS) images. Our semianalytical model of Kd(490) improves the retrieval accuracy of Kd(490) from 29.60% to 22.58% (root‐mean‐square percentage error) and 24.28 ± 18.54% (mean value ± standard deviation) to 19.29 ± 14.34% (mean absolute relative error). The Kd(490) from maps produced from MODIS images clearly present the temporal and spatial variations of Kd(490). This indicates the suitability of using the improved semianalytical model and MODIS to observe Kd(490) in turbid shallow inland waters. Long‐term MODIS observation of Kd(490) shows distinct seasonal and spatial variations. Temporally, the highest and the lowest Kd(490) were found in winter and in summer, respectively. Spatially, high Kd(490) values were primarily distributed in the southwestern and northwestern Lake Taihu, and relatively low Kd(490) values were distributed in Meiliang Bay and Gonghu Bay (the eastern part of the lake and East Bay were excluded in this study due to invalid retrieval results of Kd(490)). Many covarying points were found between Kd(490) and maximum wind speed from long‐term observations, suggesting that wind‐induced sediment resuspension is the principal driver of high Kd(490). Correlation was also observed between Kd(490) and heavy rainfall events, suggesting the regulating effect of heavy rainfall.

We used publicly available digital spatial datasets to study the area extents and their horizontal variations of two water bodies within the Danjiangkou Reservoir, China. Between 2003 and 2005, the water levels varied from 140 to 149 m above mean sea level as measured by the Geoscience Laser Altimeter System (GLAS). Detailed procedures to derive the horizontal extents and variations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) coupled with GLAS data and to verify the extents and variations spatially were provided. For the water bodies on the north and west, the surface water extents derived from four MODIS images varied between 174 and 218 km2 and from 96 to 135 km2, respectively. The extents by inundating the DEM using the GLAS data were 178–212 km2 for the water body on the north and 104–118 km2 for the water body on the west. The spatial verifications of surface water extents derived from the MODIS images versus DEM coupled with GLAS data agreed 83–93%. Within the ring areas between water/land boundaries at elevations of 140 and 147 m, and 140 and 149 m, the spatial agreement was 52–75%.

MODIS (Moderate Resolution Imaging Spectroradiometer) is a kind of new weather satellite data. Few weather satellite images obtained are all clear sky and they are always influenced by cloud more or less. Cloud is a large obstacle to remote sensing image processing and analysis all the while. In order to extract objective information more effective, cloud should be removed from the remote sensing images, which is an essential sector in the image preprocessing. Cloud detection is the most important processing before removing cloud. Taking it into account that MODIS data includes thirty-six bands, especially the infrared channels subdivided, it has realized cloud detection in MODIS images by multi-spectral synthesis method and cloud detection index in this paper. Owing to the limitation to a certainty of the above methods, an automatic cloud detection algorithm is applied based on the spatial texture analysis and neural network in this research. At last the cloud detection results gained by different ways are testified each other and analyzed by comparison. It found that the results are consistent, which shows that the cloud-contaminate pixels are detected successfully. It not only lays a good foundation for the cloud removing, but also can improve the precision of remote sensing image recognition, classification and inverse in this study.

The detection and mapping of burned areas from wildland fires is one of the most important approaches for evaluating the impacts of fire events. In this study, a novel burned area detection algorithm for rapid response applications using Moderate Resolution Imaging Spectroradiometer (MODIS) 500 m surface reflectance data was developed. Spectra from bands 5 and 6, the composite indices of the Normalized Burn Ratio, and the Normalized Difference Vegetation Index were employed as indicators to discover burned pixels. Historical statistical data were used to provide pre-fire baseline information. Differences in the current (post-fire) and historical (pre-fire) data were input into a support vector machine classifier, and the fire-affected pixels were detected and mapped by the support vector machine classification process. Compared with the existing MODIS level 3 monthly burned area product MCD45, the new algorithm is able to generate burned area maps on a daily basis when new data become available, which is more applicable to rapid response scenarios when major fire incidents occur. The algorithm was tested in three mega-fire cases that occurred in the continental USA. The experimental results were validated against the fire perimeter database generated by the Geospatial Multi-Agency Coordination Group and were compared with the MCD45 product. The validation results indicated that the algorithm was effective in detecting burned areas caused by mega-fires.

The Moderate Resolution Imaging Spectroradiometer (MODIS) is an important data set in global burned-area mapping. The MODIS global burned-area product has a coarse spatial resolution at approximately 500 m, which often introduces errors to the mapped burned areas. In this letter, a novel subpixel mapping (SPM) approach was proposed to produce burned-area maps at the fine spatial resolution similar to Landsat imagery, by exploring the spectral and spatial information provided by the second and fifth bands of MODIS. The proposed SPM approach aims to refine the estimate of burned areas, which have been detected by the MODIS global burned-area product. The performance of the proposed SPM approach was assessed with an experiment area containing six burned areas, by comparing with the MODIS burned-area product MCD45. The result shows that the average omission error decreased from 52.26% for MCD45 to 16.74% for SPM, and the average commission error decreased from 21.76% for MCD45 to 12.54% for SPM. The kappa value increased from 0.5583 for MCD45 to 0.8756 for SPM, indicating that the proposed SPM approach is effective in reducing the influence of the coarse spatial resolution of MODIS imagery in mapping a burned area and refining existing global burned-area products.

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 2 = 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.

This article presents a comprehensive subpixel water mapping algorithm to automatically produce routinely open water fraction maps in the Tibetan Plateau (TP) with the Moderate Resolution Imaging Spectroradiometer (MODIS). A multi-index threshold endmember extraction method was applied to select the endmembers from MODIS images. To incorporate endmember variability, an endmember selection strategy, called the combined use of typical and neighboring endmembers, was adopted in multiple endmember spectral mixture analysis (MESMA), which can assure a robust subpixel water fractions estimation. The accuracy of the algorithm was assessed at both the local scale and regional scale. At the local scale, a comparison using the eight pairs of MODIS/Landsat 8 Operational Land Imager (OLI) water maps demonstrated that subpixels water fractions were well retrieved with a root mean square error (RMSE) of 7.86% and determination coefficient (R2) of 0.98. At the regional scale, the MODIS water fraction map in October 2014 matches well with the TP lake data set and the Global Lake and Wetland Database (GLWD) in both latitudinal and longitudinal distribution. The lake area estimation is more consistent with the reference TP lake data set (difference of −3.15%) than the MODIS Land Water Mask (MOD44W) (difference of −6.39%).

Techniques are presented for generating thermal infrared and visible composite images from cloud-free portions of MODIS (MODerate resolution Imaging Spectroradiometer) images closely spaced in time, with a focus on studies of landfast sea ice along the East Antarctic coast. Composite image inclusion criteria are based on modified MODIS EOS cloud mask product results. The compositing process places emphasis on retaining maximum spatial resolution while minimizing computing storage space requirements. Composite images can be produced either as a regular product (e.g., on a 10-day grid), or dynamically (whenever enough information is acquired to produce a new output image). The techniques presented are applicable at any latitude, are available for all MODIS channels at their native resolution, can combine Aqua and Terra images, and can produce maps in any output projection. However, due to the polar orbit of NASA's Terra and Aqua satellites which host the MODIS instrument, more frequent coverage is produced at higher latitudes. Thus, the techniques presented are particularly applicable to polar research. Examples of summertime (visible) and wintertime (IR) composite image generation of the landfast sea ice around the Mertz Glacier region, East Antarctica, are included.

In recent years, there has been a significant increase in the use of remotely sensed data to address global issues. With the open data policy, the data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Enhanced Thematic Mapper Plus (ETM+) sensors have become a critical component of numerous applications. These two sensors have been operational for more than a decade, providing a rich archive of multispectral imagery for analysis of mutitemporal remote sensing data. This paper focuses on evaluating the radiometric calibration agreement between MODIS and ETM+ using the near-simultaneous and cloud-free image pairs over an African pseudo-invariant calibration site, Libya 4. To account for the combined uncertainties in the top-of-atmosphere (TOA) reflectance due to surface and atmospheric bidirectional reflectance distribution function (BRDF), a semiempirical BRDF model was adopted to normalize the TOA reflectance to the same illumination and viewing geometry. In addition, the spectra from the Earth Observing-1 (EO-1) Hyperion were used to compute spectral corrections between the corresponding MODIS and ETM+ spectral bands. As EO-1 Hyperion scenes were not available for all MODIS and ETM+ data pairs, MODerate resolution atmospheric TRANsmission (MODTRAN) 5.0 simulations were also used to adjust for differences due to the presence or lack of absorption features in some of the bands. A MODIS split-window algorithm provides the atmospheric water vapor column abundance during the overpasses for the MODTRAN simulations. Additionally, the column atmospheric water vapor content during the overpass was retrieved using the MODIS precipitable water vapor product. After performing these adjustments, the radiometric cross-calibration of the two sensors was consistent to within 7%. Some drifts in the response of the bands are evident, with MODIS band 3 being the largest of about 6% over 10 years, a change that will be corrected in Collection 6 MODIS processing.