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 cloud particles exhibit a range of shapes and sizes affecting a cloud’s single-scattering properties. Because they cannot be inferred from passive visible/infrared imager measurements, assumptions about the bulk single-scattering properties of ice clouds are fundamental to satellite cloud retrievals and broadband radiative flux calculations. To examine the sensitivity to ice particle model assumptions, three sets of models are used in satellite imager retrievals of ice cloud fraction, thermodynamic phase, optical depth, effective height, and particle size, and in top-of-atmosphere (TOA) and surface broadband radiative flux calculations. The three ice particle models include smooth hexagonal ice columns (SMOOTH), roughened hexagonal ice columns, and a two-habit model (THM) comprising an ensemble of hexagonal columns and 20-element aggregates. While the choice of ice particle model has a negligible impact on daytime cloud fraction and thermodynamic phase, the global mean ice cloud optical depth retrieved from THM is smaller than from SMOOTH by 2.3 (28%), and the regional root-mean-square difference (RMSD) is 2.8 (32%). Effective radii derived from THM are 3.9 μm (16%) smaller than SMOOTH values and the RMSD is 5.2 μm (21%). In contrast, the regional RMSD in TOA and surface flux between THM and SMOOTH is only 1% in the shortwave and 0.3% in the longwave when a consistent ice particle model is assumed in the cloud property retrievals and forward radiative transfer model calculations. Consequently, radiative fluxes derived using a consistent ice particle model assumption throughout provide a more robust reference for climate model evaluation compared to ice cloud property retrievals.

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.

The Multi‐angle Imaging SpectroRadiometer (MISR) provides measurements over a wider scattering‐angle range for a given location than a cross‐track scanning sensor such as the MODerate resolution Imaging Spectroradiometer (MODIS). Based on a full year (2013) of fused MISR‐MODIS datasets, we develop a variable surface roughness model for ice particles with the goal of identifying the optimal degree of roughness in the ice model for a given pixel containing single‐layer ice clouds. For the MISR‐MODIS observations over oceans, severe roughness values are often selected for a pixel when optical thickness (τ) and particle effective radius (R_eff) are large in conjunction with larger cloud heterogeneity index (H_σ) or a warmer cloud top temperature. Furthermore, τ, R_eff, and ice water path are retrieved with the optimal model and compared to operational MODIS Collection 6 (MC6) products that assume a constant roughness. In general, the retrievals based on the present optimal model lead to greater consistency with MISR measurements, and result in larger median τ by 10.1% and smaller median R_eff by 6.5% but almost identical ice water path in comparison with the MC6 counterparts. The higher average τ value is caused by a slightly larger number of large τ cases, but the smaller average R_eff value is due to the shifting of the retrieved R_eff value toward smaller values by approximately 2–4 μm in comparison to the MC6 distribution over all seasons. Both τ retrievals have similar regional and monthly variations, but a larger annual cycle of R_eff is associated with the optimal model.

The inference of ice cloud properties from remote sensing data depends on the assumed forward ice particle model, as they are used in the radiative transfer simulations that are part of the retrieval process. The Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 (MC6) ice cloud property retrievals are produced in conjunction with a single-habit ice particle model with a fixed degree of ice particle surface roughness (the MC6 model). In this study, we examine the MC6 model and five other ice models with either smoother or rougher surface textures to determine an optimal model to reproduce the angular variation of the radiation field sampled by the Multi-angle Imaging Spectroradiometer (MISR) as a function of latitude. The spherical albedo difference (SAD) method is used to infer an optimal ice particle model. The method is applied to collocated MISR and MODIS data over ocean for clouds with temperatures ≤233 K during December solstice from 2012–2015. The range of solar zenith angles covered by the MISR cameras is broader at the solstices than at other times of the year, with fewer scattering angles associated with sun glint during the December solstice than the June solstice. The results suggest a latitudinal dependence in an optimal ice particle model, and an additional dependence on the solar zenith angle (SZA) at the time of the observations. The MC6 model is one of the most optimal models on the global scale. In further analysis, the results are filtered by a cloud heterogeneity index to investigate cloudy scenarios that are less susceptible to potential 3D effects. Compared to results for global data, the consistency between measurements and a given model can be distinguished in both the tropics and extra-tropics. The SAD analysis suggests that the optimal model for thick homogeneous clouds corresponds to more roughened ice particles in the tropics than in the extra-tropics. While the MC6 model is one of the models most consistent with the global data, it may not be the most optimal model for the tropics.

Structural patterns of cloud effective radius (ER) and liquid water content (LWC) profiles are essential variables of cloud lifecycle and precipitation processes, while observing cloud profiles from passive remote-sensing sensors remains highly challenging. Understanding whether there are typical structural patterns of ER and LWC profiles in liquid clouds and how they are linked to cloud entrainment or precipitating status is critical in developing algorithms to derive cloud profiles from passive satellite sensors. This study aims to address these questions and provide a preliminary foundation for the development of liquid cloud profile retrievals for the Multi-viewing, Multi-channel and Multi-polarization Imaging (3MI) sensor aboard the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Polar System-Second Generation (EPS-SG) satellite, which is scheduled to be launched in 2025. Firstly, we simulate a large ensemble of stratocumulus cloud profiles using the Colorado State University (CSU) Regional Atmospheric Modeling System (RAMS). The empirical orthogonal function (EOF) analysis is adopted to describe the shape of simulated profiles with a limited number of elemental profile variations. Our results indicate that the first three EOFs of LWC and ER profiles can explain >90 % of LWC and ER profiles. The profiles are divided into four prominent patterns and all of these patterns can be simplified as triangle-shaped polylines. The frequency of these four patterns is found to relate to intensities of the cloud-top entrainment and precipitation. Based on these analyses, we propose a simplified triangle-shaped cloud profile parameterization scheme allowing us to represent these main patterns of LWC and ER. This simple yet physically realistic analytical model of cloud profiles is expected to facilitate the representation of cloud properties in advanced retrieval algorithms such as those developed for the 3MI/EPS-SG.