## 1. Introduction

The retrieval of microphysics of precipitating snow from Doppler radar and other remote sensing measurements, as well as snow microphysical parameterizations, requires knowledge of the form of the size distribution that allows an accurate derivation of the distribution moments that are important for descriptions of microphysical processes. Moreover, characteristics of individual snowflakes such as representative dimensional relations of mass and velocity are needed. The study of Woods et al. (2007) demonstrated an important sensitivity of the precipitation redistribution on the changes in both mass and velocity dimensional relationships for snow in bulk microphysical schemes.

In Zawadzki et al. (2010, henceforth Part I), the study of the variability of the snowflake fall velocity, the environmental parameters controlling this variability, and the uncertainties related to the velocity measurement have been presented. In this work, based on hydrodynamic theory and the results presented in Part I, we derive an approximate relation between the mass and terminal velocity of snowflake.

According to observational, laboratory, and theoretical studies, mass has a major effect on the terminal velocity of particles. Empirical power-law relations obtained by Langleben (1954) express the snowflake fall speed in terms of its melted equivalent diameter representing the mass; the coefficients are dependent on the aggregate type. The dependence on the riming intensity related to the particle density has been introduced to other empirically derived formulas (Kajikawa 1998; Barthazy and Schefold 2006). Theoretical and laboratory work on the determination of the terminal velocity as a function of the particle mass, or density, has been mainly based on hydrodynamic theory using parameterized relationships between the Reynolds number and the Best (or Davies) number, the latter expressed in terms of mass and effective cross-sectional area. By inverting this procedure, the particle mass can be estimated from the observed terminal velocity (Hanesch 1999; Schefold 2004; Lee et al. 2008). The same approach is used here as a first step to estimate the mass relation from optical spectrograph measurements of terminal velocity as a function of snowflake size. In the next step, we determine the approximate average relation between the experimentally obtained velocity–size relationship and the estimated mass–size relationship. In this way we obtain a consistent parameterization of velocity–size and mass–size relationships.

The cross-sectional area included in the calculation of the mass–velocity relation based on the Best number *X* and Reynolds number Re (*X*–Re) relationship is expected to be related the snowflake effective density as suggested by Cunningham (1978). More recently, an empirical average expression relating particle size, density, and effective area has been proposed by Heymsfield et al. (2004).

Two types of uncertainties contribute to the uncertainty of the derived relations. The first type represents fluctuations in the measured data, such as the velocity or the area ratio for a given size category, used as an important starting point for the calculations. The second type of uncertainty arises from the fact that the theoretical formulas used for the derivation of the resulting relation are not well known. These two types of uncertainties are combined when investigating the uncertainties of the derived relations.

The measurements that we use are from an optical disdrometer, which can give information not only on velocity and area for each size bin of snowflake but also on snowflake particle size distribution (PSD). Measured PSDs were used to evaluate the retrieved relation through the comparison of the expected and a measured radar reflectivity time series and also to derive some average relations between PSD bulk quantities.

We begin with the description of the experimental data used as the starting point of our calculations. In section 3, the theoretical basis and the main assumptions of our method to retrieve mass–velocity relations are described; the calculation results and uncertainties estimation are done in section 4. Section 5 is devoted to the validation of the proposed relation. Some useful applications can be found in section 6. Finally, section 7 is a summary of the results and a discussion of some of the limitations.

## 2. Experimental data used in the study

In this study we use a large dataset of snow measurements collected by a ground-based optical disdrometer [Hydrometeor Velocity and Shape Detector (HVSD)] (Barthazy et al. 2004). Measurements were taken at the Centre for Atmospheric Research Experiments (CARE) site (80 km north of Toronto) during winter 2005/06 as part of the Canadian CloudSat–Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) Validation Project (C3VP) described by Hudak et al. (2006). The HVSD measurements provide particle size and terminal fall speed for each size class [assuming negligible vertical air velocity at the height of the HVSD measurement]. The study in Part I gives a detailed description of the measurements and investigates the variability in the velocity–size power-law coefficients. It shows that three environmental parameters—surface temperature, echo top temperature and, in particular, the depth of the precipitation system—determine the fall velocity fairly well. These factors express the impact of different active microphysical processes and of the time of growth.

In this study, nine snowfall events have been selected from the same database as in Part I by inspecting the vertically pointing X-band radar (VertiX) records and retaining only the snow systems uniform in time (for more details, see Part I). An additional requirement for the present study was the availability of well-sampled time series of the particle size distributions measured by the HVSD. For each snow event and each size bin, the mean values of velocity and area ratio, as well as their standard deviations, were obtained from Part I (see Table 1 for the summary of the analyzed events). We consider that during each event we have the same dominant type of particles characterized by the single mass and area relations. Values of the measured velocity and area ratio that deviate by more than 2 standard deviations from the average value for a given size class by the obtained relations were discarded as outliers.

The retrieved mass–velocity relationship is evaluated by comparing the time series of the reflectivity factor calculated from the mass–size relationship applied to the size distribution measured by the HVSD to the reflectivity obtained from the collocated small X-band bistatic Doppler radar providing measurements a few meters above ground [Precipitation Occurrence Sensor System (POSS); Sheppard 1990; Sheppard and Joe 2000]. The HVSD and POSS data are averaged over 6-min periods. The tests of the sensitivity of the results to the time averaging of the HVSD and POSS data have been done. Only very small differences in the results have been obtained for different averaging periods. The total number of analyzed spectra is 805.

The temperature at the ground during these events varied between −17° and −2°C. The observed reflectivity factor varied between about −10 and 30 dB*Z*.

## 3. Theoretical basis and main assumptions

### a. Functional forms of the used dimensional relationships

*D*represents a reference size of snowflake.

The coefficients *a _{u}* and

*b*have been experimentally derived separately for each snow event. The range of the obtained

_{u}*a*,

_{u}*b*values underlines the variability between the events as shown in Part I. For a given snow event, the derived

_{u}*a*and

_{u}*b*, are assumed to be size independent. However, analytical relations such as those developed by Mitchell and Heymsfield (2005, henceforth MH05) or Khvorostyanov and Curry (2005, henceforth KC05) describe

_{u}*a*and

_{u}*b*for different size ranges. These variations with size are not very important for sizes that mostly contribute to the total mass for a given particle type having the same mass and area dimensional relationships. The calculated mean mass-weighted fall speed with constant coefficients used to characterize the particle sizes, where the major part of the mass is located, is only slightly different from the calculated fall speed with size-range-dependent coefficients. The calculations performed by McFarquhar and Black (2004) for different particle types also lead to the above conclusion.

_{u}The parameters *a _{m}*,

*b*in the mass–size relation (1b) are, in general, empirically derived from observations for different types of particles of different size ranges (e.g., Mitchell et al. 1990) and at different temperatures (Heymsfield et al. 2007). Our work aims at evaluating the coefficients

_{m}*a*,

_{m}*b*and relating them to the

_{m}*a*,

