Moderate Resolution Imaging Spectroradiometer Vegetation Indices Research Articles

Assessing Potential Climatic and Human Pressures in Indonesian Coastal Ecosystems Using a Spatial Data-Driven Approach

Blue carbon ecosystems are key for successful global climate change mitigation; however, they are one of the most threatened ecosystems on Earth. Thus, this study mapped the climatic and human pressures on the blue carbon ecosystems in Indonesia using multi-source spatial datasets. Data on moderate resolution imaging spectroradiometer (MODIS) ocean color standard mapped images, VIIRS (visible, infrared imaging radiometer suite) boat detection (VBD), global artificial impervious area (GAIA), MODIS surface reflectance (MOD09GA), MODIS land surface temperature (MOD11A2), and MODIS vegetation indices (MOD13A2) were combined using remote sensing and spatial analysis techniques to identify potential stresses. La Niña and El Niño phenomena caused sea surface temperature deviations to reach −0.5 to +1.2 °C. In contrast, chlorophyll-a deviations reached 22,121 to +0.5 mg m−3. Regarding fishing activities, most areas were under exploitation and relatively sustained. Concerning land activities, mangrove deforestation occurred in 560.69 km2 of the area during 2007–2016, as confirmed by a decrease of 84.9% in risk-screening environmental indicators. Overall, the potential pressures on Indonesia’s blue carbon ecosystems are varied geographically. The framework of this study can be efficiently adopted to support coastal and small islands zonation planning, conservation prioritization, and marine fisheries enhancement.

Deriving drought indices from MODIS vegetation indices (NDVI/EVI) and Land Surface Temperature (LST): Is data reconstruction necessary?

Droughts pose significant economic and ecological concerns, and considering climate change projections, timely monitoring and early warning based on satellite observations must be realized at regional to global scales. Nevertheless, whether data reconstruction is necessary to produce high-quality satellite-based time series data for drought monitoring and the data reconstruction approaches to be applied, if necessary, remain unclear. We attempted to fill this knowledge gap by investigating three widely used data reconstruction approaches, i.e., the Savitzky–Golay filter, harmonic analysis of time series (HANTS) and Whittaker Smoother, across the Lancang–Mekong river basin through the Moderate Resolution Imaging Spectroradiometer (MODIS) MOD13Q1 and MOD11A2 products for 2001–2018 and Google Earth Engine cloud platform. Several remote sensing drought indices based on the unreconstructed and reconstructed vegetation indices (VIs) and Land Surface Temperature (LST) were retrieved for basin-wide drought detection. The performance of the examined reconstruction approaches was evaluated using three statistical indices (coefficient of determination (R2), standard deviation of bias (BIAS(std)) and mean squared error (MSE)), spatial consistency with the reference dataset and capability to characterize the drought events. Data reconstruction considerably enhanced the drought index performances for drought detection; however, reconstruction was not necessary in all situations. The reconstructed drought indices exhibited higher R2 values (by 10.6–10.8%), lower BIAS(std) values (by −1.7– −12.5%), and smaller MSE values (by −5.8– −13.4%) compared to those of unreconstructed indices. For most evaluation indicators, HANTS outperformed the other methods, and Vegetation Condition Index (VCI) and Vegetation Health Index (VHI) outperformed the other drought indices. The findings highlight the importance of data reconstruction to detect and characterize drought events and the dependency of the performance of reconstruction methods on drought indices and evaluation metrics when using MODIS time series data.

MAIZE YIELD ESTIMATION IN KENYA USING MODIS

Abstract. Monitoring staple crop production can support agricultural research, business such as crop insurance, and government policy. Obtaining accurate estimates through field work is very expensive, and estimating it through remote sensing is promising. We estimated county-level maize yield for the 37 maize producing countries in Kenya from 2010 to 2017 using Moderate Resolution Imaging Spectroradiometer (MODIS) data. Support Vector Regression (SVR) and Random Forest (RF) were used to fit models with observed county level maize yield as a function of vegetation indices. The following five MODIS vegetation indices were used: green normalized difference vegetation index, normalized difference vegetation index, normalized difference moisture index, gross primary production, and fraction of photosynthetically active radiation. The models were evaluated with 5-fold leave one year out cross-validation. For SVR, R2 was 0.70, the Root Mean Square Error (RMSE) was 0.50 MT/ha and Mean Absolute Percentage Error (MAPE) was 27.6%. On the other hand for RF these were 0.69, 0.51 MT/ha and 29.3% respectively. These results are promising and should be tested in specific applications to understand if they are good enough for use.

