MODIS snow albedo bias at high solar zenith angles relative to theory and to in situ observations in Greenland

In situ measurements of snow albedo at five stations along a north-south transect in the dry-snow facies of the interior of Greenland follow the theoretically expected dependence of snow albedo with solar zenith angle (SZA). Greenland Climate Network (GC-Net) measurements from 1997 through 2007 exhibit the trend of modest surface brightening with increasing SZA on both diurnal and seasonal timescales. SZA explains up to 50% of seasonal albedo variability. The two other environmental factors considered, temperature and cloudiness, play much less significant roles in seasonal albedo variability at the five stations studied. Compared to the 10-year record of these GC-Net measurements, the five-year record of MODIS satellite-retrieved snow albedo shows a systematic negative bias for SZA larger than about 55°. Larger bias of MODIS snow albedo exists at more northerly stations. MODIS albedos successfully capture the snow albedo dependence on SZA and have a relatively good agreement with GC-Net measurements for SZA < 55°. The discrepancy of MODIS albedo with in situ albedo and with theory is determined mainly by two related factors, SZA and retrieval quality. While the spatiotemporal structure, especially zonal features, of the MODIS-retrieved albedo may be correct for large SZA, the accuracy deteriorates for SZA > 55° and often becomes physically unrealistic for SZA > 65°. This unphysical behavior biases parameterizations of surface albedo and restricts the range of usefulness of the MODIS albedo products. Seasonal-to-interannual trends in surface brightness in Greenland, and in polar (i.e., large SZA) regions in general, and model simulations of these trends, should be evaluated in light of these limitations.     

View angle effects on the discrimination of soybean varieties and on the relationships between vegetation indices and yield using off-nadir Hyperion data

Imaging spectrometry has the potential to provide improved discrimination of crop types and better estimates of crop yield. Here we investigate the potential of Hyperion to discriminate three Brazilian soybean varieties and to evaluate the relationship between grain yield and 17 narrow-band vegetation indices. Hyperion analysis focused on two datasets acquired from opposite off-nadir viewing directions but similar solar geometry: one acquired on 08 February 2005 (forward scattering) and the other on 14 January 2006 (back scattering). In 2005, the soybean canopies were observed by Hyperion at later reproductive stages than in 2006. Additional Hyperion datasets were not available due to cloud cover. To further examine the impact of viewing geometry within the same season, Hyperion data were complemented by 250 m Moderate Resolution Imaging Spectroradiometer (MODIS) images (bands 1 and 2) acquired in consecutive days (05-06 February 2005) with opposite viewing geometries (− 42° and + 44°, respectively). MODIS data analysis was used to keep reproductive stage as a constant factor while isolating the impact of viewing geometry. For discrimination purposes, multiple discriminant analysis (MDA) was applied over each dataset using surface reflectance values as input variables and a stepwise procedure for band selection. All possible Hyperion band ratios and the 17 narrow-band vegetation indices with soybean grain yield were evaluated across years through Pearson's correlation coefficients and linear regression. MODIS-derived Normalized Difference Vegetation Index (NDVI) and Simple Ratio (SR) were evaluated within the same growing season. Results showed that: (1) the three soybean varieties were discriminated with highest accuracy in the back scattering direction, as deduced from MDA classification results from Hyperion and MODIS data; (2) the highest correlation between Hyperion vegetation indices and soybean yield was observed for the Normalized Difference Water Index (NDWI) (r = + 0.74) in the back scattering direction and this result was consistent with band ratio analysis; (3) higher Hyperion correlation results were observed in the back scattering direction when compared to the forward scattering image. For the same reproductive stage, stronger shadowing effects were observed over the MODIS red band in the forward scattering direction producing lower and lesser variable reflectance for the sensor. As a result, the relationship between MODIS-derived NDVI and soybean yield improved from the forward (r of + 0.21) to the back scattering view (r of + 0.60). The same trend was observed for SR that increased from + 0.22 to + 0.58.     

Leaf area index estimation with MODIS reflectance time series and model inversion during full rotations of Eucalyptus plantations

