The role of clouds remains the largest uncertainty in climate projections. They influence solar and thermal radiative transfer and the earth's water cycle. Therefore, there is an urgent need for accurate cloud observations to validate climate models and to monitor climate change. Passive satellite imagers measuring radiation at visible to thermal infrared (IR) wavelengths provide a wealth of information on cloud properties. Among others, the cloud top height (CTH) – a crucial parameter to estimate the thermal cloud radiative forcing – can be retrieved. In this paper we investigate the skill of ten current retrieval algorithms to estimate the CTH using observations from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard Meteosat Second Generation (MSG). In the first part we compare ten SEVIRI cloud top pressure (CTP) data sets with each other. The SEVIRI algorithms catch the latitudinal variation of the CTP in a similar way. The agreement is better in the extratropics than in the tropics. In the tropics multi-layer clouds and thin cirrus layers complicate the CTP retrieval, whereas a good agreement among the algorithms is found for trade wind cumulus, marine stratocumulus and the optically thick cores of the deep convective system. In the second part of the paper the SEVIRI retrievals are compared to CTH observations from the Cloud–Aerosol LIdar with Orthogonal Polarization (CALIOP) and Cloud Profiling Radar (CPR) instruments. It is important to note that the different measurement techniques cause differences in the retrieved CTH data. SEVIRI measures a radiatively effective CTH, while the CTH of the active instruments is derived from the return time of the emitted radar or lidar signal. Therefore, some systematic differences are expected. On average the CTHs detected by the SEVIRI algorithms are 1.0 to 2.5 km lower than CALIOP observations, and the correlation coefficients between the SEVIRI and the CALIOP data sets range between 0.77 and 0.90. The average CTHs derived by the SEVIRI algorithms are closer to the CPR measurements than to CALIOP measurements. The biases between SEVIRI and CPR retrievals range from −0.8 km to 0.6 km. The correlation coefficients of CPR and SEVIRI observations vary between 0.82 and 0.89. To discuss the origin of the CTH deviation, we investigate three cloud categories: optically thin and thick single layer as well as multi-layer clouds. For optically thick clouds the correlation coefficients between the SEVIRI and the reference data sets are usually above 0.95. For optically thin single layer clouds the correlation coefficients are still above 0.92. For this cloud category the SEVIRI algorithms yield CTHs that are lower than CALIOP and similar to CPR observations. Most challenging are the multi-layer clouds, where the correlation coefficients are for most algorithms between 0.6 and 0.8. Finally, we evaluate the performance of the SEVIRI retrievals for boundary layer clouds. While the CTH retrieval for this cloud type is relatively accurate, there are still considerable differences between the algorithms. These are related to the uncertainties and limited vertical resolution of the assumed temperature profiles in combination with the presence of temperature inversions, which lead to ambiguities in the CTH retrieval. Alternative approaches for the CTH retrieval of low clouds are discussed.

The current wealth of spaceborne passive and active measurements from ultraviolet to the infrared wavelengths provides an unprecedented opportunity to construct ice cloud bulk optical property models that lead to consistent ice cloud property retrievals across multiple sensors and platforms. To infer the microphysical and radiative properties of ice clouds from these satellite measurements, the general approach is to assume an ice cloud optical property model that implicitly assumes the habit (shape) and size distributions of the ice particles in these clouds. The assumption is that this ice optical property model will be adequate for global retrievals. In this review paper, we first summarize the key optical properties of individual particles and then the bulk radiative properties of their ensemble, followed by a review of the ice cloud models developed for application to satellite remote sensing. We illustrate that the random orientation condition assumed for ice particles is arguably justified for passive remote sensing applications based on radiometric measurements. The focus of the present discussion is on the ice models used by the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Clouds and Earth’s Radiant Energy System (CERES) science teams. In addition, we briefly review the ice cloud models adopted by the Polarization and Directionality of the Earth’s Reflectance (POLDER) and the Himawari-8 Advanced Himawari Imager (AHI) for ice cloud retrievals. We find that both the MODIS Collection 6 ice model and the CERES two-habit model result in spectrally consistent retrievals.

Ice clouds are an important element in the radiative balance of the earth's climate system, but their microphysical and optical properties still are not well constrained, especially ice particle habit and the degree of particle surface roughness. In situ observations have revealed common ice particle habits and evidence for surface roughness, but these observations are limited. An alternative is to infer the ice particle shape and surface roughness from satellite observations of polarized reflectivity since they are sensitive to both particle shape and degree of surface roughness. In this study an adding–doubling radiative transfer code is used to simulate polarized reflectivity for nine different ice habits and one habit mixture, along with 17 distinct levels of the surface roughness. A lookup table (LUT) is constructed from the simulation results and used to infer shape and surface roughness from PARASOL satellite polarized reflectivity data over the ocean. Globally, the retrievals yield a compact aggregate of columns as the most commonly retrieved ice habit. Analysis of PARASOL data from the tropics results in slightly more aggregates than in midlatitude or polar regions. Some level of surface roughness is inferred in nearly 70% of PARASOL data, with mean and median roughness near σ = 0.2 and 0.15, respectively. Tropical region analyses have 20% more pixels retrieved with particle surface roughness than in midlatitude or polar regions. The global asymmetry parameter inferred at a wavelength of 0.865 μm has a mean value of 0.77 and a median value of 0.75.

The degree of surface roughness of ice particles within thick, cold ice clouds is inferred from multi-directional, multi-spectral satellite polarimetric observations over oceans, assuming a column-aggregate particle habit. An improved roughness inference scheme is employed that provides a more noise-resilient roughness estimate than the conventional best-fit approach. The improvements include the introduction of a quantitative roughness parameter based on empirical orthogonal function analysis and proper treatment of polarization due to atmospheric scattering above clouds. A global 1-month data sample supports the use of a severely roughened ice habit to simulate the polarized reflectivity associated with ice clouds over ocean. The density distribution of the roughness parameter inferred from the global 1-month data sample and further analyses of a few case studies demonstrate the significant variability of ice cloud single-scattering properties. However, the present theoretical results do not agree with observations in the tropics. In the extratropics, the roughness parameter is inferred but 74 % of the sample is out of the expected parameter range. Potential improvements are discussed to enhance the depiction of the natural variability on a global scale.