In this contribution we discuss the problems of numerical modeling of the scattered radiance field in the coupled atmosphere-ocean systems. The analysis is performed in the framework of the discrete theory of radiative transfer. For the linearized system, the solution of inverse problems is obtained