The leaf area index (LAI) of fast-growing Eucalyptus plantations is highly dynamic both seasonally and inter-annually, and is spatially variable depending on pedo-climatic conditions. LAI is very important in determining the carbon and water balance of a stand, but is difficult to measure during a complete stand rotation and at large scales. Remote-sensing methods allowing the retrieval of LAI time series with accuracy and precision are therefore necessary. Here, we tested two methods for LAI estimation from MODIS 250m resolution red and near-infrared (NIR) reflectance time series. The first method involved the inversion of a coupled model of leaf reflectance and transmittance (PROSPECT4), soil reflectance (SOILSPECT) and canopy radiative transfer (4SAIL2). Model parameters other than the LAI were either fixed to measured constant values, or allowed to vary seasonally and/or with stand age according to trends observed in field measurements. The LAI was assumed to vary throughout the rotation following a series of alternately increasing and decreasing sigmoid curves. The parameters of each sigmoid curve that allowed the best fit of simulated canopy reflectance to MODIS red and NIR reflectance data were obtained by minimization techniques. The second method was based on a linear relationship between the LAI and values of the GEneralized Soil Adjusted Vegetation Index (GESAVI), which was calibrated using destructive LAI measurements made at two seasons, on Eucalyptus stands of different ages and productivity levels. The ability of each approach to reproduce field-measured LAI values was assessed, and uncertainty on results and parameter sensitivities were examined. Both methods offered a good fit between measured and estimated LAI (R² = 0.80 and R² = 0.62 for model inversion and GESAVI-based methods, respectively), but the GESAVI-based method overestimated the LAI at young ages.     

Commentary on Wang and Zender—MODIS snow albedo bias at high solar zenith angles relative to theory and to in situ observations in Greenland

The recent paper by Wang and Zender [Wang, X., & Zender, C. S. (2010). MODIS snow albedo bias at high solar zenith angles relative to theory and to in situ observations in Greenland. Remote Sensing of Environment.] draws erroneous conclusions about solar zenith angle biases at high latitudes by not making appropriate use of the extensive quality flags available with the MODIS BRDF/Albedo. Coarse resolution MODIS white-sky albedo data are compared with actual blue-sky field albedometer measurements from the Greenland GC-Net. By utilizing large area averages of the MODIS data product that combine both high quality and poor quality data indiscriminately, the authors erroneously conclude that the accuracy deteriorates for solar zenith angle (SZA) > 55° and often becomes physically unrealistic for SZA > 65°. Once the quality flags are considered, however, the comparisons demonstrate that the MODIS product performs quite well out to the recommended limit for product use of 70° SZA. This verifies the conclusions of an earlier more rigorous evaluation performed by Stroeve et al. [Stroeve, J., Box, J. E., Gao, F., Liang, S., Nolin, A., & Schaaf, C. B. (2005). Accuracy assessment of the MODIS 16-day albedo product for snow: comparisons with Greenland in situ measurements. Remote Sensing of Environment.]. With over a decade of observations and products now available from the MODIS instrument, these data are increasingly being used to evaluate and tune climate and biogeochemical models. However, such use should take into account the documented quality and limitations of the satellite-derived product.     

Reflectance seasonality and its relation to the canopy leaf area index in an eastern Siberian larch forest: Multi-satellite data and radiative transfer analyses

