Predicting riparian evapotranspiration from MODIS vegetation indices and meteorological data

A vegetation index (VI) model for predicting evapotranspiration (ET) from data from the Moderate Resolution Imaging Spectrometer (MODIS) on the EOS-1 Terra satellite and ground meteorological data was developed for riparian vegetation along the Middle Rio Grande River in New Mexico. Ground ET measurements obtained from eddy covariance towers at four riparian sites were correlated with MODIS VIs, MODIS land surface temperatures (LSTs), and ground micrometeorological data over four years. Sites included two saltcedar (Tamarix ramosissima) and two Rio Grande cottonwood (Populus deltoides ssp. Wislizennii) dominated stands. The Enhanced Vegetation Index (EVI) was more closely correlated (r=0.76) with ET than the Normalized Difference Vegetation Index (NDVI; r=0.68) for ET data combined over sites and species. Air temperature (T_a) measured over the canopy from towers was the meteorological variable that was most closely correlated with ET (r=0.82). MODIS LST data at 1- and 5-km resolutions were too coarse to accurately measure the radiant surface temperature within the narrow riparian corridor; hence, energy balance methods for estimating ET using MODIS LSTs were not successful. On the other hand, a multivariate regression equation for predicting ET from EVI and T_a had an r²=0.82 across sites, species, and years. The equation was similar to VI-ET models developed for crop species. The finding that ET predictions did not require species-specific equations is significant, inasmuch as these are mixed vegetation zones that cannot be easily mapped at the species level.     

Modeling gross primary production of temperate deciduous broadleaf forest using satellite images and climate data

Net ecosystem exchange (NEE) of CO₂ between the atmosphere and forest ecosystems is determined by gross primary production (GPP) of vegetation and ecosystem respiration. CO₂ flux measurements at individual CO₂ eddy flux sites provide valuable information on the seasonal dynamics of GPP. In this paper, we developed and validated the satellite-based Vegetation Photosynthesis Model (VPM), using site-specific CO₂ flux and climate data from a temperate deciduous broadleaf forest at Harvard Forest, Massachusetts, USA. The VPM model is built upon the conceptual partitioning of photosynthetically active vegetation and non-photosynthetic vegetation (NPV) within the leaf and canopy. It estimates GPP, using satellite-derived Enhanced Vegetation Index (EVI), Land Surface Water Index (LSWI), air temperature and photosynthetically active radiation (PAR). Multi-year (1998-2001) data analyses have shown that EVI had a stronger linear relationship with GPP than did the Normalized Difference Vegetation Index (NDVI). Two simulations of the VPM model were conducted, using vegetation indices from the VEGETATION (VGT) sensor onboard the SPOT-4 satellite and the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard the Terra satellite. The predicted GPP values agreed reasonably well with observed GPP of the deciduous broadleaf forest at Harvard Forest, Massachusetts. This study highlighted the biophysical performance of improved vegetation indices in relation to GPP and demonstrated the potential of the VPM model for scaling-up of GPP of deciduous broadleaf forests.     

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.     

Boosted decision tree classifications of land cover over Turkey integrating MODIS,climate and topographic data

