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.     

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.     

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.     

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.     

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.     

Icy wave-cloud lunar corona and cirrus iridescence

Dual-polarization lidar data and radiosonde data are used to determine that iridescence in cirrus and a lunar corona in a thin wave cloud were caused by tiny ice crystals, not droplets of liquid water. The size of the corona diffraction rings recorded in photographs is used to estimate the mean diameter of the diffracting particles to be 14.6 μm, much smaller than conventional ice crystals. The iridescent cloud was located at the tropopause [~11-13.6 km above mean sea level (ASL)] with temperature near -70 °C, while the more optically pure corona was located at approximately 9.5 km ASL with temperature nearing -60 °C. Lidar cross-polarization ratios of 0.5 and 0.4 confirm that ice formed both the iridescence and the corona, respectively.     

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.     

Measurements of cirrus cloud backscatter color ratio with a two-wavelength lidar

We present observations of cirrus clouds from June 2006 to July 2007 performed by using a two-wavelength lidar located at Hampton University. For this time period, cirrus clouds were observed mostly in 7-13.5 km altitudes. Data analyses have been performed focusing on a color-ratio retrieval. In total, 86,369 samples from 1,689 profiles (1 min average and 15 m range resolution) containing cirrus clouds with attenuated backscatter ratio (ratio of attenuated total backscatter to the molecular backscatter) larger than 10 have been selected. The cirrus color ratio distribution shows a peak value at about 0.88 and a full width at half-maximum of 0.12.     

Improved retrievals of the optical properties of cirrus clouds by a combination of lidar methods

We focus on improvement of the retrieval of optical properties of cirrus clouds by combining two lidar methods. We retrieve the cloud's optical depth by using independently the molecular backscattering profile below and above the cloud [molecular integration (MI) method] and the backscattering profile inside the cloud with an a priori effective lidar ratio [particle integration (PI) method]. When the MI method is reliable, the combined MI-PI method allows us to retrieve the optimal effective lidar ratio. We compare these results with Raman lidar retrievals. We then use the derived optimal effective lidar ratio for retrieval with the PI method for situations in which the MI method cannot be applied.     

Calibration of the 1064 nm lidar channel using water phase and cirrus clouds

Calibration is essential to derive aerosol backscatter coefficients from elastic scattering lidar. Unlike the visible UV wavelengths where calibration is based on a molecular reference, calibration of the 1064 nm lidar channel requires other approaches, which depend on various assumptions. In this paper, we analyze two independent calibration methods which use (i) low-altitude water phase clouds and (ii) high cirrus clouds. In particular, we show that to achieve optimal performance, aerosol attenuation below the cloud base and cloud multiple scattering must be accounted for. When all important processes are considered, we find that these two independent methods can provide a consistent calibration constant with relative differences less than 15%. We apply these calibration techniques to demonstrate the stability of our lidar on a monthly scale, along with a natural reduction of the lidar efficiency on an annual scale. Furthermore, our calibration procedure allows us to derive consistent aerosol backscatter coefficients and angstrom coefficient profiles (532-1064 nm) along with column extinction-to-backscatter ratios which are in good agreement with sky radiometer inversions.     

Iterative method to determine an averaged backscatter-to-extinction ratio in cirrus clouds

An iterative method to determine an average backscatter-to-extinction ratio and extinction coefficient simultaneously in cirrus clouds is proposed. The method is based on Klett's inversion, which is constrained by the total optical depth. A signal-to-noise ratio greater than 3 at the cloud top is required for an error in the backscatter-to-extinction ratio lower than 20% to result. The method has been tested with simulated lidar signals. An application to an experimental lidar signal is discussed.     

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.     

Retrieval of cloud optical parameters from space-based backscatter lidar data

We present an approach to estimating the multiple-scattering (MS) contribution to lidar return signals from clouds recorded from space that enables us to describe in more detail the return formation at the depth where first orders of scattering dominate. Estimates made have enabled us to propose a method for correcting solutions of single-scattering lidar equations for the MS contribution. We also describe an algorithm for reconstructing the profiles of the cloud scattering coefficient and the optical thickness tau under conditions of a priori uncertainties. The approach proposed is illustrated with results for optical parameters of cirrus and stratiform clouds determined from return signals calculated by the Monte Carlo method as well as from return signals acquired with the American spaceborne lidar during the Lidar In-Space Technology Experiment (LITE).     

Determination by spaceborne backscatter lidar of the structural parameters of atmospheric scattering layers

