Shape Dependence of Falling Snow Crystals’ Microphysical Properties Using an Updated Shape Classiﬁcation

: We present ground-based in situ snow measurements in Kiruna, Sweden, using the ground-based in situ instrument Dual Ice Crystal Imager (D-ICI). D-ICI records dual high-resolution images from above and from the side of falling natural snow crystals and other hydrometeors with particle sizes ranging from 50 µ m to 4 mm. The images are from multiple snowfall seasons during the winters of 2014/2015 to 2018/2019, which span from the beginning of November to the middle of May. From our images, the microphysical properties of individual particles, such as particle size, cross-sectional area, area ratio, aspect ratio, and shape, can be determined. We present an updated classiﬁcation scheme, which comprises a total of 135 unique shapes, including 34 new snow crystal shapes. This is useful for other studies that are using previous shape classiﬁcation schemes, in particular the widely used Magono–Lee classiﬁcation. To facilitate the study of the shape dependence of the microphysical properties, we further sort these individual particle shapes into 15 different shape groups. Relationships between the microphysical properties are determined for each of these shape groups.


Introduction
The shape of ice particles is an important characteristic that affects the radiative impact of clouds. Accurate knowledge of the microphysical properties of clouds, including particle shape, is important in order to assure accurate cloud parameterizations in climate and meteorological forecast models, e.g., that presented in References [1,2]. Realistic cloud microphysical parameterizations are also essential for most cloud retrievals from satellite measurements. For instance, in order to retrieve quantities such as cloud water path or cloud effective radius, the underlying assumptions of particle shape, size, and distribution have a massive impact on the retrieval itself [3,4]. The unavoidable sensitivity of satellite retrievals to assumptions on ice particle properties, such as particle size, area, and the shape of snow crystals, is one of the dominating sources of uncertainties in cloud retrievals [5]. Although satellite observations of clouds are important validation sources for models, the role of clouds still poses one of the most substantial uncertainties in modeling the climate [6].
The shape of the snow crystals, which this paper primarily concerns, is formed by deposition, riming, aggregation, or a combination of these mechanisms [7]. Deposition refers to the process of water vapor freezing (depositing) onto ice nuclei or an existing ice crystal. If the ice crystals are in an environment with supercooled water droplets, these droplets may collide with and freeze onto the

Measurement Site
Our measurements have been carried out in Kiruna, Sweden (67.8 • N, 20.4 • E) using the D-ICI, the ground-based in situ instrument described in Kuhn and Vázquez-Martín [22]. Most of the measurements took place on the roof of the main building of the Swedish Institute of Space Physics (IRF) on the space campus at the height of about 15 m above the ground at an altitude of 425 m above mean sea level. The remaining measurements were done from the roof of the IRF Lidar laboratory located 950 m away from the space campus at the height of 5 m above ground. We consider these two locations close to each other and thus can further regard them as the same Kiruna measurement site.
Our dataset is from multiple snowfall seasons during the winters of 2014/2015 to 2018/2019. At the Kiruna site, each season lasts from about the beginning of November to the middle of May. Surface temperatures in Kiruna can reach down to −43.3 • C [24], and our data include temperatures down to −27 • C. The IRF weather database [25], based on instruments 1.7 km away from the space campus and 2.44 km from our current measurement location, provides temperature and relative humidity used in this study.

Instrument
The D-ICI recorded images of falling snow crystals and other hydrometeors. After falling through the inlet, in-focus particles are detected and then imaged twice, simultaneously from above and from the side. These dual images are taken using two high-resolution imaging systems (resolution of about 10 µm with a pixel resolution of 3.7 µm/px; for a detailed description, see [22]). The additional information provided by a second image improves the shape classification.
An inlet camera is mounted to monitor the D-ICI inlet, such that it is viewing the inlet without obstructing the path of the falling particles. The camera is used to detect inlet blockage in order to exclude, in post-processing, any data affected by this.

Image Processing
Top-view images are used to determine particle size, cross-sectional area, area ratio, and aspect ratio by the automated process presented in detail in Kuhn and Vázquez-Martín [22]. The first steps are to remove the background features from uneven illumination and to remove out-of-focus particles. The particles then have their boundaries traced. These particle boundaries are further used to determine the properties of the particle.
The maximum dimension, D max , is used to describe particle size. Here, it is determined by the smallest diameter that completely encircles the particle on the image. The cross-sectional area, A, is the area enclosed by the particle boundary on the image, determined from the number of enclosed pixels. Equation (1) defines the area ratio, A r , as follows: Equation (2) defines the aspect ratio, A s , as follows: where X is the width perpendicular to Y, the longest distance between any two points of the particle boundary.

