1. Introduction


This study has been performed in the course of development of a contrail cirrus prediction tool (CoCiP) (Schumann 2009). CoCiP predicts the IWC and the number of ice particles and the area covered by contrail cirrus. The radiative forcing by the contrail cirrus is parameterized as a function of the effective particle radius and the optical depth (Schumann et al. 2009). For this, we need to know the factor C. Knowledge about C for liquid clouds (droplets) has been derived from PSD statistics and measurements (Pontikis and Hicks 1992; Martin et al. 1994; Korolev et al. 1999). For example, Martin et al. (1994) suggested C = 0.67 for continental and 0.9 for maritime warm stratocumulus clouds. For ice clouds, some C parameterizations have been introduced for two-moment cirrus models with little theoretical or experimental support so far (Lohmann 2002; Liu et al. 2007; Ström and Gierens 2002).
The concept of the effective radius was introduced first for liquid (spherical) cloud particles by Hansen (1971) as the ratio of the third moment of the PSD to the second moment of the PSD. He used the fact that a cloud particle scatters an amount of light approximately in proportion to its area (van de Hulst 1957). Hansen and Travis (1974) discussed this radius in comparison to alternative definitions. Stephens (1978) parameterized the optical thickness of liquid clouds as the ratio between the liquid water path and the effective radius. These studies noted that the effective radius concept is valid only for short wavelength λ or large particles (i.e., for size parameters x = 2πreff/λ ≫ 1), for which the extinction efficiency Qext is close to 2, the geometric optics limit.
Cirrus cloud ice particles are composed of various types of nonspherical crystals (Weickmann 1945; Heymsfield and Iaquinta 2000; Bailey and Hallett 2009). For modeling, one often assumes prescribed particle habits (Baum et al. 2005a; Baran 2009). Recent laboratory data reveal that most ice particles at low temperatures (−40° to −70°C) have columnar shapes (with variable aspect ratio length/diameter) and that the very small crystals growing at low ice supersaturation are compact faceted polycrystals (Bailey and Hallett 2009). Often ice particles are hollow (Weickmann 1945; Schmitt and Heymsfield 2007). Ice particles may result from homogeneous freezing of initially liquid (water or solution) droplets, which then grow by deposition of water vapor. Ice particles also originate from heterogeneous freezing on preexisting ice nuclei. The different freezing processes may impact the particle habits. Once entering a subsaturated environment, ice particles will sublimate and may become more spherical again (Nelson 1998). Also, particle orientation may play a role for cloud optical properties, in particular for large ice particles (Reichardt et al. 2008; Westbrook et al. 2010).
Therefore, the concept of effective radius has been generalized to account for the nonspherical shapes (Foot 1988). Most studies now define the effective radius as proportional to the ratio between mean particle volume and mean particle projected cross-sectional area (mean shadow area) (Foot 1988; Francis et al. 1994; Fu 1996; McFarquhar and Heymsfield 1998; Wyser 1998; Hong et al. 2009). Some earlier alternative definitions (e.g., Ebert and Curry 1992) are still in use (e.g., Fusina et al. 2007).
Parameterizations of the bulk radiative properties of ice clouds for solar and terrestrial radiation are beyond the scope of this paper, but we note that they often use the concept of effective radius, both for applications in climate models and remote sensing (e.g., Ebert and Curry 1992; Fu 1996; Fu et al. 1998; Stubenrauch et al. 2004; Hong et al. 2009; Minnis et al. 1998, 2010, submitted to IEEE Trans. Geosci. Remote Sens., hereafter MIN). Various papers have discussed the usefulness and limitations of this concept (Fu 1996; Wyser 1998; Mitchell 2002; McFarquhar et al. 2003; Baran and Havemann 2004; Baran 2005; Yang et al. 2005; De Leon and Haigh 2007; Edwards et al. 2007). In many studies the optical properties are modeled as a function of ice water content, temperature, and the effective radius. Such one-moment models are independent of the volume mean radius and hence not a function of C. In principle, the effective radius can be avoided and the optical and physical properties of the clouds can be parameterized directly as a function of ice water content and temperature if crystal sizes and habits are known or can be suitably prescribed (Mitchell et al. 2008; Baran et al. 2009). In fact, many global circulation models use a prescribed effective radius, which is either given as a constant or as a function of ice water content and temperature (Liou et al. 2008). For large size parameters, the parameterizations depend also on the details of the assumed size and habit distributions. Early parameterizations were based on size distributions derived from a few field experiments. For example, Yang et al. (2000) used 30 particle size distributions as collected by Fu (1996) and Mitchell and Arnott (1994). Minnis et al. (1998) used about 10 size distributions partly from the same measurements. More recent studies (Hong et al. 2009) added further measurement cases.
Measurements of the number, area, and volume size distributions of cirrus particles as a function of particle sizes are difficult and prone to uncertainties for many reasons. Ice particles occur at sizes of about 1–8000 μm (Dowling and Radke 1990); in this paper we mainly consider the lower part of this size range. Many studies detect and size small ice particles by forward scattering particle counters (Knollenberg 1972; Baumgardner et al. 2001). Such forward scattering spectrometer probe (FSSP) instruments measure the light scattered by the particles. The signal depends on many parameters including the size and shape of the particles (Baumgardner et al. 1992). To invert the measured signals into actual particle sizes (and areas and volumes), Mie calculations or T-matrix calculations for spheres and spheroids have been used (Borrmann et al. 2000). More recent studies use the FSSP-300 for small particles, which allows identifying ice particles larger than about 0.4 μm. Early instruments underestimated the number of small ice particles. However, cirrus may contain significant concentrations of ice particles smaller than 30 μm (Baumgardner et al. 1992; Noone et al. 1993; Gayet et al. 2002; Febvre et al. 2009; Hong et al. 2009). This is the case in spite of the fact that shattering of large crystals on the probe inlet may cause an overestimate of the concentration of small ice crystals in cirrus that also contain many large ice particles (Field et al. 2003; Korolev and Isaac 2005; McFarquhar et al. 2007; Heymsfield 2007; Westbrook and Illingworth 2009). Two-dimensional optical probes are used to directly depict the shadow area of particles larger than about 25 μm (Baumgardner et al. 2001; Lawson et al. 2001). New instrument versions now allow detection of shattering by measuring the interarrival times of particles (Field et al. 2003, 2006). The particle volume for given particle sizes is estimated using empirical relationships (e.g., Baker and Lawson 2006). As an alternative to measuring the cross-sectional area of particles, one may measure the extinction coefficient β, which is directly related to the cross section of the ensemble of particles for constant extinction efficiency (Garrett et al. 2003). The polar nephelometer (Gayet et al. 1997) measures the scattering phase function of an ensemble of cloud particles for particle diameters below about 800-μm diameter. From the scattering phase function one can evaluate, among other properties, the optical extinction coefficient (Gayet et al. 2002, 2004).
Particles in contrails differ from cirrus in that they are smaller and more numerous (Schumann 1996; Gayet et al. 1996a; Baumgardner and Gandrud 1998; Sassen 1997), at least for some time. Atmospheric observations in young contrails (ages between seconds and 0.5 h) with impactors (first by Weickmann 1945) and particle size measurements (Knollenberg 1972) show ice particles composed of nearly spherical ice crystals with mean diameters in the range of 1–10 μm and narrow size distributions (Petzold et al. 1997; Poellot et al. 1999; Schröder et al. 2000), at least for cold ambient conditions (below −50°C). Contrail particles are small inside the contrail and may get larger by uptake of humidity from ambient air at the outer contrail edge (Petzold et al. 1997; Heymsfield et al. 1998; Lawson et al. 1998). The smallness of the particles causes more variability of the extinction and absorption efficiencies (Poellot et al. 1999). Febvre et al. (2009) observed optical properties of young and slightly aged (15 min) contrails and ambient cirrus. The measurements confirmed that quasispherical ice particles with diameters smaller than 5 μm control the optical properties of the plume shortly after formation. In the slightly aged contrails the optical properties are governed by larger nonspherical ice crystals. The particles with sizes smaller than 10 μm contribute about 80% to the total extinction and therefore dominate the optical properties. The shape of most of the observed ice particles of sizes smaller than 10 μm in the aged contrail strongly deviates from spherical. In aged contrail cirrus (age of 35–90 min), in particular in highly supersaturated air masses, more complex ice habits have been observed (Gayet et al. 1996a; Baumgardner and Gandrud 1998; Heymsfield et al. 1998; Lawson et al. 1998; Baumgardner et al. 2005). Heymsfield et al. (1998) and Lawson et al. (1998) find many small ice particles (diameters of 1–20 μm) but also a few large ice particles. Some of the latter could originate from mixing with ambient cirrus. Fall streaks of large ice crystals have been observed in more than 2-h aged contrail cirrus at temperatures above −40°C (Atlas et al. 2006).
Numerical simulations with the contrail–cirrus model described in Unterstrasser and Gierens (2010) show that the mean diameter evaluated over the total contrail area increases from 5 to 15–40 μm within in the first 2–3 h for a large variety of ambient conditions. Temperature has the largest impact on the crystal sizes for two reasons: the water vapor concentrations increase with temperature and more ice crystals are lost during the vortex phase in a warmer environment (Unterstrasser et al. 2008). The ambient relative humidity affects the crystal sizes to a much lesser extent, since for increasing relative humidity the two effects of increasing water vapor concentrations and of a larger number of ice crystals surviving the vortex phase nearly cancel out each other.
Often contrails show higher lidar depolarization than cirrus, indicating possibly large deviations from sphericity of the ice particles, with details depending on particle sizes, contrail age, and temperature (Freudenthaler et al. 1996; Mishchenko and Sassen 1998; Sassen and Hsueh 1998; Del Guasta and Niranjan 2001; Langford et al. 2005). Some papers list measured values of the volume mean and effective radius or diameter of cirrus and contrails (e.g., Gayet et al. 1996a,b; Schröder et al. 2000), but without discussing the importance of the ratio between the two particle scales.
So far, little is known about the transition of the habit mix from young contrails to natural cirrus. Contrail ice particles originating from freezing droplets first take a droxtal shape (Thuman and Robinson 1954; Ohtake 1970; Roth and Frohn 1998) and later may become more aspherical according to laboratory measurements (Gonda and Yamazaki 1984). The optical properties of droxtals are systematically different from those of spheres, in particular in the solar range, with less forward scattering (Yang et al. 2003; Zhang et al. 2004; Wendisch et al. 2007). After some time, contrail particles become more similar to cirrus particles. From the abovementioned observations in aged contrail cirrus, Liou et al. (1998) derived an ice crystal shape model for contrails consisting of bullet rosettes, hollow columns, and plates—hence, without near-spherical particles.
In this paper we will show that C may be larger or smaller than one and in general depends strongly on the PSD and on the particle habit. In particular, we identify upper and lower bounds for C for special cases. The value of C is deduced from measured size distributions with parameterized volume and projected area values for ice particles of given habit mixes. For contrail cirrus, a habit mix is suggested that reflects the few existing observations. Moreover, we collect measured data on number density, ice water content, and projected area or extinction coefficient from experiments in cirrus and contrails to derive a set of C values. This includes published measurements as well as new data from the recent Contrail and Cirrus Experiment (CONCERT; Voigt et al. 2010). Based on these data, the dependence of C on cirrus particle habits, volume mean radius, ambient humidity, and the age of contrails is discussed. For contrails and cirrus, a simple parameterization of the factor C suitable for two-moment models is suggested. For uncertainty estimates, approximate upper and lower bounds are provided.
2. Definitions


































