Microwave Emission Model Of Layered Snowpacks Research Articles

A novel global snow depth retrieval method considering snow metamorphism and forest influence

Accurate snow depth (SD) information at global scales is highly important for validation and calibration of hydrological model, snowmelt runoff prediction, and global water supplies estimation. Passive microwave remote sensing has been proven to be an efficient method to detect SD at regional or global scales. However, numerous heterogeneous snow property parameters and forest parameters affect the microwave radiation characteristics, which is a great challenge for SD retrieval. In this study, a novel SD retrieval algorithm that comprehensively considered the multiple sensitive factors was developed at global scales by coupling the Microwave Emission Model of Layered Snowpacks (MEMLS) model and machine learning (ML) technique. The exponential correlation length (Pec), snow density (ρ), and forest cover fraction (fforest) are three sensitive parameters for retrieving SD from passive microwave observations based on the extended Fourier amplitude sensitivity test (EFAST) algorithm. The key parameter Pec, which characterizes the snow microstructure information inside snowpacks, was optimized by minimizing the brightness temperature between satellite observations and microwave radiative transfer model simulations by integrating the MEMLS, forest impact, and atmospheric radiation attenuation. The random forest (RF) ML technique was used, and the optimal combination of independent variables was selected by correlation heatmap and importance ranking of variables. Six selected essential parameters, including Pec, vertical polarized brightness temperature differences between 18.7 and 36.5 GHz, ρ from ERA5, fforest, elevation, and longitude together with target SD output were used to train the RF model. Seventy percent of the entire data at meteorological stations were randomly chosen to train the RF model, and the remaining 30% were adopted as the validation data from 2013 to 2020 to ensure the mutual independence of training data and validation data. Compared with independent in situ measurements, the overall RMSE of SD is 13.6 cm, which was noticeably lower than those of the international SD remote sensing products and reanalysis datasets, including AMSR2 (36.2 cm), GlobSnow 3.0 (25.3 cm), ERA5 (38.2 cm), and ERA5-land (29.9 cm). The results demonstrated that the new proposed algorithm, by coupling the MEMLS model with the ML technique, could significantly enhance the overall performance in SD retrieval at global scales.

Global Sensitivity Analysis of the MEMLS Model for Retrieving Snow Water Equivalent

Sensitivity analysis (SA) of model parameters is of great importance for understanding, development, and application of models. However, the influence of snow microstructure variability on snow water equivalent retrieval from passive microwave measurements is still unclear. This article explores the parameter sensitivity of the microwave emission model of layered snowpacks (MEMLS) with improved born approximation (IBA) by using a quantitative global SA method, the extended Fourier amplitude sensitivity test (EFAST) algorithm. A deep analysis is conducted, including the sensitivity of passive microwave emission to snow parameters, the sensitivity variation analysis for different snow conditions, and the temporal properties of the parameter sensitivity. The results show the exponential correlation length, snow depth, and snow density are the three most sensitive parameters for snow without salt in the MEMLS model for the brightness temperature gradient at 18.7 and 36.5 GHz. For snow with a small salt content, the exponential correlation length, snow depth, snow temperature, and snow density are the four most sensitive parameters. Second, snow parameter variability highly affects the microwave radiation. The sensitivity values of microwave brightness temperature to snow depth gradually increase when the exponential correlation length is less than 0.25 mm and then slightly decreases with the increase of exponential correlation length and decreases along with the increase of snow density. Finally, our analysis highlights the importance to include the snow density, especially for deep snow depth, in the combination of sensitive factors in future multiparameter retrievals.

Retrieval of Snow Depth on Arctic Sea Ice from the FY3B/MWRI