_{u}*b*coefficients in (1a) in the range of snowflake sizes that dominate the total snow mass or reflectivity factor.

_{u}*A*, defined as the area of the particle image

_{r}*A*divided by that of a circle of diameter

*D*:where

*D*

_{eq}is area equivalent diameter. Two forms have been proposed for the parameterization of

*A*for the side-view particle images, depending on the snowfall event:where

_{r}*a*,

_{r}*b*,

_{r}*a*

_{re}, and

*b*

_{re}are constants determined empirically. The form selected to describe each snow event is the one with a greater correlation.

For nonspherical particles with complex structures such as ice crystals or snowflakes, the definition of size *D* is an issue in itself. This size is in general determined from particle image recorded with two-dimensional imaging probes; because it is noncircular, there are different possibilities to define the image size. Moreover, depending on the observation method, the image may refer to the side view when the image is projected on the vertical plane parallel to the flow (e.g., Barthazy et al. 2004) or the image projected on the horizontal plane normal to the flow (e.g., Mitchell et al. 1990) In investigations reported in the literature the dimension *D* has been considered in various ways—for example, as the diameter of an area-equivalent circle (i.e., the diameter that corresponds to the area as the particle shadow; Locatelli and Hobbs 1974; Francis et al. 1998), the equivalent spherical volume diameter (Brandes et al. 2007), the maximum diameter of the particle image (Kajikawa 1989; Mitchell et al. 1990), the mean of the two orthogonal extensions (Brown and Francis 1995), and the width of the enclosing box (Barthazy et al. 2004). For the mass relation a new approach is proposed by Baker and Lawson (2006), who introduce a size parameter that is a combination of the maximum dimension, width, area, and perimeter. A combination of parameters is also used by Hanesch (1999) in her study of the snowflake fall speed. In addition, the equivalent melted diameter is considered as a snowflake size parameter.

In the data collected by the HVSD the observed snowflake dimensions correspond to the two-dimensional side-view pattern. As a reference, the dimension *D* is chosen to be the maximum side-view size—that is, the maximum of the two perpendicular extensions: height of the image (vertical dimension as the snowflake falls) and width of the image (see Part I). This definition of the snowflake reference size is used in our dimensional relationships obtained from measurements of the velocity and area ratio and retrieved for mass. The dimensions normal to the flow required for the hydrodynamic calculations have to be estimated from the side-view projection.

### b. Introduction of the Reynolds and Best (Davies) numbers

*F*with the difference between the gravitational force

_{D}*mg*and the buoyancy force. However, the problem with applying this expression is the dependence of the drag coefficient

*C*on the fall speed itself. Thus, the common method for calculating the terminal velocity is from a determined relationship between the Reynolds number Re and the Best number

_{D}*X*related to the drag coefficient but having no dependence on fall speed. The two numbers are defined as (e.g., List and Schemenauer 1970)taking

*ρ*

_{snow}−

*ρ*≈

_{a}*ρ*

_{snow}, the following expression for

*X*is obtained:The index ⊥ describes the quantities related to the horizontal projection (perpendicular to the flow), while the index ‖ (introduced in Eq. 7a below) represents the projection on the vertical plane (parallel to the flow).

*ν*and

*ρ*, respectively;

_{a}*g*denotes gravitational acceleration. The expression for

*X*contains the two primary variables that determine

*u*(i.e., the mass

*m*and the effective particle area projected normal to the flow

*A*

_{⊥}, in addition to

*D*

_{*}, which denotes the chosen characteristic dimension of the particle). Introducing the area ratio normal to the flow,

*A*

_{r⊥}, and using (2),

*X*can be expressed aswhere

*D*

_{⊥}denotes the maximum diameter in the direction normal to the flow that is not measurable and therefore must be estimated. The final expression for

*X*depends on the choice of the characteristic size

*D*

_{*}. The symbols

*X*

_{‖}and

*X*

_{⊥}below describe

*X*calculated with

*D*

_{*}=

*D*(our measured maximum side-view diameter taken as reference) and

*D*

_{*}=

*D*

_{⊥}(estimated maximum diameter normal to the flow), respectively:The corresponding Re numbers are given byAssuming an idealized spheroidal shape for the snowflakes and taking

*D*

_{*}equal to the maximum diameter projected on the flow

*D*

_{⊥}, Böhm (1992) modified the general expression for

*X*(7b) by introducing the fourth root of the inverse of area ratio and the axial ratio

*ϕ*of this assumed oblate spheroid shape of snowflake:

### c. Relations between Reynolds and Best numbers

The equation relating the *X* to the Re numbers can be obtained empirically from direct measurements of particle velocity, cross-sectional area, and size (e.g., Knight and Heymsfield 1983; Heymsfield and Kajikawa 1987; Redder and Fukuta 1991). However, important uncertainties are associated with these measurements, mainly those of the cross-sectional area that is normal to the flow.

*X*derived from the boundary layer theory was developed by Böhm (1989) and Mitchell (1996) on the basis of the previous work of Abraham (1970), givingThe values of the two constants

*δ*

_{0}and

*C*

_{0}characterizing the boundary layer shape and thickness are 5.83 and 0.6, respectively, for ice particles. For aggregates, MH05 modified slightly the relation by adding an empirical term −

*a*

_{0}

*X*

^{b0}with

*a*

_{0}= 0.0017,

*b*

_{0}= 0.8. This additional term accounts for turbulent flow around larger aggregates so that there is better agreement with the measurements presented by Heymsfield et al. (2002b). KC05 proposed an alternate method to introduce this correction.

*f*(

*X*) from MH05 and KC05 by fitting log(

*X*) to an eighth-order polynomial of log(Re):with

*C*=

_{l}*C*

_{MH,l}for the MH05 relation, and

*C*=

_{l}*C*

_{KC,l}for the KC05 relation. The coefficients

*C*are given in Table 2. The

_{l}*X*–Re relations given in (11) are applied for the pairs (

*X*

_{‖}, Re

_{‖}) and (

*X*

_{⊥}, Re

_{⊥}). The differences between the values of

*X*retrieved using different

*X*–Re relationships are evident for the larger Re, which describe the snowflakes that have a larger contribution to the reflectivity values, as seen in Fig. 1, which gives a good estimate of the uncertainty associated with the relationship

*X*–Re.

### d. Assumed relations between side-view snowflake projection and the section perpendicular to the flow

*A*

_{r⊥}≈

*A*(i.e., the area ratio is independent of the angle of observation, normal or parallel to the flow) and 2) we use the relation from Schefold (2004) giving the ratio of the area projected normally to the flow to the area from the side view, given as a function of the canting angle

_{r}*α*and the side projected axial ratio

*ϵ*. The latter is the quotient of the side projected minor axis to the side projected major axis. Thus,where

*ϵ*

_{min}describes the minimal value of the axis ratio evaluated from the measurements as equal to 0.3. The last relation together with the assumption

*A*

_{r⊥}≈

*A*is also used to calculate the maximum horizontal dimension

_{r}*D*

