On intra-annual EVI variability in the dry season of tropical forest: A case study with MODIS and hyperspectral data

Considerable controversy is associated with dry season increases in the Enhanced Vegetation Index (EVI), observed using the Moderate Resolution Imaging Spectroradiometer (MODIS), compared with field-based estimates of decreasing plant productivity. Here, we investigate potential causes of intra-annual variability by comparing EVI from mature forest with field-measured Leaf Area Index (LAI) to validate space-based observations. EVI was calculated from 19 nadir and off-nadir Hyperion images in the 2005 dry season, and inspected for consistency with MODIS observations from 2004 to 2009. The objective was to evaluate the possible influence of the view-illumination geometry and of canopy foliage and leaf flush on the EVI. Spectral mixture models were used to evaluate the relationship between EVI and the shade fraction, a measure that varies with pixel brightness. MODIS LAI values were compared with LAI estimated using hemispherical photographs taken in two field campaigns in the dry season. To keep LAI and leaf flush conditions as constant variables and vary solar illumination, we used airborne Hyperspectral Mapper (Hymap) data acquired over mature forest from another region on the same day but with two distinct solar zenith angles (SZA) (29° and 53°). Results showed that intra-annual variability in MODIS and nadir Hyperion EVI in the dry season of tropical forest were driven by solar illumination effects rather than changes in LAI. The reflectance of the MODIS and Hyperion blue, red and near infrared (NIR) bands was higher at the end of the dry season because of the predominance of sunlit canopy components for the sensors due to decreasing SZA from June (44°) to September (26°). Because EVI was highly correlated with the reflectance of the NIR band used to generate it (r of + 0.98 for MODIS and + 0.88 for Hyperion), this vegetation index followed the general NIR pattern, increasing with smaller SZA towards the end of the dry season. Hyperion EVI was inversely correlated with the shade fraction (r = − 0.93). Changes in canopy foliage detected from MODIS LAI data were not consistent with LAI estimates from hemispherical photographs. Although further research is necessary to measure the impact of leaf flush on intra-annual EVI variability in the Querência region, analysis of Hymap data with fixed LAI and leaf flush conditions confirmed the influence of the illumination effects on the EVI.     

An assessment of photosynthetic light use efficiency from space: Modeling the atmospheric and directional impacts on PRI reflectance

Estimation of photosynthetic light use efficiency (ε) from satellite observations is an important component of climate change research. The photochemical reflectance index, a narrow waveband index based on the reflectance at 531 and 570 nm, allows sampling of the photosynthetic activity of leaves; upscaling of these measurements to landscape and global scales, however, remains challenging. Only a few studies have used spaceborne observations of PRI so far, and research has largely focused on the MODIS sensor. Its daily global coverage and the capacity to detect a narrow reflectance band at 531 nm make it the best available choice for sensing ε from space. Previous results however, have identified a number of key issues with MODIS-based observations of PRI. First, the differences between the footprint of eddy covariance (EC) measurements and the MODIS footprint, which is determined by the sensor's observation geometry make a direct comparison between both data sources challenging and second, the PRI reflectance bands are affected by atmospheric scattering effects confounding the existing physiological signal. In this study we introduce a new approach for upscaling EC based ε measurements to MODIS. First, EC-measured ε values were “translated” into a tower-level optical PRI signal using AMSPEC, an automated multi-angular, tower-based spectroradiometer instrument. AMSPEC enabled us to adjust tower-measured PRI values to the individual viewing geometry of each MODIS overpass. Second, MODIS data were atmospherically corrected using a Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, which uses a time series approach and an image-based rather than pixel-based processing for simultaneous retrievals of atmospheric aerosol and surface bidirectional reflectance (BRDF). Using this approach, we found a strong relationship between tower-based and spaceborne reflectance measurements (r² = 0.74, p < 0.01) throughout the vegetation period of 2006. Swath (non-gridded) observations yielded stronger correlations than gridded data (r² = 0.58, p < 0.01) both of which included forward and backscatter observations. Spaceborne PRI values were strongly related to canopy shadow fractions and varied with different levels of ε. We conclude that MAIAC-corrected MODIS observations were able to track the site-level physiological changes from space throughout the observation period.     