Seasonal Comparisons of Himawari-8 AHI and MODIS Vegetation Indices over Latitudinal Australian Grassland Sites

The Advanced Himawari Imager (AHI) on board the Himawari-8 geostationary (GEO) satellite offers comparable spectral and spatial resolutions as low earth orbiting (LEO) sensors such as the Moderate Resolution Imaging Spectroradiometer (MODIS) and Visible Infrared Imaging Radiometer Suite (VIIRS) sensors, but with hypertemporal image acquisition capability. This raises the possibility of improved monitoring of highly dynamic ecosystems, such as grasslands, including fine-scale phenology retrievals from vegetation index (VI) time series. However, identifying and understanding how GEO VI temporal profiles would be different from traditional LEO VIs need to be evaluated, especially with the new generation of geostationary satellites, with unfamiliar observation geometries not experienced with MODIS, VIIRS, or Advanced Very High Resolution Radiometer (AVHRR) VI time series data. The objectives of this study were to investigate the variations in AHI reflectances and normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and two-band EVI (EVI2) in relation to diurnal phase angle variations, and to compare AHI VI seasonal datasets with MODIS VIs (standard and sun and view angle-adjusted VIs) over a functional range of dry grassland sites in eastern Australia. Strong NDVI diurnal variations and negative NDVI hotspot effects were found due to differential red and NIR band sensitivities to diurnal phase angle changes. In contrast, EVI and EVI2 were nearly insensitive to diurnal phase angle variations and displayed nearly flat diurnal profiles without noticeable hotspot influences. At seasonal time scales, AHI NDVI values were consistently lower than MODIS NDVI values, while AHI EVI and EVI2 values were significantly higher than MODIS EVI and EVI2 values, respectively. We attributed the cross-sensor differences in VI patterns to the year-round smaller phase angles and backscatter observations from AHI, in which the sunlit canopies induced a positive EVI/ EVI2 response and negative NDVI response. BRDF adjustments of MODIS VIs to solar noon and to the oblique view zenith angle of AHI resulted in strong cross-sensor convergence of VI values (R2 > 0.94, mean absolute difference <0.02). These results highlight the importance of accounting for cross-sensor observation geometries for generating compatible AHI and MODIS annual VI time series. The strong agreement found in this study shows promise in cross-sensor applications and suggests that a denser time series can be formed through combined GEO and LEO measurement synergies.

A double peak in the seasonality of California's photosynthesis as observed from space

Abstract. Solar-induced chlorophyll fluorescence (SIF) has been shown to be a powerful proxy for photosynthesis and gross primary productivity (GPP). The recently launched TROPOspheric Monitoring Instrument (TROPOMI) features the required spectral resolution and signal-to-noise ratio to retrieve SIF from space. Here, we present a downscaling method to obtain 500 m spatial resolution SIF over California. We report daily values based on a 14 d window. TROPOMI SIF data show a strong correspondence with daily GPP estimates at AmeriFlux sites across multiple ecosystems in California. We find a linear relationship between SIF and GPP that is largely invariant across ecosystems with an intercept that is not significantly different from zero. Measurements of SIF from TROPOMI agree with MODerate Resolution Imaging Spectroradiometer (MODIS) vegetation indices – the normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and near-infrared reflectance of vegetation index (NIRv) – at annual timescales but indicate different temporal dynamics at monthly and daily timescales. TROPOMI SIF data show a double peak in the seasonality of photosynthesis, a feature that is not present in the MODIS vegetation indices. The different seasonality in the vegetation indices may be due to a clear-sky bias in the vegetation indices, whereas previous work has shown SIF to have a low sensitivity to clouds and to detect the downregulation of photosynthesis even when plants appear green. We further decompose the spatiotemporal patterns in the SIF data based on land cover. The double peak in the seasonality of California's photosynthesis is due to two processes that are out of phase: grasses, chaparral, and oak savanna ecosystems show an April maximum, while evergreen forests peak in June. An empirical orthogonal function (EOF) analysis corroborates the phase offset and spatial patterns driving the double peak. The EOF analysis further indicates that two spatiotemporal patterns explain 84 % of the variability in the SIF data. Results shown here are promising for obtaining global GPP at sub-kilometer spatial scales and identifying the processes driving carbon uptake.

