Even though the specific surface area (SSA) and the snow area index (SAI) of snow are crucial variables to determine the chemical and climatic impact of the snow cover, few data are available on the subject. We propose here a novel method to measure snow SSA and SAI. It is based on the measurement of the hemispherical infrared reflectance of snow samples using the DUFISSS instrument (DUal Frequency Integrating Sphere for Snow SSA measurement). DUFISSS uses the 1310 or 1550 nm radiation of laser diodes, an integrating sphere 15 cm in diameter, and InGaAs photodiodes. For SSA60 m² kg⁻¹, snow is usually of low density (typically 30 to 100 kg m⁻³), resulting in insufficient optical depth and 1310 nm radiation reaches the bottom of the sample, causing artifacts. The 1550 nm radiation is therefore used for SSA>60 m² kg⁻¹. Reflectance is then in the range 5 to 12% and the accuracy on SSA is 12%. We propose empirical equations to determine SSA from reflectance at both wavelengths, with that for 1310 nm taking into account the snow density. DUFISSS has been used to measure the SSA of snow and the SAI of snowpacks in polar and Alpine regions.

Snow on the ground is a complex multiphase photochemical reactor that dramatically modifies the chemical composition of the overlying atmosphere. A quantitative description of the emissions of reactive gases by snow requires knowledge of snow physical properties. This overview details our current understanding of how those physical properties relevant to snow photochemistry vary during snow metamorphism. Properties discussed are density, specific surface area, thermal conductivity, permeability, gas diffusivity and optical properties. Inasmuch as possible, equations to parameterize these properties as functions of climatic variables are proposed, based on field measurements, laboratory experiments and theory. The potential of remote sensing methods to obtain information on some snow physical variables such as grain size, liquid water content and snow depth are discussed. The possibilities for and difficulties of building a snow photochemistry model by adapting current snow physics models are explored. Elaborate snow physics models already exist, and including variables of particular interest to snow photochemistry such as light fluxes and specific surface area appears possible. On the other hand, understanding the nature and location of reactive molecules in snow seems to be the greatest difficulty modelers will have to face for lack of experimental data, and progress on this aspect will require the detailed study of natural snow samples.

The energy budget and the photochemistry of a snowpack depend greatly on the penetration of solar radiation in snow. Below the snow surface, spectral irradiance decreases exponentially with depth with a decay constant called the asymptotic flux extinction coefficient. As with the albedo of the snowpack, the asymptotic flux extinction coefficient depends on snow grain shape. While representing snow by a collection of spherical particles has been successful in the numerical computation of albedo, such a description poorly explains the decrease of irradiance in snow with depth. Here we explore the limits of the spherical representation. Under the assumption of geometric optics and weak absorption by snow, the grain shape can be simply described by two parameters: the absorption enhancement parameter B and the geometric asymmetry factor gG. Theoretical calculations show that the albedo depends on the ratio B/(1-gG) and the asymptotic flux extinction coefficient depends on the product B(1-gG). To understand the influence of grain shape, the values of B and gG are calculated for a variety of simple geometric shapes using ray tracing simulations. The results show that B and (1-gG) generally covary so that the asymptotic flux extinction coefficient exhibits larger sensitivity to the grain shape than albedo. In particular it is found that spherical grains propagate light deeper than any other investigated shape. In a second step, we developed a method to estimate B from optical measurements in snow. A multi-layer, two-stream, radiative transfer model, with explicit grain shape dependence, is used to retrieve values of the B parameter of snow by comparing the model to joint measurements of reflectance and irradiance profiles. Such measurements were performed in Antarctica and in the Alps yielding estimates of B between 0.8 and 2.0. In addition, values of B were estimated from various measurements found in the literature, leading to a wider range of values (1.0–9.9) which may be partially explained by the limited accuracy of the data. This work highlights the large variety of snow microstructure and experimentally demonstrates that spherical grains, with B = 1.25, are inappropriate to model irradiance profiles in snow, an important result that should be considered in further studies dedicated to subsurface absorption of short-wave radiation and snow photochemistry.

The Snow Microwave Radiative Transfer (SMRT) thermal emission and backscatter model was developed to determine uncertainties in forward modeling through intercomparison of different model ingredients. The model differs from established models by the high degree of flexibility in switching between different electromagnetic theories, representations of snow microstructure, and other modules involved in various calculation steps. SMRT v1.0 includes the dense media radiative transfer theory (DMRT), the improved Born approximation (IBA), and independent Rayleigh scatterers to compute the intrinsic electromagnetic properties of a snow layer. In the case of IBA, five different formulations of the autocorrelation function to describe the snow microstructure characteristics are available, including the sticky hard sphere model, for which close equivalence between the IBA and DMRT theories has been shown here. Validation is demonstrated against established theories and models. SMRT was used to identify that several former studies conducting simulations with in situ measured snow properties are now comparable and moreover appear to be quantitatively nearly equivalent. This study also proves that a third parameter is needed in addition to density and specific surface area to characterize the microstructure. The paper provides a comprehensive description of the mathematical basis of SMRT and its numerical implementation in Python. Modularity supports model extensions foreseen in future versions comprising other media (e.g., sea ice, frozen lakes), different scattering theories, rough surface models, or new microstructure models.