Satellite observations of Amazon forests show seasonal and interannual variations, but the underlying biological processes remain debated. Here we combined radiative transfer models (RTMs) with field observations of Amazon forest leaf and canopy characteristics to test three hypotheses for satellite-observed canopy reflectance seasonality: seasonal changes in leaf area index, in canopy-surface leafless crown fraction and/or in leaf demography. Canopy RTMs (PROSAIL and FLiES), driven by these three factors combined, simulated satellite-observed seasonal patterns well, explaining c. 70% of the variability in a key reflectance-based vegetation index (MAIAC EVI, which removes artifacts that would otherwise arise from clouds/aerosols and sun-sensor geometry). Leaf area index, leafless crown fraction and leaf demography independently accounted for 1, 33 and 66% of FLiES-simulated EVI seasonality, respectively. These factors also strongly influenced modeled near-infrared (NIR) reflectance, explaining why both modeled and observed EVI, which is especially sensitive to NIR, captures canopy seasonal dynamics well. Our improved analysis of canopy-scale biophysics rules out satellite artifacts as significant causes of satellite-observed seasonal patterns at this site, implying that aggregated phenology explains the larger scale remotely observed patterns. This work significantly reconciles current controversies about satellite-detected Amazon phenology, and improves our use of satellite observations to study climate-phenology relationships in the tropics.

Satellite remote sensing of vegetation at regional to global scales is undertaken at considerable variations in solar zenith angle (SZA) across space and time, yet the extent to which these SZA variations matter for the retrieval of phenology remains largely unknown. Here we examined the effect of seasonal and spatial variations in SZA on retrieving vegetation phenology from time series of the Normalized Difference Vegetation Index (NDVI) and Enhanced Vegetation Index (EVI) across a study area in southeastern Australia encompassing forest, woodland, and grassland sites. The vegetation indices (VI) data span two years and are from the Advanced Himawari Imager (AHI), which is onboard the Japanese Himawari-8 geostationary satellite. The semi-empirical RossThick-LiSparse-Reciprocal (RTLSR) bidirectional reflectance distribution function (BRDF) model was inverted for each spectral band on a daily basis using 10-minute reflectances acquired by H-8 AHI at different sun-view geometries for each site. The inverted RTLSR model was then used to forward calculate surface reflectance at three constant SZAs (20°, 40°, 60°) and one seasonally varying SZA (local solar noon), all normalised to nadir view. Time series of NDVI and EVI adjusted to different SZAs at nadir view were then computed, from which phenological metrics such as start and end of growing season were retrieved. Results showed that NDVI sensitivity to SZA was on average nearly five times greater than EVI sensitivity. VI sensitivity to SZA also varied among sites (biome types) and phenological stages, with NDVI sensitivity being higher during the minimum greenness period than during the peak greenness period. Seasonal SZA variations altered the temporal profiles of both NDVI and EVI, with more pronounced differences in magnitude among NDVI time series normalised to different SZAs. When using VI time series that allowed SZA to vary at local solar noon, the uncertainties in estimating start, peak, end, and length of growing season introduced by local solar noon varying SZA VI time series, were 7.5, 3.7, 6.5, and 11.3 days for NDVI, and 10.4, 11.9, 6.5, and 8.4 days for EVI respectively, compared to VI time series normalised to a constant SZA. Furthermore, the stronger SZA dependency of NDVI compared with EVI, resulted in up to two times higher uncertainty in estimating annual integrated VI, a commonly used remote-sensing proxy for vegetation productivity. Since commonly used satellite products are not generally normalised to a constant sun-angle across space and time, future studies to assess the sun-angle effects on satellite applications in agriculture, ecology, environment, and carbon science are urgently needed. Measurements taken by new-generation geostationary (GEO) satellites offer an important opportunity to refine this assessment at finer temporal scales. In addition, studies are needed to evaluate the suitability of different BRDF models for normalising sun-angle across a broad spectrum of vegetation structure, phenological stages and geographic locations. Only through continuous investigations on how sun-angle variations affect spatiotemporal vegetation dynamics and what is the best strategy to deal with it, can we achieve a more quantitative remote sensing of true signals of vegetation change across the entire globe and through time.

Remote sensing of phenology usually works at the regional and global scales, which imposes considerable variations in the solar zenith angle (SZA) across space and time. Variations in SZA alters the shape and profile of the surface reflectance and vegetation index (VI) time series, but this effect on remote-sensing-derived vegetation phenology has not been adequately evaluated. The objective of this study is to understand the behaviour of VIs response to SZA, and to further improve the interpretation of satellite observed vegetation dynamics, across space and time. In this study, the sensitivity of two widely used VIs—the normalised difference vegetation index (NDVI) and the enhanced vegetation index (EVI)—to SZA was investigated at four northern Australian savanna sites, over a latitudinal distance of 9.8° (~1100 km). Complete time series of surface reflectances, as acquired with different SZA configurations, were simulated using Bidirectional Reflectance Distribution Function (BRDF) parameters provided by MODerate Resolution Imaging Spectroradiometer (MODIS). The sun-angle dependency of the four phenological transition dates were assessed. Results showed that while NDVI was very sensitive to SZA, such sensitivity was nearly absent for EVI. A negative correlation was also observed between NDVI sensitivity to SZA and vegetation cover, with sensitivity declining to the same level as EVI when vegetation cover was high. Different sun-angle configurations resulted in considerable variations in the shape and magnitude of the phenological profiles. The sensitivity of VIs to SZA was generally greater during the dry season (with only active trees present) than in the wet season (with both active trees and grasses), thus, the sun-angle effect on VIs was phenophase-dependent. The sun-angle effect on NDVI time series resulted in considerable differences in the phenological metrics across different sun-angle configurations. Across four sites, the sun-angle effect caused 15.5 days, 21.6 days, and 20.5 days differences in the start, peak, and the end of the growing season derived from NDVI time series, with seasonally varying SZA at local solar noon, as compared to those metrics derived from NDVI time series with fixed SZA. In comparison, those differences in the start, peak, and end of the growing season for EVI were significantly smaller, with only 4.8 days, 4.9 days, and 3 days, respectively. Our results suggest the potential importance of considering the seasonal SZA effect on VI time series prior to the retrieval of phenological metrics.

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.