In the framework of further development of a unified computational tool for the needs of biomedical optics, we introduce an electric field Monte Carlo (MC) model for simulation of backscattering of coherent linearly polarized light from a turbid tissue-like scattering medium with a rough surface. We consider the laser speckle patterns formation and the role of surface roughness in the depolarization of linearly polarized light backscattered from the medium. The mutual phase shifts due to the photons' pathlength difference within the medium and due to reflection/refraction on the rough surface of the medium are taken into account. The validation of the model includes the creation of the phantoms of various roughness and optical properties, measurements of co- and cross-polarized components of the backscattered/reflected light, its analysis and extensive computer modeling accelerated by parallel computing on the NVIDIA graphics processing units using compute unified device architecture (CUDA). The analysis of the spatial intensity distribution is based on second-order statistics that shows a strong correlation with the surface roughness, both with the results of modeling and experiment. The results of modeling show a good agreement with the results of experimental measurements on phantoms mimicking human skin. The developed MC approach can be used for the direct simulation of light scattered by the turbid scattering medium with various roughness of the surface. (C) 2016 Society of Photo-Optical Instrumentation Engineers (SPIE)