Mapping high temperature damaged area of paddy rice along the Yangtze River using Moderate Resolution Imaging Spectroradiometer data

ABSTRACT The subtropical high temperature anomaly of 2013 in China, adversely affected paddy rice physiology, especially at the heading, flowering and grain filling stages, and left untold damage and significant rice yield reductions in that year. Considering the importance of data on crop failure induced by abnormal meteorological events, in understanding the potential impacts of climate change on agriculture and food security, and in support of crop insurance policies, the current study employs satellite remote sensing data to map and estimate the areal extent of high-temperature-induced damage on paddy rice in 2013, along the Yangtze River, a major rice producing area in China. In this regard, the daily mean and maximum air temperatures were first estimated using the Moderate Resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature (LST) products, MODIS Vegetation Indices (VI) products, altitude information from a Digital Elevation Model (DEM), and air temperature measurements at meteorological stations located within the Yangtze River Basin from 2000 to 2016, to derive fully covered mean and maximum daily air temperatures. Paddy rice cultivated areas were then identified based on the seasonal profiles of MODIS Enhanced Vegetation Index (EVI) and Land Surface Water Index (LSWI) in 2013, and map errors at provincial scales were no more than 8.7%. Finally, paddy rice growth characteristics were combined with MODIS VI products to identify heading and flowering stages, and the time series of daily mean and maximum air temperatures, paddy rice phenology and cultivated area, were integrated to map and estimate the rice crop area damaged by the anomalous high temperature. Our results show that all provinces along the Yangtze River, this study investigated, recorded a substantial damaged rice crop area, following the 2013 subtropical temperature high, and Anhui Province was the hardest hit with a damaged rice area of 63.6%. For the entire study area, a damage extent of about 5,656,000 ha, accounting for 35.9% of total paddy rice planted area, was recorded. This indicates a significant agro-meteorological disaster, and in an era of increasing climate variability and weather extremes, studies of this nature would greatly benefit global food security planning and management policies.

Detection of Vegetation Cover Change in Renewable Energy Development Zones of Southern California Using MODIS NDVI Time Series Analysis, 2000 to 2018

New solar energy facilities on public lands in the deserts of southern California are being monitored long-term to detect environmental impacts. For this purpose, we have developed a framework for detecting changes in vegetation cover region-wide using greenness index data sets from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensor. This study focused on three sites, Joshua Tree National Park (JOTR), Mojave National Preserve (MOJA), and a proximal group of solar energy Development Focus Areas (DFAs). Three MODIS vegetation indices (VIs), the normalized difference (NDVI), enhanced (EVI), and soil-adjusted (SAVI), all at 250-m spatial resolution, were evaluated using the Breaks for Additive Season and Trend (BFAST) methodology to estimate significant time series shifts (“breakpoints”) in green vegetation cover, from February 2000 to May 2018. The sample cross-correlation function with local precipitation records and comparison with timing of wildfires near the study sites for breakpoint density (proportion of area with a breakpoint) showed that NDVI had the strongest response and hence greatest sensitivity to these major disturbances compared to EVI and SAVI, supporting its use over the other VIs for subsequent analysis. Time series of NDVI breakpoint change densities for individual solar energy DFAs did not have a consistent vegetation response following construction. Bootstrap-derived 95% confidence intervals show that the DFAs have significantly larger kurtosis and standard deviation in positive NDVI breakpoint distribution than protected National Park System (NPS) sites, but no significant difference appeared in the negative distribution among all sites. The inconsistent postconstruction NDVI signal and the large number of detected breakpoints across all three sites suggested that the largest shifts in greenness are tied to seasonal and total annual precipitation amounts. Further results indicated that existing site-specific conditions are the main control on vegetation response, mostly driven by the history of human disturbances in DFAs. Although the results do not support persistent breakpoints in solar energy DFAs, future work should seek to establish links between statistical significance and physical significance through ground-based studies to provide a more robust interpretation.

