We conducted Monte Carlo simulations based on anisotropic sclera-mimicking models to examine the polarization features in Mueller matrix polar decomposition (MMPD) parameters during the refractive index matching process, which is one of the major mechanisms of optical clearing. In a preliminary attempt, by changing the parameters of the models, wavelengths, and detection geometries, we demonstrate how the depolarization coefficient and retardance vary during the refractive index matching process and explain the polarization features using the average value and standard deviation of scattering numbers of the detected photons. We also study the depth-resolved polarization features during the gradual progression of the refractive index matching process. The results above indicate that the refractive index matching process increases the depth of polarization measurements and may lead to higher contrast between tissues of different anisotropies in deeper layers. MMPD-derived polarization parameters can characterize the refractive index matching process qualitatively.