We estimate global terrestrial gross primary production (GPP) based on models that use satellite data within a simplified light-use efficiency framework that does not rely upon other meteorological inputs. Satellite-based geometry-adjusted reflectances are from the MODerate-resolution Imaging Spectroradiometer (MODIS) and provide information about vegetation structure and chlorophyll content at both high temporal (daily to monthly) and spatial (∼1 km) resolution. We use satellite-derived solar-induced fluorescence (SIF) to identify regions of high productivity crops and also evaluate the use of downscaled SIF to estimate GPP. We calibrate a set of our satellite-based models with GPP estimates from a subset of distributed eddy covariance flux towers (FLUXNET 2015). The results of the trained models are evaluated using an independent subset of FLUXNET 2015 GPP data. We show that variations in light-use efficiency (LUE) with incident PAR are important and can be easily incorporated into the models. Unlike many LUE-based models, our satellite-based GPP estimates do not use an explicit parameterization of LUE that reduces its value from the potential maximum under limiting conditions such as temperature and water stress. Even without the parameterized downward regulation, our simplified models are shown to perform as well as or better than state-of-the-art satellite data-driven products that incorporate such parameterizations. A significant fraction of both spatial and temporal variability in GPP across plant functional types can be accounted for using our satellite-based models. Our results provide an annual GPP value of ∼140 Pg C year ${}^{-1}$ for 2007 that is within the range of a compilation of observation-based, model, and hybrid results, but is higher than some previous satellite observation-based estimates.

Amazon forests experienced two severe droughts at the beginning of the 21st century: one in 2005 and the other in 2010. How Amazon forests responded to these droughts is critical for the future of the Earth’s climate system. It is only possible to assess Amazon forests’ response to the droughts in large areal extent through satellite remote sensing. Here, we used the Multi-Angle Implementation of Atmospheric Correction (MAIAC) Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation index (VI) data to assess Amazon forests’ response to droughts, and compared the results with those from the standard (Collection 5 and Collection 6) MODIS VI data. Overall, the MAIAC data reveal more realistic Amazon forests inter-annual greenness dynamics than the standard MODIS data. Our results from the MAIAC data suggest that: (1) the droughts decreased the greenness (i.e., photosynthetic activity) of Amazon forests; (2) the Amazon wet season precipitation reduction induced by El Niño events could also lead to reduced photosynthetic activity of Amazon forests; and (3) in the subsequent year after the water stresses, the greenness of Amazon forests recovered from the preceding decreases. However, as previous research shows droughts cause Amazon forests to reduce investment in tissue maintenance and defense, it is not clear whether the photosynthesis of Amazon forests will continue to recover after future water stresses, because of the accumulated damages caused by the droughts.

Leaf area index (LAI) and fraction of photosynthetically active radiation (FPAR) absorbed by vegetation are key variables in many global models of climate, hydrology, biogeochemistry, and ecology. These parameters are being operationally produced from Terra and Aqua MODIS bidirectional reflectance factor (BRF) data. The MODIS science team has developed, and plans to release, a new version of the BRF product using the multi-angle implementation of atmospheric correction (MAIAC) algorithm from Terra and Aqua MODIS observations. This paper presents analyses of LAI and FPAR retrievals generated with the MODIS LAI/FPAR operational algorithm using Terra MAIAC BRF data. Direct application of the operational algorithm to MAIAC BRF resulted in an underestimation of the MODIS Collection 6 (C6) LAI standard product by up to 10%. The difference was attributed to the disagreement between MAIAC and MODIS BRFs over the vegetation by −2% to +8% in the red spectral band, suggesting different accuracies in the BRF products. The operational LAI/FPAR algorithm was adjusted for uncertainties in the MAIAC BRF data. Its performance evaluated on a limited set of MAIAC BRF data from North and South America suggests an increase in spatial coverage of the best quality, high-precision LAI retrievals of up to 10%. Overall MAIAC LAI and FPAR are consistent with the standard C6 MODIS LAI/FPAR. The increase in spatial coverage of the best quality LAI retrievals resulted in a better agreement of MAIAC LAI with field data compared to the C6 LAI product, with the RMSE decreasing from 0.80 LAI units (C6) down to 0.67 (MAIAC) and the R² increasing from 0.69 to 0.80. The slope (intercept) of the satellite-derived vs. field-measured LAI regression line has changed from 0.89 (0.39) to 0.97 (0.25).

The Moderate Resolution Imaging Spectroradiometer (MODIS) provides daily global coverage, but the 10 km resolution of its aerosol optical depth (AOD) product is not suitable for studying spatial variability of aerosols in urban areas. Recently, a new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm was developed for MODIS which provides AOD at 1 km resolution. Using MAIAC data, the relationship between MAIAC AOD and PM_2.5 as measured by the 27 EPA ground monitoring stations was investigated. These results were also compared to conventional MODIS 10 km AOD retrievals (MOD04) for the same days and locations. The coefficients of determination for MOD04 and for MAIAC are R² =0.45 and 0.50 respectively, suggested that AOD is a reasonably good proxy for PM_2.5 ground concentrations. Finally, we studied the relationship between PM_2.5 and AOD at the intra-urban scale (⩽10 km) in Boston. The fine resolution results indicated spatial variability in particle concentration at a sub-10 kilometer scale. A local analysis for the Boston area showed that the AOD-PM_2.5 relationship does not depend on relative humidity and air temperatures below ∼7 °C. The correlation improves for temperatures above 7–16 °C. We found no dependence on the boundary layer height except when the former was in the range 250–500 m. Finally, we apply a mixed effects model approach to MAIAC aerosol optical depth (AOD) retrievals from MODIS to predict PM_2.5 concentrations within the greater Boston area. With this approach we can control for the inherent day-to-day variability in the AOD-PM_2.5 relationship, which depends on time-varying parameters such as particle optical properties, vertical and diurnal concentration profiles and ground surface reflectance. Our results show that the model-predicted PM_2.5 mass concentrations are highly correlated with the actual observations (out-of-sample R² of 0.86). Therefore, adjustment for the daily variability in the AOD-PM_2.5 relationship provides a means for obtaining spatially-resolved PM_2.5 concentrations.

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 (R² > 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.