Results from 22 airborne field campaigns, including more than 10 million high‐resolution particle images collected in cirrus formed in situ and in convective anvils, are interpreted in terms of particle shapes and their potential impact on radiative transfer. Emphasis is placed on characterizing ice particle shapes in tropical maritime and midlatitude continental anvil cirrus, as well as in cirrus formed in situ in the upper troposphere, and subvisible cirrus in the upper tropical troposphere layer. There is a distinctive difference in cirrus ice particle shapes formed in situ compared to those in anvils that are generated in close proximity to convection. More than half the mass in cirrus formed in situ are rosette shapes (polycrystals and bullet rosettes). Cirrus formed from fresh convective anvils is mostly devoid of rosette‐shaped particles. However, small frozen drops may experience regrowth downwind of an aged anvil in a regime with RH_ice > ~120% and then grow into rosette shapes. Identifiable particle shapes in tropical maritime anvils that have not been impacted by continental influences typically contain mostly single plate‐like and columnar crystals and aggregates. Midlatitude continental anvils contain single‐rimed particles, more and larger aggregates with riming, and chains of small ice particles when in a highly electrified environment. The particles in subvisible cirrus are < ~100 μm and quasi‐spherical with some plates and rare trigonal shapes. Percentages of particle shapes and power laws relating mean particle area and mass to dimension are provided to improve parameterization of remote retrievals and numerical simulations.

There is a long-standing challenge in cloud and climate models to simulate the process of ice particle riming realistically, partly due to the unrealistic parameterization of the growth of ice particle mass (m) and projected area (A) during riming. This study addresses this problem, utilizing ground-based measurements of m and ice particle maximum dimension (D) as well as theory to formulate simple expressions describing the dependence of m and A on riming. It was observed that β in the m − D power law m = α D^β appears independent of riming during the phase 1 (before the formation of graupel), with α accounting for the ice particle mass increase due to riming. This semi-empirical approach accounts for the degree of riming and renders a gradual and smooth ice particle growth process from unrimed ice particles to graupel, and thus avoids discontinuities in m and A during accretional growth. Once the graupel with quasi-spherical shape forms, D increases with an increase in m and A (phase 2 of riming). The treatment for riming is explicit, and includes the parameterization of the ice crystal–cloud droplet collision efficiency (E_c) for hexagonal columns and plates using hydrodynamic theory. In particular, E_c for cloud droplet diameters less than 10 µm are estimated, and under some conditions observed in mixed-phase clouds, these droplets can account for roughly half of the mass growth rate from riming. These physically meaningful yet simple methods can be used in models to improve the riming process.

A new satellite remote sensing method is described whereby the sensitivity of thermal infrared wave resonance absorption to small ice crystals is exploited to estimate cirrus cloud ice-particle number concentration N, effective diameter D_e and ice water content IWC. This method uses co-located observations from the Infrared Imaging Radiometer (IIR) and from the CALIOP (Cloud and Aerosol Lidar with Orthogonal Polarization) lidar aboard the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) polar orbiting satellite, employing IIR channels at 10.6 and 12.05 µm. Using particle size distributions measured over many flights of the TC4 (Tropical Composition, Cloud and Climate Coupling) and the mid-latitude SPARTICUS (Small Particles in Cirrus) field campaigns, we show for the first time that N∕IWC is tightly related to β_eff; the ratio of effective absorption optical depths at 12.05 and 10.6 µm. Relationships developed from in situ aircraft measurements are applied to β_eff derived from IIR measurements to retrieve N. This satellite remote sensing method is constrained by measurements of β_eff from the IIR and is by essence sensitive to the smallest ice crystals. Retrieval uncertainties are discussed, including uncertainties related to in situ measurement of small ice crystals (D<15 µm), which are studied through comparisons with IIR β_eff. The method is applied here to single-layered semi-transparent clouds having a visible optical depth between about 0.3 and 3, where cloud base temperature is ≤235 K. CALIPSO data taken over 2 years have been analyzed for the years 2008 and 2013, with the dependence of cirrus cloud N and D_e on altitude, temperature, latitude, season (winter vs. summer) and topography (land vs. ocean) described. The results for the mid-latitudes show a considerable dependence on season. In the high latitudes, N tends to be highest and D_e smallest, whereas the opposite is true for the tropics. The frequency of occurrence of these relatively thick cirrus clouds exhibited a strong seasonal dependence in the high latitudes, with the occurrence frequency during Arctic winter being at least twice that of any other season. Processes that could potentially explain some of these micro- and macroscopic cloud phenomena are discussed.