This study investigates the impact of using different combinations of Moderate Resolution Imaging Spectroradiometer (MODIS) and ancillary datasets on overall and per-class classification accuracies for nine land cover types modified from the classification system of the International Geosphere Biosphere Programme (IGBP). Twelve land cover maps were generated for Turkey using boosted decision trees (BDTs) based on the stepwise addition of 14 explanatory variables derived from a time series of 16-day MODIS composites between 2000 and 2006 (Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI) and four spectral bands) and ancillary climate and topographic data (minimum and maximum air temperature, precipitation, potential evapotranspiration, aspect, elevation, distance to sea and slope) at 500-m resolution. Evaluation of the 12 BDTs indicated that the BDT built as a function of all the MODIS and climate variables, aspect and elevation produced the highest degree of overall classification accuracy (79.8%) and kappa statistic (0.76) followed by the BDTs that additionally included distance to sea (DtS), and both DtS and slope. Based on an independent validation dataset derived from a pre-existing national forest map and Landsat images of Turkey, the highest overall accuracy (64.7%) and kappa coefficient (0.58) among the 12 land cover maps was achieved by using MODIS-derived NDVI time series only, followed by NDVI and EVI time series combined; NDVI, EVI and four MODIS spectral bands; and the combination of all MODIS and climate data, aspect, elevation and distance to sea, respectively. The largest improvements in producer's accuracies were observed for grasslands (+50%), barrenlands (+46%) and mixed forests (+39%) and in user's accuracies for grasslands (+53%), shrublands (+30%) and mixed forests (+28%), in relation to the lowest producer's accuracy. The results of this study indicate that BDTs can increase the accuracy of land cover classifications at the national scale.     

Satellite-based modeling of gross primary production in a seasonally moist tropical evergreen forest

A CO₂ eddy flux tower study has recently reported that an old-growth stand of seasonally moist tropical evergreen forest in Santarém, Brazil, maintained high gross primary production (GPP) during the dry seasons [Saleska, S. R., Miller, S. D., Matross, D. M., Goulden, M. L., Wofsy, S. C., da Rocha, H. R., de Camargo, P. B., Crill, P., Daube, B. C., de Freitas, H. C., Hutyra, L., Keller, M., Kirchhoff, V., Menton, M., Munger, J. W., Pyle, E. H., Rice, A. H., & Silva, H. (2003). Carbon in amazon forests: Unexpected seasonal fluxes and disturbance-induced losses. Science, 302, 1554-1557]. It was proposed that seasonally moist tropical evergreen forests have evolved two adaptive mechanisms in an environment with strong seasonal variations of light and water: deep roots system for access to water in deep soils and leaf phenology for access to light. Identifying tropical forests with these adaptive mechanisms could substantially improve our capacity of modeling the seasonal dynamics of carbon and water fluxes in the tropical zone. In this paper, we have analyzed multi-year satellite images from the VEGETATION (VGT) sensor onboard the SPOT-4 satellite (4/1998-12/2002) and the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra satellite (2000-2003). We reported temporal analyses of vegetation indices and simulations of the satellite-based vegetation photosynthesis model (VPM). The Enhanced Vegetation Index (EVI) identified subtle changes in the seasonal dynamics of leaf phenology (leaf emergence, leaf aging and leaf fall) in the forest, as suggested by the leaf litterfall data. The land surface water index (LSWI) indicated that the forest experienced no water stress in the dry seasons of 1998-2002. The VPM model, which uses EVI, LSWI and site-specific climate data (air temperature and photosynthetically active radiation, PAR) for 2001-2002, predicted high GPP in the late dry seasons, consistent with observed high evapotranspiration and estimated GPP from the CO₂ eddy flux tower.     

Development of a two-band enhanced vegetation index without a blue band

The enhanced vegetation index (EVI) was developed as a standard satellite vegetation product for the Terra and Aqua Moderate Resolution Imaging Spectroradiometers (MODIS). EVI provides improved sensitivity in high biomass regions while minimizing soil and atmosphere influences, however, is limited to sensor systems designed with a blue band, in addition to the red and near-infrared bands, making it difficult to generate long-term EVI time series as the normalized difference vegetation index (NDVI) counterpart. The purpose of this study is to develop and evaluate a 2-band EVI (EVI2), without a blue band, which has the best similarity with the 3-band EVI, particularly when atmospheric effects are insignificant and data quality is good. A linearity-adjustment factor β is proposed and coupled with the soil-adjustment factor L used in the soil-adjusted vegetation index (SAVI) to develop EVI2. A global land cover dataset of Terra MODIS data extracted over land community validation and FLUXNET test sites is used to develop the optimal parameter (L, β and G) values in EVI2 equation and achieve the best similarity between EVI and EVI2. The similarity between the two indices is evaluated and demonstrated with temporal profiles of vegetation dynamics at local and global scales. Our results demonstrate that the differences between EVI and EVI2 are insignificant (within ± 0.02) over a very large sample of snow/ice-free land cover types, phenologies, and scales when atmospheric influences are insignificant, enabling EVI2 as an acceptable and accurate substitute of EVI. EVI2 can be used for sensors without a blue band, such as the Advanced Very High Resolution Radiometer (AVHRR), and may reveal different vegetation dynamics in comparison with the current AVHRR NDVI dataset. However, cross-sensor continuity relationships for EVI2 remain to be studied.     

