Research article 23 Dec 2021
Research article  23 Dec 2021
Mass of different snow crystal shapes derived from fall speed measurements
 ^{1}Division of Space Technology, Department of Computer Science, Electrical and Space Engineering, Luleå University of Technology (LTU), 98 128, Kiruna, Sweden
 ^{2}Swedish Meteorological and Hydrological Institute (SMHI), 601 76, Norrköping, Sweden
 ^{1}Division of Space Technology, Department of Computer Science, Electrical and Space Engineering, Luleå University of Technology (LTU), 98 128, Kiruna, Sweden
 ^{2}Swedish Meteorological and Hydrological Institute (SMHI), 601 76, Norrköping, Sweden
Correspondence: Thomas Kuhn (thomas.kuhn@ltu.se)
Hide author detailsCorrespondence: Thomas Kuhn (thomas.kuhn@ltu.se)
Meteorological forecast and climate models require good knowledge of the microphysical properties of hydrometeors and the atmospheric snow and ice crystals in clouds, for instance, their size, crosssectional area, shape, mass, and fall speed. Especially shape is an important parameter in that it strongly affects the scattering properties of ice particles and consequently their response to remote sensing techniques. The fall speed and mass of ice particles are other important parameters for both numerical forecast models and the representation of snow and ice clouds in climate models. In the case of fall speed, it is responsible for the rate of removal of ice from these models. The particle mass is a key quantity that connects the cloud microphysical properties to radiative properties. Using an empirical relationship between the dimensionless Reynolds and Best numbers, fall speed and mass can be derived from each other if particle size and crosssectional area are also known.
In this study, groundbased in situ measurements of snow particle microphysical properties are used to analyse mass as a function of shape and the other properties particle size, crosssectional area, and fall speed. The measurements for this study were done in Kiruna, Sweden, during snowfall seasons of 2014 to 2019 and using the groundbased in situ Dual Ice Crystal Imager (DICI) instrument, which takes highresolution side and topview images of natural hydrometeors. From these images, particle size (maximum dimension), crosssectional area, and fall speed of individual particles are determined. The particles are shapeclassified according to the scheme presented in our previous study, in which particles sort into 15 different shape groups depending on their shape and morphology. Particle masses of individual ice particles are estimated from measured particle size, crosssectional area, and fall speed. The selected dataset covers sizes from about 0.1 to 3.2 mm, fall speeds from 0.1 to 1.6 m s^{−1}, and masses from 0.2 to 450 µg. In our previous study, the fall speed relationships between particle size and crosssectional area were studied. In this study, the same dataset is used to determine the particle mass, and consequently, the mass relationships between particle size, crosssectional area, and fall speed are studied for these 15 shape groups. Furthermore, the mass relationships presented in this study are compared with the previous studies.
For certain crystal habits, in particular columnar shapes, the maximum dimension is unsuitable for determining Reynolds number. Using a selection of columns, for which the simple geometry allows the verification of an empirical BestnumbertoReynoldsnumber relationship, we show that Reynolds number and fall speed are more closely related to the diameter of the basal facet than the maximum dimension. The agreement with the empirical relationship is further improved using a modified Best number, a function of an area ratio based on the falling particle seen in the vertical direction.
Atmospheric models need accurate knowledge of atmospheric ice crystals' and snow particles' microphysical properties to ensure realistic parameterizations (e.g. Stoelinga et al., 2003; Tao et al., 2003). These properties, including size, crosssectional area, shape, fall speed, and mass of ice particles, cannot be measured directly with remote sensing methods. Therefore, retrieval methods of cloud and snow properties also rely on good assumptions of the microphysical properties.
Particle shape is an essential parameter for retrievals of cloud properties from optical remote sensing (see e.g. Yang et al., 2008; Baum et al., 2011; Xie et al., 2011; Loeb et al., 2018). Furthermore, it can affect retrievals from active and passive microwave measurements of clouds and snowfall (e.g. Sun et al., 2011; Matrosov et al., 2012; Marchand et al., 2013; Kneifel et al., 2010; Cooper et al., 2017). Therefore, the shape dependence of the other microphysical properties is crucial to ensure accurate parameterizations. The fall speed of ice and snow crystals is a critical parameter for the modelling of the microphysical precipitation processes (Schefold et al., 2002) and the climate as it influences the lifetime of cirrus clouds, the vertical transport of water, and the snowfall rate (e.g. Mitchell et al., 2008). Ice particle mass parameterizations are required to derive ice water content (IWC). IWC is a crucial parameter to describe cloud contribution to the atmospheric models' radiation budget (Waliser et al., 2009; Thornberry et al., 2017).
Therefore, it is desirable to have datasets of falling snow particles based on simultaneous measurements of the microphysical properties maximum dimension (particle size), crosssectional area, shape, fall speed, and particle mass. If not available as measurement, particle mass or fall speed is retrievable based on all other properties. The fall speed cannot be computed directly from maximum dimension, crosssectional area, and mass because the drag force on the particle depends on the drag coefficient C_{D} that also depends on the fall speed. The dimensionless Best number X that only depends on maximum dimension, crosssectional area, and mass can eliminate this interdependency. The Best number can then help determine the Reynolds number, Re, through empirical relationships between Re and X. Finally, Re is used to calculate the fall speed, v.
For spherical particles, this Re–X relationship is well known (Abraham, 1970). Böhm (1989) suggested a modified Re–X relationship to determine v for all snow particles. Mitchell (1996) used that relationship to derive v vs. maximum dimension power laws from dimensional power laws of crosssectional area and mass. Heymsfield and Westbrook (2010) suggested a shapedependent modification of the Best number based on the area ratio. With this modified Best number, they showed that the error in fall speed determined from the Re–X relationship could be reduced for particles with open geometries, i.e. particles with low area ratio.
Instead of deriving fall speed from mass, the Re–X relationship may also be used to determine mass from measured fall speed. The Reynolds number can be derived from fall speed and then mass from X together with maximum dimension and crosssectional area. Szyrmer and Zawadzki (2010) have done this to determine average v vs. mass relationships from measurements of snow aggregates' fall speeds.
In this study, the Re–X relationship together with the modified Best number (Heymsfield and Westbrook, 2010) is used to determine masses of individual particles from measured maximum dimensions, crosssectional areas, and fall speeds given by the dataset of our previous study, VázquezMartín et al. (2021), that also includes particle shape. We analyse mass relationships as functions of maximum dimension, crosssectional area, and fall speed for different snow particle shapes. Section 2 describes the dataset used in this study. Section 3 shows the derivation of particle mass and mass relationships. Section 4 shows and discusses the resulting relationship between mass, size, crosssectional area, and fall speed. All relationships are studied separately for various particle shapes. In the same section, we also present comparisons between our mass relationships and those from previous studies. In Sect. 5, this study is summarized and concluded.
The dataset consists of 2461 highresolution dual images of falling natural snow crystals and other hydrometeors. The same dataset has been used in VázquezMartín et al. (2021). The data have been collected using DICI, the groundbased in situ instrument described in Kuhn and VázquezMartín (2020), at a site in Kiruna, Sweden (67.83^{∘} N, 20.41^{∘} E), described in VázquezMartín et al. (2020) during multiple snowfall seasons, the winters of 2014/15 to 2018/19. The images are taken when the snow particles fall into the inlet and consequently fall down the sampling tube and traverse the optical cell. In the centre of the optical cell is the sensing volume. If particles are falling through the sensing volume they are detected by the detecting optics (for a detailed description see Kuhn and VázquezMartín, 2020). Upon detection, the particles are optically imaged simultaneously from two different viewing directions. One is horizontal, recording a side view, and one is close to vertical, recording a top view. From the topview images, we can determine for each particle its maximum dimension D_{max}, which we use to describe particle size, crosssectional area A, and area ratio. From the sideview images, since they are exposed twice, we can determine fall speed. These images are highresolution (optical resolution of about 10 µm), where one pixel corresponds to 3.7 µm. The additional information dual images provide improves the shape classification carried out by looking at both top and sideview images. The particles are classified according to their shape and sorted into 15 different shape groups as described in VázquezMartín et al. (2020). A complete description of the dataset and data processing methods is given by VázquezMartín et al. (2021).
3.1 Mass derivation
The motion of hydrometeors when freefalling through the atmosphere establishes an equilibrium between two forces: the gravity and the aerodynamic drag. The resulting particle settling speed is called fall speed v. Thus, the fall speed is governed by the physical properties of the hydrometeors, including their mass and projected area, and it involves aerodynamic principles and environmental conditions. The gravitational force is proportional to the particle mass m, while the frictional or drag force is proportional to both the particle projected area, i.e. the crosssectional area A, and the square of its fall speed v. The force balance yields
where g is the gravitational acceleration, ρ_{a} the air density, and C_{D} the drag coefficient. To determine v from the particle properties m and A using this equation, the drag coefficient C_{D} has to be known as well. However, C_{D} depends on maximum dimension, shape, and on v itself. To circumvent these interdependencies, one can first determine the Best number $X={C}_{\mathrm{D}}\cdot \phantom{\rule{0.125em}{0ex}}{\mathit{\text{Re}}}^{\mathrm{2}}$ by rearranging Eq. (1) together with the Reynolds number
where η is the dynamic viscosity of air, to get an expression that does not depend on fall speed v:
Thus, X can be calculated from the particle properties D_{max}, A, and m. If the relationship between Re and X is known, one can determine Re from X. In these circumstances, Eq. (2) provides the fall speed, v. Böhm (1989) provides a relationship between Re and X for snow particles, which is shown here in the form given by Mitchell (1996)
where δ_{0} and C_{0} are unitless constants, and uses it together with the approach described above to determine v from the particle properties D_{max}, A, and m.
In a similar approach, one can determine particle mass if D_{max}, A, and v are known. For this, Re is determined from v and D_{max} using Eq. (2). Then, X is determined from Eq. (4):
Finally, m is added to the dataset using Eq. (3)
where the atmospheric conditions can be accounted for each particle by adapting η and ρ_{a} to the measured temperature and pressure.
Instead of using Eqs. (4) or (5) with one set of δ_{0} and C_{0} for all particles regardless of their shape, as proposed by Böhm (1989), Heymsfield and Westbrook (2010) suggested using a modified Best number X^{*}, replacing X in Eqs. (4) or (5), to correct for effects due to opengeometry shapes. They proposed ${X}^{*}=X\cdot {A}_{\mathrm{r}}^{\mathrm{1}/\mathrm{2}}$, where ${A}_{\mathrm{r}}=\frac{A}{\frac{\mathit{\pi}}{\mathrm{4}}\cdot {D}_{\text{max}}^{\mathrm{2}}}$ is the area ratio, which is close to 1 for compact shapes and smaller the more open the geometry is. Heymsfield and Westbrook (2010) showed that by using this approach they could reduce errors in determined fall speeds associated with opengeometry particles with low area ratios. Using our data for simple thick columns in shape group (3), we could confirm that their approach is better than the approach by Böhm (1989) without modifying X (see Appendix C). Therefore, here, we use the modified Best number X^{*}. Consequently, Eq. (6) is modified to
Note that then the Best number determined from Eq. (5) is the modified Best number X^{*}. In Eq. (5), we use δ_{0}=8.0 and C_{0}=0.35 from Heymsfield and Westbrook (2010).
3.2 Fitting relationships to data
Once mass is calculated, by fitting the following power laws to our data we can parameterize the relationships mass vs. maximum dimension, m(D_{max}); mass vs. crosssectional area, m(A); and fall speed vs. mass, v(m):
which represent straight lines on logarithmic plots. Hence, linear leastsquares fits to the logarithm of the data yield the parameters ${\stackrel{\mathrm{\u0303}}{a}}_{D}$, ${\stackrel{\mathrm{\u0303}}{b}}_{D}$, ${\stackrel{\mathrm{\u0303}}{a}}_{A}$, ${\stackrel{\mathrm{\u0303}}{b}}_{A}$, a_{m}, and b_{m}. The parameter ${\stackrel{\mathrm{\u0303}}{a}}_{D}$ corresponds to the mass at D_{max}=1 mm, ${\stackrel{\mathrm{\u0303}}{a}}_{A}$ to the mass at A=1 mm^{2}, and a_{m} to the fall speed at m=1 µg. The parameters ${\stackrel{\mathrm{\u0303}}{b}}_{D}$, ${\stackrel{\mathrm{\u0303}}{b}}_{A}$, and b_{m} are the exponents in the power laws and the slopes in the linear fits.
As seen in VázquezMartín et al. (2021), using binned data instead of individual data reduces the data spread so that fit functions based on binned data are more robust than fit functions based on individual data. Therefore, also here the data are first binned into a suitable number of bins before fitting Eqs. (8)–(10) to the data. A total of 10 mass bins (for m vs. D_{max} and m vs. A relationships) and 10 fall speed bins (for v vs. m) are used, respectively. The bins are spaced such that each bin contains as close to the same number of particles as possible. As a consequence, the bin widths are variable and specific to each shape group and thereby avoid the problem of individual bins having a disproportional effect on the fit. The binned data consist of the median values for each bin. Then, the m vs. D_{max}, m vs. A, and v vs. m relationships are fitted to the median masses vs. median maximum dimensions, median masses vs. median crosssectional areas, and median fall speeds vs. median masses, respectively. VázquezMartín et al. (2021) found that about 40 particles in a shape group (currently the lowest number in our dataset is 37) is the limit where binning can still be used. The advantages of binning become prominent only at larger numbers of particles.
3.3 Analytical derivation of relationships
These relationships may be useful for parameterizations in models and retrievals and are readily comparable to other studies. In case a suitable dataset is not available, an alternative to fitting these relationships to measured data is to derive particle mass analytically from previously determined parameterizations of crosssectional area vs. maximum dimension (A vs. D_{max}), fall speed vs. maximum dimension (v vs. D_{max}), and fall speed vs. crosssectional area (v vs. A) given by power laws
For each relationship, the inverse is also shown as the corresponding parameters are convenient for some of the derivations. The parameter a corresponds to the crosssectional area at D_{max}=1 mm, a^{′} corresponds to the maximum dimension at A=1 mm^{2}, a_{D} to the fall speed at D_{max}=1 mm, ${a}_{D}^{\prime}$ to the maximum dimension at v=1 m s^{−1}, a_{A} to the fall speed at A=1 mm^{2}, and ${a}_{A}^{\prime}$ to the crosssectional area at v=1 m s^{−1}. The parameters b, b^{′}, b_{D}, ${b}_{D}^{\prime}$, b_{A}, and ${b}_{A}^{\prime}$ are the exponents in the power laws.
The resulting power laws are
The derivation of these power laws is shown in Appendix B (Eqs. B3–B6). There, also the X vs. Re relationship is expressed as a power law instead of using Eq. (5). This can be done by approximating Eq. (5) piecewise in several regions of X with power laws Eq. (B1) (with coefficient γ and exponent δ), as done by Mitchell (1996). Note that both methods for deriving the relationships given by Eqs. (8)–(10), i.e. either the method described in Sect. 3.1 with fitting detailed in Sect. 3.2 or the alternative derivation from existing relationships described in this section, are equivalent if they are based on the same dataset. The two methods will yield the same relationships if both use the same power law approximations of X vs. Re and the same atmospheric conditions (given as constant η and ρ_{a} for the whole dataset). Thus, in this study, we have chosen to fit Eqs. (8)–(10) directly to our data (Sect. 3.2). This allows the use of environmental conditions individually for each particle and avoids the need to consider error propagation when deriving new relationships from existing ones.
4.1 Results from fitting and correlations
The particle masses have been determined from measured D_{max}, A, and v with the method described in Sect. 3.1. The m vs. D_{max}, m vs. A, and v vs. m relationships given by Eqs. (8)–(10) are then fitted to the resulting data, now consisting of D_{max}, A, v, and m, for the 15 shape groups using the fitting method based on binned data described in Sect. 3.2. Figure 1 and Table 1 show the results. For simplicity, we use short names included in Table 1 for the shape groups from here on, and Fig. 1 shows their full names. The large spread in the data represented as individual points is apparent in Figs. A1–A3 in Appendix A.
When fitting m vs. D_{max}, m vs. A, and v vs. m relationships to the binned data, we note that, in general, there is a high correlation (0.9≲R^{2}≲ 1) for most shape groups. In the following, we call the correlation coefficients ${R}_{D}^{\mathrm{2}}$, ${R}_{A}^{\mathrm{2}}$, and ${R}_{m}^{\mathrm{2}}$ to indicate which of the three relationships they belong to. For the m vs. D_{max} relationship, the only exceptions to high correlations are shape groups (1) Needles, (2) Crossed needles, and (3) Thick columns as well as (6) Stellar and (10) Spatial plates, which both have a low number of particles, all having ${R}_{D}^{\mathrm{2}}\simeq \mathrm{0.7}.$ For the m vs. A relationship, only shape groups (2) and (6), with ${R}_{A}^{\mathrm{2}}\simeq \mathrm{0.8}$, and (10), with ${R}_{A}^{\mathrm{2}}\simeq \mathrm{0.5}$, have a lower correlation. Only in the case of (10) is it uncertain if the fit function is representative of the measured data, as judged by the low ${R}_{A}^{\mathrm{2}}$. In most shape groups, the coefficients ${R}_{D}^{\mathrm{2}}$ and ${R}_{A}^{\mathrm{2}}$ are similar. Only the four groups (1), (2), (3), and (10), mentioned above with lower correlation in one of the relationships, have a distinct difference between ${R}_{D}^{\mathrm{2}}$ and ${R}_{A}^{\mathrm{2}}$. Of these, the three shape groups (1)–(3) have a better correlation for m vs. A than for m vs. D_{max}, which is consistent with a better v vs. A correlation than v vs. D_{max} for the same groups (VázquezMartín et al., 2021), given that we have derived m using measured v here.
For the v vs. m relationship, all values of ${R}_{m}^{\mathrm{2}}$ are 0.85 or higher. These high values indicate that v is better correlated to m than to D_{max} or A (see the generally lower R^{2} values reported in VázquezMartín et al., 2021). The generally very high correlations are partly also a consequence of m being derived from v rather than being an independent measurement.
4.2 Mass vs. D_{max} and A
Figure 1a and b show the m vs. D_{max} and m vs. A relationships including, for reference, the mass of liquid water spheres symbolizing rain or fog droplets given by the power laws $m=\frac{\mathit{\pi}}{\mathrm{6}}\cdot {\mathit{\rho}}_{\mathrm{w}}\cdot {D}_{\text{max}}^{\mathrm{3}}$ and $m=\frac{\mathrm{4}\cdot {\mathit{\rho}}_{\mathrm{w}}}{\mathrm{3}\cdot \sqrt{\mathit{\pi}}}\cdot {A}^{\mathrm{3}/\mathrm{2}}$, respectively, where ρ_{w}=1 g cm^{−3} is the density of liquid water. The mass of spheres is proportional to ${D}_{\text{max}}^{\mathrm{3}}$ and to ${A}^{\mathrm{3}/\mathrm{2}}$. Thus, comparing to Eqs. (8) and (9), one can see that the exponents ${\stackrel{\mathrm{\u0303}}{b}}_{D}=\mathrm{3}$ and ${\stackrel{\mathrm{\u0303}}{b}}_{A}=\mathrm{1.5}$ for spheres. The values of ${\stackrel{\mathrm{\u0303}}{a}}_{D}$ and ${\stackrel{\mathrm{\u0303}}{a}}_{A}$ for spheres are 524 µg, the mass of a droplet with 1 mm diameter, and 752 µg, the mass of a droplet with a crosssectional area A=1 mm^{2}, respectively.
4.2.1 Slopes $\stackrel{\mathrm{\u0303}}{b}$_{D} and $\stackrel{\mathrm{\u0303}}{b}$_{A}
The exponent ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ for shape groups (12) Graupel and (15) Spherical is close to the value of 3 for spheres, 2.74 and 2.84, respectively. For the same groups, ${\stackrel{\mathrm{\u0303}}{b}}_{A}$ is close to the value of 1.5 for spheres, 1.34 and 1.43 for shape groups (12) Graupel and (15) Spherical, respectively. For these shape groups, this is expected due to their spherical or roundish morphology. These exponent values, corresponding to the slopes in Fig. 1a and b, are among the highest values for all shape groups. Shape groups (6) Stellar and (11) Spatial stellar are the only other shape groups that have similarly steep m vs. D_{max} and m vs. A relationships. These two groups do not have a roundish morphology that could explain this. However, a slope similar to spherical particles may indicate that the morphology remains similar in these groups independent of size, i.e. ice particles scale equally in all three dimensions. An example for this would be hexagonal plates or columns that all have the same aspect ratio. For pristine stellar particles one may not expect such a steep slope similar to spherical particles but rather a decreasing area ratio with increasing size. Shape group (6), however, contains other shapes besides pristine stellar particles, such as rimed stellar and split stellar crystals. A particular mix of shapes may cause an apparently steep slope. Indeed, the area ratio in this shape group is approximately constant (VázquezMartín et al., 2021). Our dataset does not contain a sufficient number of stellar particles yet to analyse this further by, for example, regrouping particle shapes. Additionally, the low number of particles in this group also results in a relatively high uncertainty (${\stackrel{\mathrm{\u0303}}{b}}_{D}=\mathrm{2.61}\pm \mathrm{0.59}$ and ${\stackrel{\mathrm{\u0303}}{b}}_{A}=\mathrm{1.34}\pm \mathrm{0.29}$).
For most other shape groups, the exponent ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ varies between 1.2 and 2, and all other ${\stackrel{\mathrm{\u0303}}{b}}_{A}$ values range between 0.8 and 1.2. Three shape groups, (1) Needles, (2) Crossed Needles, and (3) Thick columns, stand out with the lowest exponents ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ of approximately 0.8 or lower. These can easily be seen in Fig. 1a as the lines with the most shallow slopes. For these groups, this is understandable due to their morphology. We have seen in VázquezMartín et al. (2021) that an increase in D_{max} (needle length) is directly proportional to A, indicating that the diameter of these needleshaped particles (needle width) remains similar when D_{max} and consequently also A are growing. Thus, D_{max} is approximately proportional to A, and predictably, both ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ and ${\stackrel{\mathrm{\u0303}}{b}}_{A}$ are close to 1 for these three shape groups. VázquezMartín et al. (2021), observing the very poor correlation between D_{max} and measured fall speed for these shape groups, argued that D_{max} is not suitable to determine the Reynolds number. Therefore, a more suitable characteristic length than D_{max} should be used to determine the Reynolds number and derive mass from it. Otherwise, the derived mass and consequently ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ are likely not useful. Jayaweera (1971) suggested a characteristic length for hexagonal crystals, for which the dimensions of the basal facet and the aspect ratio are known. We have tested this for a small subset of relatively simple columns from shape group (3). Indeed, this characteristic length, which can easily be determined from the geometric dimensions on the image, is better suited to calculate Re. If it is used instead of D_{max} in our method to determine m (see Sects. 3.1 and 3.2), then ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ is now 2.4 (for this subset of columns). Furthermore, the mass can be estimated directly from the geometric dimensions, which allows a closure study to confirm the Re–X relationship, which is central in our method (see Sect. 3.1). The details of this closure study are reported in Appendix C.
Unfortunately, for many particles the dimensions of the basal facet and the aspect ratio are not readily available (or are not defined in case of more complex particles). Consequently, the characteristic length cannot be used, and one may want to use D_{max} instead. However, as the closure study shows (see Appendix C), for shape groups (1)–(3) D_{max} should not be used as a representative size parameter instead of characteristic length. Thus, mass determined using D_{max} and relationships based on this mass should not be used.
The ratio between the exponents ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ and ${\stackrel{\mathrm{\u0303}}{b}}_{A}$ is equal to the exponent b, as can be seen from Eqs. (8), (9), and (11a). Figure 2 shows the ratios $\frac{{\stackrel{\mathrm{\u0303}}{b}}_{D}}{{\stackrel{\mathrm{\u0303}}{b}}_{A}}$ plotted vs. b, and most ratios on this plot are close to the line $\frac{{\stackrel{\mathrm{\u0303}}{b}}_{D}}{{\stackrel{\mathrm{\u0303}}{b}}_{A}}=b$ and range between 1.7 and 2. The exceptions which have ratios much below the line are the two shape groups with the lowest ${R}_{D}^{\mathrm{2}}$, groups (1) Needles and (2) Crossed needles. The ratios for shape groups (3) Thick columns, (9) Side planes, and (13) Ice particles are found slightly below the line. Of these groups, (3) and (13) are among the groups showing more uncertainty in the determined relationship, as indicated in Fig. 1a by the larger confidence regions around the fits. For group (9) Side planes, the uncertainty is smaller and cannot explain the lower ratio.
Intuitively, ${\stackrel{\mathrm{\u0303}}{b}}_{D}$, the exponent of the m–D_{max} relationship, should be larger than b, the exponent of A–D_{max}, as confirmed by the literature, such as by Mitchell (1996). For some shape groups, however, b is larger than ${\stackrel{\mathrm{\u0303}}{b}}_{D}$. Not surprisingly, groups (1)–(3) that were noticed earlier for the lowest ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ values are among these groups as well as groups (9) and (10). The latter two were noticed in VázquezMartín et al. (2021) with very poor correlations between D_{max} or A and fall speed. This problem likely indicates that for these shape groups D_{max} is not suitable as a size parameter to calculate Re. For simple thick columns, this is demonstrated in Appendix C. While suitable substitutes exist for regular shapes, such as the characteristic length suggested by Jayaweera (1971), for an arbitrary shape our current image analysis methods cannot determine a similar quantity. Thus, the modified X^{*} approach according to Heymsfield and Westbrook (2010) remains the best alternative for our study; it lessens the problem considerably for groups (1)–(3).
4.2.2 Coefficients ã_{D} and ã_{A}
All relationships but those of shape group (15) Spherical form a cluster of lines located in a smaller region in both Fig. 1a and b. The only relationship found outside this cluster is that of shape group (15), which, if extrapolated towards larger sizes or crosssectional areas, predicts larger masses than any relationship of the other shape groups. The fit coefficients ${\stackrel{\mathrm{\u0303}}{a}}_{D}$ and ${\stackrel{\mathrm{\u0303}}{a}}_{A}$ reflect this since they predict the mass at the unit reference of 1 mm for ${\stackrel{\mathrm{\u0303}}{a}}_{D}$ and 1 mm^{2} for ${\stackrel{\mathrm{\u0303}}{a}}_{A}$. These values are much larger for spheres, ${\stackrel{\mathrm{\u0303}}{a}}_{D}=\mathrm{260}$ µg and ${\stackrel{\mathrm{\u0303}}{a}}_{A}=\mathrm{404}$ µg, respectively, than for any other shape group. The secondlargest values are for shape group (12) Graupel, ${\stackrel{\mathrm{\u0303}}{a}}_{D}=\mathrm{53.9}$ µg and ${\stackrel{\mathrm{\u0303}}{a}}_{A}=\mathrm{124}$ µg, respectively; all other groups have still much smaller values. The smallest values are found for the five groups (1)–(3), (6), and (8). Of these, the three shape groups (1)–(3), which have the most shallow slopes, form the lower edge of the cluster of the m vs. D_{max} relationships for groups (1)–(14) (Fig. 1a), Similarly, these five groups (1)–(3), (6), and (8) form the lower edge of the cluster of all m vs. A relationships except for group (15) (Fig. 1b).
As can be seen in Fig. 1a and b, the power laws for (15) are close to the reference lines for liquid droplets, however, predicting somewhat lower masses. These differences may be due to several reasons. While shape group (15) Spherical may contain liquid droplets, it also contains ice particles that have a lower bulk density ρ_{ice} compared to the bulk density of liquid water ρ_{w}. Also, the small frozen rain droplets that shape group (15) contains are not perfectly spherical, which leads to overestimating mass if assuming a spherical shape. Furthermore, sizing errors cause an apparent error in fall speeds. Overestimating the size leads effectively to toolow v, which in turn yields toolow derived m.
4.3 Fall speed vs. particle mass
The exponent values b_{m}, i.e. the slopes of the v vs. m relationships in Fig. 1c, vary less than the slopes of the m vs. D_{max} and m vs. A relationships; they range only from 0.33 to 0.54. The shape groups with the highest slope values include group (15) as well as most of the groups that had the lowest slope values in the m vs. D_{max} and m vs. A relationships, ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ and ${\stackrel{\mathrm{\u0303}}{b}}_{A}$, respectively, i.e. groups (1)–(3) and (6). Rather than the slopes, different speeds at any given mass distinguish the different shapes. This can be seen with the values of a_{m}, representing the fall speed predicted by the relationships at the mass given by the reference unit of 1 µg. However, 1 µg is below the masses usually encountered for most shape groups. Therefore, it is more instructive to evaluate predicted fall speeds closer to the median of masses in the dataset. At a mass of for example 3 µg, the fall speeds vary between 0.14 and 0.53 m s^{−1}, as seen in Fig. 1c. The highest four fall speeds at this mass correspond to shape groups (15), (13), (3), and (12), in order of descending speed. These groups contain the most compact shapes. Contrarily, the group with the lowest speed at 3 µg, shape group (6), features the most open structures.
4.4 Comparison with previous studies
The mass vs. particle size (m vs. D) and fall speed vs. mass (v vs. m) relationships of the common shapes plates, dendrites, graupel, and spheres, i.e. for our shape groups (5) Plates, (6) Stellar (called dendrites in other studies), (12) Graupel, and (15) Spherical, respectively, are compared to previously published relationships based on measurements of mass of individual particles. The parameterizations of m vs. D (see Fig. 3a–c) selected for this comparison are taken from Locatelli and Hobbs (1974) L74, Heymsfield and Kajikawa (1987) H87, Kajikawa (1989) K89, Mitchell (1996) M96, and Erfani and Mitchell (2017) E17 and are listed in Table 2. For comparison with v vs. m (see Fig. 3d) of our shape groups (12) Graupel and (15) Spherical, parameterizations of measurements by L74 (see also Table 2) and measurements of Gunn and Kinzer (1949) G49 have been selected. Relationships from this study are further referred to as VM21 and are taken from Table 1. They have been determined as described in Sect. 3.2. Figure 3 shows all these relationships. For comparison, a line for speeds determined from Eq. (2) using Re calculated from the Re vs. X relationship (Eq. 4) and X given by Eq. (3) for spherical particles having a density ρ_{w}=1 g cm^{−3} is added to the v vs. m relationships in Fig. 3. This line is referred to as Re–X.
Depending on the study, the particle size D was defined somewhat differently. For VM21, as well as for H87, E17, K89, and M96, D corresponds to D_{max}. For L74, D is the diameter of an estimated circle that has the same crosssectional area as the imaged particle.
4.4.1 Plates
We note that for plates (Fig. 3a), the m and D relationship for crystal with sectorlike branches (P1b) reported by M96 is most similar to VM21. It is similar also with respect to its slope given by the value of ${\stackrel{\mathrm{\u0303}}{b}}_{D}$; all other relationships are steeper with larger values of ${\stackrel{\mathrm{\u0303}}{b}}_{D}$. Reasonably close to VM21 is also the relationship for hexagonal plates by M96, which, however, is heavier at larger sizes than about 1 mm. For those larger sizes, hexagonal plates by M96 are similar to hexagonal plates by H87, the latter having the steepest slope (with ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ even larger than 3). The relationship for P1b by H87 predicts the lightest particles below about 2 mm. At about 2 mm, it is also similar to VM21 and P1b by M96. Thick plates by H87 are heavier at larger sizes but similar to VM21 at about 0.2 mm. Our relationship VM21 for shape group (5) has a lower slope ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ than any of the other relationships from previous studies. Chen (1992) approximated hexagonal plates with spheroids and found a theoretical lower limit of 2 for ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ of plates, which the value 1.76 of VM21 seems to violate. While the selected previous studies with ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ values larger than 2 looked at particular shapes, VM21's shape group (5) represents a mixture of platelike shapes such as rimed plates, split plates, and double plates. Two of the shapes are represented with more than 40 particles, namely rimed plates (R1c) and double plates (P1o; see VázquezMartín et al., 2020), sufficient to determine their own relationships. As can be seen in Table 2 and Fig. 3a, both have steeper relationships with ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ of approximately 2.1. Double plates are composed of two plates with a small gap in between so that they almost resemble thicker plates. They are most similar to thick plates (C1h) by H87 within their size range. Most rimed plates in our dataset are thinner plates with light to moderate riming. They are most similar to hexagonal plates by M96.
4.4.2 Dendrites
For stellar particles (Fig. 3b), several m and D relationships are fairly close to VM21, for example the two relationships for P2a from H87, which has a similar slope, and K89, which has a lower slope than VM21. Relationships for another stellar particle type, P1d, are still relatively close to VM21. For example the one by M96 crosses VM21 somewhat above 1 mm and is higher for smaller particles, whereas the one by H87 is about a factor of 2 below VM21. This P1d relationship by H87 may be compared to the rimed stellars (R1d) of the same study H87. These two curves differ by about a factor of 2, with R1d having the larger mass and being very close to VM21. Riming of stellar particles adds mass without increasing their size noticeably (Erfani and Mitchell, 2017), as seen in E17, which explains the difference seen between the two mentioned relationships of H87. A similar difference is seen between the two relationships by E17 from a dataset used to study effects of riming. However, the two relationships by K89, which also feature unrimed and rimed stellar particles, respectively, do not show a significant difference. Particles included in shape group (6) Stellar of VM21 include cases of light riming. Distinguishing between unrimed and rimed stellar particles in the data of group (6) resulted in two relationships (not shown) that are both, within uncertainties, identical to the one produced from all data in shape group (6).
4.4.3 Graupel and spheres
L74 reported three m vs. D relationships for lump graupel (R4b) corresponding to three different particle densities with larger masses predicted by the relationships for higher densities. Our relationship for graupel is between L74's low and mediumdensity relationships (Fig. 3c). It is well approximated by the mass of spherical particles with a density of 0.12 g cm^{−3} (not shown in Fig. 3c), which is at the lower end of the density range reported by L74 for their medium density relationship (>0.10 to 0.25 g cm^{−3}). The relationship by H87 for lump graupel (R4b) is similar to L74's medium density. The relationship by E17 agrees also with VM21, but only around 1 mm, as their relationship has a much lower slope (${\stackrel{\mathrm{\u0303}}{b}}_{D}=\mathrm{2.16}$) than all other relationships for graupel (2.7 to 3.1). The mass of liquid water spheres $m=\frac{\mathit{\pi}}{\mathrm{6}}\cdot {\mathit{\rho}}_{\mathrm{w}}\cdot {D}^{\mathrm{3}}$ that was shown in Fig. 1a is added also to Fig. 3c as reference. Its comparison with VM21's line for shape group (15) Spherical is discussed in Sect. 4.2.2.
The v vs. m relationships from L74 (Fig. 3d) come, within their ranges, close to our relationship for shape group (12). In general, at a certain particle mass, the size and crosssectional area and thus the drag force decrease with increasing graupel particle density. This can be seen, to some extent, for the three lines by L74. However, their lines have different slopes in a way that makes them intersect with each other. Their slopes are more shallow than the relationship of VM21, and consequently they also cross that line. The slope for graupel of VM21 is more similar to that of the relationships related to spherical particles than the lines for graupel by L74. Consequently it approaches spherical particles, which represent an upper limit in speed, at a lower mass than the lines by L74.
The lines for spherical particles of G49 and Re–X are very close to each other; thus Re–X predicts well these measurements. The straight line for the shape group (15) of VM21 is at somewhat lower fall speeds below approximately 10 µg. All data but two particles in shape group (15) have m below that mass. For those two particles heavier than 10 µg the fit line VM21 overpredicts mass (see Fig. A3 in the Appendix). While VM21 represents the power law fit to our measurements of droplets and spherical and almostspherical ice particles, the two curved lines of G49 and Re–X represent only liquid droplets and, thus, an upper limit in fall speed.
This paper presents new mass and fall speed parameterizations derived from DICI measurements of natural snow, ice crystals, and other hydrometeors, covering sizes from 0.1 to 3.2 mm. Using the dataset and fall speed vs. maximum dimension and vs. crosssectional area relationships from VázquezMartín et al. (2021), where fall speeds cover ranges from 0.1 to 1.6 m s^{−1}, in this study, we have added particle masses to our dataset of measured maximum dimension, crosssectional area, and fall speed of individual particles. The calculated values of individual particle masses range from close to 0.2 to 450 µg
Mitchell (1996) presented fall speed relationships derived from power laws of crosssectional area and mass vs. maximum dimension using a relationship between Re and X. We calculate particle mass data from our measurements of maximum dimension, crosssectional area, and fall speed using the same Re–X relationship. With this new extended dataset, mass vs. maximum dimension relationships, mass vs. crosssectional area, and fall speed vs. mass, given by Eqs. (8)–(10), have been derived and studied for different particle shapes. We present the conclusions that our results led to below.

