Review Article: Global Monitoring of Snow Water Equivalent using High Frequency Radar Remote Sensing

Seasonal snow cover is the largest single component of the cryosphere in areal extent, covering an average of 46 million square km of Earth's surface (31 % of the land area) each year, and is thus an important expression of and driver of the Earth’s climate. In recent years, Northern Hemisphere spring snow cover has been declining at about the same rate (~ −13 %/decade) as Arctic summer sea ice. More than one-sixth of the world’s population relies on seasonal snowpack and glaciers for a water supply that is likely to decrease this century. Snow is also a critical component of Earth’s cold regions' ecosystems, in which wildlife, vegetation, and snow are strongly interconnected. Snow water equivalent (SWE) describes the quantity of snow stored on the land surface and is of fundamental importance to water, energy, and geochemical cycles. Quality global SWE estimates are lacking. Given the vast seasonal extent combined with the spatially variable nature of snow distribution at regional and local scales, surface observations will not be able to provide sufficient SWE information. Satellite observations presently cannot provide SWE information at the spatial and temporal resolutions required to address science and high socio-economic value applications such as water resource management and streamflow forecasting. In this paper, we review the potential contribution of X- and Ku-Band Synthetic Aperture Radar (SAR) for global monitoring of SWE. We describe radar interactions with snow-covered landscapes, characterization of snowpack properties using radar measurements, and refinement of retrieval algorithms via synergy with other microwave remote sensing approaches. SAR can image the surface during both day and night regardless of cloud cover, allowing high-frequency revisit at high spatial resolution as demonstrated by missions such as Sentinel-1. The physical basis for estimating SWE from X- and Ku-band radar measurements at local scales is volume scattering by millimetre-scale snow grains. Inference of global snow properties from SAR requires an interdisciplinary approach based on field observations of snow microstructure, physical snow modelling, electromagnetic theory, and retrieval strategies over a range of scales. New field measurement capabilities have enabled significant advances in understanding snow microstructure such as grain size, densities, and layering. We describe radar interactions with snow-covered landscapes, the characterization of snowpack properties using radar measurements, and the refinement of retrieval algorithms via synergy with other microwave remote sensing approaches. This review serves to inform the broader snow research, monitoring, and applications communities on progress made in recent decades, and sets the stage for a new era in SWE remote-sensing from SAR measurements.

Experimental and model-based investigation of the links between snow bidirectional reflectance and snow microstructure

Snow stands out from materials at the Earth’s surface owing to its unique optical properties. Snow optical properties are sensitive to the snow microstructure, triggering potent climate feedbacks. The impacts of snow microstructure on its optical properties such as reflectance are, to date, only partially understood. However, precise modelling of snow reflectance, particularly bidirectional reflectance, are required in many problems, e.g. to correctly process satellite data over snow-covered areas. This study presents a dataset that combines bidirectional reflectance measurements over 500–2500 nm and the X-ray tomography of the snow microstructure for three snow samples of two different morphological types. The dataset is used to evaluate the stereological approach from Malinka (2014) that relates snow optical properties to the chord length distribution in the snow microstructure. The mean chord length and specific surface area (SSA) retrieved with this approach from the albedo spectrum and those measured by the X-ray tomography are in excellent agreement. The analysis of the 3D images has shown that the random chords of the ice phase obey the gamma distribution with the shape parameter m taking the value approximately equal to or a little greater than 2. For weak and intermediate absorption (high and medium albedo), the simulated bidirectional reflectances reproduce the measured ones accurately but tend to slightly overestimate the anisotropy of the radiation. For such absorptions the use of the exponential law for the ice chord length distribution instead of the one measured with the X-ray tomography does not affect the simulated reflectance. In contrast, under high absorption (albedo of a few percent), snow microstructure and especially facet orientation at the surface play a significant role in the reflectance, particularly at oblique viewing and incidence.

Tomography-based investigation of concurrent snow creep and isothermal metamorphism