Classification Method
The particle shape classification is carried out manually by looking at both top-and side-view images simultaneously since these provide complementary information about the shape. The images provide enough detail so that they can be classified according to the Magono-Lee scheme. Figure 1 shows a few examples of successfully classified snow particles. In some cases, classification is relatively simple with one image only, for example, the top view of the stellar crystal of example (a). However, looking at the side view of the same example, the ice particle's shape may be mistaken for a needle type, as it looks similar to example case (b), which is of needle shape. Thus, determining the shape from one image is challenging even when using high-resolution images due to ambiguities arising in the case of unfavorable orientation in that image.
Cases (c) and (d) are more examples, which, in the side view, appear to have a similar shape. However, from the top-view images, we can distinguish two different shapes: (c) has a dendritic crystal or stellar shape, and (d) is a bullet capped with a dendrite. Similarly, particles in (e) and (f) appear to have similar shapes in the side view, i.e., hexagonal plate shape, however, from the top-view images, we can distinguish two different shapes: (e) is a capped column, and (f) is a hexagonal plate capped on one side with a thin plate. In a few cases, it is difficult to specify the shape despite the availability of top-and side-view images. For instance, we cannot distinguish if (g), (h), or (i) are the same shape, i.e., if they are graupel-like snow or densely rimed ice crystals. For such cases, we determine that they have a graupel shape.

Recent Updates
When classifying snow particles using the shape classes from the Magono-Lee classification, for some particles, a suitable shape class is missing. To these particles, we have either assigned a shape found in more recent literature or a new shape class. In both cases, we have chosen shape labels that are consistent with the Magono-Lee nomenclature. Thus, this new classification is an update of the Magono-Lee classification, which is only an extension of this original classification scheme, so that the many previous studies or datasets that have been using Magono-Lee may still be used as a reference.
Some new shapes appear in the study of Kikuchi et al. [21], which uses labels that are different from Magono-Lee. Therefore, we include these shapes with new labels consistent with Magono-Lee.
Other shapes are added but keep the same labels as in Kikuchi et al. [21], as they are not conflicting with Magono-Lee.
We also add some shapes observed by Libbrecht [23] as an extension of the Magono-Lee classification. Since shapes by Libbrecht [23] do not have any labels, only a shape name, they are added by assigning compatible labels that are consistent with the Magono-Lee nomenclature. Table 1 shows the shapes added from Libbrecht [23] and Kikuchi et al. [21].

New Shapes from Kiruna
New shapes only found in Kiruna and not yet described are also included as an extension of the Magono-Lee classification with labels compatible with Magono-Lee. These 34 new shapes are named as shown in Table 2.  Figure 2 shows example images of all the new shapes found in Kiruna. In some cases, more images of the same shape are shown in order to show the small differences that can be found within the same shape, i.e., the accepted variability within the shape classification. In other cases, side-and top-view images of the particle shape are shown for better shape description. For instance, in shape CP2d, only by looking at the side view of the first particle can we see the bullet with two plates.  Table 2 shows their completed name. The same scaling is applied to all images; a 1.0-mm scale bar is shown as reference. In some cases, more images of the same shape are shown to indicate possible differences found in the same shape. In other cases, sideand top-view images of the particle shape are shown for better shape description.
The resulting updated classification of Magono-Lee has 135 shapes, including the 34 new shapes from Kiruna. In Appendix A, Table A1 shows the complete names of all these shapes. Figures A1 and A2 show images of all the particle shapes taken by D-ICI. The 34 new shapes found in Kiruna are highlighted by thick gray frames.

