Taking into account the multiple scattering of vector radiative transfer (VRT), an iterative layering method is developed by stratifying the random media into m ulti-slabs and employing the lower-order Mueller matrix solution to derive high- order scatterings. This approach is applied to the simulation of polarized brigh tness temperature from an inhomogeneous snowpack of randomly and densely distrib uted spherical particles with nonuniform profiles of particle size, fractional v olume and physical temperature. Numerical results are compared with conventional eigen-analysis approach for a uniform scatter medium. The functional dependence of multiple scattering and emission upon the characteristic parameters such as particle size, fractional volume, physical temperature etc. are discussed, and a re also compared with the microwave remote sensing measurements of the multi-lay ers snowpack.