Long-term global data sets of vegetation Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation absorbed by vegetation (FPAR) are critical to monitoring global vegetation dynamics and for modeling exchanges of energy, mass and momentum between the land surface and planetary boundary layer. LAI and FPAR are also state variables in hydrological, ecological, biogeochemical and crop-yield models. The generation, evaluation and an example case study documenting the utility of 30-year long data sets of LAI and FPAR are described in this article. A neural network algorithm was first developed between the new improved third generation Global Inventory Modeling and Mapping Studies (GIMMS) Normalized Difference Vegetation Index (NDVI3g) and best-quality Terra Moderate Resolution Imaging Spectroradiometer (MODIS) LAI and FPAR products for the overlapping period 2000–2009. The trained neural network algorithm was then used to generate corresponding LAI3g and FPAR3g data sets with the following attributes: 15-day temporal frequency, 1/12 degree spatial resolution and temporal span of July 1981 to December 2011. The quality of these data sets for scientific research in other disciplines was assessed through (a) comparisons with field measurements scaled to the spatial resolution of the data products, (b) comparisons with broadly-used existing alternate satellite data-based products, (c) comparisons to plant growth limiting climatic variables in the northern latitudes and tropical regions, and (d) correlations of dominant modes of interannual variability with large-scale circulation anomalies such as the EI Niño-Southern Oscillation and Arctic Oscillation. These assessment efforts yielded results that attested to the suitability of these data sets for research use in other disciplines. The utility of these data sets is documented by comparing the seasonal profiles of LAI3g with profiles from 18 state-of-the-art Earth System Models: the models consistently overestimated the satellite-based estimates of leaf area and simulated delayed peak seasonal values in the northern latitudes, a result that is consistent with previous evaluations of similar models with ground-based data. The LAI3g and FPAR3g data sets can be obtained freely from the NASA Earth Exchange (NEX) website.

Phenology is an indicator of crop growth conditions, and is correlated with crop yields. In this study, a phenological approach based on a remote sensing vegetation index was explored to predict the yield in 314 counties within the US Corn Belt, divided into semi-arid and non-semi-arid regions. The Moderate Resolution Imaging Spectroradiometer (MODIS) data product MOD09Q1 was used to calculate the normalized difference vegetation index (NDVI) time series. According to the NDVI time series, we divided the corn growing season into four growth phases, calculated phenological information metrics (duration and rate) for each growth phase, and obtained the maximum correlation NDVI (Max-R²). Duration and rate represent crop growth days and rate, respectively. Max-R² is the NDVI value with the most significant correlation with corn yield in the NDVI time series. We built three groups of yield regression models, including univariate models using phenological metrics and Max-R², and multivariate models using phenological metrics, and multivariate models using phenological metrics combined with Max-R² in the whole, semi-arid, and non-semi-arid regions, respectively, and compared the performance of these models. The results show that most phenological metrics had a statistically significant (p < 0.05) relationship with corn yield (maximum R² = 0.44). Models established with phenological metrics realized yield prediction before harvest in the three regions with R² = 0.64, 0.67, and 0.72. Compared with the univariate Max-R² models, the accuracy of models built with Max-R² and phenology metrics improved. Thus, the phenology metrics obtained from MODIS-NDVI accurately reflect the corn characteristics and can be used for large-scale yield prediction. Overall, this study showed that phenology metrics derived from remote sensing vegetation indexes could be used as crop yield prediction variables and provide a reference for data organization and yield prediction with physical crop significance.

As the Earth’s population continues to grow and demand for food increases, the need for improved and timely information related to the properties and dynamics of global agricultural systems is becoming increasingly important. Global land cover maps derived from satellite data provide indispensable information regarding the geographic distribution and areal extent of global croplands. However, land use information, such as cropping intensity (defined here as the number of cropping cycles per year), is not routinely available over large areas because mapping this information from remote sensing is challenging. In this study, we present a simple but efficient algorithm for automated mapping of cropping intensity based on data from NASA’s (NASA: The National Aeronautics and Space Administration) MODerate Resolution Imaging Spectroradiometer (MODIS). The proposed algorithm first applies an adaptive Savitzky-Golay filter to smooth Enhanced Vegetation Index (EVI) time series derived from MODIS surface reflectance data. It then uses an iterative moving-window methodology to identify cropping cycles from the smoothed EVI time series. Comparison of results from our algorithm with national survey data at both the provincial and prefectural level in China show that the algorithm provides estimates of gross sown area that agree well with inventory data. Accuracy assessment comparing visually interpreted time series with algorithm results for a random sample of agricultural areas in China indicates an overall accuracy of 91.0% for three classes defined based on the number of cycles observed in EVI time series. The algorithm therefore appears to provide a straightforward and efficient method for mapping cropping intensity from MODIS time series data.