A Review of Ice Cloud Optical Property Models for Passive Satellite Remote Sensing

Correction of Thin Cirrus Absorption Effects in Landsat 8 Thermal Infrared Sensor Images Using the Operational Land Imager Cirrus Band on the Same Satellite Platform

Data from the Operational Land Imager (OLI) and the Thermal Infrared Sensor (TIRS) instruments onboard the Landsat 8 and Landsat 9 satellite platforms are subject to contamination by cloud cover, with cirrus contributions being the most difficult to detect and mask. To help address this issue, a cirrus detection channel (Band 9) centered within the 1.375-μm water vapor absorption region was implemented on OLI, with a spatial resolution of 30 m. However, this band has not yet been fully utilized in the Collection 2 Landsat 8/9 Level 2 surface temperature data products that are publicly released by U.S. Geological Survey (USGS). The temperature products are generated with a single-channel algorithm. During the surface temperature retrievals, the effects of absorption of infrared radiation originating from the warmer earth's surfaces by ice clouds, typically located in the upper portion of the troposphere and re-emitting at much lower temperatures (approximately 220 K), are not taken into consideration. Through an analysis of sample Level 1 TOA and Level 2 surface data products, we have found that thin cirrus cloud features present in the Level 1 1.375-μm band images are directly propagated down to the Level 2 surface data products. The surface temperature errors resulting from thin cirrus contamination can be 10 K or larger. Previously, we reported an empirical and effective technique for removing thin cirrus scattering effects in OLI images, making use of the correlations between the 1.375-μm band image and images of any other OLI bands located in the 0.4-2.5 μm solar spectral region. In this article, we describe a variation of this technique that can be applied to the thermal bands, using the correlations between the Level 1 1.375-μm band image and the 11-μm BT image for the effective removal of thin cirrus absorption effects. Our results from three data sets acquired over spatially uniform water surfaces and over non-uniform land/water boundary areas suggest that if the cirrus-removed TOA 11-μm band BT images are used for the retrieval of the Level 2 surface temperature (ST) data products, the errors resulting from thin cirrus contaminations in the products can be reduced to about 1 K for spatially diffused cirrus scenes.

The Impacts of Single-Scattering and Microphysical Properties of Ice Particles Smaller Than 100 µm on the Bulk Radiative Properties of Tropical Cirrus

There are large uncertainties in the single-scattering (i.e., morphologies) and microphysical (i.e., concentrations) properties of ice particles whose size are less than ~100 µm. Insufficient resolutions of the most advanced cloud probes (e.g., cloud particle imager) cannot resolve the micrometer-scale morphologies of small ice particles. Further, the shattering of large ice particles on probes’ inlets or tips causes uncertainties in the measurement of the concentrations of small ice particles. These uncertainties have large impacts on the single-scattering and microphysical properties of small ice particles that are utilized to quantify the bulk radiative properties of cirrus. In this study, the impacts of uncertainties in the morphologies and concentrations of small ice particles on the bulk radiative properties of tropical cirrus were calculated using measurements acquired during the Tropical Warm Pool-International Cloud Experiment. Five different models (i.e., budding Buckyball, Chebyshev particle, droxtal, Gaussian random sphere, and sphere) that represent the shapes of small ice particles were used to calculate the single-scattering properties. The bulk radiative properties, average phase-function (P11¯), and average asymmetry parameter (g¯) were computed by combining the measured size/habit distributions and the calculated single-scattering properties of ice particles. The impacts of the selection of varying morphologies of small particles on the bulk radiative properties were quantified. For these calculations, the possible range of the concentrations of small ice particles which depend on the degree of shattered large particles were also used. The impacts of varying the single-scattering properties of small ice particles on the bulk radiative properties were the largest in the upper parts of cirrus (T < −60 °C), while they were the smallest in the lower parts of cirrus (−45 < T < −30 °C). The impacts of uncertainties in the concentrations of small ice particles on the bulk radiative properties were largest in the lower parts of cirrus (−45 < T < −30 °C), whereas they were smallest in the upper parts of cirrus (T < −60 °C). The effect of shattering was maximum in the lower parts of cirrus, whilst it was minimum in the upper parts of cirrus. The combined impacts of uncertainties in the single-scattering (i.e., morphologies) and microphysical (i.e., concentrations) properties of small ice particles revealed variations of up to 11.2% (127.1%; 67.3%) of the integrated intensity in the forward (sideward; backward) angles in P11¯ and a corresponding change in g¯ by up to 12.61%.

