Numerical simulation of coupled vector radiative transfer equations for a multi-layer model of dense scatterers is developed. The effect from ice layer formed during melting is studied to take account of anomalous spectrum of brightness temperature at the SSM/I channels. This approach is applied to analysis of SSM/I data observed over the southern and central Greenland snow, which can show distinctly different spectral characteristics, T B19 > T B85 or T B19 < T B85. Ice layers produced in melting and ice grain sizes can significantly affect functions of thermal emission at SSM/I channels.