Scattering processes in the atmospheres of planets cause characteristic features that can be particularly well observed in polarisation. For planet Earth, both molecular and scattering by small particles imprint specific signatures in its phase curve. An unequivocal prediction of a liquid-water-loaded atmosphere is the existence of a rainbow feature at a scattering angle of around 138-144deg. Earthshine allows us to observe the primary rainbow in linear polarisation. We observed polarisation spectra of Earthshine using FORS2 at the Very Large Telescope for phase angles from 33deg to 65deg (Sun--Earth--Moon angle). The spectra were used to derive the degree of polarisation in the B, V, R, and I passbands and the phase curve from 33deg to 136deg . The new observations extend to the smallest phases that can be observed from the ground. The degree of polarisation of planet Earth is increasing for decreasing phase angles downwards of 45deg. From comparison of the phase curve observed with models of an Earth-type atmosphere we are able to determine the refractive index of water and to constrain the mean water droplet sizes to 6-7 mum. Furthermore, we can retrieve the mean cloud fraction of liquid water clouds to 0.3, and the mean optical depth of the water clouds to values between 10 and 20. Our observations allow us to discern two fundamentally different scattering mechanisms of the atmosphere of planet Earth: molecular and particle scattering. The physical and chemical properties can be retrieved with high fidelity through suitable inversion of the phase curve. Observations of polarimetric phase curves of planets beyond the Solar System shall be extremely valuable for a thorough characterisation of their atmospheres.

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.

We present a comprehensive set of spectropolarimetric observations of Earthshine as obtained by FORS2 at the VLT for phase angles from 50degree to 135degree (Sun-Earth-Moon angle), covering a spectral range from 430nm to 920nm. The degree of polarization in BVRI passbands, the differential polarization vegetation index, and the equivalent width of the O2A polarization band around 760nm are determined with absolute errors around 0.1 percent in the degree of polarization. Earthshine polarization spectra are corrected for the effect of depolarization introduced by backscattering on the lunar surface, introducing systematic errors of the order of 1 percent in the degree of polarization. Distinct viewing sceneries such as observing the Atlantic or Pacific side in Earthshine yield statistically different phase curves. The equivalent width defined for the O2A band polarization is found to vary from -5nm to +2nm. A differential polarized vegetation index is introduced and reveals a larger vegetation signal for those viewing sceneries that contain larger fractions of vegetated surface areas. We corroborate the observed correlations with theoretical models from the literature, and conclude that the Vegetation Red Edge(VRE) is a robust and sensitive signature in polarization spectra of planet Earth. The overall behaviour of polarization of planet Earth in the continuum and in the O2A band can be explained by existing models. Biosignatures such as the O2A band and the VRE are detectable in Earthshine polarization with a high degree of significance and sensitivity. An in-depth understanding of Earthshines temporal and spectral variability requires improved models of Earths biosphere, as a prerequisite to interpret possible detections of polarised biosignatures in earthlike exoplanets in the future.

Ground-based observations of the Earthshine, i.e., the light scattered by Earth to the Moon, and then reflected back to Earth, simulate space observations of our planet and represent a powerful benchmark for the studies of Earth-like planets. Earthshine spectra are strongly linearly polarized, owing to scattering by molecules and small particles in the atmosphere of the Earth and surface reflection, and may allow us to measure global atmospheric and surface properties of planet Earth. Aims. We aim to interpret already published spectropolarimetric observations of the Earthshine by comparing them with new radiative transfer model simulations including a fully realistic three-dimensional (3D) surface-atmosphere model for planet Earth. We used the highly advanced Monte Carlo radiative transfer model MYSTIC to simulate polarized radiative transfer in the atmosphere of the Earth without approximations regarding the geometry, taking into account the polarization from surface reflection and multiple scattering by molecules, aerosol particles, cloud droplets, and ice crystals. We have shown that Earth spectropolarimetry is highly sensitive to all these input parameters, and we have presented simulations of a fully realistic Earth atmosphere-surface model including 3D cloud fields and two-dimensional (2D) surface property maps. Our modeling results show that scattering in high ice water clouds and reflection from the ocean surface are crucial to explain the continuum polarization at longer wavelengths as has been reported in Earthshine observations taken at the Very Large Telescope in 2011 (3.8 % and 6.6 % at 800 nm, depending on which part of Earth was visible from the Moon at the time of the observations). We found that the relatively high degree of polarization of 6.6 % can be attributed to light reflected by the ocean surface in the sunglint region.

Ice crystal clouds in the upper troposphere can generate polarisation signals at the uK level. This signal can seriously affect very sensitive ground based searches for E- and B-mode of Cosmic Microwave Background polarisation. In this paper we estimate this effect within the ClOVER experiment observing bands (97, 150 and 220 GHz) for the selected observing site (Llano de Chajnantor, Atacama desert, Chile). The results show that the polarisation signal from the clouds can be of the order of or even bigger than the CMB expected polarisation. Climatological data suggest that this signal is fairly constant over the whole year in Antarctica. On the other hand the stronger seasonal variability in Atacama allows for a 50% of clean observations during the dry season.