The Markov chain formalism for polarized radiative transfer through a vertically inhomogeneous atmosphere is linearized comprehensively with respect to the aerosol and polarizing surface properties. For verification, numerical results are compared to those obtained by the finite difference method. We demonstrate the use of the linearized code as part of a retrieval of aerosol and surface properties for an atmosphere overlying a black and Fresnel-reflecting ocean surface.