The polarized Monte Carlo (PMC) model has been applied to study the backscattering measurement of oceanic lidar. This study proposes a PMC model for shipborne oceanic lidar simulation. This model is validated by the Rayleigh scattering experiment, lidar equation, and in-situ lidar LOOP (Lidar for Ocean Optics Profiler) returns [Opt. Express30, 8927 (2022)10.1364/OE.449554]. The relative errors of the simulated Rayleigh scattering results are less than 0.07%. The maximum mean relative error (MRE) of the simulated single scattering scalar signals and lidar equation results is 30.94%. The maximum MRE of simulated total scattering signals and LOOP returns in parallel and cross channels are 33.29% and 22.37%, respectively, and the maximal MRE of the depolarization ratio is 24.13%. The underwater light field of the laser beam is also simulated to illustrate the process of beam energy spreading. These results prove the validity of the model. Further analyses show that the measured signals of shipborne lidar LOOP are primarily from the particle single scatterings. This model is significant for analyzing the signal contributions from multiple scattering and single scattering.