_{⊥}from the maximum side dimension

*D*taken as the reference dimension:Equations (12) and (13) give the ratio

*D*/

*D*

_{⊥}decreasing with

*D*on average. For smaller sizes the mean value of this ratio is around 0.85 and decreases progressively. Since the value of 0.8 is consistent with the results from the snowflake fractal analysis by Schmitt and Heymsfield (2010), the minimal value of 0.75 is imposed on this ratio and the corresponding maximal value on

*f*as calculated in (12).

_{A}Our first assumption about the area ratio being independent on the projection plane seems reasonable for snow aggregates for which the value of area ratio has been shown to be related to particle density (e.g., Heymsfield et al. 2004). Moreover, this postulate agrees with the results of the recent theoretical work on the fractal dimensions of aggregates (Schmitt and Heymsfield 2010).

## 4. Derivation of the mass–velocity relationship

An average relation between a snowflake mass and velocity has been derived in the following steps with the first four steps calculated separately for each snow event:

- (i) Calculation of Re from the HVSD measurements of terminal velocity and area ratio from (8a) or (8b) and (13);
- (ii) Estimation of the value of
*X*from Re using a polynomial fit (11) with constants*C*_{MH,l}for the MH05 relation, and with*C*_{KC,l}for the KC05 relation; - (iii) Best estimation of mass for
*D*-size snowflake from the estimated*X*and area ratio normal to the flow; - (iv) Determination by regression of the coefficients in the power-law relation
*m*–*D*; and - (v) From all events, derivation by regression of an average relation between the power-law coefficients in the mass and velocity dimensional relationships.

### a. Best estimation of mass for D-size snowflake and its uncertainty

*X*estimated from Re, and with the help of the relations (12) and (13) and the approximation

*A*

_{r⊥}≈

*A*, the following three equations for mass are obtained from (7a), (7b), and (9):For each size bin, using (14a)–(14c) 24 different values of mass,

_{r}*m*with

_{k}*k*= 1, … , 24, are calculated from different combinations of the

*X*–Re relationships, smoothed and unsmoothed measured

*u*and

*A*, and two values of axial ratio of 0.7 and 0.8 put in (14c). Each value of mass is considered as the “measured” value of mass using different methods and interrelated through the common parameters. The geometric weighted mean is taken as the best estimate of the snowflake mass with reference diameter

_{r}*D*:The weight of each measured mass is the inverse of the square of the corresponding uncertainty (Δ log

*m*)

_{k}^{2}obtained by combining the contribution of the uncertainty related to the different quantities involved. To unify the calculation of (Δ log

*m*)

_{k}^{2}, we take the following single logarithmic form of Eqs. (14a)–(14c):Here

*C*is a parameter describing the environmental conditions, and is different by a factor of 4/

_{E}*π*in (14c) relative to (14a) and (14b). The uncertainty of

*C*will be neglected here.

_{E}*X*value can be decomposed into two terms. The first represents the uncertainty of the used

*X*–Re relationships (11) for different coefficients

*C*. The second describes the contribution to the uncertainty of the computed Re value:Assuming that the uncertainty in the

_{l}*X*–Re relation is 30% (as in Field et al. 2008), we have (Δ log

*X*)

_{Cl}

^{2}≈ 0.017. The same potential bias of log(

*X*) is assumed for

*X*calculated via (10). The second term on the rhs of (17) is evaluated from (11) asThe term in the {·} parentheses varies slightly with the value of Re but may be approximated by its mean value of about 1.8. The uncertainty in log(Re) is different for Re

_{‖}and Re

_{⊥}. Their values can be estimated from (8a) and (8b):Combining the relations (16)–(19), the uncertainty in log(

*m*) can be written asThe relations (20) and also (17) and (19) assume that for a single

_{k}*m*there is no correlation between the errors in the variables used in the calculation of

_{k}*m*.

_{k}The uncertainties in *u* and *A _{r}* can be estimated from the measurement standard deviation (Part I). Since we are interested in the medium-size and larger snowflakes that mostly contribute to the reflectivity, the estimated relative errors are equal to about 20% for both Δ

*u*/

*u*and Δ

*A*/

_{r}*A*. Note that these error values are incorrect for snowflakes smaller than about 3–4 mm. The potential systematic bias in the axial ratio and in the calculated value of

_{r}*f*is difficult to estimate; a ±30% range is chosen. The error in the measurement of

_{A}*D*is neglected, in which case Δ log

*D*≈ 0. Table 3 contains the summary of the estimated uncertainties used in (20).

With the assumed estimates for the different terms in (20) we obtained values of (Δ log*m _{k}*)

^{2}equal to 0.0321, 0.0363, and 0.0462 for the masses estimated from (14a), (14b), and (14c), respectively. This means that the relative error in our mass estimate is between 40% and 50% and is the same for all sizes

*D*since the relative uncertainties were taken as constant.

*m̃*(

*D*)] given by (15), the effect of the existing correlations between the many components of every couple of

*m*has to be taken into account. To do this, the uncertainty of log(

_{k}*m*) has to be separated into uncorrelated and fully correlated uncertainties for each combination of two values of

_{k}*m*. Then the propagation error law has to be applied. To avoid these complex calculations we only estimate the maximum uncertainty for which all

_{k}*m*are fully correlated between them. In this case, the total uncertainty is simply the linear sum of the weighted uncertainties of an individual

_{k}*m*:with

_{k}*m̃*) computed from (21) provides the average uncertainty of

*m*calculated from the relations (20) and is equal to 0.193.

_{k}### b. Determination by regression of the coefficients in the power-law relation m–D

*m̃*(

*D*) represents the best evaluation of the snowflake mass with reference diameter

*D*, we look for a dimensional relation

*m*–

*D*of the power-law form given by (1b) for a snowflake of a given size

*D*for each snow event. Taking the logarithm of (1b) yields the linear equationThe best estimate of the two coefficients is obtained by a least squares fit. Two examples are shown in Fig. 2, where the solid line shows the best-fit power law. For the nine snow events studied, the obtained coefficients

*a*and

_{m}*b*and their values are given in Table 1 and shown in Fig. 3. For comparison, the mean value obtained by Brandes et al. (2007) from a very large dataset is also superposed. However, because

_{m}*D*in their study was taken as equivalent diameter, their relation has been recalculated using the average relation

*D*

_{eq}–

*D*obtained from our data:

*D*

_{eq}= 0.69

*D*

^{0.95}in CGS units (very close to the average relation in Part I: 0.71

*D*

^{0.92}obtained from the larger dataset). Moreover, the sets of the parameters in the relations of mass and maximum size are also superposed after recalculation using our relation (13) between maximum side-view size and maximum horizontal dimension, taking the value for larger particles of 1/

For some of the cases studied, the investigation of the regression that gives the best set (*a _{m}*,

*b*) shows that a better description of the

_{m}*m*–

*D*relationship is obtained when the analysis is performed separately for two size regimes: for

*D*smaller than 2–3 mm, and for