Assessing vegetation condition in the presence of biomass burning smoke by applying the Aerosol‐free Vegetation Index (AFRI) on MODIS images

Vegetation indices (VIs) such as the Normalized Difference Vegetation Index (NDVI) are widely used for assessing vegetation cover and condition. One of the NDVI's significant disadvantages is its sensitivity to aerosols in the atmosphere, hence several atmospherically resistant VIs were formulated using the difference in the radiance between the blue and the red spectral bands. The state‐of‐the‐art atmospherically resistant VI, which is a standard Moderate Resolution Imaging Spectroradiometer (MODIS) product, together with the NDVI, is the Enhanced Vegetation Index (EVI). A different approach introduced the Aerosol‐free Vegetation Index (AFRI) that is based on the correlation between the shortwave infrared (SWIR) and the visible red bands. The AFRI main advantage is in penetrating an opaque atmosphere influenced by biomass burning smoke, without the need for explicit correction for the aerosol effect. The objective of this research was to compare the performance of these three VIs under smoke conditions. The AFRI was applied to the 2.1 µm SWIR channel of the MODIS sensor onboard the Earth Observing System (EOS) Terra and Aqua satellites in order to assess its functionality on these imaging platforms. The AFRI performance was compared with those of NDVI and EVI. All VIs were calculated on images with and without present smoke, using the surface‐reflectance MODIS product, for three case studies of fires in Arizona, California, and Zambia. The MODIS Fire Product was embedded on the images in order to identify the exact location of the active fires. Although good correlations were observed between all VIs in the absence of smoke (in the Arizona case R ² = 0.86, 0.77, 0.88 for the NDVI–EVI, AFRI–EVI, and AFRI–NDVI, respectively) under smoke conditions a high correlation was maintained between the NDVI and the EVI, while low correlations were found for the AFRI–EVI and AFRI–NDVI (0.21 and 0.16, for the Arizona case, respectively). A time series of MODIS images recorded over Zambia during the summer of 2000 was tested and showed high NDVI fluctuations during the study period due to oscillations in aerosol optical thickness values despite application of aerosol corrections on the images. In contrast, the AFRI showed smoother variations and managed to better assess the vegetation condition. It is concluded that, beneath the biomass burning smoke, the AFRI is more effective than the EVI in observing the vegetation conditions.     

Using hyperspectral remote sensing to detect and quantify southeastern pine senescence effects in red-cockaded woodpecker (Picoides borealis) habitat

