Microwave scattering coefficient of snow in MEMLS and DMRT-ML revisited: the relevance of sticky hard spheres and tomography-based estimates of stickiness

International audienceThe description of snow microstructure in microwave models is often simplified to facilitate electromagnetic calculations. Within dense media radiative transfer (DMRT), the microstructure is commonly described by sticky hard spheres (SHS). An objective mapping of real snow onto SHS is however missing which prevents measured input parameters from being used for DMRT. In contrast, the microwave emission model of layered snowpacks (MEMLS) employs a conceptually different approach, based on the two-point correlation function which is accessible by tomogra-phy. Here we show the equivalence of both electromagnetic approaches by reformulating their microstructural models in a common framework. Using analytical results for the two-point correlation function of hard spheres, we show that the scattering coefficient in both models only differs by a factor which is close to unity, weakly dependent on ice volume fraction and independent of other microstructural details. Additionally , our analysis provides an objective retrieval method for the SHS parameters (diameter and stickiness) from to-mography images. For a comprehensive data set we demonstrate the variability of stickiness and compare the SHS diameter to the optical equivalent diameter. Our results confirm the necessity of a large grain-size scaling when relating both diameters in the non-sticky case, as previously suggested by several authors

Measurement of the specific surface area of snow using infrared reflectance in an integrating sphere at 1310 and 1550 nm

International audienceEven 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 m2 kg−1, snow is usually of low density (typically 30 to 100 kg m−3), 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 m2 kg−1. 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

Simulation of Seasonal Snow Microwave TB Using Coupled Multi-Layered Snow Evolution and Microwave Emission Models

The accurate quantification of SWE has important societal benefits, including improving domestic and agricultural water planning, flood forecasting and electric power generation. However, passive-microwave SWE algorithms suffer from variations in TB due to snow metamorphism, difficult to distinguish from those due to SWE variations. Coupled snow evolution-emission models are able to predict snow metamorphism, allowing us to account for emissivity changes. They can also be used to identify weaknesses in the snow evolution model. Moreover, thoroughly evaluating coupled models is a contribution toward the assimilation of TB, which leads to a significant increase in the accuracy of SWE estimates

Surface melt on the Shackleton Ice Shelf, East Antarctica (2003–2021)

Melt on the surface of Antarctic ice shelves can potentially lead to their disintegration, accelerating the flow of grounded ice to the ocean and raising global sea levels. However, the current understanding of the processes driving surface melt is incomplete, increasing uncertainty in predictions of ice shelf stability and thus of Antarctica's contribution to sea-level rise. Previous studies of surface melt in Antarctica have usually focused on either a process-level understanding of melt through energy-balance investigations or used metrics such as the annual number of melt days to quantify spatiotemporal variability in satellite observations of surface melt. Here, we help bridge the gap between work at these two scales. Using daily passive microwave observations from the AMSR-E and AMSR-2 sensors and the machine learning approach of a self-organising map, we identify nine representative spatial distributions (“patterns”) of surface melt on the Shackleton Ice Shelf in East Antarctica from 2002/03–2020/21. Combined with output from the RACMO2.3p3 regional climate model and surface topography from the REMA digital elevation model, our results point to a significant role for surface air temperatures in controlling the interannual variability in summer melt and also reveal the influence of localised controls on melt. In particular, prolonged melt along the grounding line shows the importance of katabatic winds and surface albedo. Our approach highlights the necessity of understanding both local and large-scale controls on surface melt and demonstrates that self-organising maps can be used to investigate the variability in surface melt on Antarctic ice shelves.</p

The sensitivity of satellite microwave observations to liquid water in the Antarctic snowpack