Shape Groups
The microphysical properties (particle size, cross-sectional area, and area ratio) are determined from the particle images taken with D-ICI. After classifying the particles, one can study these properties for specific shapes or how the properties depend on them. It is challenging to study the characteristics of all 135 shapes since most particle shapes infrequently appear so that statistically significant results are hard to achieve.
Furthermore, even if we had found enough examples of every single shape to achieve statistical significance, it may be hard for the models to utilize this level of detail. Additionally, due to natural variability, the encountered variation in microphysical properties of one shape would, in many cases, be larger than shape-to-shape variations.
Therefore, to analyze the shape properties, it seemed to be more useful to use a classification based on only a small number of shape groups, instead of the 135 individual shapes. In the following section, we describe how we gather the 135 shapes into 15 shape groups, which are then used to study shape dependencies. However, these shape groups can be assembled differently into, for instance, less than 15 groups, depending on the needs of the study. In this case, 15 shape groups were defined in an attempt to simplify shape description while still maintaining some level of detail.
In a similar approach to Magono-Lee [20], we gather the shapes into 15 groups with comparable morphological characteristics. Magono-Lee's main level P = Plane crystal contains both plate-like shapes and stellar shapes; however, we wish to separate the two for this study and, hence, group them somewhat differently. Therefore, the groups used here do not exactly follow the main levels of Magono-Lee. This is not seen as incompatible to the Magono-Lee classification scheme. Rather, the detailed classification using the updated Magono-Lee shapes from Section 3 allows one to group shapes according to the needs or interests, which may differ from study to study. In this sense, the grouping used here is not seen as rigid and definitive but as a first useful attempt to study snow particle characteristics depending on shape. These 15 shape groups are listed in Table 3. Table 3. Shape groups (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)(12)(13)(14)(15).

Shape Groups (1-15)
(1) Needles and thin or long columns (2) Crossed needles and crossed columns (3) Thick columns and bullets (4) Capped columns and capped bullets (5) Plates (6) Stellar crystals (7) Bullet rosettes (8) Branches (9) Side planes (10) Spatial plates (11) Spatial stellar crystals (12) Graupel (13) Ice and melting or sublimating particles (14) Irregulars and aggregates (15) Spherical particles Thus, we sort particle shapes that look similar, such as particles with needle-like shapes and thin columns, into groups of shared morphology. Therefore, even rimed particles can be included in such a group, as long as they appear to have the same type of morphology. Graupel particles, which do not appear to have a particular shape, are grouped separately as Graupel. Tables 4-7 provide detailed descriptions of the 15 shape groups. For each group, particle images of all the shapes that define that group are shown. In Appendix B, Table A2 lists for each shape group all shapes with their labels and names. Table 4. Shape groups and their descriptions (part 1): For each of the 15 groups, images of all included shapes are shown. The complete names of these shapes are shown in Table A2 in Appendix B.

Shape Groups (1-15) and their Descriptions Images
(1) Needles and thin or long columns: Includes needles, other snow crystals with elongated features (i.e., needle-like shape), and thin or long columns, including rimed needles and rimed thin columnar crystals.
(2) Crossed needles and crossed columns: Combination of crossed needles and a combination of crossed thin or long columns.
(4) Capped columns and capped bullets: Columns and bullets with plates at one end or both ends.
(5) Plates: Hexagonal plates including solid and hollow plates, skeletal surfaces, sector plates, plates with scrolls at the end, rimed plates, split plates, double plates, concentric plates, minute hexagonal plate, plate ice crystal, and triangular and rectangular plates. Table 5. Shape groups and their descriptions (part 2): For caption, see Table 4.

Shape Groups (1-15) and their Descriptions Images
(6) Stellar crystals: Stellar and dendrite shapes including simple stars, stellar with plates or sector-like ends, dendritic crystal with plates or sector-like ends, plate with simple, sector-like or dendritic extensions, broad branch or dendritic crystal with 12 branches, stellar crystal with scrolls at ends, rimed stellar crystal, split stellar crystals, double stellar crystal, concentric plates with extensions of different forms, minute stellar crystal, dendrite ice crystal, and triangular form with extensions of different forms.
(7) Bullet rosettes: Combination of bullets or columns, bullet rosette shapes including bullet with plates or dendrites.

Shape Groups (1-15) and their Descriptions Images
(9) Side planes: Snow crystals composed of several crossed hexagonal plates including radiating plates, side planes, scale-like side planes, combination of side planes with bullets and columns, assemblage of plates, arrowhead twins, and crossed plates shape.
(10) Spatial plates: Snow crystals composed of a plate with spatially attached plates or dendrites including densely rimed plate and columns on plates.
(11) Spatial stellar crystals: Stellar crystal or dendrites with spatial plates or with dendrites, including densely rimed stellar crystal, stellar crystal with rimed spatial branches, graupel-like snow of hexagonal type or lump type, and graupel-like snow with non-rimed extensions.
(12) Graupel: Soft hail shape considering that perceptions and definitions of soft hail can vary.  Table 4.