Evaluation of satellite based primary production modelling in the semi-arid Sahel

The MODIS (Moderate Resolution Imaging Spectroradiometer) primary productivity products are evaluated against observed Above-ground Net Primary Production (AGNPP) in the semi-arid Senegal 2001. MODIS net primary productivity (NPP) modelling is a light use efficiency (LUE) based approach incorporating constraints on vegetation productivity arising from simulated radiation, water demand and temperature data from NASA's Data Assimilation Office (DAO). Annually integrated MODIS PSN (MOD17A2 net photosynthesis, Collection 4) explains more of the observed biomass variation (r² = 0.77) than MODIS fAPAR (fraction Absorbed Photosynthetically Active Radiation, Collection 4) (r² = 0.72), indicating the effect of including the canopy stress scalar (ε_s) based on DAO data combined with modelled maintenance respiration costs (of leaf and fine roots). Annual MODIS NPP (MOD17A3, Collection 4 (C4) and Collection 4.5 (C4.5)) including growth respiration and live wood maintenance respiration costs and modified DAO input (C4.5) however increases the residual unexplained observed AGNPP variance (C4 NPP; r² = 0.49) (C4.5 NPP; r² = 0.37). The overall quality of the annual NPP MODIS C4 and C4.5 products are moderate for the semi-arid Senegal because of the annual respiration cost modelling and a change in C4.5 biome-specific parameters stored in a Biome Properties Look-Up Table (BPLUT) is the main contributor to the observed discrepancy between C4 and C4.5 NPP. The dynamic range of the values of all MOD17 products was too low when compared to observed AGNPP. An estimate of canopy water stress (SIWSI; Shortwave Infrared Water Stress Index) derived from MODIS channels 2 and 6 and photosynthetically active radiation (PAR) irradiance derived from geostationary METEOSAT data were tested for primary production modelling using a stepwise linear regression analysis. PAR irradiance was combined with MODIS fAPAR into APAR (Absorbed Photosynthetically Active Radiation) explaining 79% of the observed AGNPP variation. Introducing SIWSI significantly increased the explained variance of observed AGNPP (r² = 0.89). MODIS-derived percentage tree cover was tested as a predictor based on the hypothesis that tree cover provides information on differences in respiratory costs between trees and grasses thereby accounting for variations in the LUE conversion efficiency ε. No significant reduction in residual unexplained AGNPP variance was found. Earth observation based derivation of PAR and canopy water stress from SIWSI suggest potential improvements to primary production models in semi-arid biomes that can be implemented in general NPP modelling LUE methodology.     

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).     

Generation of dense time series synthetic Landsat data through data blending with MODIS using a spatial and temporal adaptive reflectance fusion model

Landsat imagery with a 30 m spatial resolution is well suited for characterizing landscape-level forest structure and dynamics. While Landsat images have advantageous spatial and spectral characteristics for describing vegetation properties, the Landsat sensor's revisit rate, or the temporal resolution of the data, is 16 days. When considering that cloud cover may impact any given acquisition, this lengthy revisit rate often results in a dearth of imagery for a desired time interval (e.g., month, growing season, or year) especially for areas at higher latitudes with shorter growing seasons. In contrast, MODIS (MODerate-resolution Imaging Spectroradiometer) has a high temporal resolution, covering the Earth up to multiple times per day, and depending on the spectral characteristics of interest, MODIS data have spatial resolutions of 250 m, 500 m, and 1000 m. By combining Landsat and MODIS data, we are able to capitalize on the spatial detail of Landsat and the temporal regularity of MODIS acquisitions. In this research, we apply and demonstrate a data fusion approach (Spatial and Temporal Adaptive Reflectance Fusion Model, STARFM) at a mainly coniferous study area in central British Columbia, Canada. Reflectance data for selected MODIS channels, all of which were resampled to 500 m, and Landsat (at 30 m) were combined to produce 18 synthetic Landsat images encompassing the 2001 growing season (May to October). We compared, on a channel-by-channel basis, the surface reflectance values (stratified by broad land cover types) of four real Landsat images with the corresponding closest date of synthetic Landsat imagery, and found no significant difference between real (observed) and synthetic (predicted) reflectance values (mean difference in reflectance: mixed forest x? = 0.086, σ = 0.088, broadleaf x? = 0.019, σ = 0.079, coniferous x? = 0.039, σ = 0.093). Similarly, a pixel based analysis shows that predicted and observed reflectance values for the four Landsat dates were closely related (mean r² = 0.76 for the NIR band; r² = 0.54 for the red band; p < 0.01). Investigating the trend in NDVI values in synthetic Landsat values over a growing season revealed that phenological patterns were well captured; however, when seasonal differences lead to a change in land cover (i.e., disturbance, snow cover), the algorithm used to generate the synthetic Landsat images was, as expected, less effective at predicting reflectance.     