Spaceborne active lidar systems are under development to give new insight into the vertical distribution of clouds and aerosols in the atmosphere and to provide new information on variables required for improvement of forecast models and for understanding the radiative and dynamic processes that are linked to the dynamics of climate change. However, when they are operated from space, lidar systems are limited by atmospheric backscattered signals that have low signal-to-noise ratios (SNRs) on optically thin targets. Therefore specific methods of analysis have to be developed to ensure accurate determination of the geometric and optical properties of scattering layers in the atmosphere. A first approach to retrieving the geometric properties of semitransparent cloud and aerosol layers is presented as a function of false-alarm and no-detection probabilities for a given SNR. Simulations show that the geometric properties of thin cirrus clouds and the altitude of the top of the unstable atmospheric boundary layer can be retrieved with standard deviations smaller than 150 m for a vertical resolution of the lidar system in the 50-100-m range and a SNR of 3. The altitudes of the top of dense clouds are retrieved with a precision in altitude of better than 50 m, as this retrieval corresponds to a higher SNR value. Such methods have an important potential application to future spaceborne lidar missions.     

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.     

Parry arc: a polarization lidar, ray-tracing, and aircraft case study

Using simple ray-tracing simulations, the cause of the rare Parry arc has been linked historically to horizontally oriented columns that display the peculiar ability to fall with a pair of prism faces closely parallel to the ground. Although we understand the aerodynamic forces that orient the long-column axis in the horizontal plane, which gives rise to the relatively common tangent arcs of the 22 degrees halo, the mechanism leading to the Parry crystal orientation has never been resolved adequately. On 16 November 1998, at the University of Utah Facility for Atmospheric Remote Sensing, we studied a cirrus cloud producing a classic upper Parry arc using polarization lidar and an aircraft with a new high-resolution ice crystal imaging probe. Scanning lidar data, which reveal extremely high linear depolarization ratios delta a few degrees off the zenith direction, are simulated with ray-tracing theory to determine the ice crystal properties that reproduce this previously unknown behavior. It is found that a limited range of thick-plate crystal axis (length-to-diameter) ratios from ~0.75 to 0.93 generates a maximum delta approximately 2.0-5.0 for vertically polarized 0.532-mum light when the lidar is tilted 1 degrees -2 degrees off the zenith. Halo simulations based on these crystal properties also generate a Parry arc. However, although such particles are abundant in the in situ data in the height interval indicated by the lidar, one still has to invoke an aerodynamic stabilization force to produce properly oriented particles. Although we speculate on a possible mechanism, further research is needed into this new explanation for the Parry arc.     

Parameterization of cloud lidar backscattering profiles by means of asymmetrical Gaussians

A fitting procedure for cloud lidar data processing is shown that is based on the computation of the first three moments of the vertical-backscattering (or -extinction) profile. Single-peak clouds or single cloud layers are approximated to asymmetrical Gaussians. The algorithm is particularly stable with respect to noise and processing errors, and it is much faster than the equivalent least-squares approach. Multilayer clouds can easily be treated as a sum of single asymmetrical Gaussian peaks. The method is suitable for cloud-shape parametrization in noisy lidar signatures (like those expected from satellite lidars). It also permits an improvement of cloud radiative-property computations that are based on huge lidar data sets for which storage and careful examination of single lidar profiles can't be carried out.     

Simulation of multiple scattering in lidar measurements using pure rotational raman scattering (PRRS)

Abstract

Pure rotational Raman scattering signals from atmospheric gases such as nitrogen and oxygen can be used to deduce the temperature of the atmosphere. Previously, this method has been successfully implemented as a remote temperature sensing lidar system. In this paper, theoretical studies of the method have been carried out using Monte Carlo simulations for different temperature profiles from radio sonde data. The geometry of the lidar as well as the aerosol profiles of the atmosphere can be specifically defined in this method. It is important to understand whether or not multiple scattering will have a significant effect on the accuracy of temperature retrieval from the measured lidar returns. From the exact pure rotational Raman scattering matrix, we have computed the lidar returns of individual Raman lines. We have given the ratios of multiple to single scattering return signals for atmospheres without clouds, with water clouds and with cirrus clouds. The results indicate that the effect of multiple scattering does not give errors to the temperature inversion for typical atmospheric conditions.     

Parameterization of aerosol and cirrus cloud effects on reflected sunlight spectra measured from space: application of the equivalence theorem

An original methodology to account for aerosol and cirrus cloud contributions to reflected sunlight is described. This method can be applied to the problem of retrieving greenhouse gases from satellite-observed data and is based on the equivalence theorem with further parameterization of the photon path-length probability density function (PPDF). Monte Carlo simulation was used to validate this parameterization for a vertically nonhomogeneous atmosphere including an aerosol layer and cirrus clouds. Initial approximation suggests that the PPDF depends on four parameters that can be interpreted as the effective cloud height, cloud relative reflectance, and two additional factors to account for photon path-length distribution under the cloud. We demonstrate that these parameters can be efficiently retrieved from the nadir radiance measured in the oxygen A-band and from the H(2)O-saturated area of the CO(2) 2.0 microm spectral band.