Given their high albedo and low thermal conductivity, snow and sea ice are considered key reasons for amplified warming in the Arctic. Snow-covered sea ice is a more effective insulator, which greatly limits the energy and momentum exchange between the atmosphere and surface, and further controls the thermal dynamic processes of snow and ice. In this study, using the Microwave Emission Model of Layered Snowpacks (MEMLS), the sensitivities of the brightness temperatures (TBs) from the FengYun-3B/MicroWave Radiometer Imager (FY3B/MWRI) to changes in snow depth were simulated, on both first-year and multiyear ice in the Arctic. Further, the correlation coefficients between the TBs and snow depths in different atmospheric and sea ice environments were investigated. Based on the simulation results, the most sensitive factors to snow depth, including channels of MWRI and their combination form, were determined for snow depth retrieval. Finally, using the 2012–2013 Operational IceBridge (OIB) snow depth data, retrieval algorithms of snow depth were developed for the Arctic on first-year and multiyear ice, separately. Validation using the 2011 OIB data indicates that the bias and standard deviation (Std) of the algorithm are 2.89 cm and 2.6 cm on first-year ice (FYI), respectively, and 1.44 cm and 4.53 cm on multiyear ice (MYI), respectively.

Improving the Snow Volume Scattering Algorithm in a Microwave Forward Model by Using Ground-Based Remote Sensing Snow Observations

Volume scattering (VS) estimation plays a critical role in microwave emission modeling of the snowpack. However, it is challenging to obtain VS accurately for different frequencies by using the microwave emission model of layered snowpacks (MEMLS), which is one of the representative microwave emission models. This article develops a new VS method to consider frequency and exponential correlation length based on a snowfield campaign from November 2015 to April 2016 in Altay, China. Compared with the commonly used empirical and improved Born approximation (IBA) algorithms, the proposed VS algorithm exhibits better performances at both 18 and 36 GHz with a wide range of snow grain sizes. The bias of brightness temperatures at vertical polarization from the proposed algorithm against the observed brightness temperatures are 1.1 K and −0.4 K at 18 and 36 GHz, respectively; the root mean square errors (RMSEs) are 1.8 K and 2.6 K, respectively. The RMSEs decreased by 16.2 K at 18 GHz and 6.5 K at 36 GHz compared with those from the empirical methods and by 2.1 K and 22.2 K compared with those from the IBA. This work demonstrates that the VS difference between 18 and 36 GHz is larger and the dependence of VS on grain size is weaker than those represented by existing methods.

Evaluation of stereology for snow microstructure measurement and microwave emission modeling: a case study

ABSTRACT Reliable microstructure measurement of snow is a requirement for microwave radiative transfer model validation. Snow specific surface area (SSA) can be measured using stereological methods, in which snow samples are cast in the field and photographed in the laboratory. Processing stereology photographs manually by counting intersections of test cycloids with air–ice boundaries reduces the problems in binary segmentation. This paper is a case study to evaluate the repeatability of the manually stereology interpretation by two independent research groups. We further assessed how uncertainty in snow SSA influences simulated brightness temperature (TB) driven by the Microwave Emission Model of Layered Snowpacks (MEMLS), and how stereology compares to Near Infrared (NIR) camera and hand lens. Data was obtained from two alpine snow profiles from Steamboat Springs, Colorado. Results showed that stereological SSA values measured by two groups are highly consistent, and the ground radiometer measured TB at 19 and 37 GHz was successfully predicted (RMSE<3.8 K); simulations using NIR SSA and hand-lens geometric grain size (Dg ) measurements have larger errors. This conclusion was not sensitive to uncertainty in the free parameters of TB modeling.

Validation of the SNTHERM Model Applied for Snow Depth, Grain Size, and Brightness Temperature Simulation at Meteorological Stations in China