Seasonal and inter-annual variation in view angle effects on MODIS vegetation indices at three forest sites

The bi-directional reflectance distribution function (BRDF) has been widely studied across different vegetation types. However, these studies generally report values for only one point in time. We were interested in the potential for seasonal and inter-annual variation in BRDF parameters. NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on board the EOS satellites has now been collecting data for 10 years. Since BRDF parameters are reported for the individual spectral bands, these data can be used to examine intra-annual variation. However, MODIS BRDF parameters are not calculated for the various vegetation indices which are derived from the spectral bands. Our objective in this study was to use the 10 years of MODIS data now available to examine seasonal and inter-annual variation in the view angle sensitivity of three vegetation indices; the normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), and the photochemical reflectance index (PRI) at 3 flux tower sites (Harvard Forest, Howland Forest and Morgan Monroe State Forest). For these 3 sites, only EVI was significantly affected by view angle. There was also a substantial variation in the view angle sensitivity of EVI across seasons and this variation was different for backscatter vs. forward scatter data. It is possible that differences in the scattering of radiation between the spring and the fall are responsible for the seasonal difference in view angle sensitivity. There were also complimentary variations in forward and backscatter view angle sensitivity of EVI across years. The greater view angle sensitivity of EVI, as opposed to NDVI, suggests that greater care must be taken to correct for BRDF effects when using this vegetation index.     

An assessment of the black ocean pixel assumption for MODIS SWIR bands

Recent studies show that an atmospheric correction algorithm using shortwave infrared (SWIR) bands improves satellite-derived ocean color products in turbid coastal waters. In this paper, the black pixel assumption (i.e., zero water-leaving radiance contribution) over the ocean for the Moderate Resolution Imaging Spectroradiometer (MODIS) SWIR bands at 1240, 1640, and 2130 nm is assessed for various coastal ocean regions. The black pixel assumption is found to be generally valid with the MODIS SWIR bands at 1640 and 2130 nm even for extremely turbid waters. For the MODIS 1240 nm band, however, ocean radiance contribution is generally negligible in mildly turbid waters such as regions along the U.S. east coast, while some slight radiance contributions are observed in extremely turbid waters, e.g., some regions along the China east coast, the estuary of the La Plata River. Particularly, in the Hangzhou Bay, the ocean radiance contribution at the SWIR band 1240 nm results in an overcorrection of atmospheric and surface effects, leading to errors of MODIS-derived normalized water-leaving radiance at the blue reaching ~ 0.5 mW cm^− 2 μm^− 1 sr^− 1. In addition, we found that, for non-extremely turbid waters, i.e., the ocean contribution at the near-infrared (NIR) band < ~ 1.0 mW cm^− 2 μm^− 1 sr^− 1, there exists a good relationship in the regional normalized water-leaving radiances between the red and the NIR bands. Thus, for non-extremely turbid waters, such a red-NIR radiance relationship derived regionally can possibly be used for making corrections for the regional NIR ocean contributions without using the SWIR bands, e.g., for atmospheric correction of ocean color products derived from the Sea-viewing Wide Field-of-view Sensor (SeaWiFS).     

