Abstract. A majority of snow radiative transfer models (RTMs) treat snow as a collection of idealized grains rather than an organized ice–air matrix. Here we present a generalized multi-layer photon-tracking RTM that simulates light reflectance and transmittance of snow based on X-ray microtomography images, treating snow as a coherent 3D structure rather than a collection of grains. The model uses a blended approach to expand ray-tracing techniques applied to sub-1 cm3 snow samples to snowpacks of arbitrary depths. While this framework has many potential applications, this study's effort is focused on simulating reflectance and transmittance in the visible and near infrared (NIR) through thin snowpacks as this is relevant for surface energy balance and remote sensing applications. We demonstrate that this framework fits well within the context of previous work and capably reproduces many known optical properties of a snow surface, including the dependence of spectral reflectance on the snow specific surface area and incident zenith angle as well as the surface bidirectional reflectance distribution function (BRDF). To evaluate the model, we compare it against reflectance data collected with a spectroradiometer at a field site in east-central Vermont. In this experiment, painted panels were inserted at various depths beneath the snow to emulate thin snow. The model compares remarkably well against the reflectance measured with a spectroradiometer, with an average RMSE of 0.03 in the 400–1600 nm range. Sensitivity simulations using this model indicate that snow transmittance is greatest in the visible wavelengths, limiting light penetration to the top 6 cm of the snowpack for fine-grain snow but increasing to 12 cm for coarse-grain snow. These results suggest that the 5 % transmission depth in snow can vary by over 6 cm according to the snow type.