Validation of the snow process model is an important preliminary work for the snow parameter estimation. The snow grain growth is a continuous and accumulative process, which cannot be evaluated without comparing with the observations in snow season scale. In order to understand the snow properties in the Asian Water Tower region (including Xinjiang province and the Tibetan Plateau) and enhance the use of modeling tools, an extended snow experiment at the foot of the Altay Mountain was designed to validate and improve the coupled physical Snow Thermal Model (SNTHERM) and the Microwave Emission Model of Layered Snowpacks (MEMLS). By matching simultaneously the observed snow depth, geometric grain size, and observed brightness temperature (TB), with an RMSE of 1.91 cm, 0.47 mm, and 4.43 K (at 36.5 GHz, vertical polarization), respectively, we finalized the important model coefficients, which are the grain growth coefficient and the grain size to exponential correlation length conversion coefficients. When extended to 102 meteorological stations in the 2008–2009 winter, the SNTHERM predicted the daily snow depth with an accuracy of 2–4 cm RMSE, and the coupled SNTHERM-MEMLS model predicted the satellite-observed TB with an accuracy of 13.34 K RMSE at 36.5 GHz, vertical polarization, with the fractional snow cover considered.

Snow Depth Retrieval in Farmland Based on a Statistical Lookup Table from Passive Microwave Data in Northeast China

At present, passive microwave remote sensing is the most efficient method to estimate snow depth (SD) at global and regional scales. Farmland covers 46% of Northeast China and accurate SD retrieval throughout the whole snow season has great significance for the agriculture management field. Based on the results of the statistical analysis of snow properties in Northeast China from December 2017 to January 2018, conducted by the China snow investigation project, snow characteristics such as snow grain size (SGS), snow density, snow thickness, and temperature of the layered snowpack were measured and analyzed in detail. These characteristics were input to the microwave emission model of layered snowpacks (MEMLS) to simulate the brightness temperature (TB) time series of snow-covered farmland in the periods of snow accumulation, stabilization, and ablation. Considering the larger SGS of the thick depth hoar layer that resulted in a rapid decrease of simulated TBs, effective SGS was proposed to minimize the simulation errors and ensure that the MEMLS can be correctly applied to satellite data simulation. Statistical lookup tables (LUTs) for MWRI and AMSR2 data were generated to represent the relationship between SD and the brightness temperature difference (TBD) at 18 and 36 GHz. The SD retrieval results based on the LUT were compared with the actual SD and the SD retrieved by Chang’s algorithm, Foster’s algorithm, the standard MWRI algorithm, and the standard AMSR2 algorithm. The results demonstrated that the proposed algorithm based on the statistical LUT achieved better accuracy than the other algorithms due to its incorporation of the variation in snow characteristics with the age of snow cover. The average root mean squared error of the SD for the whole snow season was approximately 3.97 and 4.22 cm for MWRI and AMSR2, respectively. The research results are beneficial for monitoring SD in the farmland of Northeast China.

Evaluation of Brightness Temperature Sensitivity to Snowpack Physical Properties Using Coupled Snow Physics and Microwave Radiative Transfer Models

There are multiple existing microwave radiative transfer models (RTMs) to simulate the brightness temperature (Tb) of snowpacks. It is still challenging to have consistent Tb responses from RTMs due to individual physical formulations of the snowpack scattering process. This article examines three of the widely-used multi-layer RTMs: 1) the microwave emission model of layered snowpacks (MEMLS); 2) the dense media radiative transfer based on the quasi-crystalline approximation (QCA) of Mie scattering of densely packed sticky spheres (DMRT-QMS); and 3) the Helsinki University of Technology (HUT) model. Interestingly, these models yield slightly different Tb responses when driven by the same physical snowpack properties. Tb variations, dependent on the choice of RTMs, are then evaluated to improve the understanding of model differences in microwave emission from a snowpack. We first perform a sensitivity study of the Tb predictions from the three RTMs as a function of snow grain sizes, densities, and depths. While Tb from all three RTMs decreases with increasing snow grain sizes, it is found that a scaling factor is required to have the same amount of Tb attenuation for small grain sizes within the Rayleigh scattering regime. For larger grain sizes, however, a scaling coefficient is not enough to match the model outputs due to the different scattering assumptions of the RTMs. For a single snow layer with increasing snow depths and densities, all three RTMs exhibit Tb attenuations arising from the increase in path lengths and optical depths. Further evaluations are conducted by feeding the three RTMs with the output of a snow physics model driven by in situ weather forcing in a coupled simulation. Outputs of this coupled model include snowpack physical properties and Tbs. By using snow stratigraphy observations, another set of Tb simulation is also conducted with RTMs driven by in situ snowpit observations. The snow physics outputs from the coupled case are compared against in situ snow stratigraphy observations. And both Tb simulations are compared against ground-based microwave observations from the European Space Agency (ESA) Nordic Snow Radar Experiment (NoSREx) 2009–2012. For three consecutive years, the in situ driven Tbs have 21.0-K root-mean-squared error (RMSE) while the coupled simulations have 24.7-K RMSE. However, after isolating the dry snow period and excluding diurnal melting snow conditions, 12.2- and 6.3-K RMSEs are achieved from in situ and coupled cases, respectively, in the 2011 water year.