Shape Groups (1-15) and their Descriptions Images
(13) Ice and melting or sublimating particles: Minute ice crystals, i.e., particles that grow in cirrus clouds from minute frozen cloud droplets, such as frozen cloud particles, and ice pellets.

(14) Irregulars and aggregates:
Particles that cannot be classified into any other shape group in addition of the agglomeration of snow crystals.

Properties of Shape Groups
In this section, we analyze the properties of the shape groups defined in Section 4. More than 10,000 images were taken during snowfall events from 2014 to 2019. As mentioned in Kuhn and Vázquez-Martín [22], we only consider particles that are completely in the field of view for analysis. In total, we found 3165 particles from 67 days that fulfilled these criteria and they make up the dataset further referred to in this study. Figure 3 shows the frequency of occurrence of each shape group. We note that, in our dataset, shape group (8) Branches is the most common, occurring 20% of the time, followed by groups (14) Irregulars and aggregates, (9) Side planes, and (1) Needles and thin or long columns occurring 17%, 14%, and 11% of the time, respectively. The other shape groups have frequencies below 7.5%.  Figure 4a-c shows the particle size, D max ; cross-sectional area, A; and area ratio, A r , respectively, for each shape group. These quantities are shown by their median, represented as a point. The data spread is given by the percentiles 84.1% and 15.9%, which together correspond to ±1σ (standard deviation) if the distribution was normal (or Gaussian). The same data are shown in Table 8.

Occurrence and Properties
In Figure 4a, we note that shape group (6) Stellar crystals has the largest particle size with median D max = 1450 µm, followed by shape groups (2) Crossed needles and crossed columns, (7) Bullet rosettes, (8) Branches, (9) Side planes, and (11) Spatial stellar crystals with D max > 1000 µm. Shape groups (3) Thick columns and bullets, (5) Plates, (13) Ice and melting or sublimating particles, and (15) Spherical particles have the smallest sizes with median D max < 500 µm. Shape groups (11) Spatial stellar crystals and (14) Irregulars and aggregates have the widest spread ranging over 1000 µm, whereas shape groups (13) Ice and melting or sublimating particles and (15) Spherical particles have the smallest data spread with less than 300 µm. Figure 4b illustrates that shape group (6) Stellar crystals with the largest particle size also has the largest cross-sectional area, with median A = 6.34 × 10 −7 m 2 , followed by shape groups (7) Bullet rosettes, (8) Branches, (9) Side planes, and (11) Spatial stellar crystals with median A > 3.0 × 10 −7 m 2 . In general, groups with large D max also have large areas. However, there are exceptions, such as shape group (2) Crossed needles and crossed columns. While this shape group is among the groups with the largest D max , it is also among the shape groups with smaller areas. This is due to its shape leading to a minimal area ratio, as can be seen in Figure 4c. We note that shape groups (1) Needles and thin or long columns and (2) Crossed needles and crossed columns have the smallest area ratios with A r < 0.2. Shape group (15) Spherical particles has an area ratio close to one, as expected.  For comparison, median and the spread corresponding to all data (regardless of shape) are also shown in black dashed lines. Table 8. The total number of particles, #; particle size, D max ; cross-sectional area, A; and area ratio, A r , for the shape groups. The median and percentiles 15.9% and 84.1% are displayed. For comparison, results for all the data, i.e., all the shape groups together regardless of shape, are also shown.

