A new approach is described for calculating the mass (m) and terminal velocity (Vt) of ice particles from airborne and balloon-borne imaging probe data as well as its applications for remote sensing and modeling studies. Unlike past studies that derived these parameters from the maximum (projected) dimension (D) and habit alone, the “two-parameter approach” uses D and the particle's projected cross-sectional area (A). Expressions were developed that relate the area ratio (Ar; the projected area of an ice particle normalized by the area of a circle with diameter D) to its effective density (ρe) and to Vt.
Habit-dependent, power-law relationships between ρe and Ar were developed using analytic representations of the geometry of various types of planar and spatial ice crystals. Relationships were also derived from new or reanalyzed data for single ice particles and aggregates observed in clouds and at the ground.
The mass relationships were evaluated by comparing calculations to direct measurements of ice water content (IWC). The calculations were from Particle Measuring Systems (PMS) 2D-C and 2D-P probes of particle size distributions in ice cloud layers on 3 days during an Atmospheric Radiation Measurement (ARM) field campaign in Oklahoma; the direct measurements were from counterflow virtual impactor (CVI) observations in ice cloud layers during the field campaign. Agreement was generally to within 20%, whereas using previous mass–dimension relationship approaches usually produced larger differences. Comparison of ground-based measurements of radar reflectivity with calculations from collocated balloon-borne ice crystal measurements also showed that the new method accurately captured the vertical reflectivity structure. Improvements in the accuracy of the estimates from the earlier mass–dimension relationships were achieved by converting them to the new form. A new, more accurate mass–dimension relationship for spatial, cirrus-type crystals was deduced from the comparison.
The relationship between Vt and Ar was derived from a combination of theory and observations. A new expression accounting for the drag coefficients of large aggregates was developed from observational data. Explicit relationships for calculating Vt as a function of D for aggregates with a variety of component crystals were developed.
The development of cloud models employing sophisticated schemes for prognosing cloud ice water content (IWC), and the planned deployment of radar in space to map ice cloud properties, have provided a need to accurately characterize the properties of individual and populations of ice particles. In this study, our aim is to characterize the properties of ice particles found in midlatitude and tropical ice clouds, research important to both climate-modeling and remote sensing studies.
Climate-modeling studies require parameterizations of bulk properties—for example, mass-weighted fall speed or extinction coefficient—as a function of a prognosed variable such as the IWC. These variables have been calculated from the measured or parameterized particle size distribution by several researchers, using dimensional relationships of the form m = σDβ, A = γDκ, and Vt = aDb to represent the functional dependence, where the coefficients and exponents in the equations depended on habit. For most habits, two or three sets of coefficients were necessary to cover the full range of observed D (Heymsfield 1972; Mitchell 1996, hereafter M96).
Heymsfield and Iaquinta (2000, hereafter HI00) pointed out that there were substantial questions regarding the generality of previous m versus D and Vt versus D relationships. The relationships were subject to error because they were usually derived from different datasets. For example, A might have been taken from in situ data, whereas m had been taken from surface collections. Furthermore, most estimates of m came from observations at the ground, not from within the cloud, where the particle characteristics might have depended upon height relative to cloud top and base, temperature, mode of origin, and a host of other factors.
Given the difficulty in measuring both m and A for individual particles collected within a cirrus cloud, a geometric model was developed (HI00) to derive the three-dimensional structure of bullet rosettes and columns, which are primary cirrus crystal habits. Calculations from the model facilitated the development of a set of m versus D and Vt versus D relationships for rosettes with from one to eight bullets and for columns. The HI00 method was not general, however, in that the aspect ratios (α)—the width/length of bullets and columns—were prescribed from observations for a small number of particles observed by a balloon-borne ice crystal replicator in very cold cirrus. A more general way would have been to use the particle images more directly in some way to prescribe m and Vt.
In this study, we improve upon the HI00 study, by developing a general method for deriving m and Vt for a number of ice particle habits using both D and A. It was logical to use A and D together instead of D alone, as these two parameters are needed to derive several primary ice cloud microphysical and radiative properties: the extinction coefficient and the effective radius, among others. Furthermore, most particle spectrometers (e.g., 2D imaging probes) record both A and D. The technology to measure m from in situ measurements is not yet available; hence, some physically based method for estimating ice particle m is needed.
In the next section, we develop an approach for deriving m from two parameters, D and A, for a wide variety of ice particle shapes. The approach described uses analytic descriptions of ice particle shapes, in situ aircraft measurements, and observations at the ground and in the laboratory. In section 3, we extend this approach to the calculation of Vt, and in section 4 assess the accuracy of the two-parameter approach using in situ and radar data. Section 5 is a summary of the results and a discussion of some of the implications of the study.
2. Development of area ratio–effective density relationships
One measure that was used to quantify how much ice is contained in an ice particle is its effective density ρe. The ρe is defined as the mass divided by the volume of a circumscribed sphere. Another measure of how much ice is contained in an ice particle is the area ratio, which is A projected onto a horizontal plane, such as a photodiode array of an imaging probe, compared to the area of a circumscribed circle. Most imaging probes on aircraft have been oriented such that they measure A directly. Our goal was to examine whether increasing Ar coincided with increased ρe. In this section, we examine the relationship between Ar and ρe for a variety of ice particle types. In the first subsection, we present the results for simplified ice particle shapes; in the second subsection, the results from an analytical model; and in the third subsection, the results from laboratory and field observations.
The emphasis in this study was on the types of crystals found in cirrus clouds, including anvils. In midlatitude anvils formed by strong convection and even in deep, synoptically generated cirrus, aggregates dominate the large end of the size spectrum (Kajikawa and Heymsfield 1989). Tropical-anvil ice particle data that we amassed in the Tropical Rainfall Measuring Mission (TRMM) field programs from February to September 1999 have indicated that planar crystals (plates, stellars, and the like) and aggregates (of planar crystals and a variety of cirruslike crystals) are prevalent in anvils. Because planar crystals usually grow at temperatures primarily between −10° and −20°C, they must be transported to the upper troposphere in updrafts or regions of strong mesoscale vertical motions and then carried into the anvils. Hence, planar crystals and aggregates were also considered. In the following discussion, the Magono and Lee (1966) habit classification scheme was used. The symbols and acronyms are defined in appendix A.
a. Results for idealized particle shapes
Following are the results of calculations of ρe versus Ar relationships for three simple particle types: hexagonal columns (types C1e and C1f, solid and hollow columns, respectively); crossed columns (type C2b), representing a cluster of attached columns (simplified bullet rosettes); and planar crystals (types P1a to P1e and the like).
With reference to Fig. 1, the Ar for columns (which for sizes above 50 μm generally fall with their long axis normal to their fall direction) is
where αC is the aspect ratio for columns. The Ar for columns depends uniquely on the aspect ratio, and for typical αC of 0.25–1, Ar ranged from 0.33 to 0.64, the latter being the Ar for αC = 1. For αC > 1, the particle may have been a thick plate (type C1h), falling with the hexagonal face horizontal. In Fig. 1, the maximum projected dimension of the column is D, but most researchers have used l instead of D to represent the maximum dimension.
For this discussion, a rosette is considered to consist of four coplanar hexagonal columns, each column having a length 1/2 and width w. The assumption that the column elements are in the same plane means that none of the column elements overlapped any other element and the length of each element was the same as the projected length. The Ar was
where αB = w/(l/2). For αB = 1, ArB reached 0.76.
A planar crystal, whose dimension across its basal face was d, had an Ar given by
The Ar for a hexagonal plate was 0.83. The aspect ratio was given by αP = t/d, where t was thickness.
The ρe for a simple column was found to be
where ρb, the bulk density (m divided by crystal volume), was either 0.91 g cm−3 for solid columns or some lower value for hollow columns. The ratio 2.6/4 in Eq. (4) for a hexagon is slightly lower than the ratio π/4 for a circle. Equation (1) was used to simplify Eq. (4):
The ρe for a simplified rosette was
which, analogous to Eq. (6), was
For a planar crystal,
The ρe for the three simplified shapes are shown as a function of Ar in Fig. 2a, as derived using Eq. (5) with and without the power-law simplification. As an approximation, ρb was taken to be 0.81 g cm−3 for columns and bullet rosettes and 0.91 g cm−3 for hexagonal plates, based on Heymsfield (1972), HI00, and others. For planar crystals, αP was taken to be 0.2 and 0.05, representative values from Auer and Veal (1970).
For columns and crossed columns, ρe ≈ ∝ A2.1r.
For planar crystals, typically αP < 0.2, therefore ρeP was linearly proportional to Ar for a given αP. The ρeP were lower than ρeC and ρeB and were governed by t; t cannot be readily ascertained from imaging probe data unless the probe is oriented perpendicular to the normal sampling orientation and has high resolution (although d is not measured in that case).
We compared these calculations to observations reported in the literature of columns and bullet rosettes (M96 and HI00). Relationships between m and D and between A and D were observed from which ρe and Ar were derived. However, these comparisons were approximations because D in these studies differed somewhat from how we defined it for this paper. The resulting relationships between ρe and Ar depended on D; these are plotted in Fig. 2b. Points along each of the curves in the figure correspond to D of 200, 500, and 1000 μm; increasing ρe along each curve corresponds to decreasing particle size. The general shape of these curves compared favorably with the calculations. Mitchell's ρe curve for columns has a shallower slope than the calculations. The M96 curve for five bullets compared favorably with the crossed-column model at the lower Ar values (larger sizes) but fell below the model at the higher Ar values (smaller sizes). Although the curves may have provided reasonable estimates of ρe within the range of Ar typically observed, fixed relationships between α and D were required, and the approach may not have produced ρe versus Ar relationships that were general.
Plots of ρeP as a function of Ar based on aspect ratio expressions in Pruppacher and Klett (1978) for types P1b (crystal with sector like branches), P1c (broad branches), and P1d (stellars) are also presented in Fig. 2b. Points along the curves represent particle sizes. The ρeP decreased with increasing D, which comparison with Fig. 2a indicates was a result of the decrease in the αP with D. The Ar for type P1b particles changed only slightly with D, presumably because these crystals did not branch appreciably. Conversely, the ρeP for type P1c and P1d particles dropped more steeply with Ar than the planar crystals, shown in Fig. 2a, indicates that the αP decreased with size. The combination of the two planar crystal curves in Fig. 2b suggests that there may be a general relationship between Ar and ρeP if Ar and αP are related in some way.
The most important conclusion to be drawn from Fig. 2b is that the m versus D relationships reported previously may not be sufficiently general, because we do not know whether what is being looked at is a dependence of ρe on D or a relationship between ρe and Ar.
b. Modeling results for rosettes
Increasing evidence from observations with new instruments that provide high-resolution images of ice particles suggests that that bullet rosettes (type C2a, combination of bullets) and single bullets are a primary ice particle habit in synoptically generated cirrus clouds (HI00; Paul Lawson 2000, personal communication). Given their significance, and because their 3D geometry has been studied and estimated (e.g., HI00), it is both possible and important to study how their ρe and Ar are related.
Bullet rosettes have been found with from one to more than eight component bullets, and the bullets were oriented in a number of 3D configurations. Based on the observations of Kobayashi et al. (1976) and others, Iaquinta et al. (1995) developed an analytic model that represents the 3D rosette geometry. In this model, each bullet was seen to have 10 surfaces: six for the hexagonal body (three rectangles and three pentagons), and four for the pyramidal end (one hexagon and three pentagons). Rosettes containing from one to eight bullets, each bullet with the same dimensions, were described by one of 19 configurations. Rosettes are 3D particles and will have a larger ρe for the same Ar and D than the simplified 2D rosettes in Fig. 1 because the bullets are longer for 3D rosettes (see appendix B).
The analytic model described in HI00 was used here to examine the ρe versus Ar relationships for 3D rosettes, again in 19 configurations, with from one to eight bullets. A configuration was identified by a given number of bullets, nb, and a 3D spatial arrangement of the bullets. For each configuration, αB—ranging from 0.1 to 1, and for a wide range of sizes, were specified; the volume occupied by the bullets in each rosette was then derived. The rosette images were numerically rotated in 3D until the maximum A and associated D were found from the projections onto a horizontal plane. The corresponding Ar and ρe2 (found using the projected dimension D rather than the maximum rosette D) were derived. The ρe were related to the Ar, as shown in Fig. 3a, assuming that ρb = 0.81 g m−3 (based on HI00). The curves for all 19 configurations were plotted in the figure to show the range of scatter for all configurations. In general, ρe increased with Ar.
The modeling results were simplified to obtain ρe versus Ar relationships as a function of nb. With appropriate averaging, curves of the form
were derived. The results are plotted in Fig. 3b, and coefficients k and exponents n for nb from one to eight are shown in Table 1. For a given Ar, increased nb resulted in lower ρe. The explanation for this effect is simple: as nb increased, the αB—and hence the bullet and rosette volume and ρe—decreased to maintain the same Ar. Also implied from the figure is that increasing nb produced larger Ar for a fixed αB. The value of k ranged from 0.47 to 0.66, and a weak, inverse dependence may be seen between k and nb. The exponent increased from 1.94 to 2.52, and there was a direct relationship between n and nb. An average ρe versus Ar curve was derived for rosettes by averaging over all configurations, but this curve did not capture the dependence reliably and is not shown in Table 1. It is suggested that the three-element curve reliably captures the ρe versus Ar relationship for from two to four bullets and the six-element curve does so for from five to eight bullets.
Comparison of the 2D, four-element rosette relationship from Eq. (8) (ρe = 0.41A2.12r, with ρb = 0.81 g cm−3) and the 3D, four-element rosette relationship from Table 1 (ρe = 0.47A2r) indicates that both the k and n for the 2D rosettes are smaller than for the 3D rosettes. This is not surprising because the lengths of 3D bullets are longer than lengths of 2D bullets, even though their projected lengths are the same; thus the ρe for 3D bullets is larger for the same Ar.
Photographs of rosettes in the literature (Heymsfield 1972) gave us the impression that ρb < 0.81 g cm−3, because the hollow ends of the bullets cover a large portion of the projected area. This question was further investigated by examining images of rosettes collected and photographed at the ground in Siberia at temperatures below −30°C by Klinov (1960). Many of the bullets in the rosettes had characteristic hollow ends. For each of 180 bullets, we estimated the hollow volume (υh) and the bullet volume (υb) and derived ρb from ρb(g cm−3) = 0.91[(υb − υh)/υb], as in Heymsfield (1972). Not all of the hollows appeared to be symmetrical, nor were they all conically shaped, although most of them did appear to be. The results plotted in Fig. 3c suggest that the mean ρb was approximately 0.82 ± 0.06 g cm−3. Much of the variance was attributable to measurement uncertainties at the small sizes, although no appreciable dependence on particle size was noted. The magnitudes of ρb found in this analysis were comparable to those found by Heymsfield (1972).
1) Cirrus crystal types
A dataset containing m, A, D, and Vt for rosettes was recently collected by M. Kajikawa at Mt. Hachimantai Observatory of Akita University, 1200 m mean sea level (MSL), where the surface temperatures ranged from −9° to −11°C. Of the 45 particles observed, 39 were classified as type C2a (combination of bullets) and the remaining as type C2b (combination of columns). Of the 45, approximately 25% of rosettes had nb = 2; 25% had nb = 4; 40% had nb = 3; and 10% had nb = 5. Crystal sizes ranged from 0.038 to 0.161 cm. The A and D parameters were obtained from crystal photographs (Fig. 4), and m from the diameters of the drops resulting from melting the particles in oil. Some of the rosettes were modified during their fall to the ground, as the bullets of a few rosettes had end plates, indicating some growth in the −12° to −18°C temperature range. Nevertheless, the dataset was useful for comparison with the rosette modeling results.
The Ar values for these rosettes are plotted as a function of ρe in Fig. 3D, with curves for nb = 3 and nb = 4 from the analytical modeling results plotted in the figure for comparison. A curve fitted to the Kajikawa observatory bullet rosette data was obtained and is plotted in Fig. 3d. The k for this curve (0.46) was the same as that found for four-element rosettes (0.47), but the n was larger (2.33 vs 2.16). The ρe fell slightly below the 3D model results, but the observations supported a ρe vs Ar dependence.
Information on the ρe vs Ar relationship for rosettes was obtained using imagery from a SPEC, Inc. cloud particle imager (CPI); images from this probe have resolution of approximately 2.3 μm. The CPI was mounted on the University of North Dakota Citation aircraft while participating in an intensive operation program (IOP) near the atmospheric radiation measurement (ARM) site in Oklahoma during March 2000. The particular CPI probe used during the IOP was owned and operated by the National Center for Atmospheric Research (NCAR). A dynamically active band of cirrus, with embedded generating cells and trails (cirrus uncinus), moved over the operations area on 9 March 2000, producing almost exclusively single rosettes and aggregates of rosettes suitable for study. Temperatures at the cloud tops were approximately −50° and −30°C, respectively. The rosettes observed on this day had Ar values similar to those observed on other days during the IOP.
We compared the ρe values as estimated from the CPI images to those from the analytic model. Referring to Fig. 3f, the volume of ice in each bullet of a rosette was calculated from the w and l, which were measured from the digital images using the National Institutes of Health's image analysis software package. The ice volume of each rosette was then derived by adding together the ice volumes of all bullets in the rosette. The ρe values were derived assuming ρb = 0.81 g cm−3. Here Ar and D measurements were also taken from the images. We examined a total of 180 rosettes, with nb ranging from one to eight.
These observed ρe amounts were compared to the ρe values for the same Ar and nb as obtained from the model. The ratio of the observed to modeled ρe as a function of D is shown in Fig. 3e; plots of this ratio against other variables, such as Ar, were not as informative and are not shown. The mean ratio over the measured range of D was slightly below unity, consistent with the modeling results (see appendix B); this shows that the ρe values estimated from 2D projections of rosettes were smaller than the actual ρe by 10% to 20% because the measured bullet lengths, which represent projected lengths, are smaller than the actual bullet lengths. For reference, the Ar values for the particles considered in Fig. 3e ranged from 0.2 to 0.5. This comparison also shows that the individual bullets became increasingly asymmetric with increasing values of D (see the large rosette in Fig. 3F). Asymmetry, presumably resulting from enhanced vapor growth from ventilation for the bullets on the underside of the rosette, was not considered in the analytic model.
Information on the Ar and ρe of side planes (habit types S1–S3) was obtained from observations; 3D modeling of side-plane geometry was beyond the scope of this study. Collections of side planes—from the Kajikawa study (Fig. 4) and by Knight and Heymsfield during cold outbreaks in Boulder, Colorado, from 1984 to 1986—were used to develop a dataset containing Ar and ρe information. The Knight and Heymsfield particles were collected in a cold room, photographed in hexane, then rephotographed following melting. A portion of this dataset was presented in Detwiler et al. (1993). A total of 85 particles (15 type S1, 35 type S2, 35 type S3), which ranged in size from 0.038 to 0.168 cm, composed the dataset.
Side-plane ρe generally increased with Ar; the resulting curve fit through the data is plotted in Fig. 5a and listed in Table 1. The values for n and k, as well as for ρe, fell below those for the 3D model rosettes and the rosette observations. This was not surprising because components of side planes had elongated or linear low-density features; this was especially true for types S1 and S2.
Inspection of Figs. 3 and 5a reveal some general characteristics of rosettes and side-plane-type particles. For a given Ar, rosettes tended to have somewhat larger ρe than side planes, on average. Rosette Ar fell in the range 0.3–0.7, with a mean of 0.45 ± 0.126, whereas side-plane Ar ranged between 0.2 and 0.6, with a mean of 0.39 μm ± 0.095.
Mitchell's (1996) and Mitchell et al's (1996) relationships between A and D and between m and D for a variety of cirrus crystal types were converted into ρe versus Ar relationships for comparison; the results are presented in Figs. 3 and 5a. These studies considered side planes (380–1680 μm) and aggregates of side planes (600–4100 μm), planar polycrystals (20–450 μm), and 5-bullet rosettes (5–1000 μm). The ρe values calculated from the various relationships for these particle types are presented in Fig. 5b; Ar amounts ranged between 0.3 and 0.8, and a given Ar corresponded to a given D. The M96 relationships may provide reasonable estimates for ρe within the Ar values commonly found for particles of these types.
2) Planar crystals
We should in Eq. (9) that the critical parameter needed to derive ρeP was the αP. Wind tunnel observations by Takahashi et al. (1991, hereafter T91) and observations at the ground reported in Heymsfield and Kajikawa (1987, hereafter HK87), were used to examine the ρe versus Ar relationships for planar crystals. Takahashi et al. grew several different types of planar crystals, at water saturation and constant temperatures, which were freely supported in a vertical wind tunnel. They measured crystal diameter, area, and mass as a function of time. We are indebted to Prof. T. Takahashi for supplying us with the original wind tunnel observations. We examined the data for four growth temperatures −12.2°C, where hexagonal plates (P1a) evolved to sector plates (P1b) after 7 min; −14.4°C, where broad-branched crystals (P1c) evolved into dendrites (P1e); −16.5°C, where sector plate crystals (P1b) grew; and −18.2°C, where plates evolved into plates with simple extensions (P2e). We also examined the masses and areas of 456 unrimed and rimed planar crystals collected at the ground in Japan (HK87).
We first examined how Ar varied as a function of αP, with the hope that αP could be expressed as a function of Ar; if this proved to be so, to establish the physical basis for a ρe versus Ar relationship. The αP values from T91 and HK87 were plotted as a function of Ar (Figs. 6a and b, respectively). The planar crystal thickness was not measured in the HK87 dataset, but an average t was calculated from known parameters: t = m/[(π/4) · 0.91(ArD2)]. The wind tunnel measurements shown in Fig. 6a indicated that for Ar < 0.8, αP and Ar were related. Curves fitted through the data, plotted in Fig. 6a and listed in Table 2, showed high correlation coefficients for Ar < 0.80 but low correlation for Ar ≥ 0.80. These data suggested, but did not confirm, that a general relationship existed between αP and Ar, as there was also a direct correspondence between Ar and D (T91).
The HK87 data further supported the idea that there was a general Ar versus αP relationship (Fig. 6b), with the exceptions of thick plate crystals (type C1h) and hexagonal plates (type P1a) below a hundred or so microns. Both the wind tunnel and observational datasets had similar dependencies (see curve fits plotted in Fig. 6b and listed in Table 2), even though the HK87 dataset included both unrimed and rimed crystals. Type C1h particles did not conform to the observations, and a mean αP for these particles was derived. Because a rather wide range of particle types and sizes were considered in the HK87 observations, and they conformed to the T91 observations, Eq. (9) suggested that there is a general relationship between ρe and Ar for planar crystals with Ar < 0.80.
The ρe versus Ar data from the wind tunnel and observational datasets are plotted in Figs. 7a and 7b, along with curves fitted to the data for Ar < 0.80. The curve fits are listed in Table 1 as planar crystals. Correlation coefficients were nearly unity. In general, the curves fitted to the wind tunnel and observational datasets were comparable.
The two-parameter approach breaks down when ArP ≈ 0.83, for types P1a (hexagonal plates), P1b (sector plates), and C1h particles. For such situations, a ρe versus D relationship was developed (see Table 1, planar crystals, Ar = 0.83). Using Eq. (9), this relationship was developed with Ar = 0.83 and ρb = 0.91 g m−3; we also used the Auer and Veal (1970) αP versus D relationships for crystal types P1a and P1b. A similar relationship has not been derived for C1h particles. In general no means is available for knowing whether a crystal is a thin (P1a) or thick (C1h) plate by viewing from above. However, P1a particles are probably more common, because C1h are usually restricted to very narrow temperatures near −10° and −20°C, and P1a cover a much wider range of sizes.
Aggregate mass as a function of D was frequently studied in the 1960s through the early 1980s in Japan and in the United States (Magono and Nakamura 1965, hereafter MN65; Locatelli and Hobbs 1974; Sasyo and Matsuo 1980; Kajikawa, 1982). Video distrometers in Switzerland have been used recently to measure aggregate A from the side, Aυ (a vertical rather than horizontal cross-sectional area), and Vt (Barthazy et al. 1998). However, the A and corresponding D for aggregates have rarely been studied, so there are virtually no relationships reported between ρe, or m, and A along with the corresponding D needed to derive ρe versus Ar relationships. In the following discussion, we reanalyze data from several previous studies to extract both ρe and Ar information for aggregates, and we report on our new work, which further studied the problem for aggregates comprising several types of component crystals.
Data from the Kajikawa (1982) and MN65 studies were reanalyzed to obtain Ar and ρe information for aggregates comprising planar-type crystals. Kajikawa measured the m and D for aggregates containing from two to six crystals—at the Mt. Teine Observatory of Hokkaido University, at a height of 1025 m MSL. Each of these embryonic or “early” aggregates was made up of one of the following groups of planar crystals: platelike habits (types P1a–P1c, P2e), branched planar crystals (types P1d–P1e, P2a–P2c), and rimed dendritic crystals (types R1d). From photographs of these aggregates in oil, we derived A and rederived D for 156 aggregates; they ranged from 0.11 to 0.604 cm in size. The resulting ρe versus Ar data are plotted in Fig. 8a. For the dataset as a whole, the mean Ar = 0.40 ± 0.11, indicating that the Ar of the early aggregates were generally small. Median ρe values in Ar intervals of 0.05 are plotted with dotted lines. Values for ρe increased from approximately 0.002 g cm−3 at Ar of 0.20–0.006 g cm−3 at Ar of 0.60, with considerable scatter. Also plotted in Fig. 8a is a curve for single planar crystals; the ρe were considerably larger for the same Ar.
We have been extremely fortunate to have obtained the original data from the MN65 study (kindly furnished by T. Nakamura). In the original notes, the component crystals were referred to as aggregates of “dendritic crystals,” “powdered snow,” or “wet snowflakes.” More than 300 particles composed the dataset, but we eliminated those that may have been partially melted (where the surface temperatures were above 0°C), or that did not have A information. Mass, area, and maximum dimension were recorded for a total of 185 particles, ranging in size from 0.60 to 3.3 cm. Given the difficulty in determining the areas of aggregates with hollow or concave regions, A and Ar were probably overestimated. For the dataset as a whole, the mean Ar was 0.68 ± 0.17, indicating that these particles had relatively high Ar. The median values in Ar intervals of 0.10 were plotted with solid lines (Fig. 8a). The ρe increased from approximately 0.003 g cm−3 at Ar of 0.30–0.009 g cm−3 at Ar of 0.90. Note that MN65 reported on the bulk density of the aggregates (aggregate mass divided by the volume of a circumscribed ellipse). Our reanalysis of the MN65 data indicated that the bulk densities were larger than ρe, on average, by a factor of 1.85.
Although considerable scatter was evident in the MN65 and Kajikawa (1982) datasets, both indicated increasing ρe associated with increasing Ar. We tried numerous types of fits to the data, separately for each set and also for both sets combined. With little loss of accuracy, we fit each of the datasets to curves of the form ρe ∝ A1.5r (Table 1). The exponent 1.5 was chosen so as to be consistent with early observations of aggregates with similar types of component crystals by Sasyo and Matsuo (1980) in Japan. Sasyo and Matsuo found that Aυ and m were related by
where the range of m was 0.0001 to 0.01 g; Aυ ranged from 0.02 to 1.0 cm2; and D, although not reported, was estimated to be from 0.1 to 1 cm. Inserting the definitions of ρe and Ar into Eq. (11) yielded the equation ρe = 0.016A1.5r, which is also listed in Table 1 and plotted in Fig. 8a. Taking Aυ to be approximately 0.9A (based on MN65), the proportionality constant in the Sasyo and Matsuo ρe versus Ar relationship was, therefore, approximately the same as for the Kajikawa (1982) data (0.015 g m−3) and 60% larger than for the MN65 data.
The ρe versus Ar relationship for aggregates comprising side-plane components was obtained from collections in a cold room at NCAR in Boulder, Colorado, by Knight and Heymsfield. The m, D, A, and Ar of 77 aggregates were obtained from the digitized photographs of the particles. The ρe values are plotted as a function of Ar in Fig. 8b, with median values of ρe plotted in Ar increments of 0.1. In general, ρe increased with Ar. A power-law curve fit to the median values is listed in Table 1. These ρe values did not differ much from those for individual side planes (curve also plotted in Fig. 8b), but were about an order of magnitude higher than for aggregates of dendrites for the same Ar. The corresponding ρe values derived for aggregates of side planes reported in M96 were comparable to those observed here, but the range of Ar was not as large.
Information on ρe and Ar for more than 100 aggregates of rosettes was obtained from CPI images on 9 March 2000. The volume of each bullet in each aggregate was measured, as illustrated in Fig. 9. The ice volume of each aggregate was derived by summing over all bullets in it. This volume was increased upward by 1/0.93, because the projected lengths of bullets in a rosette were on average approximately 93% as large as the true 3D lengths (appendix B). The masses and effective densities for aggregates were obtained by assuming ρb = 0.81 g cm−3. The Ar for each aggregate was also derived from the image.
A curve fit through the ρe versus Ar data for rosette aggregates is plotted in Fig. 8c and listed in Table 1. The ρe values for rosette aggregates were equal to or fell slightly below those for single rosettes with nb = 4 over much of the Ar range and were quite close to those for nb = 6 (not shown). The corresponding ρe values derived for aggregates of S3 from M96 were comparable to those seen here, but the range of Ar was not as large.
The resulting ρe versus Ar relationships for aggregates may be summarized as follows.
The exponent in the relationship was approximately 1.5.
The relationships for the cirrus crystal types (side planes, rosettes) were comparable to each other, and the resulting ρe values for given Ar increments do not differ appreciably from those for single (unaggregated) particles.
Values of ρe for aggregates of planar crystals are much lower than for the spatial crystal types, ranging from 0.001 to 0.015 g cm−3. The bulk densities of these aggregates are higher than ρe, on average by a factor of 1.85, but are still considerably lower than those for aggregates of rosettes or side planes.
The ρe for the various types of aggregates can also be examined as a function of D to assess whether there are systematic changes in ρe with D that may be inherent in the ρe versus Ar relationships. As the symbols in Fig. 10a show, ρe generally decreased with D in much the same way as reported in M96 and other studies.
Curves of the form ρe = aDb fitted to the observations are plotted in Fig. 10a over the relevant range of D values and are listed in the middle of Table 1. The slopes of the curves are close to unity, comparable to the slopes deduced from the mass versus diameter relationships for unrimed aggregates of side planes, bullets, and columns by Locatelli and Hobbs (1974) and in M96 of −1.1 and −0.9, respectively.
The observations shown in Fig. 10a suggest to a first approximation for aggregates ρe ∝ D−1, but is modulated by Ar. Drawing upon the general dependence ρe ∝ A1.5r, we developed a hybrid relationship that embodies the dependence of ρe on both Ar and D, ρe = 0.015A1.5rD−1 (also listed in Table 1). As shown in Fig. 10b, use of this equation to predict ρe reduces the need to have explicit knowledge of the habits of the aggregates. Although this method appears promising, a consistent set of values for mass and Ar over a wide range of habits is needed to establish the generality of the approach and more accurate estimate of the coefficients.
3. Terminal velocity–area ratio relationships
In this section, we describe an approach for deriving terminal velocity using ρe and Ar discussed in section 2. The Best number (X) vs Reynolds number (Re) approach, used by a number of researchers, provided a framework for deriving the Vt, which was relatively independent of the ice particle habit. This approach eliminated the need for an explicit drag coefficient; the parameter X was related to the particle drag coefficient Cd through X = CdRe2.
The Best number is given by
where g is the gravitational acceleration constant, ρa is the air density, and ν is the kinematic viscosity of air. Because
Eq. (12) was rewritten as
The Reynolds number is given by
A number of researchers have reported relationships between Re and X of the form
Solving for Vt in Eqs. (14) and (15)
The following coefficients af and bf were recently reported (M96):
Equation (17) provided an explicit means of deriving Vt in terms of D and Ar.
The accuracy of the Vt calculation depends on the accuracy of the coefficients af and bf; values for these coefficients were based primarily on theory and have not been well established from observations, especially for the larger ice particles and for aggregates in particular. Therefore, we have used values for m and A from the various observations at the ground discussed in section 2 to examine the accuracy of the coefficients as given in Eq. (18). The data included in this analysis are from HK87 for planar crystals; the Kajikawa collections for side planes and bullet rosettes; the Kajikawa (1982) observations for aggregates; and the MN65 data for aggregates; they spanned a size range from 0.03 to 3 cm, with Re from 5 to 5000. For each of these particle types, Vt was determined directly from the measurements by combining Eqs. (12) and (15):
which is based on the power-law relationships given in Eq. (18). Also, for each particle, X and Re were calculated from Eqs. (14) and (16). The average calculated Vt values were in excellent agreement with those measured Vt for the single (unaggregated) particle types (Fig. 11b), although a few outlying, presumably bad, data points were not omitted from the dataset. For aggregates, however, Vt values were overestimated (Fig. 11a) by 15% on average. This was shown to result from an underestimate of the drag coefficients, plotted in Fig. 12. Each point represents a measurement from the data, and the heavy solid line shows Cd implied from the af and bf coefficients as given in Eq. (18). Note the deviation from the data at Re > 700. Laboratory data from List and Schemenauer (1971) and measurements for graupel from HK87 showed that Cd for Re > 700 tended to increase with increasing Re, possibly because of the onset of tumbling and irregular movements (Kajikawa 1977) or by flow through the particle. The trend was also indicated in the data used for this study.
A correction at these large sizes was found by deriving a new af coefficient and bf exponent for the largest Re category in Eq. (18):
Unfortunately, very little data are available for values between Re = 200 and Re = 600, so the drag characteristics in this region remain unaccounted for. However, a qualitative comparison of the Re vs X values to those of List and Schemenauer (1971) and in HK87 suggested that the drag coefficients were reasonable.
Several general relationships may be inferred from Eqs. (17) and the coefficients in Eqs. (18) and (20). For synoptically generated cirrus, containing particles of several hundred to one thousand microns, bf = 0.831; therefore, nominally Vt ∝ D1.49. Most studies have shown that Ar decreases with D (by dividing the M96 ,A vs D equations by π/4D2), therefore the exponent < 1.49. The exponents were recently shown to fall in the range 0.70 to 1.30 for particles of these sizes by Matrosov and Heymsfield (2000).
In this section we examine whether the two-parameter approach is both reasonably accurate and useful. The technique was first used to calculate the IWC from 2D imaging probe measurements for three ice clouds, two of which were cirrus; the IWCs in these situations were also measured directly using a counterflow virtual impactor (CVI) (Twohy et al. 1997). Comparison of the calculations and measurements allowed us to assess the accuracy of the technique. The method was applied again when our balloon-borne ice crystal replicator (Miloshevich and Heymsfield 1997) was launched from the operations site during the First ISCCP Research Experiment (FIRE)-II field program in Oklahoma and ascended into a cirrus cloud at the same time that a vertically pointing millimeter cloud radar was measuring radar reflectivity. The replicator measurements in the cloud layer were obtained about 20 km downwind of the radar measurements, but radar observations in this case indicated that the reflectivities were relatively steady in time at given heights (Orr and Kropfli 1999). Therefore, a comparison of the calculated and measured radar reflectivity provided an indirect means of assessing whether the two-parameter approach provided realistic estimates in the vertical.
Two bulk properties of particle size distributions were examined—IWC and equivalent radar reflectivity factor (Ze or dBZe). The basic parameter needed to calculate each of these parameters is m, which was derived from
a. Comparison of IWCs in ice clouds
We were extremely fortunate to have had data from cirrus clouds on 3 days to assess the performance of the two-parameter approach. The measurements were made from the University of North Dakota Citation aircraft. The NCAR–Oregon State University CVI, a probe that directly measures the condensate mass above a threshold size of about 7 μm in diameter, provided an accurate measurement of the IWC (Twohy et al. 1997). The CVI had an IWC detection limit of about 0.003 g m−3, with a 50% uncertainty at that level. The uncertainty decreased to 20% at 0.010 g m−3. High-resolution images of the ice particles from the SPEC, Inc. CPI provided the information necessary to discern habits and the appropriate k and n parameters. Finally, size distributions were measured with PMS 2D-C and 2D-P probes, for the size ranges from 33 to above 1000 μm, and 200–6000 μm, respectively. The software used to analyze the 2D-C and 2D-P data retained the features reported on by Heymsfield and Parrish (1978) but added improvements to find the true D, not the D along the array or time axis.
This improved measure for D also improved the accuracy of Ar, because Ar depends on D as well as A. The 2D-C concentration measurements were used from the minimum detectable size to approximately 1000 μm, and the 2D-P data for 1000 μm and above.4
The data for case 1 were acquired on 9 March 2000 during a Lagrangian spiral descent over the Department of Energy–ARM site in Oklahoma in a band of cirrus clouds composing a group of generating cells. Particle habits above about 100 μm consisted almost exclusively of rosettes and aggregates of rosettes (Fig. 13). The mean number of bullets measured from CPI images for 869 rosettes was nb = 5.8 ± 0.6, within a range 2 ≤ nb ≤ 12. The value nb = 6 was used for the calculations. Although rosette aggregates were observed in the larger sizes, the data, plotted in Fig. 8c, showed that rosette aggregates have approximately the same ρe over much of the observed range of Ar as do bullet rosettes with nb from 4 to 6. Therefore, the k and n for unaggregated rosettes was used for all particles.
Time series of measured and calculated IWCs averaged over 5-s intervals during the Lagrangian spiral descent are plotted in Fig. 14. The spiral commenced near −50°C and ended near −26°C (top right panel). Measured IWCs from the CVI (top left panel) increased gradually with time (height below cloud top), and peak IWCs were about 0.1 g m−3. The panels below these show the calculated IWCs from the combined 2D-C and 2D-P data using the following methods: the two-parameter approach, followed by use of various mass versus dimension relationships: the habit-independent equation of Brown and Francis (1995; hereafter BF95), based on Locatelli and Hobbs (1974); five bullet rosettes of M96 and Mitchell et al. (1996); rosettes by Heymsfield (1975), based on direct collections of cirrus particles; and the rosettes by HI00, assuming nb = 6. Coefficients in these mass versus dimension relationships are listed in Table 3.
In the second through sixth panels on the right in Fig. 14, we show the differences between the IWCs as calculated from the size distributions for each technique and those measured by the CVI, (ΔIWC), obtained for 5-s averages of each. For the two-parameter approach, typical ΔIWC values were 0.02 g m−3 (20% accuracy) or less, mostly overestimates. The ΔIWC values did not change much with time (altitude), meaning that no biases resulted from particle size, which increased on average downward (see Fig. 13). The ΔIWC values found using the BF95 approach were small in the lower half of the layer but very large and positive in the upper part of the layer, where particles were smaller. Using the HI00 and M96 techniques, we found ΔIWC that were relatively small and negative in the upper third of the spiral where the particles were relatively small, but became increasingly negative with height. The ΔIWCs calculated using the Heymsfield (1975) coefficients compared favorably with the measurements, throughout the layer.
A reliable m versus D relationship can be developed to represent the cloud layer by making use of the results from the two-parameter approach for this cirrus cloud. By taking the amount of IWC in each 2D-C or 2D-P size bin and dividing by the concentration in the bin, an average m per size bin was obtained. These m per size bin were averaged through the cloud depth. Figure 15 shows the results of this exercise, where the m have been divided by (π/6)D3 to obtain ρe. The mass versus dimension relationships used in Fig. 14 and listed in Table 3 are also converted to ρe versus D relationships and are plotted in Fig. 15. At the smaller D, ρe fell below the ρe versus D curve obtained using the BF95 coefficients, which is why the BF95 approach overpredicted the IWCs in the upper part of the layer; at the larger D, calculations obtained from the BF95 coefficients fit the “observed” values reasonably well. The ρe versus D curve using the earlier Heymsfield (1975) coefficients fit the observed ρe to within a factor of two over a wide range of D. The coefficients in the M96 and HI00 relationships yielded values for ρe that fit the observed ρe reasonably closely when D were small but departed by an increasing amount with increasing D.
The data for case 2 were acquired over a period of 2.5 h during a step-down and step-up flight pattern near the Oklahoma ARM site on 13 March 2000; temperatures ranged from −50° to −31°C (Fig. 16, top right panel). Again, habits were dominated by rosettes, with nb values similar to those found on 9 March. The IWCs were derived in 5-s intervals from the 2D-C and 2D-P data, and the n and k for nb = 6 were used. The peak IWCs were higher than on 9 March, but the mean IWC over the period was lower. Note that this case differed from the 9 March case in that generating cells with large particles were not as prevalent, and smaller particles, some below the 2D-C probe's detection threshold, may have contributed appreciable IWCs.
The results were similar to those found for 9 March. The two-parameter approach produced accurate results throughout the period, except from 2020 to 2050 UTC, where the IWCs were underestimated; we suspect that particles below the detection threshold of the 2D-C probe were contributing a modest amount of IWC to the CVI. The BF95 approach resulted in calculations of low ΔIWCs at the warmer temperatures (larger particles) but yielded relatively large ΔIWCs at the lower temperatures (smaller particles). The M96 and HI00 m versus D relationships produced accurate estimates at the lower temperatures but gave significantly underestimated IWCs at the warmer temperatures.
To examine whether m versus D relationships posited in earlier studies could be improved if they were adapted to the two-parameter technique, we took the mass and area dimensional relationships for five-bullet rosettes from M96 and combined them to obtain a ρe versus Ar relationship (see Fig. 2b). If we compare the bottom two panels in Fig. 16, we find that use of the two-parameter technique reduced the ΔIWC amounts by about 50% through most of the flight, that is,—substantially improved the IWC estimates.
The ρe versus D data from the two-parameter approach for this case was similar to that observed for the 0309 case (Fig. 15). Comparison of these data with the curves in Fig. 15 explains why the BF95 approach overestimated the IWCs at the lower temperatures (small particles) and the HI00 and M96 relationships underestimated the IWCs at the higher temperatures (large particles).
A ρe versus D relationship representing the rosettes in the cloud layers sampled on 0309 and 0313 was derived by fitting a curve to the data points in Fig. 15:
The m versus D relationship can be found by multiplying ρe by (π/6)D3.
Case 3, a deep ice cloud layer over the Oklahoma ARM site on 12 March 2000, presented a large proportion of side-plane-type S3 particles and allowed us to evaluate the performance of the two-parameter approach for habit types different than bullet rosettes. The data were acquired over a period of 1–5 h during a step-down flight pattern, with temperatures ranging from −31° to near 0°C (Fig. 17, top right panel). For T < −15°C, habits were dominated by side planes (primarily type S3) and planar polycrystals. For T ≥ −15°C, the S3 and planar polycrystal particle habits aggregated and grew branched components. Light crystal riming and liquid water (at most, a few hundredths g m−3 measured by the forward scattering spectrometer probe) were occasionally observed. The IWCs were high throughout, at times reaching or exceeding the CVI's saturation limit of approximately 1 g m−3 after about 0020 UTC.
As shown by the right panels in Fig. 17, the two-parameter technique using the side-plane coefficients in Table 1 yielded underestimates of the measured IWC (before 0020 UTC) by about 0.10 g m−3, or 20%, with perhaps a slight increase in ΔIWC with time (decreasing altitude) noted. However, the ΔIWC values obtained using the rosette coefficients were closer to the observations, with no systematic trend with time. The BF95 coefficients produced slightly larger negative ΔIWC than for side planes, and even somewhat larger negative ΔIWC were produced using the M96 relationship for side planes (which does not differ much from those for planar polycrystals, which were also observed). As the particles were primarily of type S3, which contain rosette components (see Fig. 4), it is not surprising that the rosette k and n coefficients produce reasonable IWCs.
The deduced ρe values as a function of D were examined, and were about 50% larger than the ρe values shown for rosettes in Fig. 15.
b. Comparison of reflectivities in ice clouds
On 25 November 1991, during the FIRE-II experiment, the vertically pointing National Oceanic and Atmospheric Administration (NOAA) Ka-band radar located near Coffeyville, Kansas, was probing a cirrus cloud at the same time that our balloon-borne ice crystal replicator penetrated the cloud layer. The radar observations, from 1344 to 1400 UTC, were reported by Orr and Kropfli (1999), and the replicator observations reported in HI00. The more uniform part of the radar-observed cirrus layer extended from 7.5 to 8.5 km. A ragged base and virga—indicating sublimation—extended down to 7 km, and cirrus extended up to 9 km. Shown in Fig. 18 are the individual dbZe points as a function of height over the 16-min period; the radar reflectivity pattern during the period was relatively steady, with a variation of only a few dB at a given altitude.
Our balloon-borne ice crystal replicator was launched at 1337 UTC from the same ground site as the radar and ascended at a rate of about 4 m s−1. It passed through the cloud layer between 1408 and 1417 UTC, approximately 20 km downwind of the launch site (as determined from the balloon's latitude–longitude coordinates). Given the wind velocities, this region of cirrus would have been over the ground site from approximately 1345 to 1400 UTC. Replicator-derived cloud-top and cloud-base temperatures were −55° and −30°C, respectively.
The D and Ar were measured for each replicator particle. Size distributions—that is, number concentrations (N) as a function of D—were then derived for 28 layers, each about 100 m deep, from which the dBZe values were calculated:
where Dm is the melted equivalent diameter; the equivalent radar reflectivity, which the radar measured, was given by
The summation in Eq. (23) was carried out over all particles in a given averaging layer. Rosettes dominated in the middle and lower portions of the cloud, and quasi-spherical particles—presumably rosettes in the developing stages (much like the CPI images < 100 μm shown in Fig. 13)—predominated near the top.
The calculated and measured dBZe, for nb = 2, a lower bound, and nb = 4 are compared in Fig. 18; each nb gave about the same dBZe and the plotted lines lie on top of each other in the figure. Replicator cloud top extended almost 1 km above the radar-measured top, even though the radar had a minimum detectable reflectivity of about −40 dBZe. Within the body of radar-detectable cloud, excellent qualitative and quantitative agreement was found between the calculated and measured values.
The dBZe as calculated based on previously posited m versus D relationships are plotted in Fig. 19. These various formulations—from M96, HI00, and BF95—did not compare nearly as favorably as the two-parameter approach. This finding paralleled the results for the three ARM cases. In the higher cloud levels, where the particles were relatively small, the M96 and HI00 methods underestimated dBZe by about 5 dB; in the BF95 method, dBZe values were overestimated by about 5 dB. In the mid to lower cloud levels, where the crystals were relatively large, the M96 and HI00 methods underestimated the dBZe, whereas BF95 values were comparable to the measurements. Near cloud base, where the particles were again small, M96 and HI00 approaches provided reasonable results, but the BF95 method yielded overestimates of about 5 dB.
c. Results for aggregates
In this subsection, we examine the accuracy of the ρe versus Ar and X versus Re expressions given by Eqs. (17) and (18), and we offer a simplified means of deriving aggregate Vt directly from D.
We first examine how well the two-parameter approach predicts the mass for aggregates of planar crystals. The ρe for each of the Kajikawa (1982) aggregates was calculated using k = 0.01 and n = 1.5 from Table 1 and for the aggregates in the MN65 study k = 0.015 and n = 1.5 also from Table 1. The calculated and measured ρe are compared in Fig. 20a. This comparison is less than ideal as the ρe versus Ar relationships were derived from the MN65 and Kajikawa (1982) datasets, but does offer a means of examining how well the relationships fit the data. Although a considerable amount of scatter is evident about the 1:1 line in Fig. 20a—some natural and some resulting from experimental errors—the two-parameter approach produced reasonable estimates of ρe for these aggregates. More data (from aircraft measurements using the CVI, 2D, and CPI probes in deep stratiform ice clouds) are obviously needed to refine the estimates of ρe for aggregates.
We next examine how well the coefficients in the X versus Re relationships given by Eqs. (18) and (20) predict the terminal velocity of aggregates of planar crystals when m, Ar, and D are known; this comparison is a test of the new af and bf coefficients for aggregates given by Eq. (20). As shown in Fig. 20b, the new coefficients lead to reasonable estimates of Vt. However, despite considerable effort, we have been unable to explain the differences between Kajikawa's (1982) measurements of Vt and Vt from the calculations above 50 cm s−1.
We next examine how accurately Vt can be estimated when m is unknown but can be estimated from Ar and D. These Ar values are then used to calculate m and Vt for the Kajikawa (1982) and MN65 aggregates. In general, the calculated and measured values for Vt (Fig. 20C) were in reasonable agreement; differences primarily resulted from the scatter in ρe for aggregates.
We offer the following simplified means of deriving the Vt for aggregates, information that has not been adequately available in the literature. For levels other than at the surface, Eq. (19) shows that Vt is proportional to the ratio m/A, suggesting that a Vt versus D relationship can be derived if this ratio is known as a function of D. This ratio, plotted for several different types of aggregates in Fig. 21, has never been reported for aggregates based on the individual particle data. The m/A ratios were relatively independent of size for the aggregates and side-plane aggregates examined in MN65 and Kajikawa (1982), whereas a size dependence was noted for rosette aggregates.
In Fig. 21, curves of m/A versus D were derived for aggregates with several different types of component crystals from
It is reasonable to use a mean value for Ar (0.68) in Eq. (25) for aggregates of planar crystals observed in the MN65 study as there was little dependence of Ar on D. Inserting the appropriate k and n yields the relationship m/A = 0.0056D, which applies to well developed aggregates of planar crystals. This curve fits the data points in Fig. 21a reasonably well. Likewise, using a mean Ar (0.39) for the aggregates observed in the Kajikawa (1982) study yields the relationship m/A = 0.0063D for aggregates of from 2 to 6 planar crystals. This curve also fits the data points in Fig. 21a reasonably well. Similarly, taking a mean Ar value (0.46) for aggregates of side planes yields the relationship m/A = 0.0792D; this curve also provided a reasonable fit to the data in Fig. 21b. Finally, aggregates of rosettes did show a weak dependence of Ar on D, hence it was necessary to fit a curve to the Ar versus D values. Inserting this relationship into Eq. (25) with the appropriate k and n from Table 1 resulted in the relationship, m/A = 0.0228D0.58, which fits the data quite well.
where at a temperature and pressure of −10°C and 500 hPa, a = 160 and b = 0.91 for early (Kajikawa 1982) aggregates of planar crystals in the size range 0.1 < D 0.6 cm and a = 136 and b = 0.40 for well-developed (MN65) aggregates of planar crystals in the size range 0.6 < D 3.6 cm. At a temperature and pressure of −40°C and 300 hPa appropriate for aggregates of cirrus-type particles, a = 958 and b = 0.91 for side-plane aggregates for sizes 0.05 < D 0.2 cm and a = 536 and b = 0.68 for rosette aggregates for sizes 0.05 < D 0.2 cm. The a and b coefficients for different pressure and temperature levels can be readily found through use of Eqs. (18)–(20).
Over the past 30 years, habit-dependent relationships between particle mass or terminal velocity (which depends on mass) and dimension were developed to provide a means of prescribing these properties for observational, modeling, and remote-sensing applications. Our study uses a combination of theory, numerical modeling, and observational approaches to improve upon the reliability and generality of the one-parameter (D) method. We based our approach on the idea that A and D together, rather than D alone, were more reliable estimators of m and Vt. This two-parameter technique—when used in conjunction with particle size distribution measurements from three ice clouds sampled by aircraft and a fourth cloud sampled by our balloon-borne ice crystal replicator—produced better estimates of IWC and dBZe values than using D alone. The new approach was better able to capture the changes in the ice particle properties in the vertical and horizontal. The conversion of previous mass–dimensional relationships to this new technique through the use of appropriate area to dimension relationships reported in the literature resulted in improved estimates of IWC.
The two-parameter approach relied on the development of habit-dependent relationships between Ar and ρe of the forms ρe = kAnr and m = (π/6)ρeD3. This ρe relationship worked best for single particles. The addition of D into this relationship, ρe = kAnrDα, was shown to further improve the estimates of ρe for aggregates. Efforts should be made to relate ρe to Ar for more complicated crystal shapes (side planes, rimed particles) to better quantify the n and k parameters.
Using the results from the two-parameter approach for two cases of cirrus clouds comprising mainly bullet rosette crystals, it was possible to derive an explicit ρe and m versus D relationship for rosettes. Further study is warranted, as these relationships may vary from cloud to cloud and in the vertical for a given cloud; the relationships also depend on habit. When coupled with particle size distributions, they offer the promise of providing more reliable estimates of mean and mass-weighted fall velocities, effective radii, and radar reflectivities than have been possible in the past.
It is also recommended that more measurements of the physical properties of aggregates be acquired, for three reasons. First, more measurements would further constrain the relationships between ρe, Ar, and D, using consistent rather than disparate datasets. Second, the X versus Re relationship, developed in our study for aggregates in the range 0.6 < D < 3 cm, was only an approximation because of our limited sample size; what is needed is Vt, m, and A for aggregates in this size range. Third, although we developed Vt versus D relationships for aggregates with several types of component particles, based on the ratio of m/A, more information, particularly on m and A is needed to refine the Vt versus D relationships.
The authors wish to thank T. Nakamura, T. Takahashi, N. Knight, and K. Hilburn for the use of original data for this analysis. Help from W. Hall, who developed new 2D probe processing software; J. Oppenheimer and A. Lau, who synthesized CPI images; and B. Armstrong, who edited the manuscript is greatly appreciated. Discussions with D. Mitchell, S. Matrosov, and P. Field were helpful. David Mitchell carefully reviewed this manuscript. This research was supported by funding from NASA EOS (Contract S-97894-F) and the NASA TRMM Project Office (Grant NAG5-7743).
List of Symbols
A Particle cross-sectional area; e.g., the particle's area projected onto a horizontal plane
Ar Area of an ice particle projected onto a horizontal surface divided by the area of a circumscribed circle
ArP Area ratios of simplified columns, bullet rosettes, planar crystals
Aυ Cross-sectional area of an aggregate as viewed from the side
a Coefficient in terminal velocity–diameter relationship
af Coefficient in Reynolds number–Best number relationship
b Exponent in terminal velocity–diameter relationship
bf Exponent in Reynolds number–Best number relationship
D Maximum particle dimension; the dimension of a circumscribed sphere for columns and planar crystals, the diameter of the projected image of bullet rosettes, side planes, and aggregates
Dm Melted equivalent diameter of an ice particle, in mm
d Dimension across basal face of hexagonal planar crystal
dBZe Equivalent radar reflectivity factor
g Gravitational acceleration constant
ΔIWC Ice water content, difference between calculated and measured (CVI) IWC
k Coefficient in effective density–area ratio relationship
l Length of column, bullet
m Ice particle mass
N Concentration of particles of diameter D
n Exponent in effective density–area ratio relationship
nb Number of bullets in rosette
Re Reynolds number
t Planar crystal thickness
Vt Ice particle terminal velocity
υb Volume enclosing a bullet
υb Hollow volume of bullet
υr Ratio of particle volume to the volume of a circumscribed sphere. For rosettes and aggregates, the diameter of the circumscribed sphere is taken to be the projected dimension
w Column width
x, y Coefficient, exponent in power-law type aspect ratio-area ratio relationship
X Best number
T Temperature (°C)
Ze Equivalent radar reflectivity factor
αP Aspect ratios for columns (width/length), bullet rosettes (width/length), planar crystals (thickness/length)
γ, κ Coefficient, exponent in area–diameter equation
ν Kinematic viscosity
ρa Air density
ρb Bulk density; crystal mass divided by the volume of the exterior surface of a hexagonal column or bullet
ρe Ice particle mass divided by the volume of a circumscribed sphere
ρeB, ρeP Effective density for simplified columns, bullet rosettes, planar crystals
ρw Density of liquid water
σ, β Coefficient, exponent in mass vs diameter equation
χ, δ Coefficient, exponent in density vs diameter equation
Rosette Volume Calculated from Images
The ρe values of rosettes were estimated from images from the CPI probe. However, these were underestimates because the measured dimensions of the bullets were projections and hence were smaller than the true dimensions. The projected w of a bullet in a rosette probably had nearly the same value as the actual w, but the projected l was shorter than the actual l. The error in a bullet's volume was, therefore, proportional to the error in l to a first approximation.
The analytic model was used to estimate this error. The ratio of the mean projected l of the bullets in rosettes to the actual l was calculated for the eight-rosette configurations associated with nb from five to eight. The αB ranged from 0.2 to 0.6. The calculations showed that the ratio decreased from a mean of 0.88 (a difference of 12%) at 0.2 < Ar < 0.25 to 0.74 (difference of 26%) for 0.45 Ar < 0.50 (Fig. B1).
Corresponding author address: Andrew Heymsfield, 3450 Mitchell Lane, Boulder, CO 80301. Email: email@example.com
Hereafter, subscripts denoting the type of particle associated with a ρe − (e.g., ρeB) are eliminated to reduce overly long subscripts and because it is clear which habits are being referred to.
For all particle types, D was defined by the maximum size of the projected dimension of a particle. This was the exact dimension used for rosettes by the analytic model, for side planes, for aggregates of these types, for aggregates of dendrites, and for hexagonal columns [D = (w2 + l2)0.5]. The exception was for planar crystals—D = d(1 + α2P)0.5—as αP values were usually small (Fig. 6), D ≈ d.
A number of studies have shown that the 2D-C probe concentrations are underestimates below approximately 100 μm. We chose not to add additional particles to account for this underestimate because the cirrus clouds studied had a dominance of large particles. We also noted that the 2D-C probe underestimated the Ar below 200 or 300 μm, a feature we have not considered here but we will report on in a future study.