Estimation of Snow Depth over the Qinghai-Tibetan Plateau Based on AMSR-E and MODIS Data

Snow cover over the Qinghai-Tibetan Plateau (QTP) plays an important role in climate, hydrological, and ecological systems. Currently, passive microwave remote sensing is the most efficient way to monitor snow depth on global and regional scales; however, it presents a serious overestimation of snow cover over the QTP and has difficulty describing patchy snow cover over the QTP because of its coarse spatial resolution. In this study, a new spatial dynamic method is developed by introducing ground emissivity and assimilating the snow cover fraction (SCF) and land surface temperature (LST) of the Moderate Resolution Imaging Spectroradiometer (MODIS) to derive snow depth at an enhanced spatial resolution. In this method, the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) brightness temperature and MODIS LST are used to calculate ground emissivity. Additionally, the microwave emission model of layered snowpacks (MEMLS) is applied to simulate brightness temperature with varying ground emissivities to determine the key coefficients in the snow depth retrieval algorithm. The results show that the frozen ground emissivity presents large spatial heterogeneity over the QTP, which leads to the variation of coefficients in the snow depth retrieval algorithm. The overestimation of snow depth is rectified by introducing the ground emissivity factor at 18 and 36 GHz. Compared with in situ observations, the snow cover accuracy of the new method is 93.9%, which is better than the 60.2% accuracy of the existing method (old method) which does not consider ground emissivity. The bias and root-mean-square error (RMSE) of snow depth are 1.03 cm and 7.05 cm, respectively, for the new method; these values are much lower than the values of 6.02 cm and 9.75 cm, respectively, for the old method. However, the snow cover accuracy with depths between 1 and 3 cm is below 60%, and snow depths greater than 25 cm are underestimated in Himalayan mountainous areas. In the future, the snow cover identification algorithm should be improved to identify shallow snow cover over the QTP, and topography should be considered in the snow depth retrieval algorithm to improve snow depth accuracy in mountainous areas.

Modelling the L-Band Snow-Covered Surface Emission in a Winter Canadian Prairie Environment

Detailed angular ground-based L-band brightness temperature (TB) measurements over snow covered frozen soil in a prairie environment were used to parameterize and evaluate an electromagnetic model, the Wave Approach for LOw-frequency MIcrowave emission in Snow (WALOMIS), for seasonal snow. WALOMIS, initially developed for Antarctic applications, was extended with a soil interface model. A Gaussian noise on snow layer thickness was implemented to account for natural variability and thus improve the TB simulations compared to observations. The model performance was compared with two radiative transfer models, the Dense Media Radiative Transfer-Multi Layer incoherent model (DMRT-ML) and a version of the Microwave Emission Model for Layered Snowpacks (MEMLS) adapted specifically for use at L-band in the original one-layer configuration (LS-MEMLS-1L). Angular radiometer measurements (30°, 40°, 50°, and 60°) were acquired at six snow pits. The root-mean-square error (RMSE) between simulated and measured TB at vertical and horizontal polarizations were similar for the three models, with overall RMSE between 7.2 and 10.5 K. However, WALOMIS and DMRT-ML were able to better reproduce the observed TB at higher incidence angles (50° and 60°) and at horizontal polarization. The similar results obtained between WALOMIS and DMRT-ML suggests that the interference phenomena are weak in the case of shallow seasonal snow despite the presence of visible layers with thicknesses smaller than the wavelength, and the radiative transfer model can thus be used to compute L-band brightness temperature.

