Ice Particle Habit Research Articles

Characteristics of Generating Cells in Wintertime Orographic Clouds

Abstract During the Seeded and Natural Orographic Wintertime clouds: the Idaho Experiment (SNOWIE) field campaign, cloud-top generating cells were frequently observed in the very high-resolution W-band airborne cloud radar data. This study examines multiple flight segments from three SNOWIE cases that exhibited cloud-top generating cells structures, focusing on the in situ measurements inside and outside these cells to characterize the microphysics of these cells. The observed generating cells in these three cases occurred in cloud tops of −15° to −30°C, with and without overlying cloud layers, but always with shallow layers of atmospheric instability observed at cloud top. The results also indicate that liquid water content, vertical velocity, and drizzle and ice crystal concentrations are greater inside the generating cells compared to the adjacent portions of the cloud. The generating cells were predominantly <500 m in horizontal width and frequently exhibited drizzle drops coexisting with ice. The particle imagery indicates that ice particle habits included plates, columns, and rimed and irregular crystals, likely formed via primary ice nucleation mechanisms. Understanding the sources of natural ice formation is important to understanding precipitation formation in winter orographic clouds, and is especially relevant for clouds that may be targeted for glaciogenic cloud seeding as well as to improve model representation of these clouds. Significance Statement This study presents the characteristics of cloud-top generating cells in winter orographic clouds, and documents that fine-scale generating cells are ubiquitous in clouds over complex terrain in addition to having been observed in other types of clouds. The generating cells exhibited enhanced concentrations of larger drizzle and ice particles, which suggests the environments of these fine-scale features promote ice formation and growth. The source of ice formation in winter clouds is critical to understanding and modeling the precipitation formation process. Given the ubiquity of cloud-top generating cells in many types of clouds around the world, this study further motivates the need to investigate methods for representing subgrid-scale environments to improve ice formation in numerical models.

Effects of Inhomogeneous Ice Particle Habit Distribution on Passive Microwave Radiative Transfer Simulations

Effects of Inhomogeneous Ice Particle Habit Distribution on Passive Microwave Radiative Transfer Simulations

Comparison of Cloud Properties between SGLI Aboard GCOM-C Satellite and MODIS Aboard Terra Satellite

This study presents a comprehensive comparison of Level 2.0 cloud properties between a Second-generation Global Imager (SGLI) aboard the GCOM-C satellite and a Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the Terra satellite, to better understand the qualities of cloud properties obtained from SGLI/GCOM-C launched on 23 December 2017. The cloud pixels identified as water phase by both satellite sensors are highly consistent to each other by more than 90%, although the consistency is only ~60% for ice phase cloud pixels. A comparison of cloud properties—cloud optical thickness (COT) and cloud particle effective radius (CER)—between these two satellite sensors reveals that water and ice cloud properties can have different degrees of agreement depending on underlying surface. The relative difference (RD) values of 22% (18%) and 37% (24%) for water cloud COT (CER) comparison over ocean and land surfaces and respective values of 35% (42%) and 35% (62%) for comparisons of ice cloud properties, and also other comparison metrics, suggest better agreements for water cloud properties than for ice cloud properties, and for ocean surface than for land surface. Though cloud properties differences between MODIS and SGLI can arise from inherent features of cloud retrieval algorithms, such as differences in ancillary data, surface reflectance, cloud droplet size distribution function, model for ice particle habit, etc., this study further identifies the important roles of cloud thickness and Sun and satellite positions for differences in cloud properties between SGLI and MODIS: the differences in cloud properties are found to increase for thinner clouds, higher solar zenith angle, and higher differences in viewing zenith and azimuth angles between these satellite sensors, and such differences are more distinct for water cloud properties than for ice cloud properties.

Retrieval and Evaluation of Ice Water Content from the Airborne Wyoming Cloud Radar in Orographic Wintertime Clouds during SNOWIE