*D*larger than 2–3 mm. This separation can be important in all microphysical studies where the lower-order moments of the size distribution become important.

As shown in Fig. 3, the obtained exponents *b _{m}* through the mass power law for different snow events cluster between values of 1.8 and 2.10, with the mean value close to 1.9. The value of 2.0 corresponds to an effective snowflake density decreasing with size as

*D*

^{−1}and is in good agreement with numerous observational studies of snow at the surface (Mitchell et al. 1990 from the overall observed particle types; Brandes et al. 2007) and aircraft observations (Heymsfield et al. 2002a). Moreover, most of the snowfall mass reaching the ground is generally associated with aggregates of ice crystals and theoretical studies (Westbrook et al. 2004) predict a value of 2 for the exponent in the mass–dimension power law for snow aggregates. Certainly, individual crystals with forms other than aggregates may be less well represented by this exponent; however, they contribute only a small fraction of the overall snowfall characteristics. On the other hand, an exponent equal to 2 probably does not describe correctly the snowflakes larger than about 1.5 cm because of the rather small number of the observed particles with these sizes. Contribution of these particles to the reflectivity is significantly reduced because of the non-Rayleigh effect. For example, for a snowflake with diameter of 2 cm, the reduction of backscattering cross section with respect to Rayleigh regime is about two orders of magnitude, as shown from the results of the T-matrix calculations in Matrosov et al. (2009).

*b*to 2 for the power-law relation representing the precipitating snow. The relation (22) becomeswith only the coefficient

_{m}*a*to be determined in the mass–size relationship (the symbol

_{mf}*a*denotes

_{mf}*a*with fixed

_{m}*b*= 2). All symbols used here are listed in Table 4.

_{m}On the other hand, the velocity power laws of the form (1a) obtained by Part I for different homogenous snow events show the relatively small variability of the exponent *b _{u}*. They showed that the snowflake velocity can be modeled with good accuracy with a fixed exponent

*b*= 0.18 and a varying coefficient

_{u}*a*. This fixed value of

_{u}*b*is adopted here for further calculations and the corresponding value of

_{u}*a*is represented by

_{u}*a*. Subsequently, all the calculations that lead to the best set of mass–size relationship based on the measured snowflake velocities and area ratios are repeated but with an imposed value of

_{uf}*b*= 0.18 in the fitted velocity–size relationship and

_{u}*b*= 2 in the final mass–size power law. In Table 1 the obtained values of

_{m}*a*and

_{mf}*a*are given. In Fig. 2 the dashed line show the result with imposed

_{uf}*b*= 2 using (22′), while the solid line has been obtained through regression with no restriction on the

_{m}*b*value (22). The difference is within the bounds of uncertainty.

_{m}To evaluate the impact of imposing the value of 2 for *b _{m}* we calculate, from the observed PSD for the time series of snow events, the snow ice mass content (IWC) in two ways: first, from the (

*a*,

_{m}*b*) set that was obtained by regression with

_{m}*b*varying, named IWC

_{m}_{am}, and second, using the

*a*value obtained for

_{mf}*b*= 2, named IWC

_{m}_{amf}. The root-mean-square error RMSE = 〈(IWC

_{am}− IWC

_{amf})

^{2}〉

^{1/2}= 0.004 g m

^{−3}while the root mean of fractional error RMFE = 〈[(IWC

_{am}− IWC

_{amf})/IWC

_{am}]

^{2}〉

^{1/2}= 0.07 for the entire length of the time series. The same calculation performed for the reflectivity factor gives an RMSE value of 0.35 dB; for reflectivity-weighted velocity (using fixed

*b*= 0.18), the calculated RMSE and RMFE are 0.006 m s

_{u}^{−1}and 0.007, respectively. The particles’ radar backscatter cross section is calculated based on the Mie theory assuming the mixing rule of Maxwell–Garnett for ice–air mixture [for more details, see model 5 in Fabry and Szyrmer (1999)]. We discuss the inaccuracy of the scattering calculations arising from the Mie calculations in section 5.

The estimates of the uncertainties in the coefficient *a _{mf}* and in the value of

*m*predicted by (22′) for each

*D*must take into account the effect of the correlation between calculated values of

*m̃*resulting from the contribution to their uncertainty of the correlated components of

*m*. With the assumption of a positive correlation, the magnitude of the uncertainties rises and the analytical calculations become very complex. To simplify the calculations we use the uncertainty estimate for each log

_{k}*m̃*given by the error bars in Fig. 2. By imposing a slope equal to 2, the uncertainty in

*a*is estimated by finding the range of intercept values for which the regression line covers most of the error bars representing Δ(log

_{mf}*m̃*). The obtained relative uncertainties in

*a*, denoted by Δ

_{mf}*a*/

_{mf}*a*in the last column of Table 1 for each snow event, are included in the regression calculation of the approximate relation between

_{mf}*a*and

_{mf}*a*in (23) below and shown by error bars in Fig. 4.

_{uf}### c. Derivation of average mass–velocity relationship

To reduce the number of parameters describing an individual snowflake with size *D*, we search for an approximate relation between its mass and velocity, assuming that other factors (such as the shape of the individual crystals) introduce negligible correction to the average mass–velocity relationship. This relation is obtained for the mass and velocity exponents fixed at 2.0 and 0.18, respectively.

The atmospheric conditions have some influence on the snowflake mass calculated from (14a)–(14c), and hence on the value of *a _{mf}* derived from experimental values of

*a*. This influence is mainly through the kinematic viscosity of the air

_{uf}*ν*, which is equal to the dynamic viscosity, a function of temperature, divided by the air density. The air density variation can be neglected here because all the observations were made at the ground and to a first approximation the air density changes with height (the temperature correction is evaluated at less than 2% with respect to the average temperature of −10°C). Besides the explicit dependence of

*m*on the magnitude of

*ν*given in (14), the values of the calculated Re and

*X*also change with

*ν*. Therefore, prior to the search of the mass–velocity relationship, the influence of temperature variations on the estimated value of

*a*has to be investigated.

_{mf}The snow measurements considered here were done at temperatures between −2° and −17°C. The range in dynamic viscosity corresponding to this temperature range is from 1.72 × 10^{−5} to 1.63 × 10^{−5} kg m^{−1} s^{−1}. The correction multiplication factor that normalizes the estimated mass, and therefore *a _{mf}* , was calculated with respect to calculations done for an average value of −10°C. The obtained values of the correction factor are 1.027 and 0.983 at −2° and −17°C, respectively. This range is very small and therefore suggests that the impact of the atmospheric conditions during the measurements on the derived

*a*–

_{mf}*a*relation can be neglected.

_{uf}*a*derived for

_{mf}*b*= 2 and

_{m}*b*= 0.18 are plotted as a function of the velocity coefficient

_{u}*a*obtained from the data. Each point represents one of the nine events. The analytical approximate relation between

_{uf}*a*and

_{mf}*a*is assumed to have the form log

_{uf}*a*=

_{mf}*α*+