As seen in Figs. A1–A3 in Appendix A and discussed in Sect. 4.1, the data's large spread is apparent. However, when fitting m vs. D_{max}, m vs. A, and v vs. m relationships to binned data, there are high correlation coefficients for most shape groups, with values between 0.9 and 1. The only exceptions are shape groups (1) Needles, (2) Crossed needles, (3) Thick columns, (6) Stellar, and (10) Spatial plates for the m vs. D_{max} relationship with ${R}_{D}^{\mathrm{2}}\simeq \mathrm{0.7}$ as well as for the m vs. A relationship shape groups (2) and (6) with ${R}_{A}^{\mathrm{2}}\simeq \mathrm{0.8}$ and (6) with ${R}_{A}^{\mathrm{2}}\simeq \mathrm{0.5}.$ While for all other shape groups ${R}_{D}^{\mathrm{2}}$ and ${R}_{A}^{\mathrm{2}}$ are similar, for these groups with lower R^{2}, ${R}_{D}^{\mathrm{2}}$ is lower than ${R}_{A}^{\mathrm{2}}$ for all but shape group (10), for which ${R}_{A}^{\mathrm{2}}$ is lower. For v vs. m, there is a good correlation for all 15 shape groups (see Table 1). The fact that m is derived from v contributes to a stronger correlation between both quantities.