Reliable monitoring of seasonality in the forest canopy leaf area index (LAI) in Siberian forests is required to advance the understanding of climate-forest interactions under global environmental change and to develop a forest phenology model within ecosystem modeling. Here, we compare multi-satellite (AVHRR, MODIS, and SPOT/VEGETATION) reflectance, normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and LAI with aircraft-based spectral reflectance data and field-measured forest data acquired from April to June in 2000 in a larch forest near Yakutsk, Russia. Field data in a 30 × 30-m study site and aircraft data observed around the field site were used. Larch is a dominant forest type in eastern Siberia, but comparison studies that consider multi-satellite data, aircraft-based reflectance, and field-based measurement data are rarely conducted. Three-dimensional canopy radiative transfer calculations, which are based on Antyufeev and Marshak's [Antyufeev, V.S., & Marshak, A.L. (1990). Monte Carlo method and transport equation in plant canopies, Remote Sensing of Environment, 31, 183-191] Monte Carlo photon transport method combined with North's [North, P.R. (1996). Three-dimensional forest light interaction model using a Monte Carlo method, IEEE Transactions on Geoscience and Remote Sensing, 34(4), 946-956] geometric-optical hybrid forest canopy scene, helped elucidate the relationship between canopy reflectance and forest structural parameters, including several forest floor conditions. Aircraft-based spectral measurements and the spectral response functions of all satellite sensors confirmed that biases in reflectance seasonality caused by differences in spectral response functions among sensors were small. However, some reflectance biases occur among the near infrared (NIR) reflectance data from satellite products; these biases were potentially caused by absolute calibration errors or cloud/cloud shadow contamination. In addition, reflectance seasonality in AVHRR-based NIR data was very small compared to other datasets, which was partially due to the spring-to-summer increase in the amount of atmospheric water vapor. Radiative transfer simulations suggest that bi-directional reflectance effects were small for the study site and observation period; however, changes in tree density and forest floor conditions affect the absolute value of NIR reflectance, even if the canopy leaf area condition does not change. Reliable monitoring of canopy LAI is achieved by minimizing these effects through the use of NIR reflectance difference, i.e., the difference in reflectance on the observation day from the reflectance on a snow-free/pre-foliation day. This may yield useful and robust parameters for multi-satellite monitoring of the larch canopy LAI with less error from intersensor biases and forest structure/floor differences. Further validation with field data and combined use of other index (e.g. normalized difference water index, NDWI) data will enable an extension of these findings to all Siberian deciduous forests.     

Detecting leaf phenology of seasonally moist tropical forests in South America with multi-temporal MODIS images

Leaf phenology of tropical evergreen forests affects carbon and water fluxes. In an earlier study of a seasonally moist evergreen tropical forest site in the Amazon basin, time series data of Enhanced Vegetation Index (EVI) from the VEGETATION and Moderate Resolution Imaging Spectroradiometer (MODIS) sensors showed an unexpected seasonal pattern, with higher EVI in the late dry season than in the wet season. In this study we conducted a regional-scale analysis of tropical evergreen forests in South America, using time series data of EVI from MODIS in 2002. The results show a large dynamic range and spatial variations of annual maximum EVI for evergreen forest canopies in the region. In tropical evergreen forests, maximum EVI in 2002 typically occurs during the late dry season to early wet season. This suggests that leaf phenology in tropical evergreen forests is not determined by the seasonality of precipitation. Instead, leaf phenological process may be driven by availability of solar radiation and/or avoidance of herbivory.     

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

Evaluation of MODIS land surface temperature data to estimate air temperature in different ecosystems over Africa

The estimation of near surface air temperature (Ta) is useful for a wide range of applications such as agriculture, climate related diseases and climate change studies. Air temperature is commonly obtained from synoptic measurements in weather stations. In Africa, the spatial distribution of weather stations is often limited and the dissemination of temperature data is variable, therefore limiting their use for real-time applications. Compensation for this paucity of information may be obtained by using satellite-based methods. However, the derivation of near surface air temperature (Ta), from the land surface temperature (Ts) derived from satellite is far from straight forward. Some studies have tried to derive maximum Ta from satellites through regression analysis but the accuracy obtained is quite variable according to the study. The main objective of this study was to explore the possibility of retrieving high-resolution Ta data from the Moderate Resolution Imaging Spectroradiometer (MODIS) Ts products over different ecosystems in Africa. First, comparisons between night MODIS Ts data with minimum Ta showed that MODIS nighttime products provide a good estimation of minimum Ta over different ecosystems (with (ΔTs − Ta) centered at 0 °C, a mean absolute error (MAE) = 1.73 °C and a standard deviation = 2.4 °C). Secondly, comparisons between day MODIS Ts data with maximum Ta showed that (ΔTs − Ta) strongly varies according to the seasonality, the ecosystems, the solar radiation, and cloud-cover. Two factors proposed in the literature to retrieve maximum Ta from Ts, i.e. the Normalized Difference Vegetation Index (NDVI) and the Solar Zenith Angle (SZA), were analyzed. No strong relationship between (ΔTs − Ta) and (i) NDVI and (ii) SZA was observed, therefore requiring further research on robust methods to retrieve maximum Ta.     