Conservation of threatened and endangered species requires maintenance of critical habitat. The red-cockaded woodpecker Picoides borealis (RCW) is a threatened bird species, endemic to the mixed conifer forests of the southeastern United States. RCW nests and forages preferentially in mature longleaf pine Pinus palustris, but also uses mature loblolly pine Pinus taeda and shortleaf pine Pinus echinata forests. In the last century, the extent of mature pine forests has been greatly reduced by logging. The RCW, in contrast to other woodpeckers, excavates nest cavities in living trees and senescence symptoms (year round leaf chlorosis and leaf mortality) have been observed in mature pine stands across the southeast. Widespread mortality of the mature pine forests would threaten the long-term survival of the RCW. We used airborne hyperspectral data across a portion of Ft. Benning Military Installation, Georgia, U.S.A., to determine if senescent trees can be identified and mapped and assess the likely persistence of mature pines in the RCW habitat. Univariate analysis of variance showed good separation between asymptomatic, senesced and dead physiological conditions with asymptomatic trees having significantly higher reflectance for all bands in the wavelength range between 0.719 and 1.1676 µm, senescent trees having significantly lower reflectance for bands in the range between 1.1927 and 1.3122 µm, and dead trees having significantly higher reflectance for bands in the range between 1.8151 and 1.9471 µm. Classification and Regression Tree (CART) models achieved correct classification rates and kappa statistics above 70%. CART models selected information from wavelength regions similar to those identified from the ANOVA, which likely explains their performance. Our aggregated CART model of tree senescence estimated that 141.4 ha (3%) of the study area is affected. RCW nests occurred in areas with significantly higher tree cover, and trees within foraging and home ranges did not have significantly more senescence than areas without RCW. These results indicate that although tree senescence is widespread, mortality is yet to significantly affect RCW habitat. Results and analysis of critical habitat similar to those exemplified in this study can extend our knowledge about the stressors of these important and imperiled components of biodiversity.     

Land‐cover mixing and spectral vegetation indices

Vegetation indices have been widely used as indicators of seasonal and inter‐annual variations in vegetation caused by either human activities or climate, with the overall goal of observing and documenting changes in the Earth system. While existing satellite remote sensing systems, such as NASA's Multi‐angle Imaging SpectroRadiometer (MISR) and Moderate Resolution Imaging Spectroradiometer (MODIS), are providing improved vegetation index data products through correcting for the distortions in surface reflectance caused by atmospheric particles as well as ground covers below vegetation canopy, the impact of land‐cover mixing on vegetation indices has not been fully addressed. In this study, based on real image spectral samples for two‐component mixtures of forest and common nonforest land‐cover types directly extracted from a 1.1?km MISR image by referencing a 30?m land‐cover classification, the effect of land‐cover mixing on the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI) has been quantitatively evaluated. When the areal fraction of forest was lower than 80%, both NDVI and EVI varied greatly with mixed land‐cover types, although EVI varied less than NDVI. Such a phenomenon can cause errors in applications based on use of these vegetation indices. This study suggests that methods that reduce land‐cover mixing effects should be introduced when developing new spectral vegetation indices.     

Reflectance assessment of seasonal and annual changes in biomass and CO2 uptake of a Mediterranean shrubland submitted to experimental warming and drought

We aimed to evaluate how the remote sensing vegetation indices NDVI and PRI responded to seasonal and annual changes in an early successional stage Mediterranean coastal shrubland canopy that was submitted to experimental warming and drought simulating predicted climate change for the next decades. These conditions were obtained by using a new non-intrusive methodological approach that increases the temperature and prolongs the drought period by using roofs that automatically cover the vegetation after the sunset or when it rains. On average, warming increased air temperature by 0.7 °C and soil temperature by 1.6 °C, and the drought treatment reduced soil moisture by 22%. We measured spectral reflectance at the canopy level and at the individual plant level seasonally during 4 years. Shrubland NDVI tracked the community development and activity. In control and warming treatments, NDVI increased with the years while it did not change in the drought treatment. There was a good relationship between NDVI and both community and individual plant biomass. NDVI also decreased in summer seasons when some species dry or decolour. The NDVI of E. multiflora plant individuals was lower in autumn and winter than in the other seasons, likely because of flowering. Shrubland PRI decreased only in winter, similarly to the PRI of the most dominant species, G. alypum. At this community scale, NDVI was better related than PRI to photosynthetic activity, probably because photosynthetic fluxes followed canopy seasonal greening in this complex canopy, which includes brevideciduous, annual and evergreen species and variable morphologies and canopy coverage. PRI followed the seasonal variations in photosynthetic rates in E. multiflora and detected the decreased photosynthetic rates of drought treatment. However, PRI did not track the photosynthetic rates of G. alypum plants which have lower LAIs than E. multiflora. In this community, which is in its early successional stages, NDVI was able to track biomass, and indirectly, CO₂ uptake changes, likely because LAI values did not saturate NDVI. Thus, NDVI appears as a valid tool for remote tracking of this community development. PRI was less adequate for photosynthetic assessment of this community especially for its lower LAI canopies. PRI usefulness was also species-dependent and could also be affected by flowering. These results will help to improve the interpretation of remote sensing information on the structure and physiological status of these Mediterranean shrublands, and to gain better insight on ecological and environmental controls on their ecosystem carbon dioxide exchange. They also show the possibility of assessing the impacts of climate change on shrubland communities.     