Cross-Sectional Area and Area Ratio
The cross-sectional area versus particle size (A vs. D max ) and area ratio versus particle size (A r vs. D max ) relationships for the shape groups are analyzed in this section. The A vs. D max relationship is fitted to a power law given by Kuhn and Gultepe [26]: Equation (3) yields straight lines on a logarithmic plot; therefore, the parameters C 1 and C 2 are determined from linear fits to the data expressed as log(A) versus log(D max ), so that if all particles were spherical (i.e., they would appear circular in the images), C 1 and C 2 would be equal to π 4 × 10 −12 m 2 and 2, respectively.
The area ratio is calculated according to Equation (1) but can also be expressed as a power law function: where parameters C 3 and C 4 can be obtained from Equations (3) and (1) so that C 3 = C 1 × 4 π × 10 12 m −2 and C 4 = C 2 − 2. Figure 5 shows the fitted A vs. D max (left) and A r vs. D max (right) relationships as lines for each shape group. The point on each line represents the median of the distribution of D max , the lines' endpoints correspond to the percentiles 15.9% and 84.1%, and hence, the spread of data is illustrated by the length of the line. The fit parameters C 1 and C 2 as well as the corresponding correlation coefficients, R 2 , which are used to analyze correlation, are shown in Table 9. In general, judging from the correlation coefficient, R 2 a , in all shape groups, there is a good correlation (R 2 a ≥ 0.73) in the A vs. D max relationship. The shape groups (5) Plates, (11) Spatial stellar crystals, (12) Graupel, (13) Ice and melting or sublimating particles, (14) Irregulars and aggregates, and (15) Spherical particles have R 2 a > 0.9. For most shape groups, area ratio decreases with increasing particle size, as can be seen in Figure 5. Clearly, group (15) Spherical particles is an exception, and its area ratio is almost constant, i.e., size independent. A nearly constant A r means that C 4 0 and C 2 2 (see Equations (3) and (4)), which is the case for this group. Two more groups have values of C 2 close to 2. In group (12) Graupel, many particles have a close to spheroidal shape, and this appears to be similar over all sizes encountered, so that the area ratio is again size independent. The other group with C 2 fairly close to 2 and an almost constant A r is (6) Stellar crystals. Apparently, the stellar particles in this group also remain similar at different sizes, so that A r hardly varies. While A r for group (15) Spherical particles is close to 1, A r for (6) Stellar crystals is much smaller, approximately 0.4, as one would expect.
Particles in shape groups (1) Needles and thin or long columns, (2) Crossed needles and crossed columns, and (3) Thick columns and bullets, which are highly nonspherical, have C 2 values furthest from 2. As a consequence, C 4 values are furthest from zero and these groups have the steepest size dependence of A r , as can be seen in Figure 5.  Table 9.  Table 9.  Table 9. Tables 8 and 9 as well as Figures 4 and 5 show, in addition to the results for each shape group described above, also the results of A vs. D max and A r vs. D max fits to all data (regardless of shape). These are provided as references and for situations where one might need average relationships instead of shape-dependent ones. Table 9. Area versus particle size (A vs. D max ) and area ratio versus aspect ratio (A r vs. A s ) relationships fitted to Equations (3) and (4) for each shape group and for all data, i.e., for all the shape groups regardless of shape. The parameters C 1 , C 2 , C 5 , and C 6 and the correlation coefficients R 2 a and R 2 b are shown. (1-15) A vs. D max A r vs. A s

Area Ratio and Aspect Ratio
Area ratio versus aspect ratio (A r vs. A s ) data can be fitted to the linear relationship: where C 5 and C 6 are the fit parameters. Figure 6 shows the A r vs. A s relationship for each shape group. For comparison, the equality A r = A s is displayed as a gray dashed line. We note that particles in shape groups (3) Thick columns and bullets, (5) Plates, (13) Ice and melting or sublimating particles, and (15) Spherical particles are closest to this line. As can be seen in Table 9, these groups, together with group (1) Needles and thin or long columns, have a good correlation between area and aspect ratio (R 2 b ≥ 0.7). On the other hand, shape groups (2) Crossed needles and crossed columns, (6) Stellar crystals, (7) Bullet rosettes, and (9) Side planes have a low correlation (R 2 b < 0.5). For instance, in the case of the shape group (6) Stellar crystals, this might be due to several factors: the empty space between their branches resulting in low A r values and/or significant variations in A s due to particle orientation. Shape groups (1) Needles and thin or long columns, (3) Thick columns and bullets, (5) Plates, and (15) Spherical particles have the highest correlation (R 2 b > 0.8).
D-ICI takes dual high-resolution images (side and top view) that enable a better shape classification 286 than if there were only one image per particle. From these images, we were able to determine the 287 particle size, cross-sectional area, area ratio, and aspect ratio of individual particles. 288 We separated the particles into many different shapes following largely the Magono-Lee

293
This study sought to find some characteristic microphysical properties, particle size, 294 cross-sectional area, area ratio, aspect ratio, and frequency of occurrence of similar shapes. To 295 facilitate this, we further sorted the 135 shapes into 15 shape groups based on morphological similarity.

296
Naturally, it is also possible, and for some applications more suitable, to choose fewer groups.

297
Tables 4-7 (Section 4) show a detailed description of each of these groups, and a list with the 298 complete names of the particle shapes sorted into the shape groups is given in Table B1 (Appendix B).

299
The main conclusions of this study can be summarized as follows:

