An Operational Methodology for Validating Satellite-Based Snow Albedo Measurements Using a UAV

Snow albedo is highly variable over multiple temporal and spatial scales. This variability is more pronounced in areas that experience seasonal snowpack. Satellite retrievals, physically based models and parameterizations for snow albedo all require ground-based measurements for calibration, initialization, and validation. Ground measurements are generally made using upward and downward-facing pyranometers at opportunistically located weather stations that are sparsely distributed, particularly in mountainous regions. These station-based measurements cannot capture the spatial variability of albedo across the land surface. Uncrewed Aerial Vehicles (UAVs) equipped with upward and downward-facing pyranometers provide near-surface measurements of broadband albedo that are spatially distributed across landscapes, offering improvements over in-situ sensors. At the hillslope to watershed scale albedo measurements from UAVs taken over heterogeneous terrain are a function of the spatial variability in albedo and topography within the downward-facing sensor’s field-of-view (FOV). In this research we propose methods for topographic correction of UAV snow albedo measurements and comparison to gridded satellite albedo products. These methods account for the variability of surface topography and albedo within the sensor FOV, sensor tilt, and the angular response of pyranometers. We applied the proposed methodologies to UAV snow albedo measurements collected over an alpine meadow in southwest Montana, United States (45.23°, −111.28°). Sensitivity analyses were conducted to determine the effect of altering the processing FOV (PFOV) for both topographic corrections and comparison to coincident Landsat 8-derived albedo measurements. Validation from ground-based albedo measurements showed the topographic correction to reduce albedo measurement error considerably over mildly sloping terrain. Our sensitivity analyses demonstrated that outcomes from the topographic correction and satellite comparison are highly dependent on the specified PFOV. Based on field observations and analyses of UAV albedo measurements made at different altitudes, we provide guidelines for strategizing future UAV albedo surveys. This research presents considerable advances in the standardization of UAV-based albedo measurement. We establish the foundation for future research to utilize this platform to collect near-surface validation measurements over heterogeneous terrain with high accuracy and consistency.

SNICAR-ADv3: a community tool for modeling spectral snow albedo

The Snow, Ice, and Aerosol Radiative (SNICAR) model has been used in various capacities over the last 15 years to model the spectral albedo of snow with light-absorbing constituents (LACs). Recent studies have extended the model to include an adding-doubling two-stream solver and representations of non-spherical ice particles; carbon dioxide snow; snow algae; and new types of mineral dust, volcanic ash, and brown carbon. New options also exist for ice refractive indices and solar-zenith-angle-dependent surface spectral irradiances used to derive broadband albedo. The model spectral range was also extended deeper into the ultraviolet for studies of extraterrestrial and high-altitude cryospheric surfaces. Until now, however, these improvements and capabilities have not been merged into a unified code base. Here, we document the formulation and evaluation of the publicly available SNICAR-ADv3 source code, web-based model, and accompanying library of constituent optical properties. The use of non-spherical ice grains, which scatter less strongly into the forward direction, reduces the simulated albedo perturbations from LACs by ∼9 %–31 %, depending on which of the three available non-spherical shapes are applied. The model compares very well against measurements of snow albedo from seven studies, though key properties affecting snow albedo are not fully constrained with measurements, including ice effective grain size of the top sub-millimeter of the snowpack, mixing state of LACs with respect to ice grains, and site-specific LAC optical properties. The new default ice refractive indices produce extremely high pure snow albedo (>0.99) in the blue and ultraviolet part of the spectrum, with such values only measured in Antarctica so far. More work is needed particularly in the representation of snow algae, including experimental verification of how different pigment expressions and algal cell concentrations affect snow albedo. Representations and measurements of the influence of liquid water on spectral snow albedo are also needed.

Divergence of apparent and intrinsic snow albedo over a season at a sub-alpine site with implications for remote sensing

