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.     

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.     

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.     

Laser transmission through thin cirrus clouds

A near-infrared airborne-laser transmission model for thin cirrus clouds has been developed on the basis of the successive-order-of-scattering approach to account for multiple scattering by randomly and horizontally oriented ice crystals associated with an aircraft-target system. Direct transmission and transmission due to multiple scattering are formulated specifically for this geometric system, in which scattering and absorption associated with aerosols, water vapor, and air are accounted for. A number of sensitivity experiments have been performed for investigation of the effect of aircraft-target position, cirrus cloud optical depth, and ice crystal size on laser transmission for tactical applications. We show that transmission contributions produced by orders of scattering higher than 1 are small and can be neglected. The possibility of horizontal orientation of ice crystals can enhance transmission of laser beams in the aircraft-target geometry. Transmitted energy is strongly dependent on the horizontal distance between the aircraft and the target and on the cloud optical depth as well as on whether the cloud is above or below the aircraft.     

Optical properties of contrail-induced cirrus: discussion of unusual halo phenomena

Photographs of a 120 degrees parhelion and a 22 degrees parhelion within persistent contrails are presented. These phenomena result from hexagonal plate-shaped ice crystals oriented horizontally with diameters between 300 mum and 2 mm. From our observations and reinvestigation of previous reports, we conclude that a subset of the population in persistent contrails can consist of highly regular, oriented, hexagonal plates or columns comparable to the most regular crystals in natural cirrus clouds. This is explained by measured ambient humidities below the formation conditions of natural cirrus. The resulting strong azimuthal variability of the scattering phase function impacts the radiative transfer through persistent contrails.     

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.     

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.     

Identification of Odd-Radius Halo Arcs and of 44 degrees /46 degrees Parhelia by Their Inner-Edge Polarization

Birefringence of ice causes the inner edges of refraction halos to be polarized. The direction of this polarization relates directly to the projection of the crystal main axis onto the sky. This implies that the inner-edge polarization can serve as an observational diagnostic for determining the actual nature of a halo arc if two competing explanations exist. The direction and the visibility of the inner-edge polarization of arcs and circular halos arising from usual ice crystals and from ice crystals with pyramidal ends are calculated. It is found that the observation of inner-edge polarization can be decisive for the identification of a spot that might be either a 44 degrees parhelion or a 46 degrees parhelion, of an arc that might be either a 22 degrees sunvex Parry arc or a 20 degrees Parroid arc arising from plate-oriented pyramidal crystals, and of an arc that might be either a 22 degrees suncave Parry arc or a 23 degrees Parroid arc from plate-oriented pyramidal crystals. (With a Parroid arc, a halo is that which arises from an ice wedge made up of two faces of a crystal that rotates about a vertically oriented spin axis, and the edge of the wedge is perpendicular to this spin axis.) Polarization properties of other rare arcs are discussed. Practical hints are given for observing visually the inner-edge polarization of halos.     

Approach to photorealistic halo simulations

A multiple-scattering Monte Carlo model that can produce near-photographic quality images is developed and used to simulate several dramatic halo displays. The model atmosphere contains an absorbing ozone layer plus two clear, molecular air layers with Rayleigh scattering surrounding a cloud layer and an atmospheric boundary layer with aerosol particles subject to Lorentz-Mie scattering. Halos are produced by right hexagonal or pyramidal crystals that reflect and refract according to geometric optics without diffraction, although "junk" crystals with a pronounced forward-scattering peak but no halo peaks may be included to simulate typical, faint halos. Model parameters include ozone height and content, surface and cloud pressure, cloud optical thickness, crystal shapes, orientations and abundances, atmospheric turbidity, aerosol radius, and albedo. Beams for each wavelength are sorted into small bins as halo beams if they have been scattered once only by a single crystal and otherwise as sky beams, which are smoothed and combined with the halo beams to produce images. Multiple scattering generally vitiates halos, but extremely rare halos, such as Kern's arc, can be produced if a significant fraction of crystals in optically thick clouds have identical shapes and are highly oriented. Albedo is a model by-product with potential value in climate studies.     

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.     

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.     

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.     

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.     

Image transfer through cirrus clouds. I. Ray trace analysis and wave-front reconstruction

A new technique for modeling image transfer through cirrus clouds is presented. The technique uses a ray trace to model beam propagation through a three-dimensional volume of polydisperse, hexagonal ice crystals. Beyond the cloud, the technique makes use of standard Huygens-Fresnel propagation methods. At the air-cloud interface, each wave front is resolved into a ray distribution for input to the ray trace software. Similarly, a wave front is reconstructed from the output ray distribution at the cloud-air interface. Simulation output from the ray trace program is presented and the modulation transfer function for stars imaged through cirrus clouds of varying depths is discussed.     

Halo Observations Provide Evidence of Airborne Cubic Ice in the Earth's Atmosphere

An ice crystal halo display that contains several previously unknown halo phenomena was observed in Northern Chile. Analysis of computer simulations of the halos demonstrates that most of the new halo arcs in the display can be explained by the presence of airborne and preferentially oriented crystals of cubic ice. These observations therefore provide evidence of the existence of the cubic phase of ice in the Earth's atmosphere.     

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.     

Supernumerary ice-crystal halos?

Geometric-optics singularities in the intensity profiles of refraction halos formed by randomly oriented ice crystals are softened by diffraction and decorated with fine supernumerary fringes. If the crystals have a fixed symmetry axis (as in parhelia), the geometric singularity is a square-root divergence, as in the rainbow. However, the universal curve that describes diffraction is different from the rainbow's Airy function, with weak maxima (supernumerary fringes) on the geometrically dark region inside the halo (and even fainter fringes outside); these are much smaller than their counterparts on the light side of rainbows. If the crystals have no preferred orientation (as in the 22° halo), the geometric singularity is a step. In this case the universal diffraction function has no maxima, and its supernumeraries are shoulders rather than maxima. The low contrast of the fringes is probably the main reason why supernumerary halos are rarely if ever seen.     

Symmetry in halo displays and symmetry in halo-making crystals

The relation between the symmetry in halo displays and crystal symmetry is investigated for halo displays that are generated by ensembles of crystals. It is found that, regardless of the symmetry of the constituent crystals, such displays are always left-right (L-R) symmetric if the crystals are formed from the surrounding vapor. L-R symmetry of a halo display implies here that the cross sections for formation of a halo arc on the left-hand side of the solar vertical and its right-hand side mirror image are equal. This property leaves room for two types of halo display only: a full symmetric one (mmm-symmetric), and a partial symmetric one (mm2-symmetric) in which halo constituents lack their counterparts on the other side of the parhelic circle. A partial symmetric display can occur only for point halos. Its occurrence implies that a number of symmetry elements are not present in the shape of the halo-making crystals. These elements are a center of inversion, any rotatory-inversion axis that is parallel to the crystal spin axis P, a mirror plane perpendicular to the P axis, and a twofold rotation axis perpendicular to the P axis. A simple conceptual method is presented to reconstruct possible shapes of the halo-generating crystals from the halos in the display. The method is illustrated in two examples. Halos that may occur on the Saturnian satellite Titan are discussed. The possibilities for the Huygens probe to detect these halos during its descent through the Titan clouds in 2005 are detailed.