Satellite remote sensing of tropical forest canopies and their seasonal dynamics

Seasonal patterns of tropical evergreen forest green-up in Amazonia, corresponding to drought and the dry season, have recently been detected by the Enhanced Vegetation Index (EVI) and Leaf Area Index (LAI) products of the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensor. These observations provide additional evidence for solar radiation as the primary limiting factor regulating wet-tropical ecosystem processes. However, in situ structural mechanisms for forest canopy green-up are unclear and frequently inconsistent with observations. Here, we investigate the signal of seasonal green-up at several intensively measured sites, applying a rigorous series of filters to minimize error from atmospheric contamination that is common in tropical regions. We find that, while satellite-observed forest seasonality is sensitive to data-quality filtering, statistical noise reduction and spatial averaging, the signal is robust at sites where field observations are available, and in particular for the EVI. For the sites where field data are unavailable, it appears that additional filters to those commonly used to remove cloud effects and aerosols also reduce the seasonal magnitude of the LAI. These findings imply that seasonal tropical evergreen forest green-up remains sensitive to the methodology used in removing seasonal contamination from atmospheric conditions and that further field measurements and comparisons to remote sensing are required to reduce this uncertainty.     

Satellite passive microwave remote sensing for monitoring global land surface phenology

Vegetation phenology characterizes seasonal life-cycle events that influence the carbon cycle and land-atmosphere water and energy exchange. We analyzed global phenology cycles over a six year record (2003-2008) using satellite passive microwave remote sensing based Vegetation Optical Depth (VOD) retrievals derived from daily time series brightness temperature (T_b) measurements from the Advanced Microwave Scanning Radiometer on EOS (AMSR-E) and other ancillary data inputs. The VOD parameter derives vegetation canopy attenuation at a given microwave frequency (18.7 GHz) and varies with canopy height, density, structure and water content. An error sensitivity analysis indicates that the retrieval algorithm can resolve the VOD seasonal cycle over a majority of global vegetated land areas. The VOD results corresponded favorably (p < 0.01) with vegetation indices (VIs) and leaf area index (LAI) information from satellite optical-infrared (MODIS) remote sensing, and phenology cycles determined from a simple bioclimatic growing season index (GSI) for over 82% of the global domain. Lower biomass land cover classes (e.g. savannas) show the highest correlations (R = 0.66), with reduced correspondence at higher biomass levels (0.03 < R < 0.51) and higher correlations for homogeneous land cover areas (0.41 < R < 0.83). The VOD results display a unique end-of-season signal relative to VI and LAI series, and may reflect microwave sensitivity to the timing of vegetation biomass depletion (e.g. leaf abscission) and associated changes in canopy water content (e.g. dormancy preparation). The VOD parameter is independent of and synergistic with optical-infrared remote sensing based vegetation metrics, and contributes to a more comprehensive view of land surface phenology.     

Analysis of leaf area index products from combination of MODIS Terra and Aqua data

A prototype product suite, containing the Terra 8-day, Aqua 8-day, Terra-Aqua combined 8- and 4-day products, was generated as part of testing for the next version (Collection 5) of the MODerate resolution Imaging Spectroradiometer (MODIS) leaf area index (LAI) products. These products were analyzed for consistency between Terra and Aqua retrievals over the following data subsets in North America: single 8-day composite over the whole continent and annual time series over three selected MODIS tiles (1200 × 1200 km). The potential for combining retrievals from the two sensors to derive improved products by reducing the impact of environmental conditions and temporal compositing period was also explored. The results suggest no significant discrepancies between large area (from continent to MODIS tile) averages of the Terra and Aqua 8-day LAI and surface reflectances products. The differences over smaller regions, however, can be large due to the random nature of residual atmospheric effects. High quality retrievals from the radiative transfer based algorithm can be expected in 90-95% of the pixels with mostly herbaceous cover and about 50-75% of the pixels with woody vegetation during the growing season. The quality of retrievals during the growing season is mostly restricted by aerosol contamination of the MODIS data. The Terra-Aqua combined 8-day product helps to minimize this effect and increases the number of high quality retrievals by 10-20% over woody vegetation. The combined 8-day product does not improve the number of high quality retrievals during the winter period because the extent of snow contamination of Terra and Aqua observations is similar. Likewise, cloud contamination in the single-sensor and combined products is also similar. The LAI magnitudes, seasonal profiles and retrieval quality in the combined 4-day product are comparable to those in the single-sensor 8-day products. Thus, the combined 4-day product doubles the temporal resolution of the seasonal cycle, which facilitates phenology monitoring in application studies during vegetation transition periods. Both Terra and Aqua LAI products show anomalous seasonality in boreal needle leaf forests, due to limitations of the radiative transfer algorithm to model seasonal variations of MODIS surface reflectance data with respect to solar zenith angle. Finally, this study suggests that further improvement of the MODIS LAI products is mainly restricted by the accuracy of the MODIS observations.     

