Search

Time series of the Normalized Difference Vegetation Index (NDVI) derived from satellite observations provide important information on the state of terrestrial vegetation over a wide range of spatiotemporal scales. For understanding long-term changes in terrestrial ecosystems (post-1981), data collected by the Advanced Very High Resolution Radiometer (AVHRR) on board the satellites of National Oceanic and Atmospheric Administration (NOAA) series is a unique source of information. In this paper, we describe a new processing methodology for a comprehensive AVHRR data set at 1 km spatial resolution acquired over Canada, the northern United States and Greenland post-1981. The methodology incorporates a pre-processing algorithm, Canadian AVHRR Processing System (CAPS), recently developed by the Canada Centre of Remote Sensing (CCRS), which enables highly accurate geolocation and ortho-rectification at efficiency rates of > 90%. Once image navigation is completed, our approach consists of five key steps: first, two clear-sky composites for each 10 day interval are generated from the forward or backward scattering hemisphere; second, AVHRR Channel 1 and 2 reflectances are normalized to the AVHRR/3 on board NOAA-17 to account for differences in the spectral response function among the AVHRR sensors; third, atmospheric correction is performed using the Simplified Method for Atmospheric correction (SMAC) algorithm, using standard meteorological data sets (water vapor, surface level air pressure, ozone); fourth, NDVI is calculated based on atmospherically corrected Channel 1 and 2 reflectances; and finally, the NDVI is adjusted for directional effects based on the Ross-Thick Li-Sparse Bidirectional Reflectance Distribution Function (BRDF) model. The processed NDVI data are compared to an equivalent spatially and temporally overlapping MODIS NDVI data set from 2001 to 2005 for validation. Results at continental scale indicate that time series of MODIS and AVHRR were similar for a wide range of biomes and generalized ecoregions. Analysis stratified by land cover indicated that the correlation was strongest for homogeneous land cover types, such as cropland, when compared to structurally more diverse classes, such as deciduous broadleaf forests. The comparison of the NDVI at the local scale at seven sites of the Fluxnet Canada Research Network resulted in the correlation coefficient r = 0.95. Given confidence in the processing approach, this NDVI data set can be a valuable source of information for climate and vegetation-related studies over Canada and the northern United States.

Topography and accuracy of image geometric registration significantly affect the quality of satellite data, since pixels are displaced depending on surface elevation and viewing geometry. This effect should be corrected for through the process of accurate image navigation and orthorectification in order to meet the geolocation accuracy for systematic observations specified by the Global Climate Observing System (GCOS) requirements for satellite climate data records. We investigated the impact of orthorectification on the accuracy of maximum Normalized Difference Vegetation Index (NDVI) composite data for a mountain region in north-western Canada at various spatial resolutions (1 km, 4 km, 5 km, and 8 km). Data from AVHRR on board NOAA-11 (1989 and 1990) and NOAA-16 (2001, 2002, and 2003) processed using a system called CAPS (Canadian AVHRR Processing System) for the month of August were considered. Results demonstrate the significant impact of orthorectification on the quality of composite NDVI data in mountainous terrain. Differences between orthorectified and non-orthorectified NDVI composites (ΔNDVI) adopted both large positive and negative values, with the 1% and 99% percentiles of ΔNDVI at 1 km resolution spanning values between − 0.16