For the three shape groups related to columnar or elongated shapes, i.e. shape groups (1)–(3), width rather than length or D_{max} is more closely related to a suitable characteristic length to determine Re (see Sect. 4.2.1 and Appendix C). Consequently, mass and relationships with it are not reliable. For these shape groups, ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ is close to or smaller than 1. Additionally, contrary to expectations ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ is larger than b, and the ratio of exponents ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ to ${\stackrel{\mathrm{\u0303}}{b}}_{A}$ is too low for these groups. For most other shape groups it is similar to b, as theoretically expected. Shape groups (9) and (10) (the latter with low number of particles and low correlations in relationships) show similar limitations when comparing with b. Therefore, as long as a more suitable size parameter is not available in our dataset for these shapes, mass derived from Re for these shape groups should only be used with great caution.

On a selection of 75 simple columns from shape group (3), we have done a closure study (see Appendix C) to confirm the Re–X relationship, which is central in our method (see Sect. 3.1) and used by many other studies. For this, the widths and lengths of the columns have been determined in addition to D_{max}. From these geometric dimensions, the masses of the columns have been estimated directly. Then, from each column mass, the Best number X has been determined using Eq. (3). Thus, Re and X have been determined independently and consequently compared to the X–Re relationship given by Eq. (5). This closure showed the superiority of the characteristic length L^{*} (Jayaweera, 1971) over D_{max}, confirming that D_{max} is not suitable to calculate Re and X using Eqs. (2) and (3), respectively, for columns. The closure study also showed that using the modified Best number X^{*} (Heymsfield and Westbrook, 2010) instead of the Best number X improved the agreement. The best closure for our subset of simple columns was achieved when using both characteristic length L^{*} and modified Best number X^{*} together.

