Modeling Nonspherical Hailstones

Hail research and forecasting models necessarily involve explicit or implicit}and uncertain}physical assumptions regarding hailstones’ shape, tumbling behavior, fall speed, and thermal energy transfer. Whereas most models assume spherical hailstones, we relax this assumption by using hailstone shape data from field observations to establish empirical size–shape relationships with reasonable degrees of randomness considering hailstones’ natural shape variability, capturing the observed distribution of triaxial ellipsoidal shapes. We also incorporate explicit, random tumbling of individual hailstones during their growth to simulate their free-falling behavior and the resultant changes in cross-sectional area (which affects growth by hydrometeor collection). These physical attributes are incorporated in calculating hailstones’ fall speeds, using either empirical relationships or analytical relationships based on each hailstone’s Best and Reynolds numbers. Options for drag coefficient modification are added to emulate hailstones’ rough surfaces (lobes), which then modifies their thermal energy and vapor exchange with the environment. We investigate how applying these physical assumptions about nonspherical hail to the Penn State hail growth trajectory model, coupled with Cloud Model 1 supercell simulations, impacts hail production and examine the reasons behind the resulting variability in hail statistics. The choice of hailstone size–mass relation and fall speed scheme have the strongest influence on hail sizes. Using nonspherical, tumbling hailstones reduces the number of large hailstones produced. Applying shape-specific thermal energy transfer coefficients subtly increases sizes; the effects of lobes vary depending on the fall speed scheme used. These physical assumptions, although adding complexity to modeling, can be parameterized efficiently and potentially used in microphysics schemes.