Abstract As part of the analysis following the Seeded and Natural Orographic Wintertime Storms (SNOWIE) project, the ice water content (IWC) in ice and mixed-phase clouds is retrieved from airborne Wyoming Cloud Radar (WCR) measurements aboard the University of Wyoming King Air (UWKA), which has a suite of integrated in situ IWC, optical array probes, and remote sensing measurements, and it provides a unique dataset for this algorithm development and evaluation. A sensitivity study with different idealized ice particle habits shows that the retrieved IWC with aggregate ice particle habit agrees the best with the in situ measurement, especially in ice or ice-dominated mixed-phase clouds with a correlation coefficient (rr) of 0.91 and a bias of close to 0. For mixed-phase clouds with ice fraction ratio less than 0.8, the variances of IWC estimates increase (rr = 0.76) and the retrieved mean IWC is larger than in situ IWC by a factor of 2. This is found to be related to the uncertainty of in situ measurements, the large cloud inhomogeneity, and the retrieval assumption uncertainty. The simulated reflectivity Ze and IWC relationships assuming three idealized ice particle habits and measured particle size distributions show that hexagonal columns with the same Ze have a lower IWC than aggregates, whose Ze–IWC relation is more consistent with the observed WCR Ze and in situ IWC relation in those clouds. The 2D stereo probe (2DS) images also indicate that ice particle habit transition occurs in orographic mixed-phase clouds; hence, the retrieved IWC assuming modified gamma particle size distribution (PSD) of aggregate particles tends to have a greater bias in this kind of clouds.

Sensitivity of radiative flux simulations to ice cloud parameterization over the equatorial western Pacific Ocean region

AbstractPrevious studies suggest explanations of the observed cancellation of shortwave (SW) and longwave (LW) cloud radiative effects (CREs) at the top of the atmosphere (TOA) over tropical oceans where deep convection prevails, such as interactions among cloud microphysics, radiation, and dynamics. However, simulations based on general circulation models (GCMs) show disagreement in terms of the net (SW + LW) CREs over tropical deep convective ocean regions. One of the GCM uncertainty sources is the parameterization of ice cloud bulk optical properties. In this study, a combination of active and passive satellite daytime cloud retrievals is used to study the sensitivity of radiation flux calculations to ice cloud parameterization over the equatorial western Pacific Ocean region. Three ice cloud schemes are tested. The first is a widely used scheme that assumes hexagonal column ice particles. The second scheme treats ice particles as aggregates of surface-roughened hexagonal columns. The third scheme best matches the cloud ice mass-dimension relation in the cloud microphysics scheme by assuming a mixture of two ice particle habits. The results show that the hexagonal-column-based scheme has the weakest SW CRE but strongest LW CRE among the three. In addition, cloud optical thickness and effective radius are used to cluster cold-top single-layer ice clouds into three types, which resemble thin cirrus, detrained anvil clouds, and deep convective cores, respectively. In agreement with previous studies, cloud SW heating overwhelms LW cooling in the upper portion of anvil-like clouds.

Influence of low-level blocking and turbulence on the microphysics of a mixed-phase cloud in an inner-Alpine valley

Abstract. Previous studies that investigated orographic precipitation have primarily focused on isolated mountain barriers. Here we investigate the influence of low-level blocking and shear-induced turbulence on the cloud microphysics and precipitation formation in a complex inner-Alpine valley. The analysis focuses on a mid-level cloud in a post-frontal environment and a low-level feeder cloud induced by an in-valley circulation. Observations were obtained from an extensive set of instruments including ground-based remote sensing instrumentation, in situ instrumentation on a tethered-balloon system and ground-based precipitation measurements. During this event, the boundary layer was characterized by a blocked low-level flow and enhanced turbulence in the region of strong vertical wind shear at the boundary between the blocked layer in the valley and the stronger cross-barrier flow aloft. Cloud radar observations indicated changes in the microphysical cloud properties within the turbulent shear layer including enhanced linear depolarization ratio (i.e., change in particle shape or density) and increased radar reflectivity (i.e., enhanced ice growth). Based on the ice particle habits observed at the surface, we suggest that riming, aggregation and needle growth occurred within the turbulent layer. Collisions of fragile ice crystals (e.g., dendrites, needles) and the Hallett–Mossop process might have contributed to secondary ice production. Additionally, in situ instrumentation on the tethered-balloon system observed the presence of a low-level feeder cloud above a small-scale topographic feature, which dissipated when the low-level flow turned from a blocked to an unblocked state. Our observations indicate that the low-level blocking (due to the downstream mountain barrier) created an in-valley circulation, which led to the production of local updrafts and the formation of a low-level feeder cloud. Although the feeder cloud did not enhance precipitation in this particular case (since the majority of the precipitation sublimated when falling through a subsaturated layer above), we propose that local flow effects such as low-level blocking can induce the formation of feeder clouds in mountain valleys and on the leeward slope of foothills upstream of the main mountain barrier, where they can act to enhance orographic precipitation through the seeder–feeder mechanism.