When deriving the m vs. D_{max}, m vs. A, and v vs. m relationships analytically from A vs. D_{max}, v vs. D_{max}, and v vs. A given from a suitable dataset (see Sect. 3.3), the results are equivalent to fitting to the same dataset after adding m for individual particles derived from v (see Sect. 3.1). On the one hand, fitting m vs. D_{max}, m vs. A, and v vs. m relationships to data has the advantage that the X–Re relationship from Eq. (5) can be used rather than power law approximations required for the analytical derivation of the same relationships (see B in Appendix). On the other hand, if a suitable dataset is not available, but power law relationships for A vs. D_{max}, v vs. D_{max}, and v vs. A are, the analytically derived mass relationships Eqs. (14)–(16) can be used.

The parameters ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ and ${\stackrel{\mathrm{\u0303}}{b}}_{A}$, i.e. the slopes of the m vs. D_{max} and m vs. A power laws, respectively, are highest for the shape groups (6) Stellar, (11) Spatial stellar, (12) Graupel, and (15) Spherical. For groups (12) and (15) they are close to the values expected for spheres, i.e. ${\stackrel{\mathrm{\u0303}}{b}}_{D}=\mathrm{3}$ and ${\stackrel{\mathrm{\u0303}}{b}}_{A}=\mathrm{3}/\mathrm{2}$.

