Simulating snow-covered forest bidirectional reflectance by extending hybrid geometric optical–radiative transfer model

Accurate simulations of the bidirectional reflectance distribution function (BRDF) of snow-covered forests are crucial for retrieving snow and canopy properties, which affects the analysis of net radiation, radiative forcing, and climate change. The anisotropic reflectance of snow in conifer forests is greatly influenced by the high reflectance of snow and the complex travel pathways of radiation in 3D canopy structures, which complicates BRDF simulations. Although bidirectional reflection models have been widely established for snow-free forests, few studies have attempted to simulate the BRDF in snow-covered forests. Therefore, we propose a new snow-covered forest bidirectional reflectance (SFBR) model by coupling geometric optical and radiative transfer models, which includes a snow anisotropic reflectance model (ART), needle leaf optical properties model (LIBERTY), canopy radiative transfer model (4SAIL), and geometric optical model (GOST2). The SFBR model mainly focuses on conifer forests with snow backgrounds; it can simulate anisotropic reflectance spectra from 400 nm to 2500 nm and is the only model to fully consider discontinuous canopy distribution, topographic effects, and snow anisotropy. The SFBR model was validated by simulating the reflectance over snow-covered forests, the results of which showed good consistency with simulations of 3D scene model and remote sensing observations. The proposed SFBR model has the potential to quantify radiation interactions between the snow background and forest overstory, evaluate the sensitivity of canopy reflectance to snow and forest properties, and retrieve snow-covered forest parameters from satellite-observed reflectance data. •A new bidirectional reflectance model for snow-covered forests.•It couples GO and RT models for snow, leaf, and canopy.•Canopy heterogeneity, topographic effects, and snow anisotropy are considered.•Simulated reflectance has good consistency with 3D scene model and MODIS data.•Model flexibility allows application for spectra simulation and mechanism analysis.