Improving the Radiance Assimilation Performance in Estimating Snow Water Storage across Snow and Land-Cover Types in North America

Abstract Continental-scale snow radiance assimilation (RA) experiments are conducted in order to improve snow estimates across snow and land-cover types in North America. In the experiments, the ensemble adjustment Kalman filter is applied and the Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) brightness temperature TB observations are assimilated into an RA system composed of the Community Land Model, version 4 (CLM4); radiative transfer models (RTMs); and the Data Assimilation Research Testbed (DART). The performance of two snowpack RTMs, the Dense Media Radiative Transfer–Multi-Layers model (DMRT-ML), and the Microwave Emission Model of Layered Snowpacks (MEMLS) in improving snow depth estimates through RA is compared. Continental-scale snow estimates are enhanced through RA by using AMSR-E TB at the 18.7- and 23.8-GHz channels [3% (DMRT-ML) and 2% (MEMLS) improvements compared to the cases using the 18.7- and 36.5-GHz channels] and by considering the vegetation single-scattering albedo ω [2.5% (DMRT-ML) and 4.8% (MEMLS) improvements compared to the cases neglecting ω]. The contribution of TB of the vegetation canopy to TB at the top of the atmosphere is better represented by considering ω in the RA system, and improvements in the resulting snow depth are evident for the forest land-cover type (about 5%–11%) and the taiga and alpine snow classes (about 5%–11% and 4%–8%, respectively), especially in the MEMLS case. Compared to the open-loop run (0.171-m snow depth RMSE), about 7% (DMRT-ML) and 10% (MEMLS) overall improvements of the RA performance are achieved.

Comparison of commonly-used microwave radiative transfer models for snow remote sensing

This paper reviews four commonly-used microwave radiative transfer models that take different electromagnetic approaches to simulate snow brightness temperature (TB): the Dense Media Radiative Transfer - Multi-Layer model (DMRT-ML), the Dense Media Radiative Transfer - Quasi-Crystalline Approximation Mie scattering of Sticky spheres (DMRT-QMS), the Helsinki University of Technology n-Layers model (HUT-nlayers) and the Microwave Emission Model of Layered Snowpacks (MEMLS). Using the same extensively measured physical snowpack properties, we compared the simulated TB at 11, 19 and 37GHz from these four models. The analysis focuses on the impact of using different types of measured snow microstructure metrics in the simulations. In addition to density, snow microstructure is defined for each snow layer by grain optical diameter (Do) and stickiness for DMRT-ML and DMRT-QMS, mean grain geometrical maximum extent (Dmax) for HUT n-layers and the exponential correlation length for MEMLS. These metrics were derived from either in-situ measurements of snow specific surface area (SSA) or macrophotos of grain sizes (Dmax), assuming non-sticky spheres for the DMRT models. Simulated TB sensitivity analysis using the same inputs shows relatively consistent TB behavior as a function of Do and density variations for the vertical polarization (maximum deviation of 18K and 27K, respectively), while some divergences appear in simulated variations for the polarization ratio (PR). Comparisons with ground-based radiometric measurements show that the simulations based on snow SSA measurements have to be scaled with a model-specific factor of Do in order to minimize the root mean square error (RMSE) between measured and simulated TB. Results using in-situ grain size measurements (SSA or Dmax, depending on the model) give a mean TB RMSE (19 and 37GHz) of the order of 16–26K, which is similar for all models when the snow microstructure metrics are scaled. However, the MEMLS model converges to better results when driven by the correlation length estimated from in-situ SSA measurements rather than Dmax measurements. On a practical level, this paper shows that the SSA parameter, a snow property that is easy to retrieve in-situ, appears to be the most relevant parameter for characterizing snow microstructure, despite the need for a scaling factor.

