Remote sensing of three-dimensional cirrus clouds from satellites: application to continuous-wave laser atmospheric transmission and backscattering

A satellite remote sensing methodology has been developed to retrieve 3D ice water content (IWC) and mean effective ice crystal size of cirrus clouds from satellite data on the basis of a combination of the conventional retrieval of cloud optical depth and particle size in a horizontal plane and a parameterization of the vertical cloud profile involving temperature from sounding and/or analysis. The inferred 3D cloud fields of IWC and mean effective ice crystal size associated with two impressive cirrus clouds that occurred in the vicinity of northern Oklahoma on 18 April 1997 and 9 March 2000, obtained from the Department of Energy's Atmospheric Radiation Measurement Program, have been validated against the ice crystal size distributions that were collected independently from collocated and coincident aircraft optical probe measurements. The 3D cloud results determined from satellite data have been applied to the simulation of cw laser energy propagation, and we show the significance of 3D cloud geometry and inhomogeneity and spherical atmosphere on the transmitted and backscattered laser powers. Finally, we demonstrate that the 3D cloud fields derived from satellite remote sensing can be used for the 3D laser transmission and backscattering model for tactical application.     

Laser transmission-backscattering through inhomogeneous cirrus clouds

We have developed a two-dimensional (2D) model for inhomogeneous cirrus clouds in plane-parallel and spherical geometries for the analysis of the transmission and backscattering of high-energy laser beams. The 2D extinction-coefficient and mean effective ice-crystal size fields for cirrus clouds can be determined from a combination of the remote sensing of cirrus clouds by use of the Advanced Very High Resolution Radiometer on board National Oceanic and Atmospheric Administration satellites and the vertical profiling of ice-crystal size distributions available from limited measurements. We demonstrate that satellite remote sensing of the position and the composition of high cirrus can be incorporated directly in the computer model developed for the transmission and backscattering of high-energy laser beams in realistic atmospheres. The results of laser direct transmission, forward scattering, and backscattering are analyzed carefully with respect to aircraft height, cirrus cloud optical depth, and ice-crystal size and orientation. Uncertainty in laser transmission that is due to errors in the retrieved ice-crystal size is negligible. But uncertainty of the order of 2% can be produced if the retrieved optical depth has errors of +/-0.05. With both the aircraft and the target near the cloud top, the direct transmission decreases, owing to the propagation of the laser beam through the curved portion of the cloud top. This effect becomes more pronounced as the horizontal distance between the aircraft and the target increases.     

Visible and near infrared reflectances measured from laboratory ice clouds

We present laboratory results of the 0.68 microm visible (VIS) and 1.617 microm near infrared (NIR) reflectances typically used for inferring optical depth and ice crystal size from satellite radiometers, from ice clouds generated in a temperature controlled column cloud chamber. Two types of ice crystals were produced in this experiment: small columns and dendrites with mean maximum dimensions of about 17 and 35 microm. Within experimental uncertainty, the measured reflectances from ice clouds at both wavelengths agree reasonably well with the theoretical results computed from the plane-parallel adding-doubling method for radiative transfer using the measured ice particle morphology. We demonstrate that laboratory scattering and reflectance data for thin ice clouds with optical depths less than 0.4 can be used for validation of the thin cirrus optical depth and ice crystal size that have been routinely retrieved from the satellite VIS-NIR two channel pair.     

Classification of particle effective shape ratios in cirrus clouds based on the lidar depolarization ratio

A shape classification technique for cirrus clouds that could be applied to future spaceborne lidars is presented. A ray-tracing code has been developed to simulate backscattered and depolarized lidar signals from cirrus clouds made of hexagonal-based crystals with various compositions and optical depth, taking into account multiple scattering. This code was used first to study the sensitivity of the linear depolarization rate to cloud optical and microphysical properties, then to classify particle shapes in cirrus clouds based on depolarization ratio measurements. As an example this technique has been applied to lidar measurements from 15 mid-latitude cirrus cloud cases taken in Palaiseau, France. Results show a majority of near-unity shape ratios as well as a strong correlation between shape ratios and temperature: The lowest temperatures lead to high shape ratios. The application of this technique to space-borne measurements would allow a large-scale classification of shape ratios in cirrus clouds, leading to better knowledge of the vertical variability of shapes, their dependence on temperature, and the formation processes of clouds.     