Optical Property Model for Cirrus Clouds Based on Airborne Multi-Angle Polarization Observations

We present an improved remote sensing technique to infer an optimal habit/shape model for ice particles in cirrus clouds using multi-angle polarimetric measurements at 865 nm made by the Airborne Multi-angle SpectroPolarimeter Imager (AirMSPI) instrument. The common method of ice model inference is based on intensity (total reflectivity) measurements, which is generally not applicable to optically thin ice clouds (i.e., cirrus clouds) where single scattering dominates. The new approach is able to infer an ice model in clouds with optical thicknesses smaller than 5. The improvement is made by first assuming the optical thickness retrieved using total reflectivity. Subsequently, the polarized reflectivity is calculated based on look-up tables generated from simulated polarized reflectances computed for cirrus clouds in conjunction with eight ice particle models. The ice particle model that leads to the closest fit to the measurements is regarded as the optimal ice particle model. Additionally, an alternative method is applied that does not consider polarized reflectivity. These two methods are applied to a data sample as a proof-of-concept study where AirMSPI observed a single cirrus layer. In this case study, the hexagonal column aggregate model works for most pixels both with and without considering polarized reflectivities. The retrieval cost function is high when the camera pairs with large zenith angles are included in the retrievals. This result suggests that further studies will be necessary to have a better understanding of all eight selected ice particle models at scattering angles smaller than 100°.

The Use of Superspheroids as Surrogates for Modeling Electromagnetic Wave Scattering by Ice Crystals

Electromagnetic wave scattering by ice particles is commonly modeled by defining representative habits, including droxtals, columns, plates, and aggregates, although actual particles in the atmosphere can be even much more complex. In this study, we examined a superspheroidal approximation method for modeling electromagnetic wave scattering by ice crystals. Superspheroid can be associated with a shape index (SI) defined by the particle volume and average projected area. Corresponding to realistic ice crystals, suitable superspheroid models with the same SI (that means, identical volume and average projected area) and aspect ratio can be identified as surrogates for optical property calculations. We systematically compared the optical properties of ice crystals and superspheroids at 33 microwave bands in the range of 3–640 GHz and at three representative visible or infrared wavelengths (0.66, 2.13, and 11 μm). It was found that the single-scattering properties of compact ice crystal habits and their superspheroidal model particles were quite close. For an aggregate with sparse distribution of elements, a superspheroid model produces relatively large errors because the aspect ratio may not be sufficient to describe a particle shape. However, the optical similarity of a superspheroid and an aggregate is still encouraging.

Passive Remote Sensing of Ice Cloud Properties at Terahertz Wavelengths Based on Genetic Algorithm

Ice clouds play a critical role in the balance of the earth–atmosphere radiation system, but there are some limitations in the existing remote sensing methods for ice clouds. Terahertz wave is expected to be the best waveband for retrieving ice clouds, with terahertz wavelengths in the order of the size of typical ice cloud particles. An inversion method for the remote sensing of ice clouds at terahertz wavelengths based on genetic algorithm is proposed in this paper. First, suitable channel sets in the terahertz band, which are mainly a combination of absorption lines and window regions, are determined. Then, to improve the efficiency of the generation of the retrieval database, based on the brightness temperature simulated by the atmospheric radiative transfer simulator (ARTS) for different cloud parameters, a fast forward operator is constructed using three-dimensional interpolation to simulate the brightness temperature difference between clear sky and a cloudy scene. Finally, an inversion model to retrieve the ice cloud base height, the effective particle diameter and the ice water path is established based on the genetic algorithm, and an analysis of the inversion errors is performed. The results show that the forward operator, constructed by the nearest interpolation, can accurately calculate the brightness temperature difference at a high speed. The proposed inversion method at terahertz wavelengths based on the genetic algorithm can achieve the expected scientific requirement. The absolute error of the cloud height is around 0.2 km, and the absolute error of the low ice water path (below 20 g/m2) is small, while the relative error of the high ice water path is generally maintained at around 10%, and the absolute error of the effective particle diameter is mostly around 4 μm.