Parameterizations that relate β = τ/IWP directly to the IWC and reff (Hong et al. 2009) leave the relation between reff and IWC, and hence N of ice particles, open.
3. Split into habit and size distribution factors
In this section we will split the factor C into two factors that are independent of each other, at least for an ensemble of particles with size-independent habits. This splitting will help us to understand the dependence of the value of C on the PSD and particle habits.


























For constant particle habits and aspect ratio, the two factors are independent; CPSD accounts for the PSD, while Chabit or
4. The particle size distribution factor CPSD
In this section we discuss the particle size distribution factor CPSD. We first show that 0 < CPSD ≤ 1 in general and then discuss the influence of the PSD on the value of C.
a. Lyapunov limit for CPSD for constant habits








b. CPSD for a lognormal distribution
Here we discuss the sensitivity of CPSD to the shape of the PSD np(r) for the first of three model distributions as plotted in Fig. 4, for which we can compute CPSD analytically.












In agreement with Pontikis and Hicks (1992), CPSD is directly related to the dispersion or the variance of the size distribution for the lognormal distribution. But this is not the case for other PSDs in general. Hansen and Travis (1974) and Mishchenko and Travis (1994) discuss various variance definitions, which partly scale with the fourth moment of the PSD and hence cannot generally be related to CPSD, which depends on the second and third moment of the PSD only. For general PSDs, Martin et al. (1994) find that CPSD depends on the skewness of the size distribution; CPSD gets smaller for large positive skewness (long tail at large radius). However, Martin’s result CPSD > 1 for large negative skewness contradicts the Lyapunov limit.
c. CPSD for a power-law distribution