An assessment of the MODIS collection 5 leaf area index product for a region of mixed coniferous forest

Canopy leaf area index (LAI), defined as the single-sided leaf area per unit ground area, is a quantitative measure of canopy foliar area. LAI is a controlling biophysical property of vegetation function, and quantifying LAI is thus vital for understanding energy, carbon and water fluxes between the land surface and the atmosphere. LAI is routinely available from Earth Observation (EO) instruments such as MODIS. However EO-derived estimates of LAI require validation before they are utilised by the ecosystem modelling community. Previous validation work on the MODIS collection 4 (c4) product suggested considerable error especially in forested biomes, and as a result significant modification of the MODIS LAI algorithm has been made for the most recent collection 5 (c5). As a result of these changes the current MODIS LAI product has not been widely validated. We present a validation of the MODIS c5 LAI product over a 121 km² area of mixed coniferous forest in Oregon, USA, based on detailed ground measurements which we have upscaled using high resolution EO data. Our analysis suggests that c5 shows a much more realistic temporal LAI dynamic over c4 values for the site we examined. We find improved spatial consistency between the MODIS c5 LAI product and upscaled in situ measurements. However results also suggest that the c5 LAI product underestimates the upper range of upscaled in situ LAI measurements.     

An intercomparison of three remote sensing-based energy balance models using Large Aperture Scintillometer measurements over a wheat-corn production region

This paper compares three remote sensing-based models for estimating evapotranspiration (ET), namely the Surface Energy Balance System (SEBS), the Two-Source Energy Balance (TSEB) model, and the surface temperature-vegetation index Triangle (TVT). The models used as input MODIS/TERRA products and ground measurements collected during the wheat and corn growth period in a subhumid climate at a measurement station in Yucheng, China. MODIS land surface temperature (LST) and leaf area index (LAI) products, corrected using ground-truth observations, were used in the three models. The TSEB model output of sensible (H) and latent (LE) heat fluxes were in good agreement with Large Aperture Scintillometer (LAS)-measured H and LE derived by residual (RMSD < 45 W/m²). Reasonable agreement was also obtained with the SEBS model output yielding RMSD for H of ~ 40 W/m² and LE ~ 55 W/m². However, the TVT model output resulted in poor agreement with the LAS-estimated H and LE with RMSD-values > 110 W/m². Using the uncorrected MODIS LST and LAI products resulted in a deterioration of the agreement in H and LE with LAS-estimated values for both the TSEB and SEBS models, whereas TVT performance improved marginally. These results indicate that the TSEB model yielded the closest agreement with the LAS-estimated fluxes using either the corrected or uncorrected MODIS inputs (LST and LAI). The SEBS model also computed reasonable H and LE values but was significantly more sensitive to errors in MODIS LST and LAI inputs than the TSEB model. In the TVT model, output of H and LE was unacceptable in either scenario of MODIS input which was attributable to errors in selection of the dry edge. With the TVT method, accurate determination of the dry edge end member is critical in regional ET estimation, but for humid and subhumid regions this end member may often be quite difficult to identify or encompass within a satellite scene.     

Evaluation of the JRC-TIP 0.01° products over a mid-latitude deciduous forest site