Validation of MODIS derived aerosol optical depth over the Yangtze River Delta in China

MODIS derived aerosol optical depths (AODs) at 550 nm are compared with sunphotometer CE318 measurements at 7 sites located at Yangtze River Delta (YRD) in China from July to October, 2007. The evaluation result indicates that MODIS AODs (Collection 5, C005) are in good agreement with those from CE318 in dense vegetation regions, but show more differences in those regions with complex underlying surface (such as at lake water and urban surface sites). Reasons for these differences are discussed after removing cases with significant errors caused by validation scheme. The final validation result shows that MODIS AODs are in good agreement with CE318 with a correlation coefficient of 0.85 and RMS of 0.15. 90% of MODIS cases fall in the range of Δτ = ± 0.05 ± 0.20τ, indicating MODIS aerosol retrieval algorithm, aerosol models and surface reflectance estimate are generally suitably reasonable for aerosol retrieval in YRD. However, MODIS AODs show a systemic errors with fitted line of y = 0.75x + 0.13, indicating underestimation of AOD when aerosol loadings are high. Aerosol models and surface reflectance estimations are dominant sources of MODIS aerosol retrieval errors.     

Evaluation of the MODIS LAI product using independent lidar-derived LAI: A case study in mixed conifer forest

This study presents an alternative assessment of the MODIS LAI product for a 58,000 ha evergreen needleleaf forest located in the western Rocky Mountain range in northern Idaho by using lidar data to model (R² = 0.86, RMSE = 0.76) and map LAI at higher resolution across a large number of MODIS pixels in their entirety. Moderate resolution (30 m) lidar-based LAI estimates were aggregated to the resolution of the 1-km MODIS LAI product and compared to temporally-coincident MODIS retrievals. Differences in the MODIS and lidar-derived values of LAI were grouped and analyzed by several different factors, including MODIS retrieval algorithm, sun/sensor geometry, and sub-pixel heterogeneity in both vegetation and terrain characteristics. Of particular interest is the disparity in the results when MODIS LAI was analyzed according to algorithm retrieval class. We observed relatively good agreement between lidar-derived and MODIS LAI values for pixels retrieved with the main RT algorithm without saturation for LAI LAI ≤ 4. Moreover, for the entire range of LAI values, considerable overestimation of LAI (relative to lidar-derived LAI) occurred when either the main RT with saturation or back-up algorithm retrievals were used to populate the composite product regardless of sub-pixel vegetation structural complexity or sun/sensor geometry. These results are significant because algorithm retrievals based on the main radiative transfer algorithm with or without saturation are characterized as suitable for validation and subsequent ecosystem modeling, yet the magnitude of difference appears to be specific to retrieval quality class and vegetation structural characteristics.     

MODIS spectral signals at a flux tower site: Relationships with high-resolution data, and CO2 flux and light use efficiency measurements

In this study, spectral indices were calculated from single date HyMap (3 m; 126 bands), Hyperion (30 m; 242 bands), ASTER (15/30 m; 9 bands), and a time series of MODIS nadir BRDF-adjusted reflectance (NBAR; 1 km, 7 bands) for a study area surrounding the Tumbarumba flux tower site in eastern Australia. The study involved: a) the calculation of a range of physiologically-based vegetation indices from ASTER, HyMap, Hyperion and MOD43B NBAR imagery over the flux tower site; b) comparison across scales between HyMap, Hyperion and MODIS for the normalized difference water index (NDWI) and the Red-Green ratio; c) analysis of relationships between tower-based flux and light use efficiency (LUE) measurements and seasonal and climatic constraints on growth; and d) examination of relationships between fluxes, LUE and time series of NDVI, NDWI and Red-Green ratio. Strong seasonal patterns of variation were observed in NDWI and Red/Green ratio from MODIS NBAR. Correlations between fine (3 and 30 m) and coarse (1 km) scale indices for a small region around the flux tower site were moderately good for Red/Green ratio, but poor for NDWI. Hymap NDWI values for the understorey canopy were much lower than values for the tree canopy. MODIS NDWI was negatively correlated with CO₂ fluxes during warm and cool seasons. The correlation indicated that surface reflectance, affected by a spectrally bright grassland understorey canopy, was decoupled from growth of trees with access to deep soil moisture. The application of physiologically-based indices at earth observation scale requires careful attention to applicability of band configurations, contribution of vegetation components to reflectance signals, mechanistic relationships between biochemical processes and spectral indexes, and incorporation of ancillary information into any analysis.     