Analysis of NVDI variability in response to precipitation and air temperature in different regions of Iraq, using MODIS vegetation indices

Iraq, the land of two rivers, has a history that extends back millennia and is the subject of much archaeological research. However, little environmental research has been carried out, and as such relatively little is known about the interaction between Iraq’s vegetation and climate. This research serves to fill this knowledge gap by investigating the relationship between the Normalized Difference Vegetation Index (NDVI) and two climatic factors (precipitation and air temperature) over the last decade. The precipitation and air temperature datasets are from the Water and Global Change Forcing Data ERA-Interim (WFDEI), and the NDVI dataset was extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS) at 250 m spatial resolution and 16 day temporal resolution. Three different climatic regions in Iraq, Sulaymaniyah, Wasit, and Basrah, were selected for the period of 2001–2015. This is the first study to compare these regions in Iraq, and one of only a few investigating vegetation’s relationship with multiple climatic factors, including precipitation and air temperature, particularly in a semi-arid region. The interannual, intra-annual and seasonal variability for each region is analysed to compare the different responses of vegetation growth to climatic factors. Correlations between NDVI and climatic factors are also included. Plotting annual cycles of NDVI and precipitation reveals a coherent onset, fluctuation (peak and decline), with a time lag of 4 months for Sulaymaniyah and Wasit (while for the Basrah region, high temperatures and a short rainy season was observed). The correlation coefficients between NDVI and precipitation are relatively high, especially in Sulaymaniyah, and the largest positive correlation was (0.8635) with a time lag of 4 months. The phenological transition points range between 3 and 4 month time lag; this corresponds to the duration of maturity of the vegetation. However, when correlated with air temperature, NDVI experiences an inverse relationship, although not as strong as that of NDVI and precipitation; the highest negative correlation was observed in Wasit with a time lag of 2 months (− 0.7562). The results showed that there is a similarity between temporal patterns of NDVI and precipitation. This similarity is stronger than that of NDVI and air temperature, so it can be concluded that NDVI is a sensitive indicator of the inter-annual variability of precipitation and that precipitation constitutes the primary factor in germination while the air temperature acts with a lesser effect.

Examining the short-term impacts of diverse management practices on plant phenology and carbon fluxes of Old World bluestems pasture

Burning, grazing, and baling (hay harvesting) are common management practices in grassland. To develop and adopt sustainable management practices, it is essential to better understand and quantify the impacts of management practices on plant phenology and carbon fluxes. In this study, we combined multiple data sources, including in-situ PhenoCam digital images, eddy covariance data, and satellite data (Landsat and Moderate Resolution Imaging Spectroradiometer (MODIS)) to examine the impacts of burning, baling, and grazing on canopy dynamics, plant phenology, and carbon fluxes in a pasture in El Reno, Oklahoma in 2014. Landsat images were used to assess the baling area and the trajectory of vegetation recovery. MODIS vegetation indices (VIs) were used in the Vegetation Photosynthesis Model (VPM) to estimate gross primary production (GPPVPM) at a MODIS pixel for the flux tower (baled) site. For comparison between baled and unbaled conditions, we used MODIS VIs for a neighbor MODIS pixel (unbaled) and ran VPM. Daily PhenoCam images and green chromatic coordinate (GCC) tracked canopy dynamics and plant phenology well. The grassland greened up immediately after burning in April. GCC values showed two peaks with the similar magnitude because of quick recovery of grassland after baling. Satellite-derived VIs and GPPVPM showed that the pasture recovered in one month after baling. The GPPVPM matched well (R2=0.89) with the eddy covariance-derived GPP (GPPEC). Grazing in the late growing season did not influence plant phenology (VIs and GCC) and carbon uptake (GPP) as plants were in the late growing stage. Neither did it affect GPP differently in those two conditions because of even grazing intensity. The reduction in GPP after baling was compensated by higher GPP after large rain events in late July and early September, causing little seasonal differences in GPP (-0.002gCm−2day−1) between the baled and unbaled conditions. Interactions of different management practices with climate make it complicated to understand the impacts of different management practices on carbon dynamics and plant phenology. Thus, it is necessary to further investigate the responses of pastures to different management practices under different climate regimes at multiple temporal and spatial scales.

