Polarimetric scattering from inhomogeneous random media of non-spherical scatterers under a pulse incidence is studied. The time-dependent Mueller matrix solution of vector radiative transfer for layering random media is derived. Co-polarized and cross-polarized bistatic and back-scattering are numerically simulated. The shape and intensity of polarized echoes well depict the inhomogeneous fraction profile of random scatterers. Its functional dependence upon the fraction profile, layering thickness, and other parameters are discussed. This technique is applicable to reconstruction of inhomogeneous fraction profile and inversion of the media thickness.