The polarization state of electromagnetic radiation scattered by atmospheric particles such as aerosols, cloud droplets, or ice crystals contains much more information about the optical and microphysical properties than the total intensity alone. For this reason an increasing number of polarimetric observations are performed from space, from the ground and from aircraft. Polarized radiative transfer models are required to interpret and analyze these measurements and to develop retrieval algorithms exploiting polarimetric observations. In the last years a large number of new codes have been developed, mostly for specific applications. Benchmark results are available for specific cases, but not for more sophisticated scenarios including polarized surface reflection and multi-layer atmospheres. The International Polarized Radiative Transfer (IPRT) working group of the International Radiation Commission (IRC) has initiated a model intercomparison project in order to fill this gap. This paper presents the results of the first phase A of the IPRT project which includes ten test cases, from simple setups with only one layer and Rayleigh scattering to rather sophisticated setups with a cloud embedded in a standard atmosphere above an ocean surface. All scenarios in the first phase A of the intercomparison project are for a one-dimensional plane-parallel model geometry. The commonly established benchmark results are available at the IPRT website (http://www.meteo.physik.uni-muenchen.de/~iprt).

Initially unpolarized solar radiation becomes polarized by scattering in the Earth's atmosphere. In particular molecular scattering polarizes electromagnetic radiation, but also scattering of radiation at aerosols, cloud droplets and ice crystals polarizes. Each atmospheric constituent produces a characteristic polarization signal, thus spectro-polarimetric measurements are frequently employed for remote sensing of aerosol and cloud properties. Retrieval algorithms require efficient radiative transfer models. Usually, these apply the plane-parallel approximation, assuming that the atmosphere consists of horizontally homogeneous layers. For remote sensing applications, the radiance is considered constant over the instantaneous field-of-view of the instrument and each sensor element is treated independently in plane-parallel approximation, neglecting horizontal radiation transport between adjacent pixels. In order to estimate the errors due to the IPA approximation, three-dimensional (3D) vector radiative transfer models are required. So far, only a few such models exist. Therefore, the International Polarized Radiative Transfer (IPRT) working group of the International Radiation Commission (IRC) has initiated a model intercomparison project in order to provide benchmark results for polarized radiative transfer. The group has already performed an intercomparison for one-dimensional (1D) multi-layer test cases (Emde et al., 2015). This paper presents the continuation of the intercomparison project for 2D and 3D test cases: a step cloud, a cubic cloud, and a more realistic scenario including a 3D cloud field generated by a Large Eddy Simulation model and typical background aerosols. The commonly established benchmark results for 3D polarized radiative transfer are available at the IPRT website (http://www.meteo.physik.uni-muenchen.de/~iprt).

Remote sensing instruments observe radiation being scattered and absorbed by molecules, aerosol particles, cloud droplets and ice crystals. In order to interpret and accurately model such observations, the vector radiative transfer equation needs to be solved, because scattering polarizes the initially unpolarized incoming solar radiation. A widely used approximation in radiative transfer theory is the neglect of polarization which allows to greatly simplify the radiative transfer equation. It is well known that the error caused by multiple Rayleigh scattering can be larger than 10\%, depending on wavelength and sun-observer geometry (Mishchenko et al., 1994). For homogeneous plane-parallel layers of liquid cloud droplets the error is comparatively small (below 1\%) (Hansen 1971). However, in reality clouds are not plane-parallel layers of water droplets but complex three-dimensional (3D) structures and observations of clouds usually include pixels consisting of clear and cloudy parts. In this study we revisit the question of the magnitude of error due to the neglect of polarization in radiative transfer theory for a realistic 3D cloudy atmosphere. We apply the Monte Carlo radiative transfer model MYSTIC with and without neglecting polarization and compare the results. At a phase angle of 90{\deg} and 400nm wavelength we find the maximum overestimation error of about 8% for complete clear-sky conditions. The error is reduced to about 6% in clear-sky regions surrounded by clouds due to scattering from clouds into the clear regions. Within the clouds the error is up to 4% with the highest values in cloud shadows. In backscattering direction the radiance is underestimated by about 5% in clear regions between clouds. For other sun-observer geometries, the error ranges between the two extremes. The error decreases with wavelength and in the absorption bands.

An efficient method to compute accurate polarized solar radiance spectra using the (three-dimensional) Monte Carlo model MYSTIC has been developed. Such high resolution spectra are measured by various satellite instruments for remote sensing of atmospheric trace gases. ALIS (Absorption Lines Importance Sampling) allows the calculation of spectra by tracing photons at only one wavelength: In order to take into account the spectral dependence of the absorption coefficient a spectral absorption weight is calculated for each photon path. At each scattering event the local estimate method is combined with an importance sampling method to take into account the spectral dependence of the scattering coefficient. Since each wavelength grid point is computed based on the same set of random photon paths, the statistical error is the almost same for all wavelengths and hence the simulated spectrum is not noisy. The statistical error mainly results in a small relative deviation which is independent of wavelength and can be neglected for those remote sensing applications where differential absorption features are of interest. Two example applications are presented: The simulation of shortwave-infrared polarized spectra as measured by GOSAT from which CO$_2$ is retrieved, and the simulation of the differential optical thickness in the visible spectral range which is derived from SCIAMACHY measurements to retrieve NO$_2$. The computational speed of ALIS (for one- or three-dimensional atmospheres) is of the order of or even faster than that of one-dimensional discrete ordinate methods, in particular when polarization is considered.