The expression is well behaved near n = −1, −3, −4, where the above expression is indeterminate: for n = −1,
Figure 6 shows that CPSD(n, R) is a smooth function of n, CPSD(n, R) = 1 for R = 1, and CPSD(n, R) < 1 otherwise, decreasing with growing R = R2/R1, consistent with Eq. (35). For n ≤ −5 or n ≥ 0, and R → ∞, CPSD stays finite; CPSD = (−1 − n)1/3(−4 − n)2/3/(−3 − n) > (½)1/3 = 0.7937.
Between −4 < n < −1, CPSD goes to zero for large R with a flat minimum at n = −3.
For n > −3, the large particles contribute most to volume; hence rvol → R2. For n < −2, the small particles contribute most to area; hence rarea → R1. As a consequence, there exists a subrange −3 < n < −2 where the ratio rarea/rvol is less than one and CPSD decreases strongly with R. Outside this range, in the larger subrange −4 < n < −1, the decrease is slower. The asymptotic case CPSD → 0 implies reff → 0, τ → 0 for fixed IWC. Hence, a few very large particles may dominate the mass while contributing little to the optical properties.
Obviously, all values 0 < CPSD ≤ 1 are possible. However, for the limited range R < 10, CPSD varies between about 0.7 and 1.
d. CPSD for bimodal peaks
Perhaps the most simple distribution in this connection is a bimodal distribution, consisting of two peaks, i = 1, 2, where each has a monodisperse distribution represented by Dirac functions δ(Ri, ni) and n2 = n1Rn, R = R2/R1. Hence, the peaks at r = R1 and R2 have amplitudes that differ by a factor similar to the previous power law. The size distribution np(r) is zero outside the peaks (i.e., also between the two peaks; see the full lines in Fig. 4).




Hence, this distribution is similar to that for a power law, but with n replaced by n + 1. Obviously, for n > −2, the large particles contribute most to volume; hence rvol → R2. For n < −2, the small particles contribute most to area; hence rarea → R1. As a consequence, there exists a subrange −3 < n < 0 where the ratio rarea/rvol is less than one and CPSD decreases with R. Outside this range, CPSD first decreases with R and then returns to 1 for very large values of R.
Again, all values 0 < CPSD ≤ 1 are possible. However, for the realistic range limitation R < 10, CPSD varies between about 0.4 and 1. Here, the lower limit is smaller than in the power-law case.
e. CPSD values for averaged PSDs
The factor CPSD is nonlinear in the size distributions. Hence, the value of C for the average of a set of measured PSDs is different from the average value of C for the individual PSDs. The same is true for reff and rvol. Accordingly, the mean optical properties of a statistically inhomogeneous cloud depend on the ensemble over which the mean value is defined. Moreover, the following shows that a larger cloud domain in most cases has smaller C values than a smaller domain.
For illustration, we computed C for the individual PSDs plotted in Fig. 2 and for the arithmetic average of all possible combinations of two PSDs out of this set of 30 PSDs. In none of these combinations did the PSD resulting from combining two PSDs have a larger C than the maximum of the individual C values for the two PSDs. We conjecture that this statement is generally true, although we were not able to find a formal proof. On the other hand, in nearly 40% of all combinations, the averaged C value over two PSDs is less than the minimum of the C values of the two individual PSDs. Strongest reduction of C for the average PSD is obtained in the cases where the two individual PSDs had largest C (i.e., were close to monodisperse).
In fact, monodisperse PSDs (with CPSD = 1) occur rarely. Even when all particles have the same size at a specific time and position in space but these sizes vary with time and position, its mean PSD is nonmonodisperse. The bimodal spectrum discussed above may be interpreted as the mean spectrum over two subensembles with monodisperse PSDs. The results show that the CPSD < 1 for an ensemble of clouds even when subensembles may be monodisperse.
5. The habit factor Chabit

