Ice particle mass- and projected area-dimension (m-D and A-D) power laws are commonly used in the treatment of ice cloud microphysical and optical properties and the remote sensing of ice cloud properties. Although there has long been evidence that a single m-D or A-D power law is often not valid over all ice particle sizes, few studies have addressed this fact. This study develops self-consistent m-D and A-D expressions that are not power laws but can easily be reduced to power laws for the ice particle size (maximum dimension or D) range of interest, and they are valid over a much larger D range than power laws. This was done by combining ground measurements of individual ice particle m and D formed at temperature T < −20 °C during a cloud seeding field campaign with 2-D stereo (2D-S) and cloud particle imager (CPI) probe measurements of D and A, and estimates of m, in synoptic and anvil ice clouds at similar temperatures. The resulting m-D and A-D expressions are functions of temperature and cloud type (synoptic vs. anvil), and are in good agreement with m-D power laws developed from recent field studies considering the same temperature range (−60 °C < T < −20 °C).

Observed stratocumulus to cumulus transitions (SCTs) and their sensitivity to aerosols are studied using a large‐eddy simulation (LES) model that simulates the aerosol lifecycle, including aerosol sources and sinks. To initialize, force, and evaluate the LES, we used a combination of reanalysis, satellite, and aircraft data from the 2015 Cloud System Evolution in the Trades field campaign over the Northeast Pacific. The simulations follow two Lagrangian trajectories from initially overcast stratocumulus (Sc) to the tropical shallow cumulus region near Hawaii. The first trajectory is characterized by an initially clean, well‐mixed Sc‐topped marine boundary layer (MBL), then continuous MBL deepening and precipitation onset followed by a clear SCT and a consistent reduction of aerosols that ultimately leads to an ultra‐clean layer in the upper MBL. The second trajectory is characterized by an initially polluted and decoupled MBL, weak precipitation, and a late SCT. Overall, the LES simulates the observed general MBL features. Sensitivity studies with different aerosol initial and boundary conditions reveal aerosol‐induced changes in the transition, and albedo changes are decomposed into the Twomey effect and adjustments of cloud liquid water path and cloud fraction. Impacts on precipitation play a key role in the sensitivity to aerosols: for the first case, runs with enhanced aerosols exhibit distinct changes in microphysics and macrophysics such as enhanced cloud droplet number concentration, reduced precipitation, and delayed SCT. Cloud adjustments are dominant in this case. For the second case, enhancing aerosols does not affect cloud macrophysical properties significantly, and the Twomey effect dominates.

This paper describes a new approach for representing ice microphysics in climate models. In contrast with most previous schemes, this approach does not include separate categories for cloud and precipitating ice and instead uses a single two-moment category to represent all solid hydrometeors. Thus, there is no need for an ice “autoconversion” size threshold parameter, which has a critical impact on simulated climate in the Community Atmosphere Model (CAM5) yet is poorly constrained by theory or observations. Further, in the new treatment, all ice microphysical processes and parameters, including ice effective radius and mean fall speed, are formulated self-consistently and flexibly based on empirical ice particle mass–size and projected area–size relationships. This means that the scheme can represent the physical coupling between bulk particle density, mean fall speed, and effective radius, which is not possible in current schemes. Two different methods for specifying these relationships based on observations are proposed. The new scheme is tested in global simulations using CAM5. Differences in simulations using the two methods for specifying the mass– and projected area–size relationships, particularly the cloud radiative forcing, are attributable mainly to the effects on mean ice particle fall speed, impacting sedimentation and ice water path. With some tuning of parameters involved in calculating homogeneous freezing it produces a similar climate compared to the simulations using the original CAM5 microphysics. Thus, it can produce a comparable climate while improving the physical basis and self-consistency of ice particle properties and parameters.