Effects of spectral response function on surface reflectance and NDVI measured with moderate resolution satellite sensors: Extension to AVHRR NOAA-17, 18 and METOP-A

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 global 1° by 1° NDVI data set for climate studies derived from the GIMMS continental NDVI data

A nine-year (1982–1990) global normalized difference vegetation index (NDVI) data set with a spatial resolution of 1° by 1° and a temporal resolution of one month was compiled for use in climate studies. This data set was derived from higher resolution (5–8 km) monthly continental NDVI data sets that have been processed and archived by the Global Inventory Monitoring and Modelling Studies (GIMMS) group at NASA/Goddard Space Flight Center. The continental GIMMS NDVI data sets were calculated from Global Area Coverage (GAC) data collected at daily intervals by the Advanced Very High Resolution Radiometer (AVHRR) onboard the NOAA-7, -9 and -11 satellites

The global 1° by 1° NDVI data set was produced to calculate land surface parameters for use within general circulation model J of the atmosphere (GCM). In view of this quantitative application, an evaluation is given of the representation by the NDVI data of the spectral properties of vegetation at the landsurface. Errors are defined as deviations from measurements obtained under standard conditions, i.e., conditions at the-top-of-the-atmosphere with no clouds, clear atmosphere, near-nadir viewing angles, overhead Sun, and invariant soil background. The discussion includes an assessment of (1) the data collected and processed onboard the AVHRR; (2) processing of the AVHRR data into the continental GIMMS NDVI data sets; (3) resampling of the continental data sets to a 1° by 1° data set; and (4) propagation of inconsistencies and biases from (1), (2) and (3) into the 1° by 1° global NDVI data. Examples are shown of the temporal and spatial variations in spectral properties of vegetation contained in the 1° by 1° NDVI data, and these are compared with the dynamics of biophysical parameters derived from land cover classes that were used in previous climate studies.     

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.     

Multiscale geostatistical analysis of AVHRR, SPOT-VGT, and MODIS global NDVI products

Global NDVI data are routinely derived from the AVHRR, SPOT-VGT, and MODIS/Terra earth observation records for a range of applications from terrestrial vegetation monitoring to climate change modeling. This has led to a substantial interest in the harmonization of multisensor records. Most evaluations of the internal consistency and continuity of global multisensor NDVI products have focused on time-series harmonization in the spectral domain, often neglecting the spatial domain. We fill this void by applying variogram modeling (a) to evaluate the differences in spatial variability between 8-km AVHRR, 1-km SPOT-VGT, and 1-km, 500-m, and 250-m MODIS NDVI products over eight EOS (Earth Observing System) validation sites, and (b) to characterize the decay of spatial variability as a function of pixel size (i.e. data regularization) for spatially aggregated Landsat ETM+ NDVI products and a real multisensor dataset. First, we demonstrate that the conjunctive analysis of two variogram properties - the sill and the mean length scale metric - provides a robust assessment of the differences in spatial variability between multiscale NDVI products that are due to spatial (nominal pixel size, point spread function, and view angle) and non-spatial (sensor calibration, cloud clearing, atmospheric corrections, and length of multi-day compositing period) factors. Next, we show that as the nominal pixel size increases, the decay of spatial information content follows a logarithmic relationship with stronger fit value for the spatially aggregated NDVI products (R² = 0.9321) than for the native-resolution AVHRR, SPOT-VGT, and MODIS NDVI products (R² = 0.5064). This relationship serves as a reference for evaluation of the differences in spatial variability and length scales in multiscale datasets at native or aggregated spatial resolutions. The outcomes of this study suggest that multisensor NDVI records cannot be integrated into a long-term data record without proper consideration of all factors affecting their spatial consistency. Hence, we propose an approach for selecting the spatial resolution, at which differences in spatial variability between NDVI products from multiple sensors are minimized. This approach provides practical guidance for the harmonization of long-term multisensor datasets.     

Effects of atmospheric variation on AVHRR NDVI data