Spheroids are rotationally symmetrical ellipsoids. The aspect ratio L/d is defined here as the ratio between the rotational semiaxis b = L/2 and the semiaxis a = d/2 perpendicular to this. Depending on L/d, the shape of spheroids ranges from spheres to needles (prolates) and plates (oblates). The volume of a spheroid is Vp = (⅙)πd2L, and the area Ap follows from known relationships (Mishchenko and Travis 1994; Krotkov et al. 1999). Spheroids are used to approximate real ice particle habits (e.g., in analysis of optical particle measurements; Borrmann et al. 2000). For spheroids, the value of C is one for L/d = 1 (sphere) and larger than one otherwise. Because of smaller surface to volume ratio, C for spheroids is less than that for cylinders except for extreme aspect ratios.


6. Evaluation of C for given PSD and habits
a. Habit model
This section provides parameterizations to compute the habit factor Chabit,h of ice crystals of given size for seven different habits (h = 0–6; see Table 1), including convex and concave ones.


The results imply habit factors Chabit as a function of r, as plotted in Fig. 9. As expected, Chabit equals one for spheres and is slightly larger than one for droxtals, and the largest habit factors are computed for rosettes and for plates.
b. Chabit for habit mixes
To account for cirrus particles of different habits, Baum et al. (2005a) suggested a habit mixture that prescribes the mass fractions of various habits as a function of maximum dimension Dmax. Their mixture assumed that all particles with Dmax < 60 μm are droxtals while larger particles have mixtures of more complex habits. (We note that the mix in the size range 1000 < Dmax < 2000 μm implies smaller habit factors for this size range than for the neighboring size ranges. But this is of minor importance for the present study, which concentrates on the size range below 1000 μm.) For contrails, the observations cited in the introduction are consistent with droxtals for small particles but imply a transition to complex particles already at smaller dimensions. For larger dimensions, contrail particles seem to become similar to other cirrus crystals. Therefore, we suggest the contrail–cirrus habit mix defined in Table 2. Since the models predict the volume mean radius, we express the mixture weights as a function of both r and Dmax (see Table 2).
Next we determine the effective radius and the habit factor Chabit for this contrail–cirrus habit mixture, with fractions








Figure 11 shows the resultant factor Chabit(r) (thick black curve). We see that Chabit(r) grows with volume mean radius from 1.06 for small particles (droxtals) to values of order 1.6 for large particles (mainly rosettes). The Chabit steps at r = 5, 9.5, and 23 μm are consequences of the assumed stepwise transition from pure droxtals to more complex habit mixes. For larger radius, the results are a consequence of the Baum cirrus mixture.


For other habit mixes, the results may differ considerably. For example, Fig. 12 shows, among others, this approximation together with the habit factor that would result for the habit mix of Liou et al. (1998). The Liou mix implies about a factor of 2 larger habit factor than suggested in this study. It contains 50% bullet rosettes, 30% hollow columns, and 20% plates, that is, particles with larger habit factors than for the near-spherical droxtals. For the reasons given in the introduction, and in agreement with Baum et al. (2005a), droxtals appear to be more representative for the smallest ice particles.
c. Empirical habit factors
The habit factor for the cirrus habit mixture can be compared with empirical relationships. Several authors have derived simple power laws relating the mass M (or diameter Deq of the droplet after melting the ice particle) to the projected area A of individual ice particles with certain habits, M/M0 = α(A/A0)β, with empirical coefficients α and β (depending on the assumed particle habit) and reference dimensions M0 and A0 (Cunningham 1978; Moss and Johnson 1994; Baker and Lawson 2006). These relationships imply a volume–area relationship and hence a habit factor. Relationships between mass and maximum particle dimension can be used for this purpose only when accompanied by area–dimension relationships.
For example, Baker and Lawson (2006) derived M/M0 = α(A/A0)β with α = 0.115, β = 1.218, M0 = 1 mg, and A0 = 1 mm2 for measured projected areas and volumes of ice particles. The data were derived from images of a set of ice particles collected on a dish at ground and of the resultant droplets after melting. This implies Chabit,BL = 5.67(r/r0)0.45, with r0 = 1 mm. Figure 12 shows that this habit factor increases with rvol as expected because of growing complexity of the large ice particles and may reach values up to an order of 10 for millimeter-sized ice particles. Obviously, this relationships should be applied only for sufficiently large r > 25 μm, for which the habit factor is larger than one. Gayet et al. (1996a) use a relationship Deq = D0(A/A0)γ, derived by Cunningham (1978) for bullet-rosettes with A < 0.21 mm2, γ = 0.39, D0 = 1 mm, and A0 = 1 mm2. This implies Chabit,G = 10.3(r/r0)0.564 (for 0.016 < r/r0 < 0.14). Gayet et al. (2004) use similar formulas with other coefficients. The empirical relationships derived from imaging particles provide larger habit factors than the habit mixtures of Baum et al. (2005a). Possibly the image method overestimates the average projected area because the crystals most likely come to rest on the dish in their most stable orientation with largest projected area. The large differences between the approximations from various sources suggest large uncertainties in these factors. In any case, large C values are to be expected for large rvol unless small CPSD values counteract this trend, which is probable because large particles often imply a wide size distribution with small CPSD.
d. Habit influence on C for given cirrus ice particle size distributions