Estimation of snow depth from passive microwave brightness temperature data in forest regions of northeast China

Abstract Snow depth is an important factor in water resources management in Northeast China. Forest covers 40% of Northeast China, and the presence of forests influences the accuracy of snow depth retrievals from passive microwave remote sensing data. An optimal iteration method was used to retrieve the forest transmissivities at 18 and 36 GHz based on the snow and forest microwave radiative transfer models and the snow properties measured in field experiments. The transmissivities at 18 and 36 GHz are 0.895 and 0.656 in the horizontal polarization, and 0.821 and 0.615 in the vertical polarization, respectively. Furthermore, the forest transmissivity and snow properties were input into the Microwave Emission Model of Layered Snowpacks (MEMLS) to establish a dynamic look-up table (LUT). Snow depths were retrieved from satellite passive microwave remote sensing data based on the LUT method, and these retrievals were verified by snow depth observations at 103 meteorological stations. The results showed that the bias between the retrieved and measured snow depths is very small, with root mean square errors (RMSEs) of approximately 6 cm in forest regions and 4 cm in non-forest regions. When compared with the existing snow products, the snow depth retrieved in this work presented the highest level of accuracy. The regional snow depth product in China is superior to the GlobSnow and NASA AMSR-E standard SWE products in non-forest regions, whereas the GlobSnow estimate is superior to the regional snow depth product in China and NASA AMSR-E standard SWE product estimates in forest regions. Therefore, we conclude that 1) the influence of forest on snow depth retrieval is important, and the appropriate forest parameters should be considered in the estimation of snow depth from passive microwave brightness temperature data; and 2) the snow depth retrieval algorithm based on the dynamic LUT method proved to be efficient in Northeast China.

Arctic Snow Microstructure Experiment for the development of snow emission modelling

Abstract. The Arctic Snow Microstructure Experiment (ASMEx) took place in Sodankylä, Finland in the winters of 2013–2014 and 2014–2015. Radiometric, macro-, and microstructure measurements were made under different experimental conditions of homogenous snow slabs, extracted from the natural seasonal taiga snowpack. Traditional and modern measurement techniques were used for snow macro- and microstructure observations. Radiometric measurements of the microwave emission of snow on reflector and absorber bases were made at frequencies 18.7, 21.0, 36.5, 89.0, and 150.0 GHz, for both horizontal and vertical polarizations. Two measurement configurations were used for radiometric measurements: a reflecting surface and an absorbing base beneath the snow slabs. Simulations of brightness temperatures using two microwave emission models, the Helsinki University of Technology (HUT) snow emission model and Microwave Emission Model of Layered Snowpacks (MEMLS), were compared to observed brightness temperatures. RMSE and bias were calculated; with the RMSE and bias values being smallest upon an absorbing base at vertical polarization. Simulations overestimated the brightness temperatures on absorbing base cases at horizontal polarization. With the other experimental conditions, the biases were small, with the exception of the HUT model 36.5 GHz simulation, which produced an underestimation for the reflector base cases. This experiment provides a solid framework for future research on the extinction of microwave radiation inside snow.

Differences Between the HUT Snow Emission Model and MEMLS and Their Effects on Brightness Temperature Simulation