MODIS vegetation products as proxies of photosynthetic potential along a gradient of meteorologically and biologically driven ecosystem productivity

Abstract. A direct relationship between gross ecosystem productivity (GEP) estimated by the eddy covariance (EC) method and Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation indices (VIs) has been observed in many temperate and tropical ecosystems. However, in Australian evergreen forests, and particularly sclerophyll and temperate woodlands, MODIS VIs do not capture seasonality of GEP. In this study, we re-evaluate the connection between satellite and flux tower data at four contrasting Australian ecosystems, through comparisons of GEP and four measures of photosynthetic potential, derived via parameterization of the light response curve: ecosystem light use efficiency (LUE), photosynthetic capacity (Pc), GEP at saturation (GEPsat), and quantum yield (α), with MODIS vegetation satellite products, including VIs, gross primary productivity (GPPMOD), leaf area index (LAIMOD), and fraction of photosynthetic active radiation (fPARMOD). We found that satellite-derived biophysical products constitute a measurement of ecosystem structure (e.g. leaf area index – quantity of leaves) and function (e.g. leaf level photosynthetic assimilation capacity – quality of leaves), rather than GEP. Our results show that in primarily meteorological-driven (e.g. photosynthetic active radiation, air temperature, and/or precipitation) and relatively aseasonal ecosystems (e.g. evergreen wet sclerophyll forests), there were no statistically significant relationships between GEP and satellite-derived measures of greenness. In contrast, for phenology-driven ecosystems (e.g. tropical savannas), changes in the vegetation status drove GEP, and tower-based measurements of photosynthetic activity were best represented by VIs. We observed the highest correlations between MODIS products and GEP in locations where key meteorological variables and vegetation phenology were synchronous (e.g. semi-arid Acacia woodlands) and low correlation at locations where they were asynchronous (e.g. Mediterranean ecosystems). However, we found a statistical significant relationship between the seasonal measures of photosynthetic potential (Pc and LUE) and VIs, where each ecosystem aligns along a continuum; we emphasize here that knowledge of the conditions in which flux tower measurements and VIs or other remote sensing products converge greatly advances our understanding of the mechanisms driving the carbon cycle (phenology and climate drivers) and provides an ecological basis for interpretation of satellite-derived measures of greenness.

Assessing the relationship between shire winter crop yield and seasonal variability of the MODIS NDVI and EVI images

Australian researchers have been developing robust yield estimation models, based mainly on the crop growth response to water availability during the crop season. However, knowledge of spatial distribution of yields within and across the production regions can be improved by the use of remote sensing techniques. Images of Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation indices, available since 1999, have the potential to contribute to crop yield estimation. The objective of this study was to analyse the relationship between winter crop yields and the spectral information available in MODIS vegetation index images at the shire level. The study was carried out in the Jondaryan and Pittsworth shires, Queensland , Australia . Five years (2000 to 2004) of 250m resolution, 16-day composite of MODIS Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) images were used during the winter crop season (April to November). Seasonal variability of the profiles of the vegetation index images for each crop season using different regions of interest (cropping mask) were displayed and analysed. Correlation analysis between wheat and barley yield data and MODIS image values were also conducted. The results showed high seasonal variability in the NDVI and EVI profiles, and the EVI values were consistently lower than those of the NDVI. The highest image values were observed in 2003 (in contrast to 2004), and were associated with rainfall amount and distribution. The seasonal variability of the profiles was similar in both shires, with minimum values in June and maximum values at the end of August. NDVI and EVI images showed sensitivity to seasonal variability of the vegetation and exhibited good association (e.g. r = 0.84, r = 0.77) with winter crop yields.

Rice Crop Monitoring and Yield Assessment with MODIS 250m Gridded Vegetation Products: A Case Study of Sa Kaeo Province, Thailand