*βa*. Coefficients

_{uf}*α*and

*β*have been determined by a weighted total least squares (WTLS) regression described in Krystek and Anton (2007). For each [

*a*, log(

_{uf}*a*)] pair, the corresponding uncertainties are Δ

_{mf}*a*≈ 0.2

_{uf}*a*from Part I and Δ(log

_{uf}*a*) obtained from regression (22′). These uncertainties are shown in Fig. 4. The solid line in Fig. 4 represents the least squares fitting with

_{mf}*α*= −3.02 and

*β*= 0.0065 in CGS units, givingThe uncertainties associated with the two coefficients are valid when the effect of correlation between the values of

*a*is not taken into account. The magnitudes of

_{mf}*a*obtained from (23) are shown in Table 1. The weighted average value of

_{mf}*a*for the nine analyzed events is shown in Fig. 4 by a dotted line. Its value is equal to 0.0044 g cm

_{mf}^{−2}and is very close to the average from observations of Mitchell et al. (1990) and Kingsmill et al. (2004) recalculated using (13) as in Fig. 3. The value of

*a*corresponding to this weighted average is equal to 102 cm

_{uf}^{0.82}s

^{−1}.

Using (23), snowflake velocities are calculated as a function of melted equivalent diameter, and the results are shown in Fig. 5. In the upper graph, the derived relations between snowflake velocity and melted diameter for each of the nine events are compared with the empirically obtained relations of Langleben (1954). In the lower graph we compare snowflake velocity calculated by our relation (23) for average *a _{mf}* with some other relations previously published (Locatelli and Hobbs 1974; Kajikawa 1998; Brandes et al. 2008 for velocity at −1° and −5°C and Brandes et al. 2007 for mass relation derived from observations mostly at temperatures higher than −5°).

All presented results are derived from measurements at the surface. Therefore, when applied at other pressure levels the velocity coefficients (*a _{u}*,

*a*) have to be adjusted for a change in altitude using, as a first approximation, the adjustment given in Heymsfield et al. (2007): [

_{uf}*a*(

_{uf}*p*)/

*a*(1000 hpa)] = [

_{uf}*ρ*

_{a0}/

*ρ*(

_{a}*p*)]

^{0.54}, where

*ρ*

_{a0}and

*ρ*(

_{a}*p*) are the air density at the ground level and at pressure

*p*.

## 5. Evaluation

In this study the measured reflectivity is used to evaluate the procedure for the retrieval of mass/density as a function of particle size. For each snowfall event, the mass coefficient *a _{mf}* is estimated from the velocity coefficient

*a*using (23). Knowing the mass–size relationship, the backscatter cross section of each size bin is calculated using the Mie theory from model 5 described in Fabry and Szyrmer (1999). For particles smaller than about 5 mm, the scattering at X-band is in the Rayleigh regime (i.e., that the reflectivity is proportional to mass squared, independent of the particle shape; e.g., Matrosov et al. 2009). For larger particles the non-Rayleigh effect leads to the reduced backscattering that decreases with size for particles larger than about 1.5 cm (Matrosov et al. 2009).

_{uf}The other question to consider is the effect of nonsphericity. According to Ishimoto (2008), using the finite difference time domain (FDTD) method for fractal-shaped snowflakes, the sensitivity to the particle shape become important for snowflakes with a diameter of the ice-mass-equivalent sphere greater than about 2.4 mm (which corresponds on the average to *D* = 1.2 cm). The same result has been obtained with T-matrix method in Wang et al. (2005). The inaccuracy of the Mie calculations assuming spherical form depends on the choice of mapping of particle properties into sphere parameters. By choosing in our model the area-equivalent diameter (i.e., between ice-mass-equivalent sphere diameter and maximum dimension) as diameter of the equivalent sphere, the inaccuracy of the scattering calculations arising from the Mie calculations for snowflakes not larger than 2 cm is relatively low. Moreover, as shown in Heymsfield et al. (2008), the contribution to the X-band reflectivity of the particles larger than about 1 cm (assuming the exponential size distribution with slope equal to 6 cm^{−1} that is representative for our dataset) is negligible. Furthermore, Matrosov et al. (2009) concluded that at X-band the spherical model for low-density snowflakes provides the necessary accuracy in the reflectivity calculations if the non-Rayleigh effects are not too important.

But perhaps the most direct way of considering the effect of nonsphericity on reflectivity is by noting that observed values of differential reflectivity in precipitating snow are a few decibels, indicating that the strong effect of low density of particles decreases the influence of shape on backscattering, which confirms the results of Matrosov et al. (2009).

The expected reflectivity factor is computed from the snowflake size distributions measured by the HVSD. The time series of the calculated reflectivity are compared with the reflectivity derived from the collocated POSS. The HVSD and POSS data are averaged over 6-min intervals. The scatterplot of the reflectivity calculated for all nine events versus the POSS measured reflectivity is presented in Fig. 6. The root mean standard error is equal to 2.71 dB.

Figure 7 presents an example of the reflectivity time series. The solid line gives the POSS-measured reflectivity; the dotted–dashed line is the reflectivity calculated from the mass relationship derived for the given event of 6 January 2006. This particular event is characterized by the lowest retrieved value of *a _{mf}* (see Table 1). To provide a limit on the importance of the mass–size relation on the uncertainty in the reflectivity calculations, the dashed line in Fig. 7 shows the reflectivity calculated for the PSDs of 6 January 2006 but using the mass–size relation retrieved for the event of 9 January 2006 where the estimated

*a*was the greatest (see Table 1). The relative difference in

_{mf}*a*between these two events is about 1.5. Under the Rayleigh regime, the effective backscattering cross section is proportional to the square of the mass. Then, we have Δ

_{mf}*/*

_{m}Z_{e}*Z*≈ (1 + 1.5)

_{e}^{2}− 1 for

*Z*in mm

_{e}^{6}m

^{−3}, giving the contribution of the mass uncertainty to the calculated dB

*Z*of Δ

_{m}Z_{e,dBZ}= 10 log(1 + Δ

*/*

_{m}Z_{e}*Z*) ≈ 8 dB, as can be noted in Fig. 7. For a smaller relative uncertainty in Δ

_{e}*a*/

_{mf}*a*we can write Δ

_{mf}*/*

_{m}Z_{e}*Z*≈ 2Δ

_{e}*a*/

_{mf}*a*, from which we can deduce that Δ

_{mf}

_{m}Z_{e,dBZ}= (10/ln10)Δ

*/*

_{m}Z_{e}*Z*, equal to about 8.7 Δ

_{e}*a*/

_{mf}*a*valid under Rayleigh regime.

_{mf}## 6. Applications

The retrieved relation (23) used in conjunction with the PSD information makes it possible to derive useful relations between different bulk properties of snow PSD such as equivalent reflectivity factor *Z _{e}*, snow ice water content (IWC), liquid equivalent snowfall precipitation rate

*S*, or reflectivity- and mass-weighted velocity

*U*and