The AVHRR (Advanced Very High Resolution Radiometer) series of instruments has frequently been used for vegetation studies. The 25+ year record has enabled important time-series studies. Many applications use NDVI (Normalized Difference Vegetation Index), or derivatives of it, as their operational variable. However, most AVHRR datasets have incomplete atmospheric correction, because of which there is considerable, but largely unknown, uncertainty in the significance of differences in NDVI and other short wave observations from AVHRR instruments.The purpose of this study was to gain better understanding of the impact of incomplete or lack of atmospheric correction in widely-used, publicly available processed AVHRR-NDVI long-term datasets. This was accomplished by comparison with atmospherically corrected AVHRR data at AERONET (AErosol RObotic NETwork) sunphotometer sites in 1999. The datasets included in this study are: TOA (Top Of Atmosphere) that is with no atmospheric correction; PAL (Pathfinder AVHRR Land); and an early version of the new LTDR (Long Term Data Record) NDVI. The other publicly available datasets like GIMMS (Global Inventory Modeling and Mapping studies) and GVI (Global Vegetation Index) have atmospheric error budget similar to that of TOA, because no atmospheric correction is used in either processing stream. Of the three datasets, LTDR was found to have least errors (accuracy = 0.0064 to − 0.024, precision = 0.02 to 0.037 for clear and average atmospheric conditions) followed by PAL (accuracy = − 0.145 to − 0.035, precision = 0.0606 to 0.0418), and TOA (accuracy = − 0.0791 to − 0.112, precision = 0.0613 to 0.0684). It was also observed that temporal maximum value compositing technique does not cause significant improvement of precision in regions experiencing persistently high AOT (Aerosol Optical Thickness).     

Evaluation of AVHRR PAL and GIMMS 10‐day composite NDVI time series products using SPOT‐4 vegetation data for the African continent

Global 8 km resolution AVHRR (advanced very high resolution radiometer) NDVI (normalized difference vegetation index) 10‐day composite data sets have been used for numerous local to global scale vegetation time series studies during recent years. AVHRR Pathfinder (PAL) NDVI was available from 1981 until 2001, and the new AVHRR GIMMS NDVI was available from 1981 to the present time. A number of aspects potentially introduce noise in the NDVI data set due to the AVHRR sensor design and data processing. NDVI from SPOT‐4 VGT data is considered an improvement over AVHRR, and for this reason it is important to examine how and if the differences in sensor design and processing influence continental scale NDVI composite products. In this study, the quality of these AVHRR NDVI time series are evaluated by the continental scale 1 km resolution SPOT‐4 vegetation (VGT) 10‐day composite (S10) NDVI data. Three years of AVHRR PAL (1998–2000) and seven years of GIMMS (1998–2004) have been compared to 8 km resampled SPOT‐4 VGT (1998–2004) data. The dynamic range of SPOT‐4 VGT NDVI tends to be higher than the AVHRR PAL NDVI, whereas there is an exact match between AVHRR GIMMS NDVI and SPOT‐4 VGT NDVI. Ortho‐regression analysis on annually integrated values of AVHRR PAL/GIMMS and SPOT‐4 VGT on a continental scale reveals high correlations amongst the AVHRR and the SPOT data set, with lowest RMSE (root mean square error) on the GIMMS/SPOT‐4 VGT compared to the PAL/SPOT‐4 VGT.

Analyses on decade data likewise show that a linear relation exists between Spot‐4 VGT NDVI and the two AVHRR composite products; GIMMS explaining most of the Spot‐4 VGT NDVI variance compared to PAL. These results show that the AVHRR GIMMS NDVI is more consistent with Spot‐4 VGT NDVI compared to AVHRR PAL versus Spot‐4 VGT NDVI (in terms of RMSE and dynamic range) and can therefore be considered the more accurate long time AVHRR data record. Analyses performed on monthly maximum composites and decade composite data, however, reveal intra‐annual variations in the correlation between SPOT‐4 VGT and the two AVHRR data sets, which are attributed to different cloud masking algorithms. The SPOT‐4 VGT cloud‐screening algorithm is insufficient, thereby suppressing the rainy season NDVI.     

Satellite remote sensing of rangelands in Botswana II. NOAA AVHRR and herbaceous vegetation

NOAA-7 Advanced Very High Resolution Radiometer (AVHRR) global-area coverage (GAC) data for the visible and near-infrared bands were used to investigate the relationship between the normalized difference vegetation index (NDVI) and the herbaceous vegetation in three representative rangeland types in eastern Botswana. Regressions between Landsat MSS band-7/band-5 ratios and field measurements of the cover of the live parts of herbaceous plants, above-ground biomass of live herbaceous plants and bare ground were used in conjunction with MSS data in order to interpolate the field data to 144 km² areas for comparison with blocks of nine AVHRR GAC pixels. NOAA NDVI data were formed into 10-day composites in order to remove cloud cover and extreme off-nadir viewing angles. Both individual NDVI composite data and multitemporal integrations throughout the period May 1983-April 1984 were compared with the field data.