Abstract. Billions of people in the world depend on rice as a staple food and as an income-generating crop. Asia is the leader in rice cultivation and it is necessary to maintain an up-to-date rice-related database to ensure food security as well as economic development. This study investigates general applicability of high temporal resolution Moderate Resolution Imaging Spectroradiometer (MODIS) 250m gridded vegetation product for monitoring rice crop growth, mapping rice crop acreage and analyzing crop yield, at the province-level. The MODIS 250m Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) time series data, field data and crop calendar information were utilized in this research in Sa Kaeo Province, Thailand. The following methodology was used: (1) data pre-processing and rice plant growth analysis using Vegetation Indices (VI) (2) extraction of rice acreage and start-of-season dates from VI time series data (3) accuracy assessment, and (4) yield analysis with MODIS VI. The results show a direct relationship between rice plant height and MODIS VI. The crop calendar information and the smoothed NDVI time series with Whittaker Smoother gave high rice acreage estimation (with 86% area accuracy and 75% classification accuracy). Point level yield analysis showed that the MODIS EVI is highly correlated with rice yield and yield prediction using maximum EVI in the rice cycle predicted yield with an average prediction error 4.2%. This study shows the immense potential of MODIS gridded vegetation product for keeping an up-to-date Geographic Information System of rice cultivation.

Global variability of midtropospheric carbon dioxide as measured by the Atmospheric Infrared Sounder

The Atmospheric Infrared Sounder (AIRS) on the EOS Aqua spacecraft provides accurate and consistent measurements of midtropospheric carbon dioxide (CO2) with global monthly coverage. The data are widely used for studies of vertical transport of CO2 due to large-scale dynamics (e.g., ENSO, MJO, and the Walker Circulation). The purpose of this paper is to characterize the response of CO2 in the midtroposphere, at the altitudes where AIRS is most sensitive, to geophysical changes at the surface across the globe. Our findings confirm that surface factors, as well as weather and climate patterns, impact the global variability of midtropospheric CO2 as observed by AIRS. Despite a phase lag and a reduction in the seasonal amplitude observed in AIRS CO2 relative to surface CO2 measurements in the Northern Hemisphere, a significant correlation is observed between regional variability of CO2 from AIRS and Moderate Resolution Imaging Spectroradiometer (MODIS)-derived Gross Primary Productivity at the surface, primarily in the high-latitude boreal forests during the peak of the growing season (July). A video of global AIRS CO2 and MODIS vegetation index clearly shows the seasonal drawdown of CO2 from the midtroposphere over highly vegetated areas in the northern latitudes. In the Southern Hemisphere, we see higher amplitude in the seasonal cycle, with the phase leading that of the surface. Both are indicative of interhemispheric transport.

Comparison between SMOS Vegetation Optical Depth products and MODIS vegetation indices over crop zones of the USA

The Soil Moisture and Ocean Salinity (SMOS) mission provides multi-angular, dual-polarised brightness temperatures at 1.4GHz, from which global soil moisture and vegetation optical depth (τ) products are retrieved. This paper presents a study of SMOS' τ product in 2010 and 2011 for crop zones of the USA. Retrieved τ values for 504 crop nodes were compared to optical/IR vegetation indices from the MODIS (Moderate Resolution Imaging Spectroradiometer) satellite sensor, including the Normalised Difference Vegetation Index (NDVI), Enhanced Vegetation Index (EVI), Leaf Area Index (LAI), and a Normalised Difference Water Index (NDWI) product. τ values were observed to increase during the growing season and decrease during senescence in these areas, as did MODIS vegetation indices. SMOS' τ values generally peaked later than MODIS LAI values, with an estimated time difference of about 19days. A linear regression between τ and the MODIS products was carried out for each node and values of the determination coefficient, R2, slope, b′ and intercept, b″ were found. The average R2 value varied from 0.32 to 0.35 for the different vegetation indices. The linear regression between LAI and τ produced an average slope of b′=0.06, and an average intercept of b″=0.14. The effects of crop fraction and dominant crop type were investigated and crop fraction was found to have a low effect on R2 values. R2values appeared to be lower for wheat and hay and higher for corn. b′ and b″ values had higher standard deviations for wheat but were generally close to the mean values for corn, soybean and hay.

Drought and elevation effects on MODIS vegetation indices in northern Arizona ecosystems