Lidar inelastic multiple-scattering parameters of cirrus particle ensembles determined with geometrical-optics crystal phase functions

Multiple-scattering correction factors for cirrus particle extinction coefficients measured with Raman and high spectral resolution lidars are calculated with a radiative-transfer model. Cirrus particle-ensemble phase functions are computed from single-crystal phase functions derived in a geometrical-optics approximation. Seven crystal types are considered. In cirrus clouds with height-independent particle extinction coefficients the general pattern of the multiple-scattering parameters has a steep onset at cloud base with values of 0.5-0.7 followed by a gradual and monotonic decrease to 0.1-0.2 at cloud top. The larger the scattering particles are, the more gradual is the rate of decrease. Multiple-scattering parameters of complex crystals and of imperfect hexagonal columns and plates can be well approximated by those of projected-area equivalent ice spheres, whereas perfect hexagonal crystals show values as much as 70% higher than those of spheres. The dependencies of the multiple-scattering parameters on cirrus particle spectrum, base height, and geometric depth and on the lidar parameters laser wavelength and receiver field of view, are discussed, and a set of multiple-scattering parameter profiles for the correction of extinction measurements in homogeneous cirrus is provided.     

Effect of multiple scattering on depolarization measurements with spaceborne lidars

An analytical model based on the integration of the scattering-angle and light-path manifold has been developed to quantify the effect of multiple scattering on cirrus measurements obtained with elastic polarization lidars from space. Light scattering by molecules and by a horizontally homogeneous cloud is taken into account. Lidar parameter, including laser beam divergence, can be freely chosen. Up to 3 orders of scattering are calculated. Furthermore, an inversion technique for the retrieval of cloud extinction profiles from measurements with elastic-backscatter lidars is proposed that explicitly takes multiple scattering into account. It is found that for typical lidar system parameters such as those of the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) instrument multiple scattering does not significantly affect depolarization-ratio measurements in cirrus clouds with small to moderate optical depths. For all simulated clouds, the absolute value of the difference between measured and single-scattering volume depolarization ratio is < 0.006. The particle depolarization ratio can be calculated from the measured volume depolarization ratio and the retrieved backscatter ratio without degradation of accuracy; thus characterization of the various cirrus categories in terms of the particle depolarization ratio and retrieval of cloud microphysical properties is feasible from space. The results of this study apply to polar stratospheric clouds as well.     

Lidar ratio and depolarization ratio for cirrus clouds

We report on studies of the lidar and the depolarization ratios for cirrus clouds. The optical depth and effective lidar ratio are derived from the transmission of clouds, which is determined by comparing the backscattering signals at the cloud base and cloud top. The lidar signals were fitted to a background atmospheric density profile outside the cloud region to warrant the linear response of the return signals with the scattering media. An average lidar ratio, 29 +/- 12 sr, has been found for all clouds measured in 1999 and 2000. The height and temperature dependences ofthe lidar ratio, the optical depth, and the depolarization ratio were investigated and compared with results of LITE and PROBE. Cirrus clouds detected near the tropopause are usually optically thin and mostly subvisual. Clouds with the largest optical depths were found near 12 km with a temperature of approximately -55 degrees C. The multiple-scattering effect is considered for clouds with high optical depths, and this effect lowers the lidar ratios compared with a single-scattering condition. Lidar ratios are in the 20-40 range for clouds at heights of 12.5-15 km and are smaller than approximately 30 in height above 15 km. Clouds are usually optically thin for temperatures below approximately -65 degrees C, and in this region the optical depth tends to decrease with height. The depolarization ratio is found to increase with a height at 11-15 km and smaller than 0.3 above 16 km. The variation in the depolarization ratio with the lidar ratio was also reported. The lidar and depolarization ratios were discussed in terms of the types of hexagonal ice crystals.     