In multiple linear regressions, the cover and biomass of live herbaceous plants and bare ground measurements accounted for 42, 56 and 19 per cent respectively of the variation in NDVI. When factors were included in I he regression models to specify the site and date of acquisition of the data, between 93 and 99 per cent of the variation in NDVI was accounted for. The total herbaceous biomass at the end of the season was positively related to integrated NDVI, up lo the maximum biomass observed in a 12km × 12km area (590kgha^?1)- These results give a different regression of herbaceous biomass values on integrated AVHRR NDVI to that reported by Tucker et at. (1985 b) for Senegalese grasslands. The effect of the higher cover of the tree canopy in Botswana on this relationship and on the detection of forage available to livestock is discussed.     

Rainfall and vegetation monitoring in the Savanna Zone of the Democratic Republic of Sudan using the NOAA Advanced Very High Resolution Radiometer

NOAA-6 and NOAA-7 Advanced Very High Resolution Radiometer (AVHRR) global-area coverage (4?km ground resolution) data were obtained at three-day intervals throughout each of the four-month periods covering the 1980, 1983 and 1984 growing seasons, between latitudes 10° and 22° North in the Democratic Republic of Sudan. Daily rainfall data for twelve meteorological stations spanning the Savanna Zone were analysed. Rainfall in Sudan during 1980 was below normal, but in 1983 and 1984 there were moderate and severe droughts. The satellite data were used to calculate normalized difference vegetation index (NDVI) values from the visible and near-infrared bands of the satellite data. These were processed into ten-day composite data sets using the AVHRR thermal-infrared channel as a cloud screen and a temporal compositing procedure that reduces cloud contamination and selects viewing angles closest to nadir.

The ten-day composite NDVI values and the integrals of NDVI for each growing season were found to be closely correlated with rainfall. The constants of regressions between NDVI and rainfall were lower in 1983 and 1984 than in 1980, which suggests there was reduced water-use efficiency by the rangeland vegetation in drought years. It was found that July and August NDVI values were closely related to the integrated NDVI values; hence early- and mid-season NDVI data could be used to predict annual primary production. Images showing the geographical distribution of values of NDVI prepared for the three years clearly illustrate the effects of the 1983 and 1984 droughts, compared with the higher rainfall of 1980. The precision of the relationship between rainfall and the vegetation indices for the meteorological stations encourages the view that NOAA AVHRR GAC composite NDVI values can be used to monitor effective rainfall in the Savanna Zone of the Democratic Republic of Sudan     

Spatial and temporal accuracy of coarse resolution products of NOAA-AVHRR NDVI data

Coarse resolution products of maximum value composite NDVI data, derived from the AVHRR on board the NOAA-7, -9 and -11 series of polar-orbiting satellites, are becoming increasingly available. One such product is the ARTEMIS African decadal NDVT data set, available on CD-ROM for the years 1981 to 1991. These data have inherent spatial and temporal errors which arise as a result of the sampling procedures involved in their generation. This paper describes the way in which the ARTEMIS NDVI data are produced. It then describes the spatial and temporal accuracy of the data, estimated for a study area in East Africa, in terms of: 1. their spatial resolution, 2. pixel locational errors, 3. sampling biases and systematic errors introduced during the production of the data, and 4. the presence of signal noise in multi-temporal profiles of the data. The paper concludes that these errors may prove a considerable limitation to the usefulness of the data for distinguishing specific habitats on the ground over small study areas, and that in the further development of applications for these coarse resolution data sets, these errors must be taken into consideration.     

Global land monitoring from AVHRR: potential and limitations