The results for C for individual habit types (see various black symbols in Fig. 11) scatter between 0.14 and 1.52. This large scatter is no longer unexpected because of highly different size distributions and particle habits. The largest values C are computed for the size distribution number 25, which stays flat until a kink near a particle size d of 40 μm and then decreases rapidly (see Fig. 13). In this case, both rvol and reff are close to 20 μm, and C is near 1 therefore. The smallest C value is computed for the size distribution number 1, which decreases smoothly following closely a power law with exponent −3. In view of our results for power laws (Fig. 7), the small value of C was to be expected for this case.
Far less scatter results when evaluating the factor C for a given habit mixture. The blue symbols in Fig. 11 show the C result for the individual PSDs using the contrail cirrus habit mix. The C values show a decreasing trend with particle size for small particles (because of increasing width of the PSDs) and then a slight increase, presumably due to growing contributions from large irregular ice particles. These results will be used below to relate C to rvol.
7. C from experimental data
As a guide for proper C parameterizations, this section presents measured values. We first discuss results for cirrus without obvious anthropogenic influences, and then for contrails. The experimental data used are described in the appendix.
a. Experimental data for C in cirrus
The results derived from measurements in cirrus are plotted in Fig. 14. The data labeled Baum et al. (2005a) were derived similarly to those shown in Fig. 11: instead of using the contrail–cirrus habit mix, these data were computed with habit factors of Yang et al. (2000) and PSDs of Key et al. (2002) as before, but using the habit mix of Baum et al. (2005a). The differences of the results for the two habit mixes (cirrus and contrail cirrus) are small because of small number concentrations in the size range below 60 μm.
The results show C values scattering mostly between 0.4 and 1.1. Most of the data support the scatter and trends derived from cirrus mix of Baum et al. (2005a). These data show mean (min–max) values for C = 0.7(0.4−1.1) for rvol/μm = 22(5−55). An exception is the data of Gayet et al. (2006) for large radius, which show smaller values down to 0.25. PSDs with large dispersion would explain these low values, but underestimates of habit factors and particle shattering may have also contributed to these small values. In fact, the value of C gets smaller when the number N of small particles increases [see Eq. (11)]. However, even for constant volume, an increase of N may be partly balanced by an increase in A, so shattering does not always imply smaller values of C. However, shattering is most important in the presence of many large ice particles, and such cases are likely connected with small negative power-law PSD exponents n. Therefore, we expect that shattering causes a steeper PSD with smaller n. For small negative n, C gets smaller for decreasing n (see Fig. 6)—hence shattering of large particles is likely to cause a decrease of C. In the cirrus dataset of de Reus et al. (2009), incidences of shattering from the CIP instrument were removed based on particle interarrival times (Field et al. 2006), and this dataset gives results for C consistent with the other datasets. Although we cannot exclude particle shattering, we do not expect that shattering has large effects on the values of C on average, at least for the range of rvol discussed in this paper.
b. Experimental data for C in contrails
Because contrail ice crystals are of different origin, have small mean size with narrow PSD, and are nearly spherical, the value of C may be different for contrail cirrus compared to other cirrus, particularly at small ages. Only a few measurements are available to evaluate C for contrails. In the appendix we collect the existing information and add new results from recent in situ measurements in contrails during CONCERT (Voigt et al. 2010). The results are plotted versus rvol in Fig. 15 and versus plume age in Fig. 16. Both figures use logarithmic abscissa to cover the highly variable sizes and ages.
As a whole, Fig. 15 shows C values with considerable scatter from 0.3 to 1.2. The values are closer to one for small particles and the data generally indicate a decrease of C with rvol, similar to the cirrus cases. The smallest volume mean radius in these data is near 0.33–0.5 μm (Petzold et al. 1997; Schröder et al. 2000), causing significant deviations from Qext = 2. For small contrail ice particles, the C values are close to one because of narrow size distributions (Petzold et al. 1997). Smaller C values arise because of widening size distributions. The C values of Schröder et al. (2000) and Febvre et al. (2009) in contrails are smaller than those obtained by the same teams in ambient cirrus, presumably because of larger habit factors for the larger cirrus particles. For larger sizes, some of the large C values in contrails may be caused by nearly monodisperse size distributions and strongly nonspherical particles (e.g., in warm and humid atmospheres) (Gayet et al. 1996a). However, large C values may be also a consequence of a large habit factor assumed for large particles (see Fig. 12). The C values show no systematic dependence on the aircraft type (Voigt et al. 2010) or the fuel sulfur content (Petzold et al. 1997). In general, the C values are within the range found in Fig. 14 for cirrus.
Figure 16 shows the dependency on plume age. For fresh contrails (<10 s), C is in between 0.7 and 0.9. Such rather large values are expected for small spectral dispersion and near-spherical particle shape (droxtals). For larger plume ages, the data exhibit large scatter (around 0.7 ± 0.4) without obvious systematic trends. With growing plume age, some ice particles grow in size and the spectral dispersion and the habits may approach common cirrus-cloud values (Schröder et al. 2000). Some of the contrails were reported as “evaporating” (Schröder et al. 2000). In these cases, sublimation may cause a rounding of ice particles (Nelson 1998) and hence a reduction of C.
c. Variation of C in cirrus and contrails with relative humidity
For a few cases (Schröder et al. 2000; de Reus et al. 2009; Voigt et al. 2010) we have simultaneous measurements of PSDs together with relative humidity over ice RHi. RHi was measured in situ with different frost point hygrometers (Schröder et al. 2000; de Reus et al. 2009) or a fluorescence hygrometer combined with temperature measurements (Voigt et al. 2010). The cited publications report estimated uncertainties of RHi of about 15% (Voigt et al. 2010; for details, see the references cited). Hence, we can investigate the impact of humidity on C (see Fig. 17). Of course, the measured RHi is the humidity at the time of ice particle measurement, not the humidity during cloud formation. This may contribute to the scatter of the data. The data of de Reus et al. (2009) were obtained mostly at very low temperatures and therefore may be well comparable to contrail cases.
Contrails form by condensing water mainly emitted by the aircraft engines; they later take up water from ambient air in ice supersaturated air masses. In subsaturated air masses the particles sublimate. Contrails formed by heavy four-engine aircrafts can survive for several minutes even for considerable subsaturation since large amounts of water vapor are emitted. Contrails of smaller two-engine aircrafts can evaporate after a few seconds (Sussmann and Gierens 2001). For supersaturated air masses, the results indicate a wide range of C, without clear indications for a dependence on humidity. For subsaturated cases, the CONCERT C values seem to increase with decreasing humidity. The humidity measurements by Schröder et al. (2000) show rather low humidity (down to 50% for contrail ages up to 120 s, and 60% in ambient cirrus) and large scatter for C.
The data suggest that for high ambient humidity, the PSD widens (causing reduced PSD factors). In subsaturated air masses the measured PSD contains less particle volume in the larger size bins, causing smaller volume mean radius. The effective radius gets also reduced, perhaps even more quickly, because of more efficient sublimation of particles with large surface. In addition we expect a rounding of the shape of the sublimating particles, although this cannot be detected with the instruments used. The variability of humidity obviously explains a large part of the scatter of the data seen in Figs. 14 and 15.
8. Parameterization of C
a. Parameterizations used in previous models
The need for parameterization of C has been acknowledged by all teams who use two-moment cirrus models. Models that assume a lognormal PSD with fixed dispersion (e.g., Spichtinger and Gierens 2009) effectively prescribe the value of C with this assumption [see Eq. (41)]. Other models use parameterizations as plotted in Fig. 18.
Lohmann (2002) includes a prognostic equation for the ice particle number concentration in a global climate model and prescribes C as a function of a given temperature-dependent effective radius. Lohmann et al. (2008) essentially use the same relationship but with C as a function of the volume mean radius that follows from the computed IWC and ice particle number concentration (see Fig. 18). This function implies a monotonic decrease of C with rvol that may be reasonable for spherical particles but does not account for more complex particle habits. Hence, this model may underestimate the optical depth of thick natural cirrus clouds. The underestimate could be even stronger when shattering causes an underestimate of C in the data.
Ström and Gierens (2002) study the impact of aircraft particle emissions on the radiative properties of contrails. They use a prognostic model for IWC and the number density of ice crystals and relate the effective radius to the mean ice particle mass, which can be expressed as a function of rvol. For large rvol, their result is similar to the Lohmann model. At small rvol, the tendency to large C values is unrealistic, as we now know. It overemphasizes the climate impact of very small ice particles.
Liu et al. (2007, 2009) and Penner et al. (2009) parameterize C as a function of temperature and IWC based on data derived by Wyser (1998). We have plotted their C value [with the coefficients as given in Penner et al. (2009)] for a median midlatitude IWC–temperature correlation (Schiller et al. 2008). Observed IWC values show large scatter relative to this median. Therefore, we also show curves for the same IWC multiplied with 0.1 or 10. The C parameterization gives reasonable values in the range as observed but tends to overshoot above one for small IWC values and low temperatures. This could cause an overestimate of the radiative impact of thin and subvisible cirrus and an underestimate of the impact of thick cirrus. In any case, these models do not account for changes in the particle habits.
b. New parameterization




