Surface albedo is an essential variable to determine the Earth's surface energy budget, in particular for snow-covered areas where it is involved in one of the most powerful positive feedback loops of the climate system. In situ measurements of broadband and spectral albedo are therefore common. However they are subject to several artefacts. Here we investigate the sensitivity of spectral albedo measurements to surface slope, and we propose simple correction algorithms to retrieve the intrinsic albedo of a slope from measurements, as if it were flat. For this, we first derive the analytical equations relating albedo measured on a slope to intrinsic direct and diffuse albedo, the apportionment between diffuse and direct incoming radiation, and slope inclination and aspect. The theory accounts for two main slope effects. First, the slope affects the proportion of solar radiation intercepted by the surface relative to that intercepted by the upward-looking, horizontal, sensor. Second, the upward- and downward-looking sensors receive reduced radiation from the sky and the surface respectively and increased radiation from neighbouring terrain. Using this theory, we show that (i) slope has a significant effect on albedo (over 0.01) from as little as a ≈1∘ inclination, causing distortions of the albedo spectral shape; (ii) the first-order slope effect is sufficient to fully explain measured albedo up to ≈15∘, which we designate “small-slope approximation”; and (iii) for larger slopes, the theory depends on the neighbouring slope geometry and land cover, leading to much more complex equations. Next, we derive four correction methods from the small-slope approximation, to be used depending on whether (1) the slope inclination and orientation are known or not, (2) the snow surface is free of impurities or dirty, and (3) a single or a time series of albedo measurements is available. The methods applied to observations taken in the Alps on terrain with up to nearly 20∘ slopes prove the ability to recover intrinsic albedo with a typical accuracy of 0.03 or better. From this study, we derive two main recommendations for future field campaigns: first, sloping terrain requires more attention because it reduces the measurement accuracy of albedo even for almost invisible slopes (1–2∘). Second, while the correction of the slope effect is possible, it requires additional information such as the spectral diffuse and direction partitioning and if possible the actual slope inclination and aspect, especially when the absence of impurities can not be assumed.

The presence of light-absorbing particles (LAPs) in snow leads to a decrease in short-wave albedo affecting the surface energy budget. However, the understanding of the impacts of LAPs is hampered by the lack of dedicated datasets, as well as the scarcity of models able to represent the interactions between LAPs and snow metamorphism. The present study aims to address both these limitations by introducing a survey of LAP concentrations over two snow seasons in the French Alps and an estimation of their impacts based on the Crocus snowpack model that represents the complex interplays between LAP dynamics and snow metamorphism. First, a unique dataset collected at Col du Lautaret (2058 m a.s.l., above sea level, French Alps) for the two snow seasons 2016–2017 and 2017–2018 is presented. This dataset consists of spectral albedo measurements, vertical profiles of snow specific surface area (SSA), density and concentrations of different LAP species. Spectral albedos are processed to estimate SSA and LAP absorption-equivalent concentrations near the surface of the snowpack. These estimates are then compared to chemical measurements of LAP concentrations and SSA measurements. Our dataset highlights, among others, large discrepancies between two measurement techniques of black carbon (BC) concentrations in snow (namely thermal-optical and laser-induced incandescence). Second, we present ensemble snowpack simulations of the multi-physics version of the detailed snowpack model Crocus, forced with in situ meteorological data, as well as dust and BC deposition fluxes from an atmospheric model. The temporal variations of near-surface LAP concentrations and SSA are most of the time correctly simulated. The simulated seasonal radiative forcing of LAPs is 33 % higher for the 2017–2018 snow season than for the 2016–2017 one, highlighting a strong variability between these two seasons. However, the shortening of the snow season caused by LAPs is similar with 10 ± 5 and 11 ± 1 d for the first and the second snow seasons, respectively. This counter-intuitive result is attributed to two small snowfalls at the end of the first season and highlights the importance in accounting for meteorological conditions to correctly predict the impact of LAPs. The strong variability of season shortening caused by LAPs in the multi-physics ensemble for the first season (10 ± 5 d) also points out the sensitivity of model-based estimations of LAP impact on modelling uncertainties of other processes. Finally, the indirect impact of LAPs (i.e. the enhancement of energy absorption due to the acceleration of the metamorphism by LAPs) is negligible for the 2 years considered here, which is contrary to what was found in previous studies for other sites.