Microwave emission models are a critical component of snow water equivalent retrieval algorithms applied to passive microwave measurements. Several such emission models exist, but their differences need to be systematically compared. This paper compares the basic theories of two models: the multiple-layer Helsinki University of Technology (HUT) model and the microwave emission model of layered snowpacks (MEMLS). By comparing the mathematical formulation side by side, three major differences were identified: 1) by assuming that the scattered intensity is mostly (96%) in the forward direction, the HUT model simplifies the radiative transfer equation in $4\pi$ space into two one-flux equations, whereas MEMLS uses a two-flux theory; 2) the HUT scattering coefficient is much larger than the one of MEMLS; and 3) MEMLS considers the trapped radiation inside snow due to internal reflection by a six-flux model, which is not included in HUT. Simulation experiments indicate that the large scattering coefficient of the HUT model compensates for its large forward scattering ratio to some extent, but the effects of one-flux simplification and the trapped radiation still result in different $T_{B}$ simulations between the HUT model and MEMLS. The models were compared with observations of natural snow cover at Sodankyla, Finland; Churchill, Canada; and Colorado, USA. No optimization of the snow grain size was performed. It shows that the HUT model tends to underestimate $T_{B}$ for deep snow. MEMLS with the physically based improved Born approximation performed best among the models, with a bias of −1.4 K and a root-mean-square error of 11.0 K.

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

Abstract. The 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 tomography. 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 tomography 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.

Error Characterization of Coupled Land Surface-Radiative Transfer Models for Snow Microwave Radiance Assimilation