Intrinsic albedo is the bihemispherical reflectance of a substance with a smooth surface. Conversely, the apparent albedo is the bihemispherical reflectance of the same substance with a rough surface. For snow, the surface is often rough, and these two optical quantities have different uses: intrinsic albedo is used in scattering equations whereas apparent albedo should be used in energy balance models. Complementing numerous studies devoted to surface roughness and its effect on snow reflectance, this work analyzes a timeseries of intrinsic and apparent snow albedos over a season at a sub-alpine site using an automated terrestrial laser scanner to map the snow surface topography. An updated albedo model accounts for shade, and in situ albedo measurements from a field spectrometer are compared to those from a spaceborne multispectral sensor. A spectral unmixing approach using a shade endmember (to address the common problem of unknown surface topography) produces grain size and impurity solutions; the modeled shade fraction is compared to the intrinsic and apparent albedo difference. As expected and consistent with other studies, the results show that intrinsic albedo is consistently greater than apparent albedo. Both albedos decrease rapidly as ablation hollows form during melt, combining effects of impurities on the surface and increasing roughness. Intrinsic broadband albedos average 7 % greater than apparent albedos, with the difference being about 6 % in the near-infrared or 3–4 % if the average (planar) topography is known and corrected. Field measurements of spectral surface reflectance confirm that multispectral sensors see the apparent albedo but lack the spectral resolution to distinguish between darkening from ablation hollows versus low concentrations of impurities. In contrast, measurements from the field spectrometer have sufficient resolution to discern darkening from the two sources. Based on these results, conclusions are: 1) impurity estimates from multispectral sensors are only reliable for relatively dirty snow with high snow fraction; 2) a shade endmember must be used in spectral mixture models, even for in situ spectroscopic measurements; and 3) snow albedo models should produce apparent albedos by accounting for the shade fraction. The conclusion re-iterates that albedo is the most practical snow reflectance quantity for remote sensing.

Review article: Parameterizations of snow-related physical processes in land surface models

Snow on land surface plays a vital role in the interaction between land and atmosphere in the state-of-the-art land surface models (LSMs) and the real world. Since the snow cover affects the snow albedo and the ground and soil heat fluxes, it is crucial to detect snow cover changes accurately. It is challenging to acquire observation data for snow cover, snow albedo, and snow depth; thus, an excellent alternative is to use the simulation data produced by the LSMs that calculate the snow-related physical processes. The LSMs show significant differences in the complexities of the snow parameterizations in terms of variables and processes considered. Thus, the synthetic intercomparisons of the snow physics in the LSMs will help the improvement of each LSM. This study revealed and discussed the differences in the parameterizations among LSMs related to snow cover fraction, snow albedo, and snow density. We selected the most popular and well-documented LSMs embedded in the Earth System Model or operational forecasting systems. We examined single layer schemes, including the Unified Noah Land Surface Model (Noah LSM), the Hydrology Tiled ECMWF Scheme of Surface Exchanges over Land (HTESSEL), the Biosphere-Atmosphere Transfer Scheme (BATS), the Canadian Land Surface Scheme (CLASS), and multilayer schemes of intermediate complexity including the Community Noah Land Surface Model with Multi-Parameterization Options (Noah-MP), the Community Land Model version 5 (CLM 5), the Joint UK Land Environment Simulator (JULES), and the Interaction Soil-Biosphere-Atmosphere (ISBA). First, we identified that BATS, Noah-MP, JULES, and ISBA reflect the snow depth and roughness length to parameterize snow cover fraction, and CLM 5 accounts for the standard deviation of the elevation value for the snow cover decay function. Second, CLM 5 and BATS are relatively complex, so that they explicitly take into account the solar zenith angle, black carbon, mineral dust, organic carbon, and ice grain size for the determinations of snow albedo. Besides, JULES and ISBA are also complicated model which concerns ice grain size, solar zenith angle, new snow depth, fresh snowfall rate, and surface temperature for the albedo scheme. Third, HTESSEL, CLM 5, and ISBA considered the effects of both wind and temperature in the determinations of the new snow density. Especially, ISBA and JULES considered internal snow characteristics such as snow viscosity, snow temperature, and vertical stress for parameterizing new snow density. The future outlook discussed geomorphic and vegetation-related variables for the further improvement of the LSMs. Previous studies clearly show that spatio-temporal variation of snow is due to the influence of altitude, slope, and vegetation condition. Therefore, we recommended applying geomorphic and vegetation factors such as elevation, slope, time-varying roughness length, vegetation indexes, or optimized parameters according to the land surface type to parameterize snow-related physical processes.