Evaluation of optical satellite remote sensing for rice paddy phenology in monsoon Asia using a continuous in situ dataset

In monsoon Asia, optical satellite remote sensing for rice paddy phenology suffers from atmospheric contaminations mainly due to frequent cloud cover. We evaluated the quality of satellite remote sensing of paddy phenology: (1) through continuous in situ observations of a paddy field in Japan for 1.5 years, we investigated phenological signals in the reflectance spectrum of the paddy field; (2) we tested daily satellite data taken by Terra/Aqua MODIS (MOD09 and L1B products) with regard to the agreement with the in situ data and the influence of cloud contamination. As a result, the in situ spectral characteristics evidently indicated some phenological changes in the rice paddy field, such as irrigation start, padding, heading, harvest and ploughing. The Enhanced Vegetation Index (EVI) was the best vegetation index in terms of agreement with the in situ data. More than 65% of MODIS observations were contaminated with clouds in this region. However, the combined use of Terra and Aqua decreased the rate of cloud contamination of the daily data to 43%. In conclusion, the most robust dataset for monitoring rice paddy phenology in monsoon Asia would be daily EVI derived from a combination of Terra/MODIS and Aqua/MODIS.     

Monitoring of forage conditions with MODIS imagery in the Xilingol steppe,Inner Mongolia

A study was conducted to determine the potential suitability of Terra/MODIS imagery for monitoring short‐term phenological changes in forage conditions in a semi‐arid region. The study sites included four meadow steppes and six typical steppes in the Xilingol steppe in central Inner Mongolia, China. The live biomass, dead standing biomass, total biomass, crude protein (CP) concentration and standing CP were estimated from early April to late October using the Enhanced Vegetation Index (EVI) values from Terra imagery (500?m?pixels). Applying regression models, the EVI accounted for 80% of the variation in live biomass, 42% of the dead biomass, 77% of the total biomass, 11% of the CP concentration and 74% of the standing CP. MODIS/EVI is superior to AVHRR/NDVI when estimating forage quantity. Applying these results, the seasonal changes in live biomass and the standing CP could be described in the selected four sites with different degrees of grazing intensity. Generally, the increase in grazing intensity tended to decrease live biomass and standing CP. It was suggested that the EVI obtained from Terra imagery was an available predictor of the forage condition as measured by live biomass and standing CP. The MODIS/EVI values could provide information on the suitable timing of cutting for hay‐making and nutritive value to range managers.     

Sensitivity of vegetation indices to atmospheric aerosols: continental-scale observations in Northern Asia

Satellite observations play an important role in characterization of the interannual variation of vegetation. Here, we report anomalies of two vegetation indices for Northern Asia (40°N-75°N, and 45°E-179°E), using images from the SPOT-4 VEGETATION (VGT) sensor over the period of April 1, 1998 to November 20, 2001. The Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI), which are correlated to a number of vegetation properties (e.g., net primary production, leaf area index), were compared. The results show that there is a large disagreement between NDVI and EVI anomalies in 1998 and 1999 for Northern Asia. The NDVI anomaly in 1998 was largely affected by atmospheric contamination, predominantly aerosols from extensive forest fires in that year. The EVI anomaly in 1998 was less sensitive to residual atmospheric contamination, as it is designed to be, and thus EVI is a useful alternative vegetation index for the large-scale study of vegetation. The EVI anomaly also suggests that potential vegetation productivity in Northern Asia was highest in 1998 but declined substantially in 2001, consistent with precipitation data from 1998-2001.     