300
• The dual images provides more information about the particle shape.

301
• We found 34 snow particles that are not yet described and that warrant their own shape in the 302 updated shape classification.  Table 9.

Summary and Conclusions
This paper presents data and images of falling snow and other hydrometeors in the size range from 50 µm to 4 mm by the D-ICI during snowfall events in Kiruna, Sweden from 2014 to 2019. The D-ICI takes dual high-resolution images (side and top view) that enable a better shape classification than if there was only one image per particle. From these images, we were able to determine the particle size, cross-sectional area, area ratio, and aspect ratio of individual particles.
We separated the particles into many different shapes following largely the Magono-Lee classification scheme. This resulted in an updated classification with 135 shapes, including 34 new shapes found in Kiruna during this period, which are shown in Figures A1 and A2. The new shapes are highlighted by gray frames. This study sought to find some characteristic microphysical properties, particle size, cross-sectional area, area ratio, aspect ratio, and frequency of occurrence of similar shapes. To facilitate this, we further sorted the 135 shapes into 15 shape groups based on morphological similarity. Naturally, it is also possible and, for some applications, more suitable to choose fewer groups.
Tables 4-7 (Section 4) show detailed descriptions of each of these groups, and a list with the complete names of the particle shapes sorted into the shape groups is given in Table A2 (Appendix B).
The main conclusions of this study can be summarized as follows: • The dual images provides more information about the particle shape.
• We found 34 snow particles that are not yet described and that warrant their own shape in the updated shape classification.
• In our dataset, shape groups (1) Needles and thin or long columns, (8) Branches, (9) Side planes, and (14) Irregulars and aggregates are the most common, occurring more often than 10% of the time.
• In general, groups with large particle sizes also have large areas. However, there are exceptional shape groups with the largest particle sizes but with smaller areas, as it is in the case of shape group (2) Crossed needles and crossed columns. We can see in Figure 4c, in which we note that shape groups (1) Needles and thin or long columns, (2) Crossed needles and crossed columns, and (3) Thick columns and bullets have the smallest area ratios, A r < 0.4.

•
In general, there is a good correlation of the A vs. D max relationship for the shape groups, with correlation coefficients varying from 0.73 for shape group (2) Crossed needles and crossed columns to over 0.9 for shape groups (5) Plates, (11) Spatial stellar crystals, (12) Graupel, (13) Ice and melting or sublimating particles, (14) Irregulars and aggregates, and (15) Spherical particles (see Table 9).
• In most shape groups, area ratio decreases with increasing particle size. This is strongest for the shape groups with the lowest area ratios: (1) Needles and thin or long columns, (2) Crossed needles and crossed columns, and (3) Thick columns and bullets, which consequently also have the lowest values of C 2 . The shape groups (15) Spherical particles, (14) Graupel, and (6) Stellar crystals are exceptions and have area ratios that are almost constant.
• Shape groups (1) Needles and thin or long columns, (3) Thick columns and bullets, (5) Plates, and (15) Spherical particles have the highest correlation in the A r vs. A s relationship, with R 2 b larger than 0.8 (see Table 9). On the other hand, shape groups (2) Crossed needles and crossed columns, (6) Stellar crystals, (7) Bullet rosettes, and (9) Side planes have correlation coefficients lower than 0.5. Shape groups (6) Stellar crystals and (7) Bullet rosettes have a particularly low A r vs. A s correlation of 0.25 despite a A vs. D max correlation larger than 0.82. For instance, in the case of the shape group (6) Stellar crystals, this might be due to several factors; the empty space between their branches resulting in low A r values and/or variations in A s due to particle orientation.
This shows that, with the D-ICI and the classification in shapes and shape groups, particle characteristics can be studied for groups of similar shapes. Since fall speed can be measured by analyzing side-view images taken by the D-ICI [22], the shape dependence of fall speed can also be studied. The resulting parameterizations of the snow microphysical properties may be useful for improving our understanding of precipitation in cold climates in addition to helping improve the cloud microphysical parameterizations in the climate and forecast models.

Conflicts of Interest:
The authors declare no conflict of interest. Table A1 shows the updated classification of natural snow crystals and other hydrometeors. Original particle shapes are from the Magono-Lee classification. Added shapes from Libbrecht [23] are marked as [Li] and from Kikuchi et al. [21] are marked as [Ki]. Gray background and [KRN] mark the new particle shapes found in Kiruna.