Ground measurements for the validation of land surface temperatures derived from AATSR and MODIS data

An experimental site was set up in a large, flat and homogeneous area of rice crops for the validation of satellite derived land surface temperature (LST). Experimental campaigns were held in the summers of 2002-2004, when rice crops show full vegetation cover. LSTs were measured radiometrically along transects covering an area of 1 km². A total number of four thermal radiometers were used, which were calibrated and inter-compared through the campaigns. Radiometric temperatures were corrected for emissivity effects using field emissivity and downwelling sky radiance measurements. A database of ground-based LSTs corresponding to morning, cloud-free overpasses of Envisat/Advanced Along-Track Scanning Radiometer (AATSR) and Terra/Moderate Resolution Imaging Spectroradiometer (MODIS) is presented. Ground LSTs ranged from 25 to 32 °C, with uncertainties between ± 0.5 and ± 0.9 °C. The largest part of these uncertainties was due to the spatial variability of surface temperature. The database was used for the validation of LSTs derived from the operational AATSR and MODIS split-window algorithms, which are currently used to generate the LST product in the L2 level data. A quadratic, emissivity dependent split-window equation applicable to both AATSR and MODIS data was checked as well. Although the number of cases analyzed is limited (five concurrences for AATSR and eleven for MODIS), it can be concluded that the split-window algorithms work well, provided that the characteristics of the area are adequately prescribed, either through the classification of the land cover type and the vegetation cover, or with the surface emissivity. In this case, the AATSR LSTs yielded an average error or bias of − 0.9 °C (ground minus algorithm), with a standard deviation of 0.9 °C. The MODIS LST product agreed well with the ground LSTs, with differences comparable or smaller than the uncertainties of the ground measurements for most of the days (bias of + 0.1 °C and standard deviation of 0.6 °C, for cloud-free cases and viewing angles smaller than 60°). The quadratic split-window algorithm resulted in small average errors (+ 0.3 °C for AATSR and 0.0 °C for MODIS), with differences not exceeding ± 1.0 °C for most of the days (standard deviation of 0.9 °C for AATSR and 0.5 °C for MODIS).     

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.     

Can a satellite-derived estimate of the fraction of PAR absorbed by chlorophyll (FAPARchl) improve predictions of light-use efficiency and ecosystem photosynthesis for a boreal aspen forest?