9. Conclusions
We have discussed the relationship C = (rarea/rvol)2 = rvol/reff between the volume-specific radius and the effective radius. For ice particles with constant extinction efficiency, the optical depth is determined when knowing the IWP, the number concentration of ice particles, and the factor C. For small ice particles, as in young contrails, or for thermal wavelengths, the deviation of the effective extinction efficiency from the geometric optics limit has to be taken into account (Key et al. 2002; Yang et al. 2005; Hong et al. 2009). The ratio C is a purely geometrical parameter and depends in general strongly on the PSD and on the particle habits; C can be expressed as the product of a habit and a PSD-related factor—C = CPSDChabit. We have shown strictly that 0 < CPSD ≤ 1 for arbitrary PSD in general and equal to 1 for a monodisperse PSD. Moreover, Chabit ≥ 1 for convex and concave particles with random orientation in general and is equal to 1 for spherical particles.
In general, C may get very small for PSDs that decay with size similar to a power law with exponent n between −4 and −1 (−3 and 0 for bimodal spectra). This has been derived analytically for simple model PSDs. For such wide and slowly decaying PSDs, most of the particle volume (and ice water content) is controlled by a few large particles, while most of the cross-sectional area (and optical depth) is controlled by the many small particles. The same is true for lognormal size distributions with large dispersion. Since C is a nonlinear function of the size distributions, the mean optical properties of a statistically inhomogeneous cloud depend also on the ensemble over which the mean value is defined. A larger cloud domain tends to have smaller average C values than a smaller domain.
We provide an analysis of PSD and habit data of Yang et al. (2000) and Key et al. (2002). The analysis provides specific values of C as a function of particle sizes and habits. In addition, we analyze the value of C from previous and recent in situ measurements. In the various measured cases, for rvol below about 100 μm, CPSD varies typically between values of order 0.4 and 1, and Chabit ranges between 1 and about 2. Because of growing aspect ratios, Chabit grows with the particle size. The modeled habit factor is somewhat smaller than assumed in data analysis of particle measurements.
For contrail cirrus, a new particle habit mix is suggested. This habit mix includes near-spherical (droxtal)-shaped particles for the very small particles and allows for more complex shapes already at rather small maximum dimension. In the size range of contrails crystals, the habit factor for this mix is about a factor of 2 smaller than for the habit mix of Liou et al. (1998).
For contrails and cirrus we suggest a new parameterization to estimate C in two-moment models as a function of IWC, number of ice particles, and possibly relative humidity. For parameter studies reasonable upper and lower bounds are estimated. For fresh contrails, C = 0.7 ± 0.3 because of near-spherical particles and PSDs with small dispersion. Any uncertainty in C has direct impact on the accuracy of the computed optical depth of the cloud.
In the new parameterization, C is decreasing less quickly with rvol than in older parameterizations (Lohmann et al. 2008; Penner et al. 2009). This increases the relative contribution of large particles to optical depth. As a consequence, the climate impact of freshly nucleated small ice crystals may be smaller than estimated before.
For further progress, one would need more basic understanding of the processes and parameters controlling the ice habits and particle size distributions.
Acknowledgments
This work was performed within the DLR project “Climate-Compatible Air Transport System” (CATS). The CONCERT campaign was organized by the HGF junior research group “Impact of Aircraft Emissions on the Heterogeneous Chemistry of the Tropopause Region” (AEROTROP). For computation of the extinction efficiency of spheroids, T-matrix results from Josef Gasteiger were used. Philip Kotter helped in the formulation of the mathematical properties of convex and concave bodies. Stephan Borrmann made the data from de Reus et al. (2009) available to us. Discussions with Pat Minnis and Joyce Penner during the “Aviation Climate Change Research Initiative” (ACCRI) helped our understanding.
REFERENCES
Atlas, D., Z. Wang, and D. P. Duda, 2006: Contrails to cirrus—Morphology, microphysics, and radiative properties. J. Appl. Meteor. Climatol., 45 , 5–19.
Bailey, M. P., and J. Hallett, 2009: A comprehensive habit diagram for atmospheric ice crystals: Confirmation from the laboratory, AIRS II, and other field studies. J. Atmos. Sci., 66 , 2888–2899.
Baker, B., and R. P. Lawson, 2006: Improvement in determination of ice water content from two-dimensional particle imagery. Part I: Image-to-mass relationships. J. Appl. Meteor. Climatol., 45 , 1282–1290.
Baran, A. J., 2005: The dependence of cirrus infrared radiative properties on ice crystal geometry and shape of the size-distribution function. Quart. J. Roy. Meteor. Soc., 131 , 1129–1142.
Baran, A. J., 2009: A review of the light scattering properties of cirrus. J. Quant. Spectrosc. Radiat. Transfer, 110 , 1239–1260.
Baran, A. J., and S. Havemann, 2004: The dependence of retrieved cirrus ice-crystal effective dimension on assumed ice-crystal geometry and size-distribution function at solar wavelengths. Quart. J. Roy. Meteor. Soc., 130 , 2153–2167.
Baran, A. J., P. J. Connolly, and C. Lee, 2009: Testing an ensemble model of cirrus ice crystals using midlatitude in situ estimates of ice water content, volume extinction coefficient and the total solar optical depth. J. Quant. Spectrosc. Radiat. Transfer, 110 , 1579–1598.
Baum, B. A., A. J. Heymsfield, P. Yang, and S. T. Bedka, 2005a: Bulk scattering properties for the remote sensing of ice clouds. Part I: Microphysical data and models. J. Appl. Meteor., 44 , 1885–1895.
Baum, B. A., P. Yang, A. J. Heymsfield, S. Platnick, M. D. King, Y-X. Hu, and S. T. Bedka, 2005b: Bulk scattering properties for the remote sensing of ice clouds. Part II: Narrowband models. J. Appl. Meteor., 44 , 1896–1911.
Baumgardner, D., and B. E. Gandrud, 1998: A comparison of the microphysical and optical properties of particles in an aircraft contrail and mountain wave cloud. Geophys. Res. Lett., 25 , 1129–1132.
Baumgardner, D., J. E. Dye, B. W. Gandrud, and R. G. Knollenberg, 1992: Interpretation of measurements made by the Forward Scattering Spectrometer Probe (FSSP-300) during the Airborne Arctic Stratospheric Expedition. J. Geophys. Res., 97 , 8035–8046.
Baumgardner, D., H. Jonsson, W. Dawson, D. O’Connor, and R. Newton, 2001: The cloud, aerosol and precipitation spectrometer: A new instrument for cloud investigations. Atmos. Res., 59–60 , 251–264.
Baumgardner, D., H. Chepfer, G. B. Raga, and G. L. Kok, 2005: The shapes of very small cirrus particles derived from in situ measurements. Geophys. Res. Lett., 32 , L01806. doi:10.1029/2004GL021300.
Borrmann, S., B. Luo, and M. Mishchenko, 2000: Application of the T-matrix method to the measurement of aspherical (ellipsoidal) particles with forward scattering optical particle counters. J. Aerosol Sci., 31 , 789–799.
Cauchy, A., 1908: Oeuvres Complètes d’Augustin Cauchy. Series 1, Vol. 2, Gauthier-Villars, 167–177.
Chylek, P., 1977: Extinction cross sections of arbitrarily shaped randomly oriented nonspherical particles. J. Opt. Soc. Amer., 67 , 1348–1350.
Cunningham, R. M., 1978: Analysis of particle spectral data from optical array (PMS) 1D and 2D sensors. Preprints, Fourth Symp. on Meteorological Observations and Instrumentation, Denver, CO, Amer. Meteor. Soc., 345–350.
De Leon, R. R., and J. D. Haigh, 2007: Infrared properties of cirrus clouds in climate models. Quart. J. Roy. Meteor. Soc., 133 , 273–282.
Del Guasta, M., and K. Niranjan, 2001: Observation of low-depolarization contrails at Florence (Italy) using a 532–1064-nm polarization lidar. Geophys. Res. Lett., 28 , 4067–4070.
de Reus, M., and Coauthors, 2009: Evidence for ice particles in the tropical stratosphere from in-situ measurements. Atmos. Chem. Phys., 9 , 6775–6792.
Dowling, D. R., and L. F. Radke, 1990: A summary of the physical properties of cirrus clouds. J. Appl. Meteor., 29 , 970–978.
Ebert, E. E., and J. A. Curry, 1992: A parameterization of ice cloud optical properties for climate models. J. Geophys. Res., 97 , 3831–3836.
Edwards, J., S. Havemann, J-C. Thelen, and A. Baran, 2007: A new parametrization for the radiative properties of ice crystals: Comparison with existing schemes and impact in a GCM. Atmos. Res., 83 , 19–35.
Febvre, G., and Coauthors, 2009: On optical and microphysical characteristics of contrails and cirrus. J. Geophys. Res., 114 , D02204. doi:10.1029/2008JD010184.
Field, P. R., R. Wood, P. R. A. Brown, P. H. Kaye, E. Hirst, R. Greenway, and J. A. Smith, 2003: Ice particle interarrival times measured with a fast FSSP. J. Atmos. Oceanic Technol., 20 , 249–261.
Field, P. R., A. J. Heymsfield, and A. Bansemer, 2006: Shattering and particle interarrival times measured by optical array probes in ice clouds. J. Atmos. Oceanic Technol., 23 , 1357–1371.
Fisz, M., 1970: Wahrscheinlichkeitsrechnung und mathematische Statistik. VEB Deutscher Verlag der Wissenschaften, 528 pp.
Foot, J. S., 1988: Some observations of the optical properties of clouds. II: Cirrus. Quart. J. Roy. Meteor. Soc., 114 , 145–164.
Francis, P. N., A. Jones, R. W. Saunders, K. P. Shine, A. Slingo, and Z. Sun, 1994: An observational and theoretical study of the radiative properties of cirrus: Some results from ICE’89. Quart. J. Roy. Meteor. Soc., 120 , 809–848.
Freudenthaler, V., F. Homburg, and H. Jäger, 1996: Optical parameters of contrails from lidar measurements: Linear depolarization. Geophys. Res. Lett., 23 , 3715–3718.
Fu, Q., 1996: An accurate parameterization of the solar radiative properties of cirrus clouds for climate models. J. Climate, 9 , 2058–2082.
Fu, Q., and K. N. Liou, 1993: Parameterization of the radiative properties of cirrus clouds. J. Atmos. Sci., 50 , 2008–2025.
Fu, Q., P. Yang, and W. B. Sun, 1998: An accurate parameterization of the infrared radiative properties of cirrus clouds for climate models. J. Climate, 11 , 2223–2237.
Fusina, F., P. Spichtinger, and U. Lohmann, 2007: Impact of ice supersaturated regions and thin cirrus on radiation in the midlatitudes. J. Geophys. Res., 112 , D24S14. doi:10.1029/2007JD008449.
Garrett, T. J., H. Gerber, D. G. Baumgardner, C. H. Twohy, and E. M. Weinstock, 2003: Small, highly reflective ice crystals in low-latitude cirrus. Geophys. Res. Lett., 30 , 2132. doi:10.1029/2003GL018153.
Gayet, J-F., G. Febvre, G. Brogniez, H. Chepfer, W. Renger, and P. Wendling, 1996a: Microphysical and optical properties of cirrus and contrails: Cloud field study on 13 October 1989. J. Atmos. Sci., 53 , 126–138.
Gayet, J-F., G. Febvre, and H. Larsen, 1996b: The reliability of the PMS FSSP in the presence of small ice crystals. J. Atmos. Oceanic Technol., 13 , 1300–1310.
Gayet, J-F., O. Crépel, J-F. Fournol, and S. Oshchepkov, 1997: A new airborne polar nephelometer for the measurements of optical and microphysical cloud properties. Part I: Theoretical design. Ann. Geophys., 15 , 451–459.