Global Vegetation Index ( GVI) time series of visible, near-IR and thermal IR Advanced Very High Resolution Radiometer (AVHRR)weekly composite data with a 015^° spatial resolution collected from NOAA-9 and -11 satellites have been used to develop a prototype global land monitoring system. The system is based on standardized anomalies of the Normalized Difference Vegetation Index (NDVI) and channel 4 brightness temperature ( T₄ )for the period April 1985-September 1994. Processing included: post-launch updated calibration; cloud screening; filling in the cloud induced data gaps by monthly averaging and spatial interpolation; suppressing residual noise by smoothing; calculating 5-year monthly means and standard deviations of NDVI and T₄and their standardized anomalies. The derived anomalies show potential for detecting and interpreting the seasonal cycle and statistically significant interannual variability. Yet, discontinuities and residua! trends can be traced in time series of NDVI and T4. Discontinuities result from the switch from NOAA-9 to NOAA-11 in 1988, and the Mount Pinatubo eruption in 1991. Trends are a combined effect of satellite orbit drift and a possible persistent error in post-launch calibration of solar channels. The orbit drift affects the solar and thermal IR channels through systematic variation of illumination geometry and diurnal heating/cooling of the surface and atmosphere, respectively. Examples are given to illustrate the magnitude of these effects, which reduce the ability to monitor small year-to-year surface changes. The present system yields more accurate results in geographic regions, where atmospheric, angular and diurnal variability effects have a lesser impact on the derived anomalies, i.e. over vegetated areas outside the tropics during local summers. For global-scale monitoring, angular, atmospheric and diurnal variability corrections must be incorporated in the present system.     

Multi-platform comparisons of MODIS and AVHRR normalized difference vegetation index data

The relationship between AVHRR-derived normalized difference vegetation index (NDVI) values and those of future sensors is critical to continued long-term monitoring of land surface properties. The follow-on operational sensor to the AVHRR, the Visible/Infrared Imager/Radiometer Suite (VIIRS), will be very similar to the NASA Earth Observing System's Moderate Resolution Imaging Spectroradiometer (MODIS) sensor. NDVI data derived from visible and near-infrared data acquired by the MODIS (Terra and Aqua platforms) and AVHRR (NOAA-16 and NOAA-17) sensors were compared over the same time periods and a variety of land cover classes within the conterminous United States. The results indicate that the 16-day composite NDVI values are quite similar over the composite intervals of 2002 and 2003, and linear relationships exist between the NDVI values from the various sensors. The composite AVHRR NDVI data included water and cloud masks and adjustments for water vapor as did the MODIS NDVI data. When analyzed over a variety of land cover types and composite intervals, the AVHRR derived NDVI data were associated with 89% or more of the variation in the MODIS NDVI values. The results suggest that it may be possible to successfully reprocess historical AVHRR data sets to provide continuity of NDVI products through future sensor systems.     

The effect of the temporal resolution of NDVI data on season onset dates and trends across Canadian broadleaf forests

Satellite remote sensing has the potential to contribute to plant phenology monitoring at spatial and temporal scales relevant for regional and global scale studies. Historically, temporal composites of satellite data, ranging from 8 days to 16 days, have been used as a starting point for satellite-derived phenology data sets. In this study we assess how the temporal resolution of such composites affects the estimation of the start of season (SOS) by: 1) calibrating a relationship between satellite derived SOS with in situ leaf unfolding (LU) of trembling aspen (Populus tremuloides) across Canada and 2) quantifying the sensitivity of calibrated satellite SOS estimates and trends, over Canadian broadleaf forests, to the temporal resolution of NDVI data. SOS estimates and trends derived from daily NDVI data were compared to SOS estimates and trends derived from multiday NDVI composites that retain the exact date of the maximum NDVI value or that assume the midpoint of the multiday interval as the observation date. In situ observations of LU dates were acquired from the PlantWatch Canada network. A new Canadian database of cloud and snow screened daily 1-km resolution National Oceanic and Atmospheric Administration advanced very high resolution radiometer surface reflectance images was used as input satellite data. The mean absolute errors of SOS dates with respect to in situ LU dates ranged between 13 and 40 days. SOS estimates from NDVI composites that retain the exact date of the maximum NDVI value had smaller errors (~ 13 to 20 days). The sensitivity analysis reinforced these findings: SOS estimates from NDVI composites that use the exact date had smaller absolute deviations from the LU date (0 to − 5 days) than the SOS estimates from NDVI composites that use the midpoint (− 2 to − 27 days). The SOS trends between 1985 and 2007 were not sensitive to the temporal resolution or compositing methods. However, SOS trends at individual ecozones showed significant differences with the SOS trends from daily NDVI data (Taiga plains and the Pacific maritime zones). Overall, our results suggest that satellite based estimates of vegetation green-up dates should preferably use sub-sampled NDVI composites that include the exact observation date of the maximum NDVI to minimize errors in both, SOS estimates and SOS trend analyses. For trend analyses alone, any of the compositing methods could be used, preferably with composite intervals of less than 28 days. This is an important finding, as it suggests that existing long-term 10-day or 15-day NDVI composites could be used for SOS trend analyses over broadleaf forests in Canada or similar areas. Future studies will take advantage of the growing in situ phenology networks to improve the validation of satellite derived green-up dates.     