We used daily MODerate resolution Imaging Spectroradiometer (MODIS) imagery obtained over a five-year period to analyze the seasonal and inter-annual variability of the fraction of absorbed photosynthetically active radiation (FAPAR) and photosynthetic light use efficiency (LUE) for the Southern Old Aspen (SOA) flux tower site located near the southern limit of the boreal forest in Saskatchewan, Canada. To obtain the spectral characteristics of a standardized land area to compare with tower measurements, we scaled up the nominal 500 m MODIS products to a 2.5 km × 2.5 km area (5 × 5 MODIS 500 m grid cells). We then used the scaled-up MODIS products in a coupled canopy-leaf radiative transfer model, PROSAIL-2, to estimate the fraction of absorbed photosynthetically active radiation (APAR) by the part of the canopy dominated by chlorophyll (FAPAR_chl) versus that by the whole canopy (FAPAR_canopy). Using the additional information provided by flux tower-based measurements of gross ecosystem production (GEP) and incident PAR, we determined 90-minute averages for APAR and LUE (slope of GEP:APAR) for both the physiologically active foliage (APAR_chl, LUE_chl) and for the entire canopy (APAR_canopy, LUE_canopy).The flux tower measurements of GEP were strongly related to the MODIS-derived estimates of APAR_chl (r² = 0.78) but only weakly related to APAR_canopy (r² = 0.33). Gross LUE between 2001 and 2005 for LUE_chl was 0.0241 µmol C µmol^− 1 PPFD whereas LUE_canopy was 36% lower. Time series of the 5-year normalized difference vegetation index (NDVI) were used to estimate the average length of the core growing season as days of year 152-259. Inter-annual variability in the core growing season LUE_chl (µmol C µmol^− 1 PPFD) ranged from 0.0225 in 2003 to 0.0310 in 2004. The five-year time series of LUE_chl corresponded well with both the seasonal phase and amplitude of LUE from the tower measurements but this was not the case for LUE_canopy. We conclude that LUE_chl derived from MODIS observations could provide a more physiologically realistic parameter than the more commonly used LUE_canopy as an input to large-scale photosynthesis models.     

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.     

Mapping paddy rice agriculture in South and Southeast Asia using multi-temporal MODIS images

In this paper, we developed a new geospatial database of paddy rice agriculture for 13 countries in South and Southeast Asia. These countries have ∼ 30% of the world population and ∼ 2/3 of the total rice land area in the world. We used 8-day composite images (500-m spatial resolution) in 2002 from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard the NASA EOS Terra satellite. Paddy rice fields are characterized by an initial period of flooding and transplanting, during which period a mixture of surface water and rice seedlings exists. We applied a paddy rice mapping algorithm that uses a time series of MODIS-derived vegetation indices to identify the initial period of flooding and transplanting in paddy rice fields, based on the increased surface moisture. The resultant MODIS-derived paddy rice map was compared to national agricultural statistical data at national and subnational levels. Area estimates of paddy rice were highly correlated at the national level and positively correlated at the subnational levels, although the agreement at the national level was much stronger. Discrepancies in rice area between the MODIS-derived and statistical datasets in some countries can be largely attributed to: (1) the statistical dataset is a sown area estimate (includes multiple cropping practices); (2) failure of the 500-m resolution MODIS-based algorithm in identifying small patches of paddy rice fields, primarily in areas where topography restricts field sizes; and (3) contamination by cloud. While further testing is needed, these results demonstrate the potential of the MODIS-based algorithm to generate updated datasets of paddy rice agriculture on a timely basis. The resultant geospatial database on the area and spatial distribution of paddy rice is useful for irrigation, food security, and trace gas emission estimates in those countries.     

Improved algorithm for MODIS satellite retrievals of aerosol optical thickness over land in dusty atmosphere: Implications for air quality monitoring in China

A new algorithm, using the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite reflectance and aerosol single scattering properties simulated from a chemistry transport model (GEOS-Chem), is developed to retrieve aerosol optical thickness (AOT) over land in China during the spring dust season. The algorithm first uses a “dynamic lower envelope” approach to sample the MODIS dark-pixel reflectance data in low AOT conditions, to derive the local surface visible (0.65 μm)/near infrared (NIR, 2.1 μm) reflectance ratio. Joint retrievals of AOT at 0.65 μm and surface reflectance at 2.1 μm are then performed, based on the time, location, and spectral-dependent single scattering properties of the dusty atmosphere as simulated by the GEOS-Chem. A linearized vector radiative transfer model (VLIDORT) that simultaneously computes the top-of-atmosphere reflectance and its Jacobian with respect to AOT, is used in the forward component of the inversion of MODIS reflectance to AOT. Comparison of retrieved AOT results in April and May of 2008 with AERONET observations shows a strong correlation (R = 0.83), with small bias (0.01), and small RMSE (0.17); the figures are a substantial improvement over corresponding values obtained with the MODIS Collection 5 AOT algorithm for the same study region and time period. The small bias is partially due to the consideration of dust effect at 2.1 μm channel, without which the bias is − 0.05. The surface PM₁₀ (particulate matter with diameter less than 10 μm) concentrations derived using this improved AOT retrieval show better agreement with ground observations than those derived from GEOS-Chem simulations alone, or those inferred from the MODIS Collection 5 AOT. This study underscores the value of using satellite reflectance to improve the air quality modeling and monitoring.     