Light-absorbing particles (LAPs) such as black carbon or mineral dust are some of the main drivers of snow radiative transfer. Small amounts of LAPs significantly increase snowpack absorption in the visible wavelengths where ice absorption is particularly weak, impacting the surface energy budget of snow-covered areas. However, linking measurements of LAP concentration in snow to their actual radiative impact is a challenging issue which is not fully resolved. In the present paper, we point out a new method based on spectral irradiance profile (SIP) measurements which makes it possible to identify the radiative impact of LAPs on visible light extinction in homogeneous layers of the snowpack. From this impact on light extinction it is possible to infer LAP concentrations present in each layer using radiative transfer theory. This study relies on a unique dataset composed of 26 spectral irradiance profile measurements in the wavelength range 350–950 nm with concomitant profile measurements of snow physical properties and LAP concentrations, collected in the Alps over two snow seasons in winter and spring conditions. For 55 homogeneous snow layers identified in our dataset, the concentrations retrieved from SIP measurements are compared to chemical measurements of LAP concentrations. A good correlation is observed for measured concentrations higher than 5 ng g−1 (r2=0.81) despite a clear positive bias. The potential causes of this bias are discussed, underlining a strong sensitivity of our method to LAP optical properties and to the relationship between snow microstructure and snow optical properties used in the theory. Additional uncertainties such as artefacts in the measurement technique for SIP and chemical contents along with LAP absorption efficiency may explain part of this bias. In addition, spectral information on LAP absorption can be retrieved from SIP measurements. We show that for layers containing a unique absorber, this absorber can be identified in some cases (e.g. mineral dust vs. black carbon). We also observe an enhancement of light absorption between 350 and 650 nm in the presence of liquid water in the snowpack, which is discussed but not fully elucidated. A single SIP acquisition lasts approximately 1 min and is hence much faster than collecting a profile of chemical measurements. With the recent advances in modelling LAP–snow interactions, our method could become an attractive alternative to estimate vertical profiles of LAP concentrations in snow.

Most models simulating snow albedo assume a flat and smooth surface, neglecting surface roughness. However, the presence of macroscopic roughness leads to a systematic decrease in albedo due to two effects: (1) photons are trapped in concavities (multiple reflection effect) and (2) when the sun is low, the roughness sides facing the sun experience an overall decrease in the local incidence angle relative to a smooth surface, promoting higher absorption, whilst the other sides have weak contributions because of the increased incidence angle or because they are shadowed (called the effective-angle effect here). This paper aims to quantify the impact of surface roughness on albedo and to assess the respective role of these two effects, with (1) observations over varying amounts of surface roughness and (2) simulations using the new rough surface ray-tracing (RSRT) model, based on a Monte Carlo method for photon transport calculation. The observations include spectral albedo (400–1050 nm) over manually created roughness surfaces with multiple geometrical characteristics. Measurements highlight that even a low fraction of surface roughness features (7 % of the surface) causes an albedo decrease of 0.02 at 1000 nm when the solar zenith angle (θ_s) is larger than 50^∘. For higher fractions (13 %, 27 % and 63 %), and when the roughness orientation is perpendicular to the sun, the decrease is of 0.03–0.04 at 700 nm and of 0.06–0.10 at 1000 nm. The impact is 20 % lower when roughness orientation is parallel to the sun. The observations are subsequently compared to RSRT simulations. Accounting for surface roughness improves the model observation agreement by a factor of 2 at 700 and 1000 nm (errors of 0.03 and 0.04, respectively) compared to simulations considering a flat smooth surface. The model is used to explore the albedo sensitivity to surface roughness with varying snow properties and illumination conditions. Both multiple reflections and the effective-angle effect have a greater impact with low specific surface area (SSA; <10 m² kg⁻¹). The effective-angle effect also increases rapidly with θ_s at large θ_s. This latter effect is larger when the overall slope of the surface is facing away from the sun and has a roughness orientation perpendicular to the sun. For a snowpack where artificial surface roughness features were created, we showed that a broadband albedo decrease of 0.05 may cause an increase in the net shortwave radiation of 80 % (from 15 to 27 W m⁻²). This paper highlights the necessity of considering surface roughness in the estimation of the surface energy budget and opens the way for considering natural rough surfaces in snow modelling.