Development of a daily long term record of NOAA-14 AVHRR land surface temperature over Africa

We developed a new 6-year daily, daytime and nighttime, NOAA-14 AVHRR based land surface temperature (LST) dataset over continental Africa for the period 1995 through 2000. The processing chain was developed within the Global Inventory Modeling and Mapping System (GIMMS) at NASA's Goddard Space Flight Center. This paper describes the processing methodology used to convert the Global Area Coverage Level-1b data into LST and collateral data layers, such as sun and view geometries, cloud mask, local time of observation, and latitude and longitude. We used the Ulivieri et al. [Ulivieri, C., M.M. Castronuovo, R. Francioni, and A. Cardillo (1994), A split window algorithm for estimating land surface temperature from satellites, Adv. Space Research, 14(3):59-65.] split window algorithm to determine LST values. This algorithm requires as input values of surface emissivity in AVHRR channels 4 and 5. Thus, we developed continental maps of emissivity using an ensemble approach that combines laboratory emissivity spectra, MODIS-derived maps of herbaceous and woody fractional cover, and the UNESCO FAO soil map. A preliminary evaluation of the resulting LST product over a savanna woodland in South Africa showed a bias of < 0.3 K and an uncertainty of < 1.3 K for daytime retrievals (< 2.5 K for night). More extensive validation is required before statistically significant uncertainties can be determined. The LST production chain described here could be adapted for any wide field of view sensor (e.g., MODIS, VIIRS), and the LST product may be suitable for monitoring spatial and temporal temperature trends, or as input to many process models (e.g., hydrological, ecosystem).     

Tracking the calibration stability and consistency of the 3.7, 11.0 and 12.0 μm channels of the NOAA-KLM AVHRR with MODIS

The NOAA-KLM satellites (NOAA-15 to 18) are the current polar-orbiting operational environmental satellites (POES) that carry the Advanced Very High Resolution Radiometer (AVHRR). This study examines the calibration stability and consistency of all three infrared channels (3.7, 11.0 and 12.0 μm) of AVHRR onboard NOAA-15 to 18. The short-term stability is examined from variations of the scan-by-scan gain response, while the long-term stability and calibration consistency are examined by tracking the trends of gain response and measured scene brightness temperatures. The relative differences of observed scene brightness temperatures among NOAA-15 to 18 AVHRR are determined using MODIS as a transfer radiometer based on observations from simultaneous nadir overpasses (SNO). Results show that variations of the scan-to-scan gain responses are within 0.10% under normal operational conditions, while long-term gain changes over six years from 2001 to 2006 vary from 2 to 4% depending on channel. Long-term trending results show that total six-year drifts in observed brightness temperature from NOAA-15 to 18 AVHRR are less than 0.5 K for a given scene temperature in the 250 to 270 K range for the 3.7, 11.0 and 12.0 μm channels, respectively. The calibration consistency is examined for a scene temperature range of 220 to 290 K. The temperature biases among NOAA-16 to 18 AVHRR are within ±0.5 K for the 11.0 and 12.0 μm channels. For NOAA-15 AVHRR, biases of –2.0 K at 11.0 μm and –1.5 K at 12.0 μm are found in comparison with others at the low end of the temperature range. For the 3.7 μm channel, relative biases up to a few degrees among NOAA-15 to 18 could be found at low brightness temperatures.     

Validation of clear-sky radiances over oceans simulated with MODTRAN4.2 and global NCEP GDAS fields against nighttime NOAA15-18 and MetOp-A AVHRR data