The Joint Research Centre Two-stream Inversion Package (JRC-TIP) makes use of white sky albedo products—derived from MODIS and MISR observations in the visible and near-infrared domain—to deliver consistent sets of information about the terrestrial environments that gave rise to these data. The baseline version of the JRC-TIP operates at a spatial resolution of 0.01° and yields estimates of the Probability Distribution Functions (PDFs) of the effective canopy Leaf Area Index (LAI), the canopy background albedo, the vegetation scattering properties, as well as, the absorbed, reflected and transmitted fluxes of the vegetation canopy. In this contribution the evaluation efforts of the JRC-TIP products are extended to the deciduous forest site of Hainich (Germany) where multiannual datasets of in-situ estimates of canopy transmission—derived from LAI-2000 observations—are available. As a Fluxnet site, Hainich offers access to camera acquisitions from fixed locations in and above the canopy that are being used in phenological studies. These images qualitatively confirm the seasonal patterns of the effective LAI, canopy transmission and canopy absorption products (in the visible range of the solar spectrum) derived with the JRC-TIP. Making use of the LAI-2000 observations it is found that 3/4 of the JRC-TIP products lie within a ± 0.15 interval around the in-situ estimates of canopy transmission and absorption. The largest discrepancies occur at the end of the senescence phase when the scattering properties of the vegetation (evidenced by the pictures) and the images qualitatively confirm the seasonal patterns of the effective LAI, canopy transmission and canopy absorption products (in the visible range of the solar spectrum) derived with the JRC-TIP. Making use of the LAI-2000 observations it is found that 3/4 of the JRC-TIP products lie within a ± 0.15 interval around the in-situ estimates of canopy transmission and absorption. The largest discrepancies occur at the end of the senescence phase when the scattering properties of the vegetation (evidenced by the pictures) and the effective LAI (also derived from LAI-2000 measurements) are experiencing large simultaneous changes. It was also found that the seasonal pattern of vegetation scattering properties derived from MISR observations in the near-infrared varies together with the Excess Green index computed from the various channels of the camera data acquired at the top of the canopy.     

Evaluation of Collections 4 and 5 of the MODIS Gross Primary Productivity product and algorithm improvement at a tropical savanna site in northern Australia

In this study, we assessed the accuracy of the MODIS (Moderate Resolution Imaging Spectroradiometer) GPP (gross primary productivity) Collections 4.5, 4.8 and 5 along with Leaf Area Index (LAI), fraction of absorbed Photosynthetically Active Radiation (fPAR), light use efficiency (LUE) and meteorological variables that are used to estimate GPP for a northern Australian savanna site. Results of this study indicated that the MODIS products captured the seasonal variation in GPP, LAI and fPAR well. Using the index of agreement (IOA), it was found that Collections 4.5 and 4.8 (IOA 0.89 respectively) agreed reasonably well with flux tower measurements between 2001 and 2006. It was also found that MODIS Collection 4.5 predicted the dry season GPP well (Relative Predictive Error (RPE) 4.17%, IOA 0.72 and Root Mean Square Error (RMSE) of 1.05 g C m^− 2 day^− 1), whilst Collection 4.8 performed better in capturing wet season dynamics (RPE 1.11%, IOA 0.80 and RMSE of 0.91 g C m^− 2 day^− 1). Although the wet season magnitude of GPP was predicted well by Collection 4.8, an examination of the inputs to the GPP algorithm revealed that MODIS fPAR was too high, but this was compensated by PAR and LUE that was too low. Although LAI and fPAR estimated by Collection 5 were more accurate, GPP for this Collection resulted in a much lower value (RPE 25%) due to errors in other factors. Recalculation of MODIS GPP using site specific input parameters indicated that MODIS fPAR was the main reason for the differences between MODIS and tower derived GPP followed by LUE and meteorological inputs. GPP calculated using all site specific values agreed very well with tower data on an annual basis (IOA 0.94, RPE 6.06% and RMSE 0.83 g C m^− 2 day^− 1) but the early initiation of the growing season calculated by the MODIS algorithm was improved when the vapor pressure deficit (VPD) function was replaced with a soil water deficit function. The results of this study however, reinforce previous findings in water limited regions, like Australia, and incorporation of soil moisture in a LUE model is needed to accurately estimate the productivity.     

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.     

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.     

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.