Northern Arizona ecosystems are particularly sensitive to plant-available moisture and have experienced a severe drought with considerable impacts on ecosystems from desert shrub and grasslands to pinyon-juniper and conifer forests. Long-term time-series from the Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) are used to monitor recent regional vegetation activity and temporal patterns across various ecosystems. Surface air temperature, solar radiation and precipitation are used to represent meteorological anomalies and to investigate associated impacts on vegetation greenness. Vegetation index anomalies in the northern Arizona ecosystem have a decreasing trend with increasing surface air temperature and decreasing precipitation. MODIS NDVI and EVI anomalies are likely sensitive to the amount of rainfall for northern Arizona ecosystem conditions, whereas inter-annual variability of surface air temperature accounts for MODIS NDVI anomaly variation. The higher elevation area shows the slow vegetation recovery through trend analysis from MODIS vegetation indices for 2000–2011 within the study domain and along elevation.

Evaluation of optical remote sensing to estimate actual evapotranspiration and canopy conductance

We compared estimates of actual evapotranspiration (ET) produced with six different vegetation measures derived from the MODerate resolution Imaging Spectroradiometer (MODIS) and three contrasting estimation approaches using measurements from eddy covariance flux towers at 16 FLUXNET sites located over six different land cover types. The aim was to assess optimal approaches in using optical remote sensing to estimate ET. The first two approaches directly regressed various MODIS vegetation indices (VIs) and products such as leaf area index (LAI) and fraction of photosynthetically active radiation (fPAR) with ET and evaporative fraction (EF). In the third approach, the Penman–Monteith (PM) equation was inverted to obtain surface conductance (Gs), for dry plant canopies. The Gs values were then regressed against the MODIS data products and used to parameterize the PM equation for retrievals of ET. Jack-Knife cross-validation was used to evaluate the various regression models against observed ET. The PM-Gs approach provided the lowest root mean square error (RMSE), and highest determination coefficients (R2) across all sites, with an average RMSE=38Wm−2 and R2=0.72. Direct regressions of observed ET against the VIs resulted in an average RMSE=60Wm−2 and R2=0.22, while the EF regressions an average RMSE=42Wm−2 and R2=0.64. The MODIS LAI and fPAR product produced the poorest estimates of ET (RMSE>44Wm−2 and R2<0.6); while the VIs each performed best for some of the land cover types. The enhanced vegetation index (EVI) produced the best ET estimates for evergreen needleleaf forest (RMSE=28.4Wm−2, R2=0.66). The normalized difference vegetation index (NDVI) best estimated ET in grassland (RMSE=23.8Wm−2 and R2=0.68), cropland (RMSE=29.2Wm−2 and R2=0.86) and woody savannas (RMSE=25.4Wm−2 and R2=0.82), while the VI-based crop coefficient (Kc) yielded the best estimates for evergreen and deciduous broadleaf forests (RMSE=27Wm−2 and R2=0.7 in both cases). Using the ensemble-average of ET as estimated using NDVI, EVI and Kc we computed global grids of dry canopy conductance (Gc) from which annual statistics were extracted to characterise different functional types. The resulting Gc values can be used to parameterize land surface models.

Winter wheat mapping using temporal signatures of MODIS vegetation index data

Because most land-cover types have distinct seasonal changes and corresponding reflectance characteristics in remotely sensed images, the signatures in time-series data are useful for discriminating different land covers. Although temporal signatures have been used to classify different land-cover types, they have not been fully exploited to classify specific crops, and the influence of low resolution should be evaluated. The aims of this study were to seek an effective method to classify specific crops using the temporal signatures in coarse time-series data and to examine the applicability of the data for crop classification as well. A winter wheat-producing region in China was selected for this case study. Moderate-Resolution Imaging Spectroradiometer (MODIS) 8-day composite land surface reflectance product (MOD09Q1) data with a 250 m spatial resolution were used to calculate the vegetation index data, which was applied to detect the properties of live green plants. The noise in the time series was filtered to minimize the classification uncertainties. The curve shape in the time-series vegetation index profile was used as the major metric to classify winter wheat, and other phenological metrics extracted from the data were used conjunctly as auxiliary functions to improve the separability. The metrics for winter wheat classification were quantified in the large fields with relatively pure pixels. Winter wheat was successfully extracted from the MODIS vegetation index data, and the MODIS-derived result was validated with a fine-resolution (19.5 m) thematic map derived from images collected by the charge-coupled device sensor on board the China–Brazil Earth Resources Satellite (CBERS). It showed that the MODIS-derived result had inevitable low-resolution bias, and the errors of commission and omission were 32.3 and 33.8%, respectively. The overall classification effect of the MODIS-derived result relied upon the distribution of pixel purity in the study area.