Midlatitude cirrus cloud climatology from the Facility for Atmospheric Remote Sensing. IV. Optical displays

In this fourth of a series of papers that describe long-term cloud research at the Facility for Atmospheric Remote Sensing at Salt Lake City, Utah, an approximately 10-year record of polarization lidar and photographic observations is analyzed to characterize the occurrence of optical displays in our local varieties of midlatitude cirrus clouds. The frequencies of occurrence of various types of halo, arc, and corona displays are evaluated according to their appearance and longevity over nominal 1-h observation periods and to the meteorological source of the cirrus. We find that complex halo-arc displays are rare at our locale and that even the so-called common 22 degree halo occurs infrequently as a complete long-lived ring. For example, only approximately 6% of the 1561-h daytime cirrus periods have bright and prolonged 22 degree halos, although a total of 37.3% have some indications of this halo, even if they are brief and fragmentary. Other fairly frequent features are the 22 degree upper tangent arc (8.6%), 22 degree parhelia (8.5%), and solar corona (7.2%). Of the optical displays observed, 83.6% are refraction based, only 1.9% are due to reflection phenomena, and a surprising 15.4% are caused by diffraction. Complex halo-arc displays are disproportionally associated with cirrus formed in tropical or subtropical airflow and also contain more horizontally oriented planar ice crystals. Lidar linear depolarization ratios from a subset of vivid displays show significant differences between halo- and the corona-producing cirrus, reflecting the effects of particle shape. Halos are associated with relatively warm cirrus that contain randomly and horizontally oriented planar ice crystals, whereas the colder corona cirrus produce much stronger depolarization from crystals too small to be uniformly oriented. Comparisons are made with available information from other locales, and we attempt to explain the geographical differences in terms of basic cirrus cloud processes.     

Halos in cirrus clouds: why are classic displays so rare?

Upper tropospheric cirrus clouds consist of hexagonal ice crystals, which geometrical ray-tracing-theory predicts should regularly produce a variety of optical phenomena such as vivid 22 degrees and 46 degrees halos. Yet, cirrus inconsistently generate such optical displays, while a class of more exotic displays are reported, albeit rarely. I review current knowledge of the cirrus cloud microphysical factors that control ice crystal shape, and hence halo/arc formation, but also appeal to halo enthusiasts to help investigate the causes of unusually complex, brilliant, or rare optical displays. Currently, a wealth of meteorological information can be tapped from the Internet to help advance our knowledge of the basic meteorological factors leading to these rare events.     

Corona-producing ice clouds: a case study of a cold mid-latitude cirrus layer

A high (14.0-km), cold (-71.0 degrees C) cirrus cloud was studied by ground-based polarization lidar and millimeter radar and aircraft probes on the night of 19 April 1994 from the Cloud and Radiation Testbed site in northern Oklahoma. A rare cirrus cloud lunar corona was generated by this 1-2-km-deep cloud, thus providing an opportunity to measure the composition in situ, which had previously been assumed only on the basis of lidar depolarization data and simple diffraction theory for spheres. In this case, corona ring analysis indicated an effective particle diameter of ~22 mum. A variety of in situ data corroborates the approximate ice-particle size derived from the passive retrieval method, especially near the cloud top, where impacted cloud samples show simple solid crystals. The homogeneous freezing of sulfuric acid droplets of stratospheric origin is assumed to be the dominant ice-particle nucleation mode acting in corona-producing cirrus clouds. It is speculated that this process results in a previously unrecognized mode of acid-contaminated ice-particle growth and that such small-particle cold cirrus clouds are potentially a radiatively distinct type of cloud.     

Simulation of infrared scattering from ice aggregates by use of a size-shape distribution of circular ice cylinders