Do flowers affect biomass estimate accuracy from NDVI and EVI?

The Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) are vegetation indices widely used in remote sensing of above-ground biomass. Because both indexes are based on spectral features of plant canopy, NDVI and EVI may suffer reduced accuracy in estimating above-ground biomass when flower signals are mixed in the plant canopy. This paper addresses how flowers influence the estimation of above-ground biomass using NDVI and EVI for an alpine meadow with mixed yellow flowers of Halerpestes tricuspis (Ranunculaceae). Field spectral measurements were used in combination with simulated reflectance spectra with precisely controlled flower coverage by applying a linear spectral mixture model. Using the reflectance spectrum for the in-situ canopy with H. tricuspis flowers, we found no significant correlation between above-ground biomass and EVI (p?=?0.17) or between above-ground biomass and NDVI (p?=?0.78). However, both NDVI and EVI showed very good prediction of above-ground biomass with low root mean square errors (RMSE?=?43 g m^?2 for NDVI and RMSE?=?43 g m^?2 for EVI, both p < 0.01) when all the flowers were removed from the canopies. Simulation analysis based on the in-situ measurements further indicated that high variation in flower coverage among different quadrats could produce more noise in the relationship between above-ground biomass and NDVI, or EVI, which results in an evident decline in the accuracy of above-ground biomass estimation. Therefore, the study suggests that attention should be paid both to the flower fraction and the heterogeneity of flower distribution in the above-ground biomass estimation via NDVI and EVI.     

Analysis of time-series MODIS 250 m vegetation index data for crop classification in the U.S. Central Great Plains

The global environmental change research community requires improved and up-to-date land use/land cover (LULC) datasets at regional to global scales to support a variety of science and policy applications. Considerable strides have been made to improve large-area LULC datasets, but little emphasis has been placed on thematically detailed crop mapping, despite the considerable influence of management activities in the cropland sector on various environmental processes and the economy. Time-series MODIS 250 m Vegetation Index (VI) datasets hold considerable promise for large-area crop mapping in an agriculturally intensive region such as the U.S. Central Great Plains, given their global coverage, intermediate spatial resolution, high temporal resolution (16-day composite period), and cost-free status. However, the specific spectral-temporal information contained in these data has yet to be thoroughly explored and their applicability for large-area crop-related LULC classification is relatively unknown. The objective of this research was to investigate the general applicability of the time-series MODIS 250 m Enhanced Vegetation Index (EVI) and Normalized Difference Vegetation Index (NDVI) datasets for crop-related LULC classification in this region. A combination of graphical and statistical analyses were performed on a 12-month time-series of MODIS EVI and NDVI data from more than 2000 cropped field sites across the U.S. state of Kansas. Both MODIS VI datasets were found to have sufficient spatial, spectral, and temporal resolutions to detect unique multi-temporal signatures for each of the region's major crop types (alfalfa, corn, sorghum, soybeans, and winter wheat) and management practices (double crop, fallow, and irrigation). Each crop's multi-temporal VI signature was consistent with its general phenological characteristics and most crop classes were spectrally separable at some point during the growing season. Regional intra-class VI signature variations were found for some crops across Kansas that reflected the state's climate and planting time differences. The multi-temporal EVI and NDVI data tracked similar seasonal responses for all crops and were highly correlated across the growing season. However, differences between EVI and NDVI responses were most pronounced during the senescence phase of the growing season.     

A simple and effective radiometric correction method to improve landscape change detection across sensors and across time