Interpretation of lidar ratio and depolarization ratio of ice clouds using spaceborne high-spectral-resolution polarization lidar.

The backscattering coefficient (β), lidar ratio (S), and depolarization ratio (δ) of ice particles were estimated over a wide range of effective radii to interpret spaceborne 355-nm high-spectral-resolution lidar data from the ATLID sensor onboard the EarthCARE satellite. Five randomly oriented ice particle shapes (3D ice) and two quasi-horizontally oriented particle types (2D ice) were analyzed using five effective angles. The size dependence of β, S, and δ was examined using physical optics and geometrical optics integral equation methods. Differences in β for the same effective radius and ice water content among particle types exceeded one order of magnitude. S-δ relations are useful for inferring ice particle habit and orientation using ATLID data from EarthCARE.

Ice Cloud Properties From Himawari-8/AHI Next-Generation Geostationary Satellite: Capability of the AHI to Monitor the DC Cloud Generation Process

The Japan Meteorological Agency (JMA) successfully launched the Himawari-8 (H-8) new-generation geostationary meteorological satellite with the Advanced Himawari Imager (AHI) sensor on October 7, 2014. The H-8/AHI level-2 (L2) operational cloud property products were released by the Japan Aerospace Exploration Agency during September 2016. The Voronoi light scattering model, which is a fractal ice particle habit, was utilized to develop the retrieval algorithm called “Comprehensive Analysis Program for Cloud Optical Measurement” (CAPCOM-INV)-ice for the AHI ice cloud product. In this paper, we describe the CAPCOM-INV-ice algorithm for ice cloud products from AHI data. To investigate its retrieval performance, retrieval results were compared with 2000 samples of the ice cloud optical thickness and effective particle radius values. Furthermore, AHI ice cloud products are evaluated by comparing them with the MODIS collection-6 (C6) products. As an experiment, cloud property retrievals from AHI measurements, with an observation interval time of 2.5 min and ground-based rainfall observation radar data (the latter of which is supplied by the JMA, with a 1-km grid mesh), are used to investigate the generation processes of deep convective (DC) cloud in the vicinity of the Kyushu island, Japan. It revealed that AHI measurements have the capability of monitoring the growth processes, including variation of the cloud properties and the precipitation in the DC cloud.

A Versatile Method for Ice Particle Habit Classification Using Airborne Imaging Probe Data

AbstractA versatile method to automatically classify ice particle habit from various airborne optical array probes is presented. The classification is achieved using a multinomial logistic regression model. For each airborne probe, the model determines the particle habit (among six classes) based on a large set of geometrical and textural descriptors extracted from the two‐dimensional image of a particle. The technique is applied and evaluated using three probes with significantly different specifications: the high volume precipitation spectrometer, the two‐dimensional stereo probe, and the cloud particle imager. Performance and robustness of the method are assessed using standard machine learning tools on the basis of thousands of images manually labeled for each of the considered probes. The three classifiers show good performance characterized by overall accuracies and Heidke skill scores above 90%. Depending on the application and user preferences, the classification scheme can be easily adapted. For a more precise output, intraclass subclassification can be achieved in a nested fashion, illustrated here with columnar crystals and aggregates. A comparative study of the classification output obtained with the three probes is presented for two aircraft flight periods selected when the three probes were operating together. Results are globally consistent in term of proportions of habit identified (once blurry and partial images have been automatically discarded). A perfect agreement is not expected as the three considered probes are sensitive to different particle size range.

Simulation of the brightness temperatures observed by the visible infrared imaging radiometer suite instrument (Erratum)