Satellite-based and mesoscale regression modeling of monthly air and soil temperatures over complex terrain in Turkey

Simple regression algorithms were developed to quantify spatio-temporal dynamics of minimum and maximum air temperatures (T_min and T_max, respectively) and soil temperature for a depth of 0-5 cm (T_soil-5cm) across complex terrain in Turkey using Moderate Resolution Imaging Spectroradiometer (MODIS) data at a 500-m resolution. A total of 762 16-day MODIS composites (127 images × 6 bands) between 2000 and 2005 were averaged over a monthly basis to temporally match monthly T_min, T_max, and T_soil-5cm from 83 meteorological stations. A total of 60 (28 temporally averaged plus 32 time series-based) linear regression models of T_min, T_max, and T_soil-5cm were developed using best subsets procedure as a function of a combination of 12 explanatory variables: six MODIS bands of blue, red, near infrared (NIR), middle infrared (MIR), normalized difference vegetation index (NDVI), and enhanced vegetation index (EVI); four geographical variables of latitude, longitude, altitude, and distance to sea (DtS); and two temporal variables of month, and year. The best multiple linear regression models elucidated 65% (RMSE = 5.9 °C), 65% (RMSE = 5.1 °C), and 57% (RMSE = 6.9 °C) of variations in T_min, T_max, and T_soil-5cm, respectively, under a wide range of T_min (−34 to 25 °C), T_max (0.2-47 °C) and T_soil-5cm (−9 to 40 °C) observed at the 83 stations.     

Monitoring fall foliage coloration dynamics using time-series satellite data

Fall foliage coloration is a phenomenon that occurs in many deciduous trees and shrubs worldwide. Measuring the phenology of fall foliage development is of great interest for climate change, the carbon cycle, ecology, and the tourist industry; but little effort has been devoted to monitoring the regional fall foliage status using remotely-sensed data. This study developed an innovative approach to monitoring fall foliage status by means of temporally-normalized brownness derived from MODIS (Moderate Resolution Imaging Spectroradiometer) data. Specifically, the time series of the MODIS Normalized Difference Vegetation Index (NDVI) was smoothed and functionalized using a sigmoidal model to depict the continuous dynamics of vegetation growth. The modeled temporal NDVI trajectory during the senescent phase was further combined with the mixture modeling to deduce the temporally-normalized brownness index which was independent of the surface background, vegetation abundance, and species composition. This brownness index was quantitatively linked with the fraction of colored and fallen leaves in order to model the fall foliage coloration status. This algorithm was tested by monitoring the fall foliage coloration phase using MODIS data in northeastern North America from 2001 to 2004. The MODIS-derived timing of foliage coloration phases was compared with in-situ measurements, which showed an overall absolute mean difference of less than 5 days for all foliage coloration phases and about 3 days for near peak coloration and peak coloration. This suggested that the fall foliage coloration phase retrieved from the temporally-normalized brownness index was qualitatively realistic and repeatable.     

Estimating fractional cover of photosynthetic vegetation, non-photosynthetic vegetation and bare soil in the Australian tropical savanna region upscaling the EO-1 Hyperion and MODIS sensors