Satellite data offer unrivaled utility in monitoring and quantifying large scale land cover change over time. Radiometric consistency among collocated multi-temporal imagery is difficult to maintain, however, due to variations in sensor characteristics, atmospheric conditions, solar angle, and sensor view angle that can obscure surface change detection. To detect accurate landscape change using multi-temporal images, we developed a variation of the pseudoinvariant feature (PIF) normalization scheme: the temporally invariant cluster (TIC) method. Image data were acquired on June 9, 1990 (Landsat 4), June 20, 2000 (Landsat 7), and August 26, 2001 (Landsat 7) to analyze boreal forests near the Siberian city of Krasnoyarsk using the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and reduced simple ratio (RSR). The temporally invariant cluster (TIC) centers were identified via a point density map of collocated pixel VIs from the base image and the target image, and a normalization regression line was created to intersect all TIC centers. Target image VI values were then recalculated using the regression function so that these two images could be compared using the resulting common radiometric scale. We found that EVI was very indicative of vegetation structure because of its sensitivity to shadowing effects and could thus be used to separate conifer forests from deciduous forests and grass/crop lands. Conversely, because NDVI reduced the radiometric influence of shadow, it did not allow for distinctions among these vegetation types. After normalization, correlations of NDVI and EVI with forest leaf area index (LAI) field measurements combined for 2000 and 2001 were significantly improved; the r² values in these regressions rose from 0.49 to 0.69 and from 0.46 to 0.61, respectively. An EVI “cancellation effect” where EVI was positively related to understory greenness but negatively related to forest canopy coverage was evident across a post fire chronosequence with normalized data. These findings indicate that the TIC method provides a simple, effective and repeatable method to create radiometrically comparable data sets for remote detection of landscape change. Compared to some previous relative radiometric normalization methods, this new method does not require high level programming and statistical skills, yet remains sensitive to landscape changes occurring over seasonal and inter-annual time scales. In addition, the TIC method maintains sensitivity to subtle changes in vegetation phenology and enables normalization even when invariant features are rare. While this normalization method allowed detection of a range of land use, land cover, and phenological/biophysical changes in the Siberian boreal forest region studied here, it is necessary to further examine images representing a wide variety of ecoregions to thoroughly evaluate the TIC method against other normalization schemes.     

Spectral sensing of foliar water conditions in two co-occurring conifer species: Pinus edulis and Juniperus monosperma

Many fundamental ecosystem properties and dynamics are determined by plant water stress, particularly in dryland ecosystems where water is usually limiting. Indeed, under severe drought, plant water stress and associated insect infestations can produce landscape-scale mortality. Despite the fundamental importance of plant water stress in determining properties and dynamics at ecosystem and landscape scales, approaches for remotely sensing plant water stress are largely lacking, particularly for conifers. We evaluated the remotely sensed detection of foliar drought stress in two conifer species, Pinus edulis and Juniperus monosperma, which are co-dominants of extensive-juniper woodlands in North America, the first of which experienced extensive mortality in association with a recent drought. Needle spectra were made on these species in the field using an integrating sphere and portable spectrometer. Two indices of foliar water condition, plant water content (% of dry mass) and plant water potential, were compared to five spectral analyses: continuum removal of the 970 and 1200 nm water absorption features, the Normalized Difference Water Index (NDWI), the Normalized Difference Vegetation Index (NDVI), and the red edge wavelength position. For P. edulis, plant water content was significantly correlated with four of the five indices: NDVI (R²=0.71) and NDWI (R²=0.68) which exhibited stronger relationships than 970 nm continuum removal (R²=0.57) or red edge position (R²=0.45). All five indices were significantly correlated with P. edulis water content when trees undergoing mortality were included in analyses (R²=0.60-0.93). Although the correlations were weaker than for plant water content, plant water potential was significantly correlated with NDWI (R²=0.49), 970 nm (R²=0.44), NDVI (R²=0.35), and red edge (R²=0.34); again all five indices had significant relationships when trees undergoing mortality were included (R²=0.51-0.86). The relationships were weaker for J. monosperma: water content was significantly related to 970 nm (R²=0.50) and 1200 nm (R²=0.37) continuums and NDVI (R²=0.33), while water potential was related only to 1200 nm (R²=0.40). Our results demonstrate a critical link between plant physiological characteristics tied to water stress and associated spectral signatures for two extensive co-occurring conifer species.     

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.     