The exponent values b_{m}, i.e. the slopes of v vs. m, range from 0.33 to 0.55. These b_{m} slopes do not distinguish the different shapes as seen by the ${\stackrel{\mathrm{\u0303}}{b}}_{D}$ slopes for m vs. D_{max}. Instead, different speeds at any given mass are characteristic for the different shapes, with the highest fall speed for (15) Spherical and the lowest for (6) Stellar, which has shapes with open structures.

We compared our m vs. D_{max} and v vs. m relationships with other mass relationships given by previous studies. The shape groups compared in this study are (5) Plates, (6) Stellar, (12) Graupel, and (15) Spherical. Our results agree reasonably well with the references used.

For graupel and spheres (Sect. 4.4.3), Locatelli and Hobbs (1974) L74 reported m vs. D relationships for lump graupel (R4b) with different particle densities (high, medium, and low). Our relationship for graupel is between L74's low and mediumdensity relationships, and it is well approximated by the mass of spherical particles with a density of 0.12 g cm^{−3} (not shown in Fig. 3c).

Looking at v vs. m, the two lines for spherical particles of G49 and Re–X, corresponding to a line for speeds determined from Eq. (2) using Re from Eqs. (3) and (4) for spherical particles with density ρ_{w}=1 g cm^{−3}, are very close to each other. We report somewhat lower speeds for the shape group (15) Spherical VM21. This difference may be due to shape group (15) in VM21 consisting of any spherical or almostspherical particle, including ice, whereas the two lines of G49 and Re–X are exclusively for liquid droplets.
These resulting parameterizations may improve our understanding of precipitation in cold climates and improve the microphysical parameterizations in the climate and forecast models. Through these relationships, we can determine particle masses based on fall speed and particle sizes.
The particle mass relationships are derived analytically from a relationship between the Reynolds and Best numbers, in addition to A vs. D_{max}, v vs. D_{max}, and v vs. A power laws given by Eqs. (11a)–(13a). Section 3 has briefly presented this approach of deriving the particle mass analytically. The m vs. D_{max}, m vs. A, and v vs. m relationships given by this approach are equivalent to fitting to individual data. Indeed we get identical results in the ${\stackrel{\mathrm{\u0303}}{a}}_{D}$, ${\stackrel{\mathrm{\u0303}}{b}}_{D}$, ${\stackrel{\mathrm{\u0303}}{a}}_{A}$, ${\stackrel{\mathrm{\u0303}}{b}}_{D}$, a_{m}, and b_{m} parameters if using X^{*} vs. Re as power law
where γ and δ are the parameters in the power law. We determine these parameters by fitting the power law to Eq. (5) over ranges of Re corresponding to each shape group. For this, we first calculate Re for all particles in a shape group and determine X^{*} using Eq. (5) for this set of Re values. Then, we do a linear fit to the logarithm of X^{*} vs. logarithm of Re. Consequently, for each shape group, we get one set of γ and δ.
We express Re as a power law in D_{max} using Eq. (2) and replacing v with the power law given by Eq. (12a):
Now we can determine the particle mass m using Eq. (7) and express it as a function of particle size D_{max}, area A, or fall speed v. Consequently, the mass relationship as a function of particle size D_{max} given by Eq. (8) can be derived as follows (using Eqs. B1, B2, 11a, and the area ratio ${A}_{\mathrm{r}}=\frac{A}{\frac{\mathit{\pi}}{\mathrm{4}}\cdot {D}_{\text{max}}^{\mathrm{2}}}$):
The mass relationship as a function of A given by Eq. (9) can be derived as follows (using Eqs. B1, 11b, and 13a and expressing Re as a power law in A):
The mass relationship as a function of v given by Eq. (10) can be derived as follows (using Eqs. B1, 12b, and 13b and expressing Re as a power law in v):
From Eq. (B5), we can determine v(m) as follows:
Selecting a simple shape with area ratio noticeably below 1, we can test if the modified Best number approach by Heymsfield and Westbrook (2010) yields better results than using Best numbers and the approach by Böhm (1989) (see Sect. 3.1 for details about these approaches). If Reynolds and Best numbers (Re and X, respectively) can be determined independently, i.e. without using any X–Re relationship, then they can be compared with the X–Re relationship by Böhm (1989) given by Eq. (5). Thus, it represents a closure study that can confirm the X–Re relationship, which this and other studies rely on. The following explains the method and results when applied to simple columns, a small subset of our data.
For a simplegeometry shape we can calculate the particle mass from the geometrical dimensions and, thus, determine both X and Re independently. Then X vs. Re, or alternatively X^{*} vs. Re, can be compared to the empirical relationship given by Eqs. (4) or (5). Needles or columns would be suitable shapes as they have low area ratios and a simple geometry. Looking at particles in the shape group (1) Needles reveals that it contains many bundles of needles and only few pristine needles. Shape group (3) Thick columns, on the other hand, contains many simple columns. Therefore, we have selected 75 columns from shape group (3) for this comparison study. Figure C1 shows examples of the selected columns.
Most columns fall horizontally so that width and length can be easily determined from the topview images. We estimate that the length may be underestimated on the order of up 15 % due to deviations from alignment of the column axis in the image plane. On the other hand, the geometrically determined mass, m_{geom}, may be overestimated for part of the columns that show signs of cavities or hollowing of faces (see Fig. C1).
For columns, D_{max}, which is similar to the column's length, is not a suitable representative size parameter to determine Re, as we discuss in Sect. 4.2.1 and VázquezMartín et al. (2021). A characteristic length ${L}^{*}={A}_{\mathrm{t}}/P$, where A_{t} is the total surface area and P the perimeter of the particle projected to the flow (see Eq. 1381 in Pruppacher and Klett, 2010), can be used instead. For columns, L^{*} can be determined from width and length (Jayaweera, 1971) and is more closely related to the width. Now, Re can be determined from measured fall speed and L^{*}. The Best number, according to Eq. (3), can be determined from measured crosssectional area A and D_{max}. Note that here, D_{max} represents the same size parameter best suited to calculate Re as in Eq. (2). Thus, not only Re should be determined from L^{*} (instead of D_{max}) but also X. Then, X can be determined from measured A in addition to calculated m_{geom} and L^{*}.
Consequently, X vs. Re can be plotted and compared to the X–Re relationship (Eq. 5). Figure C2 shows X vs. Re determined using either D_{max} or L^{*}. The points related to D_{max} (blue triangles) do not match well the empirical relationship X–Re (with δ_{0}=5.83 and C_{0}=0.6) by Böhm (1989) based on a theoretical treatment by Abraham (1970). This confirms that, as argued above, D_{max} is not suitable to determine Re or X for this shape. The points X vs. Re determined using L^{*}, on the other hand, are much closer to the empirical relationship. Thus, this closure experiment comparing independently determined Re and X to the X–Re relationship demonstrates the superiority of characteristic length L^{*} over D_{max} as a particle size parameter when dealing with particle mass m.
The points X vs. Re can be transformed into X^{*} vs. Re, where ${X}^{*}=X\cdot {A}_{\mathrm{r}}^{\mathrm{1}/\mathrm{2}}$ is the modified Best number suggested by Heymsfield and Westbrook (2010). The resulting points (using L^{*}) are also shown in Fig. C2 (green “x”) and provide an even better closure to the empirical X–Re relationship. Heymsfield and Westbrook (2010) used D_{max} and not characteristic length L^{*} (they focused on shapes with open geometries for which characteristic length is difficult to determine). The closure for our columns using X^{*} and D_{max} (cyan crosses in Fig. C2) is not as good as using X^{*} and L^{*} but still better than using the unmodified Best number X and D_{max}. Thus, for columns we can conclude that the modified Best number represents an improvement over the Best number. In addition, the superiority of characteristic length L^{*} over D_{max} for columns is given also when working with the modified Best number X^{*}.
In addition to the empirical relationship X–Re by Böhm (1989), also the relationship by Heymsfield and Westbrook (2010) (δ_{0}=8.0 and C_{0}=0.35) for their the modified Best number approach, used in our study, is shown in Fig. C2. The two lines are relatively close to each other. Thus, the above conclusions of superiority of characteristic length L^{*} over D_{max} and improvement when using a modified Best number rather than the Best number remain valid regardless of which relationship is used as comparison.
Latin letters  
A  Crosssectional area 
a  Coefficient in the A vs. D_{max} relationship 
a^{′}  Coefficient in the D_{max} vs. A relationship 
a_{A}  Coefficient in the v vs. A relationship 
${a}_{A}^{\prime}$  Coefficient in the A vs. v relationship 
${\stackrel{\mathrm{\u0303}}{a}}_{A}$  Coefficient in the m vs. A relationship 
a_{D}  Coefficient in the v vs. D_{max} relationship 
${a}_{D}^{\prime}$  Coefficient in the D_{max} vs. v relationship 
${\stackrel{\mathrm{\u0303}}{a}}_{D}$  Coefficient in the m vs. D_{max} relationship 
a_{m}  Coefficient in the v vs. m relationship 
A_{t}  Total surface area 
b  Exponent in the A vs. D_{max} relationship 
b^{′}  Exponent in the D_{max} vs. A relationship 
b_{A}  Exponent in the v vs. A relationship 
${b}_{A}^{\prime}$  Exponent in the A vs. v relationship 
${\stackrel{\mathrm{\u0303}}{b}}_{A}$  Exponent in the m vs. A relationship 
b_{D}  Exponent in the v vs. D_{max} relationship 
${b}_{D}^{\prime}$  Exponent in the D_{max} vs. v relationship 
${\stackrel{\mathrm{\u0303}}{b}}_{D}$  Exponent in the m vs. D_{max} relationship 
b_{m}  Exponent in the v vs. m relationship 
C_{0}  Unitless constant in the Re vs. X relationship 
C_{D}  Drag coefficient 
g  Acceleration of gravity 
L^{*}  Characteristic length 
m  Particle mass 
P  Perimeter of projected particle image 
${R}_{A}^{\mathrm{2}}$  Correlation coefficient in the m vs. A relationship 
${R}_{D}^{\mathrm{2}}$  Correlation coefficient in the m vs. D_{max} relationship 
${R}_{m}^{\mathrm{2}}$  Correlation coefficient in the v vs. m relationship 
Re  Reynolds number 
v  Fall speed 
X  Best number 
Greek letters  
δ  Exponent in the X^{*} vs. Re relationship 
δ_{0}  Unitless constant in Re vs. X relationship 
η  Dynamic viscosity of air 
γ  Coefficient in the X^{*} vs. Re relationship 
ρ_{a}  Air density 
ρ_{w}  Density of liquid water 
The presented data will be available at the Swedish National Data Service (SNDID: 2021125 Version 2; https://doi.org/10.5878/2dnq5x15, Kuhn and VázquezMartín, 2021).
TK and SVM performed the conceptualization; TK prepared the resources and the instrumentation; SVM and TK performed the experiments and data collection; SVM and TK prepared the formal analysis; SVM and TK carried out the data curation; SVM prepared the original draft; SVM, TK, and SE contributed to changes and writing during review and revisions; SVM prepared the visualization that includes tables and figures; TK carried out the closure study in Appendix C and prepared Figs. C1 and C2; TK and SE carried out the supervision of the research project.
The contact author has declared that neither they nor their coauthors have any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We would like to thank the Swedish Institute of Space Physics (IRF) at Kiruna for offering its facilities for our instrument. In addition we thank David L. Mitchell and Carl G. Schmitt for providing valuable feedback and suggestions that helped us to improve our publication.
This research has been supported by the Graduate School of Space Technology at the Luleå University of Technology.
This paper was edited by Timothy Garrett and reviewed by David Mitchell and Carl Schmitt.
Abraham, F. F.: Functional Dependence of Drag Coefficient of a Sphere on Reynolds Number, Phys. Fluids, 13, 2194–2195, https://doi.org/10.1063/1.1693218, 1970. a, b
Baum, B. A., Yang, P., Heymsfield, A. J., Schmitt, C. G., Xie, Y., Bansemer, A., Hu, Y.X., and Zhang, Z.: Improvements in Shortwave Bulk Scattering and Absorption Models for the Remote Sensing of Ice Clouds, J. Appl. Meteorol. Clim., 50, 1037–1056, https://doi.org/10.1175/2010JAMC2608.1, 2011. a
Böhm, H. P.: A General Equation for the Terminal Fall Speed of Solid Hydrometeors, J. Atmos. Sci., 46, 2419–2427, https://doi.org/10.1175/15200469(1989)046<2419:AGEFTT>2.0.CO;2, 1989. a, b, c, d, e, f, g, h, i
Chen, J.P.: Numerical Simulations on the Redistribution of Atmospheric Trace Chemicals Through Cloud Processes, PhD thesis, The Pennsylvania State University, United States, 1992. a
Cooper, S. J., Wood, N. B., and L'Ecuyer, T. S.: A variational technique to estimate snowfall rate from coincident radar, snowflake, and fallspeed observations, Atmos. Meas. Tech., 10, 2557–2571, https://doi.org/10.5194/amt1025572017, 2017. a
Erfani, E. and Mitchell, D. L.: Growth of ice particle mass and projected area during riming, Atmos. Chem. Phys., 17, 1241–1257, https://doi.org/10.5194/acp1712412017, 2017. a, b, c
Gunn, R. and Kinzer, G. D.: The terminal velocity of fall for water droplets in stagnant air, J. Meteorol., 6, 243–248, https://doi.org/10.1175/15200469(1949)006<0243:TTVOFF>2.0.CO;2, 1949. a, b
Heymsfield, A. J. and Kajikawa, M.: An Improved Approach to Calculating Terminal Velocities of Platelike Crystals and Graupel, J. Atmos. Sci., 44, 1088–1099, https://doi.org/10.1175/15200469(1987)044<1088:AIATCT>2.0.CO;2, 1987. a, b, c
Heymsfield, A. J. and Westbrook, C. D.: Advances in the Estimation of Ice Particle Fall Speeds Using Laboratory and Field Measurements, J. Atmos. Sci., 67, 2469–2482, https://doi.org/10.1175/2010JAS3379.1, 2010. a, b, c, d, e, f, g, h, i, j, k, l
Jayaweera, K. O. L. F.: Calculations of Ice Crystal Growth, J. Atmos. Sci., 28, 728–736, https://doi.org/10.1175/15200469(1971)028<0728:COICG>2.0.CO;2, 1971. a, b, c, d
Kajikawa, M.: Observation of the Falling Motion of Early Snowflakes, J. Meteorol. Soc. Jpn. Ser. II, 67, 731–738, https://doi.org/10.2151/jmsj1965.67.5_731, 1989. a, b, c
Kneifel, S., Löhnert, U., Battaglia, A., Crewell, S., and Siebler, D.: Snow scattering signals in groundbased passive microwave radiometer measurements, J. Geophys. Res.Atmos., 115, D16214, https://doi.org/10.1029/2010JD013856, 2010. a
Kuhn, T. and VázquezMartín, S.: Microphysical properties and fall speed measurements of snow ice crystals using the Dual Ice Crystal Imager (DICI), Atmos. Meas. Tech., 13, 1273–1285, https://doi.org/10.5194/amt1312732020, 2020. a, b
Kuhn, T. and VázquezMartín, S.: Snow ice particle microphysical properties and fall speed from particle images taken in Kiruna (Sweden) 2014–2018, Swedish National Data Service [data set], Version 2.0, https://doi.org/10.5878/2dnq5x15, 2021.
Locatelli, J. D. and Hobbs, P. V.: Fall speeds and masses of solid precipitation particles, J. Geophys. Res., 79, 2185–2197, https://doi.org/10.1029/JC079i015p02185, 1974. a, b, c, d
Loeb, N. G., Yang, P., Rose, F. G., Hong, G., SunMack, S., Minnis, P., Kato, S., Ham, S.H., Smith, W. L., Hioki, S., and Tang, G.: Impact of Ice Cloud Microphysics on Satellite Cloud Retrievals and Broadband Flux Radiative Transfer Model Calculations, J. Climate, 31, 1851–1864, https://doi.org/10.1175/JCLID170426.1, 2018. a
Magono, C. and Lee, C. W.: Meteorological classification of natural snow crystals, J. Fac. Sci. Hokkaido Univ., Vol. II, 321–335, available at: http://hdl.handle.net/2115/8672 (last access: 20 October 2020), 1966. a
Marchand, R., Mace, G., Hallar, A., McCubbin, I., Matrosov, S., and Shupe, M.: Enhanced Radar Backscattering due to Oriented Ice Particles at 95 GHz during StormVEx, J. Atmos. Ocean. Tech., 30, 2336–2351, https://doi.org/10.1175/JTECHD1300005.1, 2013. a
Matrosov, S., Mace, G., Marchand, R., Shupe, M., Hallar, A., and McCubbin, I.: Observations of Ice Crystal Habits with a Scanning Polarimetric WBand Radar at Slant Linear Depolarization Ratio Mode, J. Atmos. Ocean. Tech., 29, 989–1008, https://doi.org/10.1175/JTECHD1100131.1, 2012. a
Mitchell, D. L.: Use of Mass and AreaDimensional Power Laws for Determining Precipitation Particle Terminal Velocities, J. Atmos. Sci., 53, 1710–1723, https://doi.org/10.1175/15200469(1996)053<1710:UOMAAD>2.0.CO;2, 1996. a, b, c, d, e, f, g, h, i
Mitchell, D. L., Rasch, P., Ivanova, D., McFarquhar, G., and Nousiainen, T.: Impact of small ice crystal assumptions on ice sedimentation rates in cirrus clouds and GCM simulations, Geophys. Res. Lett., 35, L09806, https://doi.org/10.1029/2008GL033552, 2008. a
Pruppacher, H. R. and Klett, J. D.: Microphysics of clouds and precipitation, vol. 18 of Atmospheric and Oceanographic Sciences Library, 2nd edn., Springer, Dordrecht, Heidelberg, London, New York, https://doi.org/10.1007/9780306481000, 2010. a
Schefold, R., Baschek, B., Uest, M., and Barthazy, E.: Fall velocity and axial ratio of snowflakes, Proceedings of ERAD, Copernicus GmbH, 84–89, 2002. a
Stoelinga, M. T., Hobbs, P. V., Mass, C. F., Locatelli, J. D., Colle, B. A., Houze, R. A., Rangno, A. L., Bond, N. A., Smull, B. F., Rasmussen, R. M., Thompson, G., and Colman, B. R.: Improvement of Microphysical Parameterization through Observational Verification Experiment, B. Am. Meteorol. Soc., 84, 1807–1826, https://doi.org/10.1175/bams84121807, 2003. a
Sun, W., Hu, Y., Lin, B., Liu, Z., and Videen, G.: The impact of ice cloud particle microphysics on the uncertainty of ice water content retrievals, J. Quant. Spectrosc. Ra., 112, 189–196, https://doi.org/10.1016/j.jqsrt.2010.04.003, 2011. a
Szyrmer, W. and Zawadzki, I.: Snow Studies. Part II: Average Relationship between Mass of Snowflakes and Their Terminal Fall Velocity, J. Atmos. Sci., 67, 3319–3335, https://doi.org/10.1175/2010JAS3390.1, 2010. a
Tao, W.K., Simpson, J., Baker, D., Braun, S., Chou, M.D., Ferrier, B., Johnson, D., Khain, A., Lang, S., Lynn, B., Shie, C.L., Starr, D., Sui, C.H., Wang, Y., and Wetzel, P.: Microphysics, radiation and surface processes in the Goddard Cumulus Ensemble (GCE) model, Meteorol. Atmos. Phys., 82, 97–137, 2003. a
Thornberry, T. D., Rollins, A. W., Avery, M. A., Woods, S., Lawson, R. P., Bui, T. V., and Gao, R.S.: Ice water contentextinction relationships and effective diameter for TTL cirrus derived from in situ measurements during ATTREX 2014, J. Geophys. Res.Atmos., 122, 4494–4507, https://doi.org/10.1002/2016JD025948, 2017. a
VázquezMartín, S., Kuhn, T., and Eliasson, S.: Shape Dependence of Falling Snow Crystals' Microphysical Properties Using an Updated Shape Classification, Appl. Sci., 10, 1163, https://doi.org/10.3390/app10031163, 2020. a, b, c
VázquezMartín, S., Kuhn, T., and Eliasson, S.: Shape dependence of snow crystal fall speed, Atmos. Chem. Phys., 21, 7545–7565, https://doi.org/10.5194/acp2175452021, 2021. a, b, c, d, e, f, g, h, i, j, k, l, m
Waliser, D. E., Li, J.L. F., Woods, C. P., Austin, R. T., Bacmeister, J., Chern, J., Del Genio, A., Jiang, J. H., Kuang, Z., Meng, H., Minnis, P., Platnick, S., Rossow, W. B., Stephens, G. L., SunMack, S., Tao, W.K., Tompkins, A. M., Vane, D. G., Walker, C., and Wu, D.: Cloud ice: A climate model challenge with signs and expectations of progress, J. Geophys. Res.Atmos., 114, D00A21, https://doi.org/10.1029/2008JD010015, 2009. a
Xie, Y., Yang, P., Kattawar, G., Baum, B., and Hu, Y.: Simulation of the optical properties of plate aggregates for application to the remote sensing of cirrus clouds, Appl. Optics, 50, 1065–1081, https://doi.org/10.1364/AO.50.001065, 2011. a
Yang, P., Zhang, Z., Kattawar, G. W., Warren, S. G., Baum, B. A., Huang, H.L., Hu, Y. X., Winker, D., and Iaquinta, J.: Effect of Cavities on the Optical Properties of Bullet Rosettes: Implications for Active and Passive Remote Sensing of Ice Cloud Properties, J. Appl. Meteorol. Clim., 47, 2311–2330, https://doi.org/10.1175/2008JAMC1905.1, 2008. a
 Abstract
 Introduction
 Dataset
 Methods
 Results and discussions
 Summary and conclusions
 Appendix A: Mass relationships for the shape groups
 Appendix B: Mass derivation using power laws
 Appendix C: Closure study – Reynolds and Best numbers for simple thick columns
 Appendix D: Nomenclature
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Dataset
 Methods
 Results and discussions
 Summary and conclusions
 Appendix A: Mass relationships for the shape groups
 Appendix B: Mass derivation using power laws
 Appendix C: Closure study – Reynolds and Best numbers for simple thick columns
 Appendix D: Nomenclature
 Data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References