This work extends the previous study of Trishchenko et al. [Trishchenko, A. P., Cihlar, J., & Li, Z. (2002). Effects of spectral response function on surface reflectance and NDVI measured with moderate resolution satellite sensors. Remote Sensing of Environment 81 (1), 1–18] that analyzed the spectral response function (SRF) effect for the Advanced Very High Resolution Radiometer (AVHRR) onboard the NOAA satellites NOAA-6 to NOAA-16 as well as the Moderate Resolution Imaging Spectroradiometer (MODIS), the VEGETATION sensor (VGT) and the Global Imager (GLI). The developed approach is now applied to cover three new AVHRR sensors launched in recent years on NOAA-17, 18, and METOP-A platforms. As in the previous study, the results are provided relative to the reference sensor AVHRR NOAA-9. The differences in reflectance among these three radiometers relative to the AVHRR NOAA-9 are similar to each other and range from − 0.015 to 0.015 (− 20% to + 2% relative) for visible (red) channel, and from − 0.03 to 0.02 (− 5% to 5%) for the near infrared (NIR) channel. The absolute change in the Normalized Difference Vegetation Index (NDVI) ranged from − 0.03 to + 0.06. Due to systematic biases of the visible channels toward smaller values and the NIR channels toward slightly larger values, the overall systematic biases for NDVI are positive. The polynomial approximations are provided for the bulk spectral correction with respect to the AVHRR NOAA-9 for consistency with previous study. Analysis was also conducted for the SRF effect only among the AVHRR-3 type of radiometer on NOAA-15, 16, 17, 18 and METOP-A using AVHRR NOAA-18 as a reference. The results show more consistency between sensors with typical correction being under 5% (or 0.01 in absolute values). The AVHRR METOP-A reveals the most different behavior among the AVHRR-3 group with generally positive bias for visible channel (up to + 5%, relative), slightly negative bias for the NIR channel (1%–2% relative), and negative NDVI bias (− 0.02 to + 0.005). Polynomial corrections are also suggested for normalization of AVHRR on NOAA-15, 16, 17 and METOP-A to AVHRR NOAA-18.

A new technology was developed at the Canada Centre for Remote Sensing (CCRS) for generating Canada-wide and North America continental scale clear-sky composites at 250 m spatial resolution for all seven MODIS land spectral bands (B1–B7). The MODIS Level 1B (MOD02) swath level data are used as input to circumvent the problems with image distortion in the mid latitude and polar regions inherent to the global sinusoidal (SIN) projection utilized for the standard MODIS data products. The MODIS 500 m land bands B3 to B7 are first downscaled to 250 m resolution using an adaptive regression and normalization scheme for compatibility with the 250 m bands B1 and B2. A new method has been developed to produce the mask of clear-sky, cloud and cloud shadow at 250 m resolution. It shows substantial advantages in comparison with the MODIS 250 m standard cloud masks. The testing of new cloud mask showed that it is in reasonable agreement with the MODIS 1-km standard product once it is aggregated to 1-km scale, while the cloud shadow detection looks more reliable with the new methodology. Nevertheless, more quantitative analyses of the presented scene identification technique are required to understand its performance over the range of input scenes in various seasons. The new clear-sky compositing scheme employs a scene-dependent decision matrix. It is demonstrated that this new scheme provides better results than any others based on a single compositing criterion, such as maximum NDVI or minimum visible reflectance. To account for surface bi-directional properties, two clear-sky composites for the same time period are produced by separating backward scattering and forward scattering geometries, which separate pixels with the sun-satellite relative azimuth angles within 90°–270° and outside of this range. Comparison with Landsat imagery and with MODIS standard composite products demonstrated the advantage of the new technique for screening cloud and cloud shadow, and generating high spatial resolution MODIS clear-sky composites. The new data products are mapped in the Lambert Conformal Conic (LCC) projection for Canada and the Lambert Azimuthal Equal-Area (LAEA) projection for North America. Presently this activity is limited to MODIS/TERRA due to known problems with band-to-band registration and noisy SWIR channels on MODIS/AQUA.