The scalar optical properties (extinction coefficient, mass extinction coefficient, single-scattering albedo, and asymmetry parameter) of a distribution of randomly oriented ice aggregates are simulated generally to well within 4% accuracy by use of a size-shape distribution of randomly oriented circular ice cylinders at wavelengths in the terrestrial window region. The single-scattering properties of the ice aggregates are calculated over the whole size distribution function by the finite-difference time-domain and improved geometric optics methods. The single-scattering properties of the size-shape distribution of circular ice cylinders are calculated by the T-matrix method supplemented by scattering solutions obtained from complex-angular-momentum theory. Moreover, radiative-transfer studies demonstrate that the maximum error in brightness temperature space when the size-shape distribution of circular ice cylinders is used to represent scattering from ice aggregates is only approximately 0.4 K The methodology presented should find wide applicability in remote sensing of ice cloud and parameterization of cirrus cloud scalar optical properties in climate models.     

Cirrus cloud iridescence: a rare case study

On the evening of 25 November 1998, a cirrus cloud revealing the pastel colors of the iridescence phenomenon was photographed and studied by a polarization lidar system at the University of Utah Facility for Atmospheric Remote Sensing (FARS). The diffraction of sunlight falling on relatively minute cloud particles, which display spatial gradients in size, is the cause of iridescence. According to the 14-year study of midlatitude cirrus clouds at FARS, cirrus rarely produce even poor iridescent patches, making this particularly long-lived and vivid occurrence unique. In this unusually high (13.2-14.4-km) and cold (-69.7 degrees to -75.5 degrees) tropopause-topped cirrus cloud, iridescence was noted from approximately 6.0 degrees to approximately 13.5 degrees from the Sun. On the basis of simple diffraction theory, this indicates the presence of particles of 2.5-5.5-microm effective diameter. The linear depolarization ratios of delta = 0.5 measured by the lidar verify that the cloud particles were nonspherical ice crystals. The demonstration that ice clouds can generate iridescence has led to the conclusion that iridescence is rarely seen in midlatitude cirrus clouds because populations of such small particles do not exist for long in the presence of the relatively high water-vapor supersaturations needed for ice-particle nucleation.     

Effects of ice-crystal structure on halo formation: cirrus cloud experimental and ray-tracing modeling studies

During the 1986 Project FIRE (First International Satellite Cloud Climatology Project Regional Experiment) field campaign, four 22° halo-producing cirrus clouds were studied jointly from a groundbased polarization lidar and an instrumented aircraft. The lidar data show the vertical cloud structure and the relative position of the aircraft, which collected a total of 84 slides by impaction, preserving the ice crystals for later microscopic examination. Although many particles were too fragile to survive impaction intact, a large fraction of the identifiable crystals were columns and radial bullet rosettes, with both displaying internal cavitations, and radial plate-column combinations. Particles that were solid or displayed only a slight amount of internal structure were relatively rare, which shows that the usual model postulated by halo theorists, i.e., the randomly oriented, solid hexagonal crystal, is inappropriate for typical cirrus clouds. With the aid of new ray-tracing simulations for hexagonal hollow ended column and bullet-rosette models, we evaluate the effects of more realistic ice-crystal structures on halo formation and lidar depolarization and consider why the common halo is not more common in cirrus clouds.     

Error analysis of Raman differential absorption lidar ozone measurements in ice clouds

A formalism for the error treatment of lidar ozone measurements with the Raman differential absorption lidar technique is presented. In the presence of clouds wavelength-dependent multiple scattering and cloud-particle extinction are the main sources of systematic errors in ozone measurements and necessitate a correction of the measured ozone profiles. Model calculations are performed to describe the influence of cirrus and polar stratospheric clouds on the ozone. It is found that it is sufficient to account for cloud-particle scattering and Rayleigh scattering in and above the cloud; boundary-layer aerosols and the atmospheric column below the cloud can be neglected for the ozone correction. Furthermore, if the extinction coefficient of the cloud is ?0.1 km(-1), the effect in the cloud is proportional to the effective particle extinction and to a particle correction function determined in the limit of negligible molecular scattering. The particle correction function depends on the scattering behavior of the cloud particles, the cloud geometric structure, and the lidar system parameters. Because of the differential extinction of light that has undergone one or more small-angle scattering processes within the cloud, the cloud effect on ozone extends to altitudes above the cloud. The various influencing parameters imply that the particle-related ozone correction has to be calculated for each individual measurement. Examples of ozone measurements in cirrus clouds are discussed.     