A versatile method for computing optimized snow albedo from spectrally fixed radiative variables: VALHALLA v1.0

In climate models, the snow albedo scheme generally calculates only a narrowband or broadband albedo, which leads to significant uncertainties. Here, we present the Versatile ALbedo calculation metHod based on spectrALLy fixed radiative vAriables (VALHALLA version 1.0) to optimize spectral snow albedo calculation. For this optimization, the energy absorbed by the snowpack is calculated by the spectral albedo model Two-streAm Radiative TransfEr in Snow (TARTES) and the spectral irradiance model Santa Barbara DISORT Atmospheric Radiative Transfer (SBDART). This calculation takes into account the spectral characteristics of the incident radiation and the optical properties of the snow based on an analytical approximation of the radiative transfer of snow. For this method, 30 wavelengths, called tie points (TPs), and 16 reference irradiance profiles are calculated to incorporate the absorbed energy and the reference irradiance. The absorbed energy is then interpolated for each wavelength between two TPs with adequate kernel functions derived from radiative transfer theory for snow and the atmosphere. We show that the accuracy of the absorbed energy calculation primarily depends on the adaptation of the irradiance of the reference profile to that of the simulation (absolute difference <1 W m−2 for broadband absorbed energy and absolute difference <0.005 for broadband albedo). In addition to the performance in terms of accuracy and calculation time, the method is adaptable to any atmospheric input (broadband, narrowband) and is easily adaptable for integration into a radiative scheme of a global or regional climate model.

Impact of dust-cloud-radiation interactions on surface albedo: a case study of "Tiramisu" snow in Urumqi, China

Dust-cloud-surface radiation interactions (DCRI) is a complex nonlinear relation referring to the influences of both atmospheric dust and dust-on-snow on surface albedo. A “Tiramisu” snow event occurred on December 1st, 2018, in Urumqi, China, providing an excellent testbed for exploring the comprehensive effect induced by atmospheric dust and those deposited atop fresh snowpack on surface radiation. A detailed analysis indicates that the decrease of snow albedo by 0.17–0.26 (22–34%) is contributed by the effects both the dust-cloud interactions and dust-on-snow at synoptic scale in this case. In particular, dust well mixed with ice clouds at altitudes of 2.5–5.5 km disrupted the “seeder-feeder” structure of clouds and heterogeneous ice nucleation. Dust-induced changes in the low layer of ice cloud (3.3–5.5 km) under a low temperature of –20 °C resulted in a 31.8% increase in the ice particle radius and 84.6% in the ice water path, which acted to indirectly buffer the incident solar radiation reaching the surface. Dust particles deposited on the snow surface further caused snow darkening since the snow albedo was found to decrease by 11.8–23.3%. These findings underscore the importance of considering the comprehensive effect of dust-cloud-radiation interactions in the future.

Retrieval of Snow Albedo and Total Ozone Column from Single-View MSI/S-2 Spectral Reflectance Measurements over Antarctica

We proposed a simple algorithm to retrieve the total ozone column and snow properties (spectral albedo and effective light absorption path) using the high spatial resolution single–view MSI/S-2 measurements over Antarctica. In addition, the algorithm allows the retrieval of the snow grain size on a scale of 10–20 m. This algorithm should be useful for the understanding of intra-pixel total ozone and snow albedo variability in complement to satellite observations performed on a much coarser spatial resolution scale (0.3–1 km and even larger spatial scales).