A novel ocean color index to detect floating algae in the global oceans

Various types of floating algae have been reported in open oceans and coastal waters, yet accurate and timely detection of these relatively small surface features using traditional satellite data and algorithms has been difficult or even impossible due to lack of spatial resolution, coverage, revisit frequency, or due to inherent algorithm limitations. Here, a simple ocean color index, namely the Floating Algae Index (FAI), is developed and used to detect floating algae in open ocean environments using the medium-resolution (250- and 500-m) data from operational MODIS (Moderate Resolution Imaging Spectroradiometer) instruments. FAI is defined as the difference between reflectance at 859 nm (vegetation “red edge”) and a linear baseline between the red band (645 nm) and short-wave infrared band (1240 or 1640 nm). Through data comparison and model simulations, FAI has shown advantages over the traditional NDVI (Normalized Difference Vegetation Index) or EVI (Enhanced Vegetation Index) because FAI is less sensitive to changes in environmental and observing conditions (aerosol type and thickness, solar/viewing geometry, and sun glint) and can “see” through thin clouds. The baseline subtraction method provides a simple yet effective means for atmospheric correction, through which floating algae can be easily recognized and delineated in various ocean waters, including the North Atlantic Ocean, Gulf of Mexico, Yellow Sea, and East China Sea. Because similar spectral bands are available on many existing and planned satellite sensors such as Landsat TM/ETM+ and VIIRS (Visible Infrared Imager/Radiometer Suite), the FAI concept is extendable to establish a long-term record of these ecologically important ocean plants.     

Using lidar and effective LAI data to evaluate IKONOS and Landsat 7 ETM+ vegetation cover estimates in a ponderosa pine forest

Structural and functional analyses of ecosystems benefit when high accuracy vegetation coverages can be derived over large areas. In this study, we utilize IKONOS, Landsat 7 ETM+, and airborne scanning light detection and ranging (lidar) to quantify coniferous forest and understory grass coverages in a ponderosa pine (Pinus ponderosa) dominated ecosystem in the Black Hills of South Dakota. Linear spectral mixture analyses of IKONOS and ETM+ data were used to isolate spectral endmembers (bare soil, understory grass, and tree/shade) and calculate their subpixel fractional coverages. We then compared these endmember cover estimates to similar cover estimates derived from lidar data and field measures. The IKONOS-derived tree/shade fraction was significantly correlated with the field-measured canopy effective leaf area index (LAI_e) (r²=0.55, p<0.001) and with the lidar-derived estimate of tree occurrence (r²=0.79, p<0.001). The enhanced vegetation index (EVI) calculated from IKONOS imagery showed a negative correlation with the field measured tree canopy effective LAI and lidar tree cover response (r²=0.30, r=−0.55 and r²=0.41, r=−0.64, respectively; p<0.001) and further analyses indicate a strong linear relationship between EVI and the IKONOS-derived grass fraction (r²=0.99, p<0.001). We also found that using EVI resulted in better agreement with the subpixel vegetation fractions in this ecosystem than using normalized difference of vegetation index (NDVI). Coarsening the IKONOS data to 30 m resolution imagery revealed a stronger relationship with lidar tree measures (r²=0.77, p<0.001) than at 4 m resolution (r²=0.58, p<0.001). Unmixed tree/shade fractions derived from 30 m resolution ETM+ imagery also showed a significant correlation with the lidar data (r²=0.66, p<0.001). These results demonstrate the power of using high resolution lidar data to validate spectral unmixing results of satellite imagery, and indicate that IKONOS data and Landsat 7 ETM+ data both can serve to make the important distinction between tree/shade coverage and exposed understory grass coverage during peak summertime greenness in a ponderosa pine forest ecosystem.