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 (DMRTQMS); 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.