Clouds play a large role in the Earth’s global energy budget, but the impact of cirrus clouds is still widely questioned and researched. Cirrus clouds reside high in the atmosphere and due to cold temperatures are comprised of ice crystals. Gaining a better understanding of ice cloud optical properties and the distribution of cirrus clouds provides an explanation for the contribution of cirrus clouds to the global energy budget. Using radiative transfer models (RTMs), accurate simulations of cirrus clouds can enhance the understanding of the global energy budget as well as improve the use of global climate models. A newer, faster RTM such as the visible infrared imaging radiometer suite (VIIRS) fast radiative transfer model (VFRTM) is compared to a rigorous RTM such as the line-by-line radiative transfer model plus the discrete ordinates radiative transfer program. By comparing brightness temperature (BT) simulations from both models, the accuracy of the VFRTM can be obtained. This study shows root-mean-square error <0.2 K for BT difference using reanalysis data for atmospheric profiles and updated ice particle habit information from the moderate-resolution imaging spectroradiometer collection 6. At a higher resolution, the simulated results of the VFRTM are compared to the observations of VIIRS resulting in a <1.5 % error from the VFRTM for all cases. The VFRTM is validated and is an appropriate RTM to use for global cloud retrievals.

Simulation of the brightness temperatures observed by the visible infrared imaging radiometer suite instrument

Clouds play a large role in the Earth’s global energy budget, but the impact of cirrus clouds is still widely questioned and researched. Cirrus clouds reside high in the atmosphere and due to cold temperatures are comprised of ice crystals. Gaining a better understanding of ice cloud optical properties and the distribution of cirrus clouds provides an explanation for the contribution of cirrus clouds to the global energy budget. Using radiative transfer models (RTMs), accurate simulations of cirrus clouds can enhance the understanding of the global energy budget as well as improve the use of global climate models. A newer, faster RTM such as the visible infrared imaging radiometer suite (VIIRS) fast radiative transfer model (VFRTM) is compared to a rigorous RTM such as the line-by-line radiative transfer model plus the discrete ordinates radiative transfer program. By comparing brightness temperature (BT) simulations from both models, the accuracy of the VFRTM can be obtained. This study shows root-mean-square error <0.2 K for BT difference using reanalysis data for atmospheric profiles and updated ice particle habit information from the moderate-resolution imaging spectroradiometer collection 6. At a higher resolution, the simulated results of the VFRTM are compared to the observations of VIIRS resulting in a <1.5 % error from the VFRTM for all cases. The VFRTM is validated and is an appropriate RTM to use for global cloud retrievals.

Investigation of ice particle habits to be used for ice cloud remote sensing for the GCOM-C satellite mission

Abstract. In this study, various ice particle habits are investigated in conjunction with inferring the optical properties of ice clouds for use in the Global Change Observation Mission-Climate (GCOM-C) satellite programme. We develop a database of the single-scattering properties of five ice habit models: plates, columns, droxtals, bullet rosettes, and Voronoi. The database is based on the specification of the Second Generation Global Imager (SGLI) sensor on board the GCOM-C satellite, which is scheduled to be launched in 2017 by the Japan Aerospace Exploration Agency. A combination of the finite-difference time-domain method, the geometric optics integral equation technique, and the geometric optics method is applied to compute the single-scattering properties of the selected ice particle habits at 36 wavelengths, from the visible to the infrared spectral regions. This covers the SGLI channels for the size parameter, which is defined as a single-particle radius of an equivalent volume sphere, ranging between 6 and 9000 µm. The database includes the extinction efficiency, absorption efficiency, average geometrical cross section, single-scattering albedo, asymmetry factor, size parameter of a volume-equivalent sphere, maximum distance from the centre of mass, particle volume, and six nonzero elements of the scattering phase matrix. The characteristics of calculated extinction efficiency, single-scattering albedo, and asymmetry factor of the five ice particle habits are compared. Furthermore, size-integrated bulk scattering properties for the five ice particle habit models are calculated from the single-scattering database and microphysical data. Using the five ice particle habit models, the optical thickness and spherical albedo of ice clouds are retrieved from the Polarization and Directionality of the Earth's Reflectances-3 (POLDER-3) measurements, recorded on board the Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar (PARASOL) satellite. The optimal ice particle habit for retrieving the SGLI ice cloud properties is investigated by adopting the spherical albedo difference (SAD) method. It is found that the SAD is distributed stably due to the scattering angle increases for bullet rosettes with an effective diameter (Deff) of 10 µm and Voronoi particles with Deff values of 10, 60, and 100 µm. It is confirmed that the SAD of small bullet-rosette particles and all sizes of Voronoi particles has a low angular dependence, indicating that a combination of the bullet-rosette and Voronoi models is sufficient for retrieval of the ice cloud's spherical albedo and optical thickness as effective habit models for the SGLI sensor. Finally, SAD analysis based on the Voronoi habit model with moderate particle size (Deff = 60 µm) is compared with the conventional general habit mixture model, inhomogeneous hexagonal monocrystal model, five-plate aggregate model, and ensemble ice particle model. The Voronoi habit model is found to have an effect similar to that found in some conventional models for the retrieval of ice cloud properties from space-borne radiometric observations.