Determination of cloud effective particle size from the multiple-scattering effect on lidar integration-method temperature measurements

A method is presented that permits the determination of the cloud effective particle size from Raman- or Rayleigh-integration temperature measurements that exploits the dependence of the multiple-scattering contributions to the lidar signals from heights above the cloud on the particle size of the cloud. Independent temperature information is needed for the determination of size. By use of Raman-integration temperatures, the technique is applied to cirrus measurements. The magnitude of the multiple-scattering effect and the above-cloud lidar signal strength limit the method's range of applicability to cirrus optical depths from 0.1 to 0.5. Our work implies that records of stratosphere temperature obtained with lidar may be affected by multiple scattering in clouds up to heights of 30 km and beyond.     

Influence of atmospheric and system parameters on multiple scattering in spaceborne backscatter lidar measurements

We have simulated backscatter signals of spaceborne lidar systems with the help of a Monte Carlo model. Calculations were performed for various combinations of system parameters. As typical examples of atmospheric observation targets, two kinds of cirrus cloud and two kinds of aerosol were considered. Both total multiple scattering and the significance of individual higher scattering orders are discussed. For all cases, an approximate multiple scattering factor F was calculated that can be used to correct the single-scattering lidar equation to account also for multiple scattering.     

Multiple scattering from clear atmosphere obscured by transparent crystal clouds in satellite-borne lidar sensing

The contribution of multiple scattering to a spaceborne lidar return from clear molecular atmosphere obscured by transparent upper-level crystal clouds is assessed by the use of the variance-reduction Monte Carlo technique. High anisotropy of scattering in the forward direction by polydispersions of ice crystals is the basis of a significant effect of multiple scattering for small values of the lidar receiver field of view. Because of scattering by large nonspherical crystal particles, the lidar signal backscattered from the molecular atmosphere under the cloud increases significantly compared with the single-scattering return. The ratio of the multiple-to-single-scattering contributions from the clear atmosphere hidden by the clouds is greater than from the crystal clouds themselves, and it is proportional to the values of cloud optical thickness.     

Airborne digital holographic system for cloud particle measurements

An in-line holographic system for in situ detection of atmospheric cloud particles [Holographic Detector for Clouds (HOLODEC)] has been developed and flown on the National Center for Atmospheric Research C-130 research aircraft. Clear holograms are obtained in daylight conditions at typical aircraft speeds of 100 m s(-1). The instrument is fully digital and is interfaced to a control and data-acquisition system in the aircraft via optical fiber. It is operable at temperatures of less than -30 degrees C and at typical cloud humidities. Preliminary data from the experiment show its utility for studies of the three-dimensional spatial distribution of cloud particles and ice crystal shapes.     

Laboratory measurements of spectral reflection from ice clouds of various habits

Spectral light from 550 to 650 nm reflected from the surface of an ice cloud produced in a temperature-controlled column is measured at seven different angles between 16.7 degrees and 29.9 degrees . Cloud optical depth (tau) is determined from the extinction of a 670 nm laser and is corrected for forward scattering using a Monte Carlo ray-tracing algorithm. Reflection measurements are compared to expectations from a plane-parallel radiative transfer model with input parameters based on the measured tau and a phase function for the observed ice crystal types. The plane-parallel radiative transfer model can be used to interpret the measured reflection for tau less than about 0.4 for this particular experiment, ideal for providing a validation data set to assist with the development of satellite bidirectional remote sensing.