In this thesis a new three–dimensional (3D) vector radiative transfer model, the Solver for Polarized Atmospheric Radiative Transfer Applications (SPARTA) is introduced, validated against benchmark results, and applied to scientific problems. SPARTA employs the statistical forward Monte Carlo technique for efficient column–response pixel–based radiance calculations comprising polarization for 3D inhomogeneous cloudless and cloudy atmospheres. By means of SPARTA, two scientific issues in the field of radiative transfer are investigated. A sensitivity study has been conducted to illustrate the errors introduced by neglecting the effects of polarization in radiation simulations. Two atmospheric scenarios have been considered: a pure one–dimensional (1D) Rayleigh atmosphere and two–dimensional (2D) realistic inhomogeneous dust fields. In case of a purely molecular atmosphere, these errors strongly depend on molecular anisotropy, isotropic reflection, and more importantly, on single scattering albedo and optical thickness (saturation occurs for τ close to 1). Overall errors in the reflected field range up to about 10.5%. On the other hand, for rather high optical thickness, the bias induced by ignoring polarization for realistic inhomogeneous atmospheres is negligible (less than 1%). In addition, solar radiative transfer simulations for LIDAR–measured fields of optical properties of Saharan dust have been performed in order to quantify the effects induced by neglecting the horizontal photon transport and internal inhomogeneities (3D radiative effects) in radiance simulations including polarization. Results are presented for two exemplary mineral dust fields constructed from LIDAR observations. For each case, three radiative calculations are investigated: a 1D calculation according to the plane–parallel (1D mode); an Independent Pixel Approximation (IPA mode); and the 2D mode. The differences in domain–averaged normalized radiances of reflection and transmission are insignificant between the 1D or IPA and 2D calculation modes. However, local differences were observed since extinction is hinge on horizontal spatial variability. In the areas with large spatial variability in optical thickness, the radiance fields of the 2D mode differ about ±20% for the first and second Stokes elements (I, Q) from the fields of the 1D mode. This work points to a brand–new field: the quantification of the sensitivity of polarization to 3D radiative effects