PHIPS–HALO: the airborne Particle Habit Imaging and Polar Scattering probe – Part 1: Design and operation

Abstract. The number and shape of ice crystals present in mixed-phase and ice clouds influence the radiation properties, precipitation occurrence and lifetime of these clouds. Since clouds play a major role in the climate system, influencing the energy budget by scattering sunlight and absorbing heat radiation from the earth, it is necessary to investigate the optical and microphysical properties of cloud particles particularly in situ. The relationship between the microphysics and the single scattering properties of cloud particles is usually obtained by modelling the optical scattering properties from in situ measurements of ice crystal size distributions. The measured size distribution and the assumed particle shape might be erroneous in case of non-spherical ice particles. There is a demand to obtain both information correspondently and simultaneously for individual cloud particles in their natural environment. For evaluating the average scattering phase function as a function of ice particle habit and crystal complexity, in situ measurements are required. To this end we have developed a novel airborne optical sensor (PHIPS-HALO) to measure the optical properties and the corresponding microphysical parameters of individual cloud particles simultaneously. PHIPS-HALO has been tested in the AIDA cloud simulation chamber and deployed in mountain stations as well as research aircraft (HALO and Polar 6). It is a successive version of the laboratory prototype instrument PHIPS-AIDA. In this paper we present the detailed design of PHIPS-HALO, including the detection mechanism, optical design, mechanical construction and aerodynamic characterization.

Ice microphysics retrieval in the convective systems of the Indian Ocean during the CINDY–DYNAMO campaign

The last decades have shown the relevance of polarimetric radar measurements in various domains of Earth science to characterize different types of particles in natural media such as snow covers, vegetation canopy and clouds. In atmospheric physics, the need to accurately define the frozen hydrometeor properties (type, size, density) in precipitating systems is directly related to our understanding of the convective processes. These microphysical properties are also particularly critical in setting ice parameterization in models.Ground-based polarimetric radars can provide information on the most likely type of ice particles present in a sampled volume. These classifications can be used to build ice particle habit statistics once they are validated.The SPolKa classification scheme (PID, Particle IDentificator) for various types of frozen hydrometeors is compared to in-situ measurements collected during the CINDY–DYNAMO campaign in the Indian Ocean. The French Falcon-20 flew inside the polarimetric radar area measuring the particle habits with several in-situ microphysical probes (FSSP, 2DS, CPI, PIP, 2DP, Nevzorov). These in-situ data are organized as catalog of images and matched to the classifications proposed by Magono and Lee (1966) and Kikuchi et al. (2013). The PID classification and the Precipitation Imaging Probe (PIP) data images are compared over two sequences representing two different atmospheric situations: a large stratiform area sampled on November 27th, 2011 and some convective activity embedded in stratiform precipitation on December 8th, 2011.The general agreement is very good for most species. Since the PIP only provides us with a 2D image of the particles' shadow, some species are difficult to identify unambiguously and could be matched with more than one PID. This is particularly critical for species with complex history.

Dynamical and Microphysical Evolution during Mixed-Phase Cloud Glaciation Simulated Using the Bulk Adaptive Habit Prediction Model