Winter wheat area estimation from MODIS-EVI time series data using the Crop Proportion Phenology Index

The global distribution of croplands is of critical interest to a wide group of end-users. Different crops have their own representative phenological stages during their growing seasons, which differ considerably from other natural vegetation types. During the last decade, the Moderate Resolution Imaging Spectroradiometer (MODIS) has become a key tool for vegetation monitoring because of its high temporal resolution, extensive scope, and rapid availability of various products. However, mixed pixels caused by the moderate spatial resolution produce significant errors in crop area estimation. Here we propose a Crop Proportion Phenology Index (CPPI) to express the quantitative relationship between the MODIS vegetation index (VI) time series and winter wheat crop area. The utility of this index was tested in two experimental areas in China: one around Tongzhou and the other around Shuyang, as representative districts around a metropolis and a rural area, respectively. The CPPI performed well in these two regions, with the root mean square error (RMSE) in fractional crop area predictions ranging roughly from 15% in the individual pixels to 5% above 6.25km2. The training samples containing mixtures of crop types mitigated the challenges of pure end-member selection in a spectral mixture analysis. A small number of training samples are sufficient to generate the CPPI, which is adaptable to other crop types and larger regions. Estimating the CPPI parameters across larger spatial scales helped improve the stability of the model.

Seasonal variations in litter production and its relation with MODIS vegetation indices in a semi-deciduous forest of Mato Grosso

Litter production is related to canopy processes including the timing and amount of leaf development, reproduction and net primary production. However, quantifying spatial and temporal patterns in litter production is complicated in Amazonian semi-deciduous forests because of high spatial heterogeneity and seasonal variation in rainfall. Here, we use monthly measurements of litter production and Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) composites over a 6-year period to assess whether MODIS can be used to quantify litter production of tropical semi-deciduous forests. Original MODIS NDVI and EVI values were poorly related to the seasonal periodicity in litter production, but after using singular spectrum analysis (SSA) signal extraction techniques, clear relationships between litter production and the NDVI and EVI emerged. These results indicate that MODIS NDVI and EVI data are useful for detecting temporal patterns in litter production for Amazonian semi-deciduous forests if signal extraction analyses such as SSA are conducted.

Assessing viewing and illumination geometry effects on the MODIS vegetation index (MOD13Q1) time series: implications for monitoring phenology and disturbances in forest communities in Queensland, Australia

Time series analysis of satellite data can be used to monitor temporal dynamics of forested environments, thus providing spatial data for a range of forest science, monitoring and management issues. The moderate resolution imaging spectroradiometer (MODIS) vegetation index (MOD13Q1) product has potential for monitoring forest dynamics and disturbances. However, the suitability of the product to accurately measure temporal changes due to phenology and disturbances is questionable as the effects of variable solar and viewing geometry have not been removed from these data. This study aimed to examine the impact that viewing and illumination geometry differences had on MOD13Q1 vegetation index values, and their subsequent ability to map changes arising from phenology and disturbances in a number of forest communities in Queensland, Australia. MOD13Q1 normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) were compared to normalized NDVI and EVI (NDVInormalized and EVInormalized), which were derived from the reflectance modelled from a bidirectional reflectance distribution function (BRDF)/albedo parameters product (MCD43A1) using fixed viewing and illumination geometry. Time series plots of the vegetation index values from a number of pixels representing different forest types and known disturbances showed that the NDVInormalized time series was more effective at capturing the changes in vegetation than the NDVI. MOD13Q1 NDVI showed higher seasonal amplitude, but was less accurate at capturing phenology and disturbances compared to the NDVInormalized. The EVI was less affected by variable viewing and illumination geometry in terms of amplitude, but was affected in terms of time shift in periodicities providing erroneous information on phenology. More studies in different environments with appropriate vegetation phenology reference data will be needed to confirm these observations.