A time series of burned land areas was generated for a 23 year period (1984–2006) using 10-day composites of AVHRR data. The study area covers 1.6 million km2 of boreal forest in western Canada. The algorithm was intended to be consistent throughout the study period and region, and to avoid commission errors, so as to obtain a reliable sample of temporal trends in burned area in the region. The algorithm relies on temporal comparisons of several spectral indices (GEMI, BAI), as well as near infrared reflectance. It emphasizes the stability of the post-fire signal, to avoid false detections associated with cloud, cloud shadows, missed data and radiometric or geometric calibration between AVHRR sensors. Final results show a very consistent temporal adjustment to official statistics and fire perimeters, with very low commission error (

For quantitative studies of vegetation dynamics, satellite data need to be corrected for spurious effects. In this study, we have applied several changes to an earlier advanced very high resolution radiometer (AVHRR) processing methodology (ABC3; [Remote Sens. Environ. 60 (1997) 35; J. Geophys. Res.-Atmos. 102 (1997) 29625; Can. J. Remote Sens. 23 (1997) 163]), to better represent the various physical processes causing contamination of the AVHRR measurements. These included published recent estimates of the NOAA-11 and NOAA-14 AVHRR calibration trajectories for channels 1 and 2; the best available estimates for the water vapour, aerosol and ozone amounts at the time of AVHRR data acquisition; an improved bidirectional reflectance algorithm that also takes into consideration surface topography; and an improved image screening algorithm for contaminated pixels. Unlike the previous study that compared the composite images to a single-date AVHRR image, we employed coincident TM images to approximate the AVHRR pixel field of view during the data acquisition. Compared to ABC3, the modified procedure ABC3V2 was found to improve the accuracy of AVHRR pixel reflectance estimates, both in the sensitivity (slope) of the regression and in r2. The improvements were especially significant in AVHRR channel 1. In comparison with reference values derived from two full TM scenes, the corrected AVHRR surface reflectance estimates had average standard errors values of ±0.009 for AVHRR C1, ±0.019 for C2, and ±0.04 for NDVI; the corresponding r2 values were 0.55, 0.80, and 0.50, respectively. The changes in ABC3V2 were not able to completely remove interannual variability for land cover types with little or no vegetation cover, which would be expected to remain stable over time, and they increased the interannual variability of mixed forest and grassland. These results are attributed to a combination of increased sensitivity to interannual dynamics on one hand, and the inability to remove all sources of noise for barren or sparsely vegetated northern land cover types on the other.

We report the results of a modeling study on the sensitivity of normalized difference vegetation index (NDVI) and surface reflectance to differences in instrument spectral response functions (SRF) for various Advanced Very High Resolution Radiometers (AVHRR) onboard the National Oceanic and Atmospheric Administration's (NOAA) satellites NOAA-6–16 as well as the Moderate Resolution Imaging Spectroradiometer (MODIS), the Vegetation sensor (VGT), and the Global Imager (GLI). Modeling results were validated against real satellite observations employing AVHRR/NOAA-14 and -15 and MODIS, with a very good agreement. It is shown that for identical atmospheric state and similar surface spectral reflectance, the NDVI and spectral reflectances are sensitive to the sensor's SRF. Relative to a reference SRF for AVHRR/NOAA-9, the differences in reflectance among the AVHRRs range from −25% to 12% for visible channel (red) and from −2% to 4% for near-infrared (NIR) channel. Absolute change in NDVI among various AVHRRs ranged from −0.02 to 0.06. The most significant difference was observed for the AVHRR/3. Consistent results were obtained with the AVHRR sensors aboard the following afternoon satellites: NOAA-9, -11, and -12, whereas important discrepancies were found for other AVHRRs aboard NOAA-6 and -10 and especially those launched more recently (NOAA-15 and -16). Reflectance and NDVI measured by MODIS channels 1 and 2 also exhibit significant differences (up to 30–40%) relative to AVHRR. GLI and VGT have some specific features that should be taken into account when intercomparing surface or top of the atmosphere (TOA) reflectance as well as NDVI. Sensitivity of the SRF effect to variable atmospheric state (water vapor, aerosol, and ozone) was also investigated. Polynomial approximations are provided for bulk spectral correction with respect to AVHRR/NOAA-9.