Impact of BRDF Spatiotemporal Smoothing on Land Surface Albedo Estimation
Authors: Yang Jian , Shuai Yanmin , Duan Junbo , Xie Donghui , Zhang Qingling , Zhao Ruishan 
Published: April 2022

Surface albedo, as a key parameter determining the partition of solar radiation at the Earth’s surface, has been developed into a satellite-based product from various Earth observation systems to serve numerous global or regional applications. Studies point out that apparent uncertainty can be introduced into albedo retrieval without consideration of surface anisotropy, which is a challenge to albedo estimation especially from observations with fewer angular samplings. Researchers have begun to introduce smoothed anisotropy prior knowledge into albedo estimation to improve the inversion efficiency, or for the scenario of observations with signal or poor angular sampling. Thus, it is necessary to further understand the potential influence of smoothed anisotropy features adopted in albedo estimation. We investigated the albedo variation induced by BRDF smoothing at both temporal and spatial scales over six typical landscapes in North America using MODIS standard anisotropy products with high quality BRDF inversed from multi-angle observations in 500 m and 5.6 km spatial resolutions. Components of selected typical landscapes were assessed with the confidence of the MCD12 land cover product and 30 m CDL (cropland data layer) classification maps followed by an evaluation of spatial heterogeneity in 30 m scale through the semi-variogram model. High quality BRDF of MODIS standard anisotropy products were smoothed in multi-temporal scales of 8 days, 16 days, and 32 days, and in multi-spatial scales from 500 m to 5.6 km. The induced relative and absolute albedo differences were estimated using the RossThick-LiSparseR model and BRDFs smoothed before and after spatiotemporal smoothing. Our results show that albedo estimated using BRDFs smoothed temporally from daily to monthly over each scenario exhibits relative differences of 11.3%, 12.5%, and 27.2% and detectable absolute differences of 0.025, 0.012, and 0.013, respectively, in MODIS near-infrared (0.7–5.0 µm), short-wave (0.3–5.0 µm), and visible (0.3–0.7 µm) broad bands. When BRDFs of investigated landscapes are smoothed from 500 m to 5.6 km, variations of estimated albedo can achieve up to 36.5%, 37.1%, and 94.7% on relative difference and absolute difference of 0.037, 0.024, and 0.018, respectively, in near-infrared (0.7–5.0 µm), short wave (0.3–5.0 µm), and visible (0.3–0.7 µm) broad bands. In addition, albedo differences caused by temporal smoothing show apparent seasonal characteristic that the differences are significantly higher in spring and summer than those in autumn and winter, while albedo differences induced by spatial smoothing exhibit a noticeable relationship with sill values of a fitted semi-variogram marked by a correlation coefficient of 0.8876. Both relative and absolute albedo differences induced by BRDF smoothing are demonstrated to be captured, thus, it is necessary to avoid the smoothing process in quantitative remote sensing communities, especially when immediate anisotropy retrievals are available at the required spatiotemporal scale.