An accurate and globally representative forward radiative transfer model (RTM) is needed to explore improvements in sea surface temperature (SST) retrievals from spaceborne infrared observations. This study evaluates the biases in top-of-atmosphere (TOA) brightness temperatures (BT) modeled with the moderate resolution transmission (MODTRAN4.2) band RTM, bounded by a Fresnel's reflective flat sea surface. This model is used to simulate global clear-sky Advanced Very High Resolution Radiometer (AVHRR) nighttime BTs from NOAA-15 through 18 and MetOp-A platforms for one full day of 18 February 2007. Inputs to RTM (SST fields and vertical profiles of atmospheric relative humidity, temperature, pressure, and geopotential height) are specified from the National Centers for Environmental Prediction's (NCEP) Global Data Assimilation System (GDAS) data. Model BTs in AVHRR channels 3B (3.7 μm), 4 (11 μm), and 5 (12 μm) are then compared with their respective measured counterparts, available in the NESDIS operational SST files. Ideally, the RTM should match the observations, but in fact, the modeled BTs are biased high with respect to the AVHRR BTs. The “Model minus Observation” (M − O) bias ranges from about 0 to 2 K, depending upon spectral band, view zenith angle, and sea and atmosphere state at the retrieval point. The bias asymptotically decreases towards confidently clear-sky conditions, but it never vanishes and invariably shows channel-specific dependencies on view zenith angle and geophysical conditions (e.g., column water vapor and sea-air temperature difference). Fuller exploration of the potential of the current RTM (e.g., adding global vertical aerosol profiles) or improvements to its input (NCEP SST and atmospheric profiles) may reduce this bias, but they cannot fully reconcile its spectral and angular structure. The fact that the M − O biases are closely reproducible for five AVHRR sensors flown onboard different platforms adds confidence in the validation approach employed in this study. We emphasize the need for establishing a globally adequate forward RTM for the use in SST modeling and retrievals. A first test of the RTM adequacy is its ability, when used in conjunction with the global fields from the numerical weather prediction models, to reproduce the TOA clear-sky radiances measured by satellite sensors.     

Evaluating different NDVI composite techniques using NOAA-14 AVHRR data

Normalized Difference Vegetation Index (NDVI) data derived from Advanced Very High Resolution Radiometer (AVHRR) data are influenced by cloud contamination, which is common in individual AVHRR scenes. Maximum value compositing (MVC) of NDVI data has been employed to minimize cloud contamination. Two types of weekly NDVI composites were built for crop seasons in summer: one from all available AVHRR data (named the traditional NDVI composite) and the other from solely cloud-free AVHRR data (named the conditional NDVI composite). The MVC method was applied to both composites. The main objective of this study was to compare the two types of NDVI composites using Texas data. The NDVI seasonal profiles produced from the conditional NDVI composites agreed with the field measured leaf area index (LAI) data, reaching maximum values at similar times. However, the traditional NDVI composites showed irregular patterns, primarily due to cloud contamination. These study results suggest that cloud detection for individual AVHRR scenes should be strongly recommended before producing weekly NDVI composites. Appropriate AVHRR data pre-processing is important for composite products to be used for short-term vegetation condition and biomass studies, where the traditional NDVI composite data do not eliminate cloud-contaminated pixels. In addition, this study showed that atmosphere composition affected near-infrared reflectance more than visible reflectance. The near-infrared reflectance was increasingly adjusted through atmospheric correction.     

Global evaluation of four AVHRR-NDVI data sets: Intercomparison and assessment against Landsat imagery