_{Z}*U*. The time series of the snowflake PSDs are provided for each analyzed snowfall case by the optical imager HVSD.

_{M}*U*and

_{Z}*Z*in dB

_{e}*Z*and log(

*a*):Note that

_{mf}*U*and

_{Z}*a*are in CGS units. Introducing a weighted mean value for

_{mf}*a*equal to 0.0044 in CGS units as shown in Fig. 4, the average

_{mf}*U*–

_{Z}*Z*relation becomesfor precipitating low-density snow. The last term in (24) introduces the estimated sensitivity of

_{e}*U*to the error in mass–dimensional relationship. The uncertainty in

_{Z}*U*resulting from the uncertainty in the value of

_{Z}*a*may be estimated as Δ

_{mf}*U*(cm s

_{z}^{−1}) ≈ 50Δ

*a*/

_{mf}*a*. The relation (24) when

_{mf}*a*is given is expected to provide a more accurate

_{mf}*U*–

_{Z}*Z*relation compared to (24′). For higher altitudes, the air density correction for

_{e}*U*has to be introduced as given at the end of section 4. The mean ratio of calculated reflectivity-weighted velocity to mass-weighted velocity is equal to 0.93.

_{Z}The second useful relationship between two bulk quantities is the power law relating the equivalent reflectivity factor *Z _{e}* and the snowfall precipitation rate

*S*:

*Z*=

_{e}*χS*, with

^{ω}*Z*in mm

_{e}^{6}m

^{−3}and

*S*in mm h

^{−1}. In Fig. 9, the calculated snow precipitation rate

*S*is plotted as a function of the calculated

*Z*in dB

_{e}*Z*. The solid line shows the best fit obtained by the WTLS method applied to the logarithmic form of the

*Z*–

_{e}*S*power law. The uncertainties taken into account for

*Z*and

_{e}*S*in WTLS calculations describe the contributions from the uncertainties in snowflake fall speed and mass represented by the uncertainties in

*a*and

_{uf}*a*, and the correlation between them. The obtained values of the constants are

_{mf}*χ*= 494 and

*ω*= 1.44.

Our *Z _{e}*–

*S*relationship is in a general agreement with some of relations previously derived empirically, semiempirically, or theoretically that are drawn in Fig. 9 [Ohtake and Henmi (1970) for snowflake consisting of spatial dendrites; Fujiyoshi et al. (1990) for 1-min and 30-min average; Matrosov (1992) for snow density of 0.02 g cm

^{−3}; Zawadzki et al. (1993) recalculated using our average relation

*S*–IWC:

*S*= 3.3IWC

^{1.05}, with IWC in g m

^{−3}; and seven relations retrieved by Huang et al. (2010)]. However, compared with the recent work of Matrosov et al. (2009), our value of

*χ*is much larger for similar values of the exponent

*ω*. This difference may be partly explained by the applied snow velocity relationship. Our measured velocity for low-density snowflakes as shown in Part I is in general lower than that obtained from the relation used in Matrosov’s work.

The increase of snowflake density (i.e., the value of *a _{mf}* in our parameterization) is related to the increasing velocity described by

*a*. For the same snowflake spectrum, this increase results in the increase of both

_{uf}*Z*(proportional to

_{e}*S*(proportional to the product

*a*×

_{mf}*a*). To estimate the impact of these changes on the coefficients in the

_{uf}*Z*–

_{e}*S*relation, we present in Fig. 10 the calculations of

*Z*and

_{e}*S*repeated 3 times for the time series of PSDs on 9 January 2006 using 1)

*a*and

_{mf}*a*obtained for this event (solid circle), 2) these two coefficients retrieved for a different event, 6 January 2006 (open circle), and 3) weighted average values of these coefficients as shown in Fig. 4 (symbol “X”). The calculation results presented in Fig. 10 show that the coefficients of the power law

_{uf}*Z*–

_{e}*S*for low-density snowfall do not exhibit a significant dependence on the assumed mass–velocity relationship if (23) is satisfied. Matrosov et al. (2009) show the slight dependence of the coefficient

*χ*on the assumed mass dimensional relationship. However, they used the density independent expression for snowflake velocity. Our result suggests that for the unrimed snowfall of our study, the observed spread in the

*Z*–

_{e}*S*power law is mainly due to the PSDs differences, as shown in Matrosov et al. (2009), and to a lesser degree to the assumed mass relation.

## 7. Summary and discussion

Both modeling and remote sensing studies require the parameterizations of bulk quantities. Generally for snow the accuracy of these calculated variables and the relations between them are dependent on the parameterization of the PSDs used and on the mass and velocity dimensional relationships as well. The variability of snowflake velocity discussed in Part I between different snowfall events is related here to the variability in the mass dimensional relationship based on the idea that the reliable estimators of snowflake velocity are its density together with its size. Even if only light riming (or no riming) cases were used here, we see appreciable differences in (23) and all the derived relationships.

Radar observables provide an alternative verification of microphysical description of snow. Our ultimate goal is to derive the microphysical relations that are consistent with radar measurements and suitable for modeling and optimized radar data assimilation. We have limited here our interest to precipitating unrimed or lightly rimed snowflake aggregates for which we derived the mass–size–velocity relationship that is a simpler expression of the general relations of mass, area, size, and velocity developed in MH05 and KC05. In our relationship the area is eliminated: the cross-sectional area, related to the details of the crystal types composing the aggregates, is closely related to the particle mass, and therefore its direct impact on the snowflake fall speed is reduced. In this way, in our formulation the variability in the fall speed for snowflakes with the same size is only attributed to particle density (or mass).

The main goal of this work was to retrieve an average relation between unrimed or slightly rimed snowflake velocity and mass using the observational data and theoretical calculations. The retrieval is through the following five steps: 1) calculation of the Re number from HVSD measurements of terminal velocity and area ratio for each snowfall event; 2) estimation of the *X* number from Re; 3) retrieval of mass for each size bin from the *X*–Re relationship and the estimated area ratio normal to the flow; 4) determination by regression of the coefficients in the power-law relation *m*–*D* for each event; and 5) from all snow events, derivation by regression of an average relation between the power-law coefficients in the mass and velocity dimensional relationships (setting the exponents to the fixed values of 2 and 0.18). This relation is representative for the snowflake sizes dominating the snow bulk quantities such as mass content or reflectivity factor. Given all other uncertainties, the variation of the mass and velocity coefficients with particle size is not taken into account. These coefficients and the derived relation are considered to be representative for the range of sizes that dominates the radar signal and where the majority of snowfall mass is located. For very light precipitating snow with an important contribution to IWC of particles smaller than about 1–2 mm, the derived relation is not valid. If applying the relations obtained in this study, it is important to account for possible differences in the definition of snowflake reference size, taken here as maximum side-view dimension.

The evaluation of the results is made by comparing the time series of the reflectivity factor calculated for a derived mass–size relationship for an individual snowflake and applied for the size distribution measured by the HVSD, with reflectivity obtained from the zeroth moment of the spectrum measured by the collocated POSS.

