We have implemented Raman scattering in a vector radiative transfer model for coupled atmosphere and ocean systems. A sensitivity study shows that the Raman scattering contribution is greatest in clear waters and at longer wavelengths. The Raman scattering contribution may surpass the elastic scattering contribution by several orders of magnitude at depth. The degree of linear polarization in water is smaller when Raman scattering is included. The orientation of the polarization ellipse shows similar patterns for both elastic and inelastic scattering contributions. As polarimeters and multipolarization-state lidars are planned for future Earth observing missions, our model can serve as a valuable tool for the simulation and interpretation of these planned observations.