Advanced Very High Resolution Radiometer (AVHRR) data with their long-term (1981-current) global coverage at frequent intervals provide unique opportunities to explore vegetation dynamics related to climate variability, climate change, and land-use driven changes of land cover. Several AVHRR-derived Normalized Difference Vegetation Index (NDVI) data sets exist, each based on the AVHRR Global Area Coverage archive but differing in their processing to correct for sensor and atmospheric effects. This paper presents a global comparative analysis for the land surface involving four AVHRR-derived NDVI data sets: (1) Pathfinder AVHRR Land (PAL); (2) Global Inventory Modeling and Mapping Studies (GIMMS); (3) Land Long Term Data Record (LTDR) version 3 (V3); and (4) Fourier-Adjustment, Solar zenith angle corrected, Interpolated Reconstructed (FASIR). Our aims are two-fold: (1) to assess the level of agreement of the medians, trends, and variances, as well as the correlation between the four AVHRR-NDVI data sets from 1982 to 1999; and (2) to independently assess the performance of each AVHRR-NDVI data set, and that of Moderate Resolution Imaging Spectroradiometer (MODIS) NDVI, using 11,764 Landsat samples of 20 × 20 km² located globally covering every major land-cover type. For the AVHRR-NDVI intercomparison equal medians, variance, and trends, and no correlation between all the respective AVHRR-NDVI data sets were found for 9.9%, 45.5%, 48.1% and 38.4% of the total land surface, respectively (p ≥ 0.05). For the four AVHRR-NDVI data sets we found: (1) consistent trends for the tundra and particularly Australia; (2) inconsistent trends for Europe, Africa, and the Sahel; and (3) moderately consistent trends for the rest of the terrestrial land surface including North America and China. The PAL and LTDR V3 data sets lack calibration, as evidenced by the presence of apparent trends in desert areas. In the Landsat-NDVI vs. AVHRR-NDVI comparison of absolute values the LTDR V3 data set performed best, whereas in the comparison of temporal-change values the GIMMS data set performed best. In both analyses MODIS-NDVI performed better than any AVHRR-NDVI data set. The simple average of the four AVHRR-NDVI data sets produced better results than either AVHRR-NDVI data set alone, indicating that the errors between the data sets are at least partially unrelated. This research emphasizes the implications of AVHRR-NDVI data set choice for studies assessing the vegetation response to climate change and modeling of the terrestrial carbon balance.     

Vegetation greening in Spain detected from long term data (1981–2015)

ABSTRACT

This study describes a newly developed high-resolution (1.1 km) Normalized Difference Vegetation Index dataset for the peninsular Spain and the Balearic Islands (Sp_1km_NDVI). This dataset is developed based on National Oceanic and Atmospheric Administration–Advanced Very High Resolution Radiometer (NOAA–AVHRR) afternoon images, spanning the past three decades (1981–2015). After a careful pre-processing procedure, including calibration with post-launch calibration coefficients, geometric and topographic corrections, cloud removal, temporal filtering, and bi-weekly composites by maximum NDVI-value, we assessed changes in vegetation greening over the study domain using Mann-Kendall and Theil-Sen statistics. Our trend results were compared with those derived from some widely recognized global NDVI datasets [e.g. the Global Inventory Modelling and Mapping Studies 3rd generation (GIMMS3g), Smoothed NDVI (SMN) and Moderate-Resolution Imaging Spectroradiometer (MODIS)]. Results demonstrate that there is a good agreement between the annual trends based on Sp_1km_NDVI product and other datasets. Nonetheless, we found some differences in the spatial patterns of the NDVI trends at the seasonal scale. Overall, in comparison to the available global NDVI datasets, Sp_1km_NDVI allows for characterizing changes in vegetation greening at a more-detailed spatial and temporal scale. In specific, our dataset provides relatively long-term corrected satellite time series (>30 years), which are crucial to understand the response of vegetation to climate change and human-induced activities. Also, given the complex spatial structure of NDVI changes over the study domain, particularly due to the rapid land intensification processes, the spatial resolution (1.1 km) of our dataset can provide detailed spatial information on the inter-annual variability of vegetation greening in this Mediterranean region and assess its links to climate change and variability.     

Monitoring the phenology of Tunisian grazing lands

Abstract

Normalized difference vegetation index (NDVI) data obtained from the Advanced Very High Resolution Radiometer (AVHRR) on board NOAA-9 have been analysed to assess their utility for monitoring the vegetation of Tunisian grazing lands. Preliminary analysis shows that the NDVI provides a sensitive indicator of monthly variations in biomass which correlate with spatial and temporal changes in growing conditions. Investigations suggest that the percentage contribution of the soil background to total recorded reflectance, provides an important limiting factor to the sensitivity of the NDVI, creating a threshold beyond which the accuracy of this index becomes less reliable.     

Navigation of AVHRR images using ARGOS or TBUS bulletins

Abstract

An automatic AVHRR navigation method has been tested using ARGOS or TBUS orbital elements as input data. Satisfactory results have been obtained on 16 NOAA-9 and NOAA-10 orbits between 19 June 1988 and 26 September 1988: a mean error of 3·5km and a standard deviation of this error of 1·2km with ARGOS, a mean error of 2·3km and a standard deviation of 0·6km with TBUS. The positioning error consists of a pixelerror which is systematic for each satellite (?1·6 pixels for NOAA-9 and 2·0 pixels for NOAA-10) and a variable line error which is mainly due to the positioning error of the satellite on its trajectory.