Surface melting on the Antarctic Ice Sheet has been monitored by satellite microwave radiometry for over 40 years. Despite this long perspective, our understanding of the microwave emission from wet snow is still limited, preventing the full exploitation of these observations to study supraglacial hydrology. Using the Snow Microwave Radiative Transfer (SMRT) model, this study investigates the sensitivity of microwave brightness temperature to snow liquid water content at frequencies from 1.4 to 37 GHz. We first determine the snowpack properties for eight selected coastal sites by retrieving profiles of density, grain size and ice layers from microwave observations when the snowpack is dry during wintertime. Second, a series of brightness temperature simulations is run with added water. The results show that (i) a small quantity of liquid water (≈0.5 kg m−2) can be detected, but the actual quantity cannot be retrieved out of the full range of possible water quantities; (ii) the detection of a buried wet layer is possible up to a maximum depth of 1 to 6 m depending on the frequency (6–37 GHz) and on the snow properties (grain size, density) at each site; (iii) surface ponds and water-saturated areas may prevent melt detection, but the current coverage of these waterbodies in the large satellite field of view is presently too small in Antarctica to have noticeable effects; and (iv) at 1.4 GHz, while the simulations are less reliable, we found a weaker sensitivity to liquid water and the maximal depth of detection is relatively shallow (&lt;10 m) compared to the typical radiation penetration depth in dry firn (≈1000 m) at this low frequency. These numerical results pave the way for the development of improved multi-frequency algorithms to detect melt intensity and the depth of liquid water below the surface in the Antarctic snowpack.</p

On the assimilation of optical reflectances and snow depth observations into a detailed snowpack model

International audienceThis paper examines the ability of optical re-flectance data assimilation to improve snow depth and snow water equivalent simulations from a chain of models with the SAFRAN meteorological model driving the detailed multi-layer snowpack model Crocus now including a two-stream radiative transfer model for snow, TARTES. The direct use of reflectance data, allowed by TARTES, instead of higher level snow products, mitigates uncertainties due to commonly used retrieval algorithms. Data assimilation is performed with an ensemble-based method, the Sequential Importance Resampling Particle filter , to represent simulation uncertainties. In snowpack mod-eling, uncertainties of simulations are primarily assigned to meteorological forcings. Here, a method of stochastic perturbation based on an autoregressive model is implemented to explicitly simulate the consequences of these uncertainties on the snowpack estimates. Through twin experiments, the assimilation of synthetic spectral reflectances matching the MODerate resolution Imaging Spectroradiometer (MODIS) spectral bands is examined over five seasons at the Col du Lautaret, located in the French Alps. Overall, the assimilation of MODIS-like data reduces by 45 % the root mean square errors (RMSE) on snow depth and snow water equivalent. At this study site, the lack of MODIS data on cloudy days does not affect the assimilation performance significantly. The combined assimilation of MODIS-like reflectances and a few snow depth measurements throughout the 2010/2011 season further reduces RMSEs by roughly 70 %. This work suggests that the assimilation of optical reflectances has the potential to become an essential component of spatialized snowpack simulation and forecast systems. The assimilation of real MODIS data will be investigated in future works

Aquarius Brightness Temperature Variations at Dome C and Snow Metamorphism at the Surface

The Antarctic Plateau is a promising site to monitor microwave radiometers' drift, and to inter-calibrate microwave radiometers, especially 1.4 GHz (L-band) radiometers on board the Soil Moisture and Ocean Salinity (SMOS), and AquariusSAC-D missions. The Plateau is a thick ice cover, thermally stable in depth, with large dimensions, and relatively low heterogeneities. In addition, its high latitude location in the Southern Hemisphere enables frequent observations by polar-orbiting satellites, and no contaminations by radio frequency interference. At Dome C (75S, 123E), on the Antarctic Plateau, the substantial amount of in-situ snow measurements available allows us to interpret variations in space-borne microwave brightness temperature (TB) (e.g. Macelloni et al., 2007, 2013, Brucker et al., 2011, Champollion et al., 2013). However, to analyze the observations from the Aquarius radiometers, whose sensitivity is 0.15 K, the stability of the snow layers near the surface that are most susceptible to rapidly change needs to be precisely assessed. This study focuses on the spatial and temporal variations of the Aquarius TB over the Antarctic Plateau, and at Dome C in particular, to highlight the impact of snow surface metamorphism on the TB observations at L-band