Using the measurements of snow velocity and retrieved snowflake mass together with the observed PSDs, different bulk quantities have been calculated and approximate average relations between them are derived. The sensitivity to the parameter *a _{mf}* or

*a*, related by (23), is significant for calculated

_{uf}*Z*,

_{e}*S*, IWC,

*U*, and

_{Z}*U*and for the derived expression

_{M}*U*–

_{Z}*Z*. However, the coefficients in the power law

_{e}*Z*–

_{e}*S*are not affected by the changes in the interdependent

*a*and

_{mf}*a*and are in a generally good agreement with the relations obtained from the literature for the same meteorological conditions. Taking into account the fact that the lower size detected by the HVSD is about 0.3 mm, the observed PSDs are truncated mainly at the lower size limit and the calculated relations may be not applicable to the light snow precipitation for which the effect of this truncation can be important.

_{uf}Because of this interdependence of snowflake mass and velocity, the coefficients in the mass and velocity dimensional relationships selected in the microphysics scheme have to be related to each other to ensure consistency between calculated quantities. The requirement of the self-consistency between snowflake mass and velocity may be partly satisfied by incorporating our retrieved relation (23) into microphysical calculations.

Heavy rimed snow is excluded from our analysis, and consequently the average relationships derived in this study may not be valid for this type of snow, as in the cases of reflectivity higher than about 25–30 dB*Z*. Our relationships may as well not be representative of other geographical conditions. The sensitivity of the proposed relations to the PSDs parameterization has to be taken into account when applying to higher levels in the atmosphere.

In Part III of our study, we will present the third important element of the bulk snow microphysics parameterization: the PSD representation. Further study will be devoted to relating the developed snow parameterization to the environmental conditions (see Part I).

The authors are grateful to EunSil Jung for providing data used in this study. The editing of the manuscript by Aldo Bellon is greatly appreciated. The authors also acknowledge the comments of Andy Heymsfield and anonymous reviewer in improving and clarifying certain portions of the manuscript.

## REFERENCES

Abraham, F. F., 1970: Functional dependence of drag coefficient of a sphere on Reynolds number.

,*Phys. Fluids***13****,**2194–2195.Baker, B., , and R. P. Lawson, 2006: Improvement in determination of ice water content from two-dimensional particle imagery. Part I: Image-to-mass relationships.

,*J. Appl. Meteor. Climatol.***45****,**1282–1290.Barthazy, E., , and R. Schefold, 2006: Fall velocity of snowflakes of different riming degree and crystal types.

,*Atmos. Res.***82****,**391–398.Barthazy, E., , S. Göke, , R. Schefold, , and D. Högl, 2004: An optical array instrument for shape and fall velocity measurements of hydrometeors.

,*J. Atmos. Oceanic Technol.***21****,**1400–1416.Böhm, H., 1989: A general equation for the terminal fall speed of solid hydrometeors.

,*J. Atmos. Sci.***46****,**2419–2427.Böhm, J. P., 1992: A general hydrodynamic theory for mixed-phase microphysics. Part I: Drag and fall speed of hydrometeors.

,*Atmos. Res.***27****,**253–274.Brandes, E. A., , K. Ikeda, , G. Zhang, , M. Schönhuber, , and R. Rasmussen, 2007: A statistical and physical description of hydrometeor distributions in Colorado snow storms using a video disdrometer.

,*J. Appl. Meteor. Climatol.***46****,**634–650.Brandes, E. A., , K. Ikeda, , G. Thompson, , and M. Schönhuber, 2008: Aggregate terminal velocity/temperature relations.

,*J. Appl. Meteor. Climatol.***47****,**2729–2736.Brown, P. R. A., , and P. N. Francis, 1995: Improved measurements of the ice water content in cirrus using a total-water probe.

,*J. Atmos. Oceanic Technol.***12****,**410–414.Cunningham, M. R., 1978: Analysis of particle spectral data from optical array (PMS) 1D and 2D sensors. Preprints,

*Fourth Symp. on Meteorological Observations and Instrumentation,*Denver, CO, Amer. Meteor. Soc., 345–349.Fabry, F., , and W. Szyrmer, 1999: Modeling of the melting layer. Part II: Electromagnetics.

,*J. Atmos. Sci.***56****,**3593–3600.Field, P. R., , J. Heymsfield, , A. Bansemer, , and C. H. Twohy, 2008: Determination of the combined ventilation factor and capacitance for ice crystal aggregates from airborne observations in a tropical anvil cloud.

,*J. Atmos. Sci.***65****,**376–391.Francis, P. N., , P. Hignett, , and A. Macke, 1998: The retrieval of cirrus cloud properties from aircraft multi-spectral reflectance measurements during EUCREX’93.

,*Quart. J. Roy. Meteor. Soc.***124****,**1273–1291.Fujiyoshi, Y., , T. Endoh, , T. Yamada, , K. Tsuboki, , Y. Tachibana, , and G. Wakahama, 1990: Determination of a

*Z*–*R*relationship for snowfall using a radar and high sensitivity snow gauges.,*J. Appl. Meteor.***29****,**147–152.Hanesch, M., 1999: Fall velocity and shape of snowflakes. Ph.D. thesis, Swiss Federal Institute of Technology, 117 pp.

Heymsfield, A. J., , and M. Kajikawa, 1987: An improved approach to calculating terminal velocities of plate-like crystals and graupel.

,*J. Atmos. Sci.***44****,**1088–1099.Heymsfield, A. J., , A. Bansemer, , P. R. Field, , S. L. Durden, , J. L. Stith, , J. E. Dye, , W. Hall, , and C. A. Grainger, 2002a: Observations and parameterizations of particle size distributions in deep tropical cirrus and stratiform precipitating clouds: Results from in situ observations in TRMM field campaigns.

,*J. Atmos. Sci.***59****,**3457–3491.Heymsfield, A. J., , S. Lewis, , A. Bansemer, , J. Iaquinta, , L. M. Miloshevich, , M. Kajikawa, , C. Twohy, , and M. R. Poellot, 2002b: A general approach for deriving the properties of cirrus and stratiform ice cloud particles.

,*J. Atmos. Sci.***59****,**3–29.Heymsfield, A. J., , A. Bansemer, , C. Schmitt, , C. Twohy, , and M. R. Poellot, 2004: Effective ice particle densities derived from aircraft data.

,*J. Atmos. Sci.***61****,**982–1003.Heymsfield, A. J., , A. Bansemer, , and C. Twohy, 2007: Refinements to ice particle mass dimensional and terminal velocity relationships for ice clouds. Part I: Temperature dependence.

,*J. Atmos. Sci.***64****,**1047–1067.Heymsfield, A. J., , A. Bansemer, , S. Matrosov, , and L. Tian, 2008: The 94-GHz radar dim band: Relevance to ice cloud properties and CloudSat.