Snow microwave radiance assimilation (RA) or brightness temperature data assimilation (DA) has shown promise for improving snow water equivalent (SWE) estimation. A successful RA study requires, however, an analysis of the error characteristics of coupled land surface-radiative transfer models (LSM/RTMs). This paper focuses on the Community Land Model version 4 (CLM4) as the land-surface model and on the microwave emission model for layered snowpacks (MEMLS) and the dense media radiative transfer multilayer (DMRT-ML) model as RTMs. Using the National Aeronautics and Space Administration Cold Land Processes Field Experiment (CLPX) data sets and through synthetic experiments, the errors of the coupled CLM4/DMRT-ML and CLM4/MEMLS are characterized by: 1) evaluating the CLM4 snowpack state simulations; 2) assessing the performance of RTMs in simulating the brightness temperature $(T_{B} )$ ; and 3) analyzing the correlations between the SWE error $(\varepsilon\_\text{SWE})$ and the $T_{B}$ error $(\varepsilon\_T_{B} ) $ from the RA perspective. The results using the CLPX data sets show that, given a large error of the snow grain radius $(\varepsilon\_r_{e} )$ under dry snowpack conditions (along with a small error of the snow temperature $(\varepsilon\_T_\mathrm{snow})$ ), the correlations between $\varepsilon\_\text{SWE}$ and $\varepsilon\_T_{B} $ are mainly determined by the relationship between $\varepsilon\_r_{e}$ and the snow depth error $(\varepsilon\_d_{\mathrm{snow}})$ or the snow density error $(\varepsilon\_\rho_{\mathrm{snow}} ) $ . The synthetic experiments were carried out for the CLPX region (shallow snowpack conditions and the Rocky Mountains (deep snowpack conditions using the atmospheric ensemble reanalysis produced by the coupled DA Research Testbed/Community Atmospheric Model (CAM4. The synthetic experiments support the results from the CLPX data sets and show that the errors of soil (the water content and the temperature, snow wetness, and snow temperature mostly result in positive correlations between $\varepsilon\_\text{SWE}$ and $\varepsilon\_T_{B} $ . CLM4/DMRT-ML and CLM4/MEMLS tend to produce varying RA performance, with more positive and negative correlations between $\varepsilon\_\text{SWE}$ and $\varepsilon\_T_{B} $ , respectively. These results suggest the necessity of using multiple snowpack RTMs in RA to improve the SWE estimation at the continental scale. The results in this paper also show that the magnitude of $\varepsilon\_r_{e}$ and its relationship to $\varepsilon\_\text{SWE} $ are important for the RA performance. Most of the SWE estimations in RA are improved when $\varepsilon\_\text{SWE}$ and $\varepsilon\_r_{e} $ show a high positive correlation (greater than 0.5).

MEMLS3&amp;a: Microwave Emission Model of Layered Snowpacks adapted to include backscattering

Abstract. The Microwave Emission Model of Layered Snowpacks (MEMLS) was originally developed for microwave emissions of snowpacks in the frequency range 5–100 GHz. It is based on six-flux theory to describe radiative transfer in snow including absorption, multiple volume scattering, radiation trapping due to internal reflection and a combination of coherent and incoherent superposition of reflections between horizontal layer interfaces. Here we introduce MEMLS3&amp;amp;a, an extension of MEMLS, which includes a backscatter model for active microwave remote sensing of snow. The reflectivity is decomposed into diffuse and specular components. Slight undulations of the snow surface are taken into account. The treatment of like- and cross-polarization is accomplished by an empirical splitting parameter q. MEMLS3&amp;amp;a (as well as MEMLS) is set up in a way that snow input parameters can be derived by objective measurement methods which avoid fitting procedures of the scattering efficiency of snow, required by several other models. For the validation of the model we have used a combination of active and passive measurements from the NoSREx (Nordic Snow Radar Experiment) campaign in Sodankylä, Finland. We find a reasonable agreement between the measurements and simulations, subject to uncertainties in hitherto unmeasured input parameters of the backscatter model. The model is written in Matlab and the code is publicly available for download through the following website: http://www.iapmw.unibe.ch/research/projects/snowtools/memls.html.

Snow thickness retrieval from L-band brightness temperatures: a model comparison

AbstractThe Soil Moisture and Ocean Salinity (SMOS) satellite’s L-band (1.4 GHz) measurements have been used to retrieve Snow thickness over thick sea Ice in a previous study. Here we consider brightness temperature simulations for 2.5–4.5m thick Arctic multi-year Ice and compare the results of the relatively simple emission model (M2013) used previously for the retrieval with simulations from a more complex model (T2011) that combines a sea-Ice version of the Microwave Emission Model for Layered Snowpacks (MEMLS) with a thermodynamic model. We find that L-band brightness temperature is mainly determined by Ice temperature. In the M2013 model, Ice temperature in turn is mainly determined by surface temperature and Snow thickness, and this dependence has been used previously to explain the potential for a Snow thickness retrieval. Our comparisons suggest that the M2013 retrieval model may benefit from a more sophisticated thermodynamic calculation of the Ice temperature or from using independent temperature data (e.g. from 6 GHz channels). In both models, horizontally polarized brightness temperatures increase with Snow thickness while holding surface temperature, Ice thickness and Snow density near constant. The increase in the T2011 model is steeper than in M2013, suggesting a higher sensitivity to Snow thickness than found earlier.

The microwave emissivity variability of snow covered first-year sea ice from late winter to early summer: a model study

Abstract. Satellite observations of microwave brightness temperatures between 19 GHz and 85 GHz are the main data sources for operational sea-ice monitoring and retrieval of ice concentrations. However, microwave brightness temperatures depend on the emissivity of snow and ice, which is subject to pronounced seasonal variations and shows significant hemispheric contrasts. These mainly arise from differences in the rate and strength of snow metamorphism and melt. We here use the thermodynamic snow model SNTHERM forced by European Re-Analysis (ERA) interim data and the Microwave Emission Model of Layered Snowpacks (MEMLS), to calculate the sea-ice surface emissivity and to identify the contribution of regional patterns in atmospheric conditions to its variability in the Arctic and Antarctic. The computed emissivities reveal a pronounced seasonal cycle with large regional variability. The emissivity variability increases from winter to early summer and is more pronounced in the Antarctic. In the pre-melt period (January–May, July–November) the standard deviations in surface microwave emissivity due to diurnal, regional and inter-annual variability of atmospheric forcing reach up to Δε = 0.034, 0.043, and 0.097 for 19 GHz, 37 GHz and 85 GHz channels, respectively. Between 2000 and 2009, small but significant positive emissivity trends were observed in the Weddell Sea during November and December as well as in Fram Strait during February, potentially related to earlier melt onset in these regions. The obtained results contribute to a better understanding of the uncertainty and variability of sea-ice concentration and snow-depth retrievals in regions of high sea-ice concentrations.