Abstract A bulk microphysics scheme predicting ice particle habit evolution has been implemented in the Weather Research and Forecasting Model. Large-eddy simulations are analyzed to study the effects of ice habit and number concentration on the bulk ice and liquid masses, dynamics, and lifetime of Arctic mixed-phase boundary layer clouds. The microphysical and dynamical evolution simulated using the adaptive habit scheme is compared with that assuming spherical particles with a density of bulk ice or a reduced density and with mass–dimensional parameterizations. It is found that the adaptive habit method returns an increased (decreased) ice (liquid) mass as compared to spheres and provides a more accurate simulation as compared to dendrite mass–size relations. Using the adaptive habit method, simulations are then completed to understand the microphysical and dynamical interactions within a single-layer mixed-phase stratocumulus cloud observed during flight 31 of the Indirect and Semi-Direct Aerosol Campaign. With cloud-top longwave radiative cooling as a function of liquid mass acting as the primary dynamic driver of turbulent eddies within these clouds, the consumption of liquid at the expense of ice growth and subsequent sedimentation holds a strong control on the cloud lifetime. Ice concentrations ≥ 4 L−1 collapse the liquid layer without any external maintaining sources. Layer maintenance is possible at 4 L−1 when a constant cloud-top cooling rate or the water mass lost due to sedimentation is supplied. Larger concentrations require a more substantial source of latent or sensible heat for mixed-phase persistence.

The Three-Dimensional Morphology of Simulated and Observed Convective Storms over Southern England

AbstractA set of high-resolution radar observations of convective storms has been collected to evaluate such storms in the Met Office Unified Model during the Dynamical and Microphysical Evolution of Convective Storms (DYMECS) project. The 3-GHz Chilbolton Advanced Meteorological Radar was set up with a scan-scheduling algorithm to automatically track convective storms identified in real time from the operational rainfall radar network. More than 1000 storm observations gathered over 15 days in 2011 and 2012 are used to evaluate the model under various synoptic conditions supporting convection. In terms of the detailed three-dimensional morphology, storms in the 1500-m grid length simulations are shown to produce horizontal structures a factor of 1.5–2 wider compared to radar observations. A set of nested model runs at grid lengths down to 100 m show that the models converge in terms of storm width, but the storm structures in the simulations with the smallest grid lengths are too narrow and too intense compared to the radar observations. The modeled storms were surrounded by a region of drizzle without ice reflectivities above 0 dBZ aloft, which was related to the dominance of ice crystals and was improved by allowing only aggregates as an ice particle habit. Simulations with graupel outperformed the standard configuration for heavy-rain profiles, but the storm structures were a factor of 2 too wide and the convective cores 2 km too deep.

Ice particle habit and surface roughness derived from PARASOL polarization measurements

Abstract. Ice clouds are an important element in the radiative balance of the earth's climate system, but their microphysical and optical properties still are not well constrained, especially ice particle habit and the degree of particle surface roughness. In situ observations have revealed common ice particle habits and evidence for surface roughness, but these observations are limited. An alternative is to infer the ice particle shape and surface roughness from satellite observations of polarized reflectivity since they are sensitive to both particle shape and degree of surface roughness. In this study an adding–doubling radiative transfer code is used to simulate polarized reflectivity for nine different ice habits and one habit mixture, along with 17 distinct levels of the surface roughness. A lookup table (LUT) is constructed from the simulation results and used to infer shape and surface roughness from PARASOL satellite polarized reflectivity data over the ocean. Globally, the retrievals yield a compact aggregate of columns as the most commonly retrieved ice habit. Analysis of PARASOL data from the tropics results in slightly more aggregates than in midlatitude or polar regions. Some level of surface roughness is inferred in nearly 70% of PARASOL data, with mean and median roughness near σ = 0.2 and 0.15, respectively. Tropical region analyses have 20% more pixels retrieved with particle surface roughness than in midlatitude or polar regions. The global asymmetry parameter inferred at a wavelength of 0.865 μm has a mean value of 0.77 and a median value of 0.75.

Aircraft Observations of Ice Particle Properties in Stratiform Precipitating Clouds

This study presented airborne measurements of ice particle properties in three stratiform precipitating clouds over northern China. By using horizontal observations at selected altitudes, the distributions of ice water content (IWC), particle habits, and particle size spectrum parameters were investigated. The cloud cases were characterized by high IWC values due to the existence of embedded convective cells. Liquid water contents were rather low with the maxima of less than 0.3 g m−3and the general values of less than 0.1 g m−3. The occurrence of large dendritic crystals as well as rimed capped columns and branched crystals indicated that ice seeding from the above cloud layer (6 km altitude or above) contributed significantly to both high ice crystal number concentrations and IWCs. Horizontal observations at selected levels suggested the general decreasing trend of IWC with decreasing temperature only in part of the cloud layers but not throughout the cold layer of the multilayered stratiform clouds. Both exponential and gamma functions were used to characterize the particle size spectrum parameters. The slope parameter values of exponential distributions were primarily in the range of 103–104 m−1. In comparison, slope values of the gamma distribution fits spanned more and a relationship was found between the dispersion and slope values.