<p>Once on the ground, the microstructure of snow, i.e. the three-dimensional arrangement of ice and pores, quickly evolves with metamorphism and deforms under the overburden of the overlaying snow. Understanding these concurrent processes is important to predict the evolution of the physical and mechanical properties of snow which are crucial for many applications, such as avalanche forecasting. To this end, we monitored oedometric creep tests of snow under isothermal conditions at -8.6°C for about one week with X-ray tomography. We investigated the evolution of recent snow under a constant load of around 4 kPa, where both ice matrix creep and metamorphism are active. Our time-series comprises one of the most highly-resolved images of snow microstructure evolution, with a temporal resolution of 3 h and spatial resolution of 8.5 microns and thousands of images. Interestingly, we observed distinct effects of the overburden and of the vapor transport on the microstructure evolution. In particular, the quantification of the ice bond network through the Euler characteristic and the min-cut surface shows that metamorphism progressively increases the bond size almost independently of the applied overburden, while the application of an overburden yields a rapid increase of the bond coordination number. These distinct impacts exhibit the difficulty to accurately reproduce the time evolution of recent snow by snow cover models, whose snow microstructure representation with density and snow type remains too coarse. </p>

Experimental and model-based investigation of the links between snow bidirectional reflectance and snow microstructure

Snow stands out from materials at the Earth's surface owing to its unique optical properties. Snow optical properties are sensitive to the snow microstructure, triggering potent climate feedbacks. The impacts of snow microstructure on its optical properties such as reflectance are, to date, only partially understood. However, precise modelling of snow reflectance, particularly bidirectional one, are required in many problems, e.g. to process correctly satellite data over snow-covered areas. This study presents a dataset that combines bidirectional reflectance measurements over 500–2500 nm and the X-ray tomography of the snow microstructure for three snow samples of two different morphological types. The dataset is used to evaluate the stereological approach from Malinka (2014) that relates snow optical properties to the chord length distribution in the snow microstructure. The mean chord length and SSA retrieved with this approach from the albedo spectrum and those measured by the X-ray tomography are in excellent agreement. The analysis of the 3D images has shown that the random chords of the ice phase obey the gamma distribution with the shape parameter m taking the value approximately equal or a little greater than 2. For weak and intermediate absorption (high and medium albedo), the simulated bidirectional reflectances reproduce the measured ones accurately but tend to slightly overestimate the anisotropy of the radiation. For such absorptions the use of the exponential law for the ice chord length distribution instead of the one measured with the X-ray tomography does not affect the simulated reflectance. In contrast, under high absorption (albedo of a few percent), snow microstructure and especially facet orientation at the surface, plays a significant role for the reflectance, particularly at oblique viewing and incidence.

Deep ice layer formation in an alpine snowpack: monitoring and modeling

Ice layers may form deep in the snowpack due to preferential water flow, with impacts on the snowpack mechanical, hydrological and thermodynamical properties. This detailed study at a high-altitude alpine site aims to monitor their formation and evolution thanks to the combined use of a comprehensive observation dataset at a daily frequency and state-of-the-art snow-cover modeling with improved ice formation representation. In particular, daily SnowMicroPen penetration resistance profiles enabled us to better identify ice layer temporal and spatial heterogeneity when associated with traditional snowpack profiles and measurements, while upward-looking ground penetrating radar measurements enabled us to detect the water front and better describe the snowpack wetting when associated with lysimeter runoff measurements. A new ice reservoir was implemented in the one-dimensional SNOWPACK model, which enabled us to successfully represent the formation of some ice layers when using Richards equation and preferential flow domain parameterization during winter 2017. The simulation of unobserved melt-freeze crusts was also reduced. These improved results were confirmed over 17 winters. Detailed snowpack simulations with snow microstructure representation associated with a high-resolution comprehensive observation dataset were shown to be relevant for studying and modeling such complex phenomena despite limitations inherent to one-dimensional modeling.