Comprehensive Analysis of Snow BRDF Variations by Assessing the Improved Kernel-Driven BRDF Model

Understanding the variations in the bidirectional reflectance distribution function (BRDF) and albedo over snow surface under various conditions is important for interpreting the surface–atmosphere processes of the cryosphere, and the kernel-driven model is among the most popular methods to obtain this information for a comprehensive analysis. Recently, the RossThick-LiSparseReciprocal-Snow (RTLSRS) model was developed to better characterize the anisotropic reflectance of snow and shows strong potential for integration into operational remote sensing algorithms for snow BRDF/albedo retrieval. To comprehensively test the ability of the RTLSRS model to reproduce snow reflectance, the fitting accuracy to different multi-angular data derived from ground, tower, aircraft, and satellite platforms across the full optical wavelength range were demonstrated in this study. Special attention in this study was directed to analyzing the model performance under extreme illumination observation geometries, particularly with respect to the retrieval accuracy and stability under large Solar Zenith Angles (SZAs) and different Relative Azimuth Angles (RAAs). The model performance for silt-polluted snow surface with different concentrations is also assessed to provide necessary supplementation, relative to “pure” snow surface in the previous study. The main findings of this study are summarized as follows: (1) The RTLSRS model exhibits strong robustness under various SZAs; even when the SZA exceeds 80°, the model maintains high accuracy in BRDF reconstruction, with root mean square error (RMSE) values below 0.05. (2) The model also demonstrates satisfactory inversion capability when observations deviate from the principal plane (PP); the model can achieve fitting accuracy with R2 approaching 0.5 and RMSE below 0.05 for MODIS data. (3) In the spectral range below 1300 nm, the RTLSRS model effectively reconstructs the scattering characteristics of snow surfaces with light impurity levels (<20 g/0.5 m2). (4) The spectral shape of snow reflectance remains consistent across different view zenith angles (VZAs) in general. However, the variations caused by different SZAs can be as high as 38.49% and such SZA-induced difference can result in WSA estimation discrepancy of up to 63.43%. This comprehensive assessment further affirms and demonstrates the applicability of the RTLSRS model for the first time in fitting observations across different platforms with various optical wavelengths and geometries, and provides an improved understanding to analyze BRDF variations for the user community.