The evolution of microphysical and optical properties of an A380 contrail in the vortex phase

Abstract. A contrail from a large-body A380 aircraft at cruise in the humid upper troposphere has been probed with in-situ instruments onboard the DLR research aircraft Falcon. The contrail was sampled during 700 s measurement time at contrail ages of about 1–4 min. The contrail was in the vortex regime during which the primary wake vortices were sinking 270 m below the A380 flight level while the secondary wake remained above. Contrail properties were sampled separately in the primary wake at 90 and 115 s contrail age and nearly continously in the secondary wake at contrail ages from 70 s to 220 s. The scattering phase functions of the contrail particles were measured with a polar nephelometer. The asymmetry parameter derived from these data is used to distinguish between quasi-spherical and aspherical ice particles. In the primary wake, quasi-spherical ice particles were found with concentrations up to 160 cm−3, mean effective diameter Deff of 3.7 μm, maximum extinction of 7.0 km−1, and ice water content (IWC) of 3 mg m−3 at slightly ice-subsaturated conditions. The secondary and primary wakes were separated by an almost particle-free wake vortex gap. The secondary wake contained clearly aspherical contrail ice particles with mean Deff of 4.8 μm, mean (maximum) concentration, extinction, and IWC of 80 (350) cm−3, 1.6 (5.0) km−1, and 2.5 (10) mg m−3, respectively, at conditions apparently above ice-saturation. The asymmetry parameter in the secondary wake decreased with contrail age from 0.87 to 0.80 on average indicating a preferential aspherical ice crystal growth. A retrieval of ice particle habit and size with an inversion code shows that the number fraction of aspherical ice crystals increased from 2% initially to 56% at 4 min contrail age. The observed crystal size and habit differences in the primary and secondary wakes of an up to 4 min old contrail are of interest for understanding ice crystal growth in contrails and their climate impact. Aspherical contrail ice particles cause less radiative forcing than spherical ones.

The impact of ice particle roughness on the scattering phase matrix

The goal of this study is to explore the influence of ice particle habit (or shape) and surface roughness on the scattering phase matrix. As an example, reported here are the results for two wavelengths: 0.67 and 1.61 μm. For this effort, a database of single-scattering properties has been computed for a set of habits including hexagonal plates, hollow and solid columns, hollow and solid 3D bullet rosettes, droxtals, aggregates of solid columns, and aggregates of plates. The database provides properties for each of the habits at 101 wavelengths between 0.45 and 2.24 μm for smooth, moderately roughened, and severely roughened particles. At each wavelength, the scattering properties are provided at 233 discrete particle diameters ranging from 2 to 10,000 μm. A single particle size distribution from a very cold ice cloud sampled during the CRYSTAL-FACE field campaign ( T cld=–76 °C) is used to illustrate the influence of habit and roughness on the phase matrix. In all, four different habit mixtures are evaluated. The nonzero elements of the phase matrix are shown to be quite sensitive to the assumed habit, particularly in the case of − P 12/ P 11 that is associated with the degree of linear polarization of scattered radiation. Surface roughness is shown to smooth out maxima in the scattering phase function and in the other elements of the phase matrix, consistent with other studies. To compare with the theoretical simulations of the phase matrix for smooth and roughened particles, a full year of cloud-aerosol lidar with orthogonal polarization (CALIOP) data from 2008 is analyzed to provide global statistics on the values of P 11 and P 22/ P 11 in the backscattering direction. In a comparison of two of the habit mixtures (one used for MODIS Collection 5 and another that incorporates new habits including hollow bullet rosettes and aggregates of plates) with the CALIOP data, the values for P 11 are higher regardless of the degree of particle surface roughness, and the values for P 22/ P 11 are lower than those for CALIOP. Further investigation is warranted to better understand this discrepancy.