,*Geophys. Res. Lett.***35****,**L03802. doi:10.1029/2007GL031361.Huang, G. J., , V. N. Bringi, , R. Cifelli, , D. Hudak, , and W. A. Petersen, 2010: A methodology to derive radar reflectivity–liquid equivalent snow rate relations using C-band radar and a 2D video disdrometer.

,*J. Atmos. Oceanic Technol.***27****,**637–651.Hudak, D., , H. Barker, , P. Rodriguez, , and D. Donovan, 2006: The Canadian CloudSat Validation Project.

*Proc. Fourth European Conf. on Radar in Hydrology and Meteorology,*Barcelona, Spain, ERAD, 609–612.Ishimoto, H., 2008: Radar backscattering computations for fractal-shaped snowflakes.

,*J. Meteor. Soc. Japan***86****,**459–469.Kajikawa, M., 1989: Observation of the falling motion of early snowflakes. Part II: On the variation of falling velocity.

,*J. Meteor. Soc. Japan***67****,**731–738.Kajikawa, M., 1998: Influence of riming on the fall velocity of dendritic snow crystals.

,*J. Fac. Sci. Hokkaido Univ. Ser. 7***11****,**169–174.Khvorostyanov, V. I., , and J. A. Curry, 2005: Fall velocities of hydrometeors in the atmosphere: Refinements to a continuous power law.

,*J. Atmos. Sci.***62****,**4343–4357.Kingsmill, D., and Coauthors, 2004: TRMM common microphysics products: A tool for evaluating spaceborne precipitation retrieval algorithms.

,*J. Appl. Meteor.***43****,**1598–1618.Knight, N. C., , and A. Heymsfield, 1983: Measurements and interpretation of hailstone density and terminal velocity.

,*J. Atmos. Sci.***40****,**1510–1516.Krystek, M., , and M. Anton, 2007: A weighted total least-squares algorithm for fitting a straight line.

,*Meas. Sci. Technol.***18****,**3438–3442.Langleben, M. P., 1954: The terminal velocity of snowflakes.

,*Quart. J. Roy. Meteor. Soc.***80****,**174–181.Lee, G. W., and Coauthors, 2008: Snow microphysical processes and variation of effective density–diameter relationships.

*Proc. Fifth European Conf. on Radar in Meteorology and Hydrology,*Helsinki, Finland, ERAD, 7.2.List, R., , and R. S. Schemenauer, 1970: Free-fall behavior of planar snow crystals, conical graupel, and small hail.

,*J. Atmos. Sci.***28****,**110–115.Locatelli, J. D., , and P. V. Hobbs, 1974: Fall speeds and masses of solid precipitation particles.

,*J. Geophys. Res.***79****,**2185–2197.Magono, C., , and T. Nakamura, 1965: Aerodynamic studies of falling snowflakes.

,*J. Meteor. Soc. Japan***43****,**139–147.Matrosov, S. Y., 1992: Radar reflectivity in snowfall.

,*IEEE Trans. Geosci. Remote Sens.***30****,**454–461.Matrosov, S. Y., , and A. J. Heymsfield, 2008: Estimating ice content and extinction in precipitating cloud systems from CloudSat radar measurements.

,*J. Geophys. Res.***113****,**D00A05. doi:10.1029/2007JD009633.Matrosov, S. Y., , C. Campbell, , D. Kingsmill, , and E. Sukovich, 2009: Assessing snowfall rates from X-band radar reflectivity measurements.

,*J. Atmos. Oceanic Technol.***26****,**2324–2339.McFarquhar, G. M., , and R. A. Black, 2004: Observations of particle size and phase in tropical cyclones: Implications for mesoscale modeling of microphysical processes.

,*J. Atmos. Sci.***61****,**422–439.Mitchell, D. L., 1996: Use of mass– and area–dimensional power laws for determining precipitation particle terminal velocities.

,*J. Atmos. Sci.***53****,**1710–1723.Mitchell, D. L., , and A. J. Heymsfield, 2005: Refinements in the treatment of ice particle terminal velocities, highlighting aggregates.

,*J. Atmos. Sci.***62****,**1637–1644.Mitchell, D. L., , R. Zhang, , and R. L. Pitter, 1990: Mass-dimensional relationships for ice particles and the influence of riming on snowfall rates.

,*J. Appl. Meteor.***29****,**153–163.Ohtake, T., , and T. Henmi, 1970: Radar reflectivity of aggregated snowflakes. Preprints,

*14th Radar Meteorology Conf.,*Tucson, AZ, Amer. Meteor. Soc., 209–210.Redder, C. R., , and N. Fukuta, 1991: Empirical equations of ice crystal growth microphysics for modeling and analysis. II: Fall velocity.

,*Atmos. Res.***26****,**489–507.Schefold, R., 2004: Messungen von Schneeflocken: Die Fallgeschwindigkeit und eine Abschätzung weiterer Grössen. Ph.D. thesis, Swiss Federal Institute of Technology (ETH Diss. 15431), 182 pp.

Schmitt, C. G., , and A. J. Heymsfield, 2010: The dimensional characteristics of ice crystal aggregates from fractal geometry.

,*J. Atmos. Sci.***67****,**1605–1616.Sheppard, B. E., 1990: The measurement of raindrop size distributions using a small Doppler radar.

,*J. Atmos. Oceanic Technol.***7****,**255–268.Sheppard, B. E., , and P. I. Joe, 2000: Automated precipitation detection and typing in winter: A two-year study.

,*J. Atmos. Oceanic Technol.***17****,**1493–1507.Wang, Z., , G. M. Heymsfield, , L. Li, , and A. J. Heymsfield, 2005: Retrieving optically thick ice cloud microphysical properties by using airborne dual-wavelength radar measurements.

,*J. Geophys. Res.***110****,**D19201. doi:10.1029/2005JD005969.Westbrook, C. D., , R. C. Ball, , P. R. Field, , and A. J. Heymsfield, 2004: Universality in snowflake aggregation.

,*Geophys. Res. Lett.***31****,**L15104. doi:10.1029/2004GL020363.Woods, C. P., , M. T. Stoelinga, , and J. D. Locatelli, 2007: The IMPROVE-1 storm of 1–2 February 2001. Part III: Sensitivity of a mesoscale model simulation to the representation of snow particle types and testing of a bulk microphysical scheme with snow habit prediction.

,*J. Atmos. Sci.***64****,**3927–3948.Zawadzki, I., , P. Zwack, , and A. Frigon, 1993: A study of a CASP storm: Analysis of radar data.

,*Atmos.–Ocean***31****,**175–199.Zawadzki, I., , E. Jung, , and G. W. Lee, 2010: Snow studies. Part I: A study of natural variability of snow terminal velocity.

,*J. Atmos. Sci.***67****,**1591–1604.

Summary of the characteristics of the analyzed events. All parameters are given in CGS units; times are in UTC. Solid symbols indicate events for which surface air temperature was above −10°C. The associated symbols are for further reference.

Coefficients in (11).

Assumed range of uncertainty (%) for the variables used to estimate the *m*–*D* relationship in (20).

List of symbols.