Abstract: In order to study the polarization properties of scattered light from far ship wakes, the polarized Monte Carlo program for light propagation was used to simulate the 2-dimensional spatial distributions of degree of polarization(DOP) of ship wakes with different polarized incident lights. Radial distributions of DOP of scattered light from far ship wakes were analyzed with different bubble radius, scattering coefficients and thicknesses of bubbles film. The results show that spatial distributions of DOP have the properties of azimuth selecting which decreases with the increase of bubble radius, and the difference could be observed with different polarization states of the incident light. Besides, DOP of forward scattering light of ship wakes decreases with the increase of the scattering coefficients and thicknesses of bubble film.