Quantitative estimation of fractional cover of photosynthetic vegetation (f_PV), non-photosynthetic vegetation (f_NPV) and bare soil (f_BS) is critical for natural resource management and for modeling carbon dynamics. Accurate estimation of fractional cover is especially important for monitoring and modeling savanna systems, subject to highly seasonal rainfall and drought, grazing by domestic and native animals, and frequent burning. This paper describes a method for resolving f_PV, f_NPV and f_BS across the ~ 2 million km² Australian tropical savanna zone with hyperspectral and multispectral imagery. A spectral library compiled from field campaigns in 2005 and 2006, together with three EO-1 Hyperion scenes acquired during the 2005 growing season were used to explore the spectral response space for f_PV, f_NPV and f_BS. A linear unmixing approach was developed using the Normalized Difference Vegetation Index (NDVI) and the Cellulose Absorption Index (CAI). Translation of this approach to MODerate resolution Imaging Spectroradiometer (MODIS) scale was assessed by comparing multiple linear regression models of NDVI and CAI with a range of indices based on the seven MODIS bands in the visible and shortwave infrared region (SWIR) using synthesized MODIS surface reflectance data on the same dates as the Hyperion acquisitions. The best resulting model, which used NDVI and the simple ratio of MODIS bands 7 and 6 (SWIR3/SWIR2), was used to generate a time series of fractional cover from 16 day MODIS nadir bidirectional reflectance distribution function-adjusted reflectance (NBAR) data from 2000-2006. The results obtained with MODIS NBAR were validated against grass curing measurement at ten sites with good agreement at six sites, but some underestimation of f_NPV proportions at four other sites due to substantial sub-pixel heterogeneity. The model was also compared with remote sensing measurements of fire scars and showed a good matching in the spatio-temporal patterns of grass senescence and posterior burning. The fractional cover profiles for major grassland cover types showed significant differences in relative proportions of f_PV, f_NPV and f_BS, as well as large intra-annual seasonal variation in response to monsoonal rainfall gradients and soil type. The methodology proposed here can be applied to other mixed tree-grass ecosystems across the world.     

Regional mapping of gross light-use efficiency using MODIS spectral indices

Direct estimation of photosynthetic light-use efficiency (LUE) from space would be of significant benefit to LUE-based models which use inputs from remote sensing to estimate terrestrial productivity. The Photochemical Reflectance Index (PRI) has shown promise in tracking LUE at the leaf- to small canopy levels, but its use at regional to global scales still remains a challenge. In this study, we used different formulations of PRI calculated from the MODIS ocean band centered at 531 nm and a set of alternative reference bands at 488, 551, and 678 nm to explore the relationship between PRI and LUE where LUE was measured at eight eddy covariance flux towers located in the boreal forest of Saskatchewan, Canada. The magnitude and variability of LUE was significantly lower at the times when useful MODIS ocean band images were available (i.e. around midday under clear-sky conditions) relative to the rest of the growing season. PRI₆₇₈ (reference band at 678 nm) showed the strongest relationship (r² = 0.70) with LUE_90a (i.e. 90-minute mean LUE calculated using Absorbed Photosynthetically Active Radiation, APAR), but only when all sites were combined. Overall, the relationships between the MODIS PRIs and LUE_90a were always stronger for observations closer to the backscatter direction and there were no significant differences in the strength of the correlations whether LUE was calculated based on incident PAR or on APAR. Predictions of ecosystem photosynthesis at the time of the MODIS overpasses were significantly improved by multiplying either PAR or APAR by MODIS PRI (r² improved from 0.09 to 0.44 and 0.54 depending on the PRI formulation).We used our PRI-LUE model to create a regional LUE_90a map for the three cover types covering 47,500 km² around the flux sites. The MODIS PRI-derived LUE_90a map appeared to capture more realistic spatial heterogeneity of LUE across the landscape compared to a daily LUE map derived using the look-up table in the MODIS GPP (MOD17) algorithm. While our LUE map is only a snapshot of minimum regional LUE_90a values, with appropriate gap-filling methods it could be used to improve regional-scale monitoring of GPP. Moreover, the strong relationship between midday and daily LUE on clear days (r² = 0.93) indicates that instantaneous MODIS observations of LUE_90a could be used to estimate daily LUE. Finally, pixel shadow fraction from the 5-Scale geometric-optical model was closely related to both MODIS PRI and tower-derived LUE suggesting that differences in stand leaf area and in diffuse illumination among flux sites play a role in the relationship we observed between LUE and MODIS PRI.