the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# On the ice-nucleating potential of warm hydrometeors in mixed-phase clouds

### Agathe Chouippe

### Markus Uhlmann

### Jan Dušek

### Thomas Leisner

The question as to whether or not the presence of warm hydrometeors in clouds may play a significant role in the nucleation of new ice particles has been debated for several decades. While the early works of Fukuta and Lee (1986) and Baker (1991) indicated that it might be irrelevant, the more recent study of Prabhakaran et al. (2020) suggested otherwise. In this work, we attempt to quantify the ice-nucleating potential using high-fidelity flow simulation techniques around a single hydrometeor and use favorable considerations to upscale the effects to a collective of ice particles in clouds. While we find that ice nucleation may be significantly enhanced in the vicinity of a warm hydrometeor and that the affected volume of air is much larger than previously estimated, it is unlikely that this effect alone causes the rapid enhancement of ice nucleation observed in some types of clouds, mainly due to the low initial volumetric ice concentration. Furthermore, it is demonstrated that the excess nucleation rate does not primarily depend on the rate at which cloud volume is sampled by the meteors' wakes but is rather limited by the exposure time of ice-nucleating particles to the wake, which is estimated to be of the order of few microseconds. It is suggested to further investigate this phenomenon by tracking the trajectories of ice-nucleating particles in order to obtain a parametrization which can be implemented into existing cloud models to investigate second-order effects such as ice enhancement after the onset of glaciation.

The formation of hydrometeors in clouds is of great importance for the prediction of weather and cloud electrification, as well as
for the hydrological cycle and, thus, eventually for the evolution of climate. However, despite its relevance and great research effort over
the past decades, many aspects remain poorly understood. One such puzzle is the discrepancy between the concentration of
ice particles and that of available
ice nuclei (IN) in airborne observations by several orders of magnitude (Pruppacher and Klett, 2010), which has been observed
for various cloud types
(Koenig, 1963; Auer et al., 1969; Hobbs, 1969; Hobbs and Rangno, 1985; Mossop, 1985; Hogan et al., 2002).
This phenomenon has been termed *ice enhancement*, and various mechanisms which amplify primary ice nucleation
have been proposed to explain the observed surplus in ice particles (Field et al., 2017) .

The most promising class of enhancement mechanisms is the so-called *secondary ice production* (SIP), whereby new ice is formed
from preexisting ice particles. Commonly accepted SIP includes rime splintering (Hallett and Mossop, 1974), fragmentation of ice
(Vardiman, 1978; Takahashi et al., 1995; Bacon et al., 1998) and freezing drops (Hobbs and Alkezweeny, 1968).
Most of these mechanisms have been implemented into cloud models with explicit ice microphysics
(see Field et al., 2017, for a recent overview), which, however, are not capable of satisfactorily explaining the
large amount of ice particles in observations. Further mechanisms have been proposed in the past
whose relative importance is still to be evaluated.

While studying the effect of supersaturation on primary ice formation, Gagin (1972) proposed several mechanisms which locally produce high values of supersaturation and, thus, regions of significantly enhanced nucleation activity. He suggested that freezing hydrometeors, which attain higher temperatures than the surrounding air due to the release of latent heat, may cause transient supersaturation by simultaneous evaporation and heat transfer and linked it to the observation of satellite drops which have previously been observed experimentally during the freezing of supercooled droplets by Dye and Hobbs (1968). This hypothesis was corroborated by Nix and Fukuta (1974), who investigated numerically the transient freezing process of an isolated droplet.

Another phenomenon which leads to localized supersaturation around hydrometeors is the riming process of ice particles in mixed-phase clouds with high liquid-water content (Gagin, 1972). Here, sufficiently large ice particles collect supercooled droplets from their surroundings, which then accumulate on its surface and subsequently freeze, leading to similar non-equilibrium conditions as observed for the freezing drop. The supersaturation around a riming ice particle has been investigated numerically by Fukuta and Lee (1986), assuming steady conditions. Indeed, they found that the air in the vicinity of a warm ice particle can be highly supersaturated, while the magnitude and spatial extent strongly increase with increasing difference in surface to ambient temperature. However, even at low cloud temperatures, the supersaturated regions do not extend far away from the hydrometeor in their simulations, which led to the conclusion of Rangno and Hobbs (1991) that the overall ice enhancement is likely to be negligible.

While the previously mentioned studies focused on the quantification of the supersaturation field, Baker (1991) attempted to quantify the actual effect on ice enhancement around warm hydrometeors. He concluded that, although significant ice enhancement factors appear to be possible, the affected air mass seems to be too small to substantially contribute to the explanation of the discrepancy between ice crystal and ice nucleus concentrations.

Recently, the potential of hail and rain to nucleate droplets gained new attention when Prabhakaran et al. (2017) observed condensation in the wake of cold droplets in a moist convection apparatus. In their experiment, a pressurized mixture of sulfur hexafluoride / helium was used, and operating conditions were chosen such that slight supersaturations led to homogeneous condensation, which could indeed be observed in the wake of larger falling droplets. However, their results can only be qualitatively transferred to Earth's atmosphere, mainly due to the difference in primary nucleation mechanism. Moreover, the mechanism to create supersaturation is different from the previously discussed one in the sense that the droplets are colder than the surrounding air.

To overcome these shortcomings, Prabhakaran et al. (2020) conducted a similar laboratory experiment using moist air which has been seeded with aerosol particles (APs) as an operating medium to investigate heterogeneous nucleation in the wake of a warm falling droplet. Again, significant nucleation was observed in the wake of the falling hydrometeor. While the extent of wake-induced nucleation has not been quantified, the affected volume of air appears to be significantly larger than the predictions that Fukuta and Lee (1986) and Baker (1991) were suggesting. This discrepancy might be related to the strong assumptions on the flow made in these early works, which neglect important features of the wake of falling objects that are nowadays accessible to numerical simulations (Johnson and Patel, 1999; Bouchet et al., 2006; Zhou and Dušek, 2015). In particular, it has been shown that both unsteadiness and vortical structures of the flow strongly affect the temperature and vapor concentration in the vicinity and far downstream of a falling sphere (Bagchi et al., 2001; de Stadler et al., 2014) and thus strongly affect the distribution of supersaturation around hydrometeors (Chouippe et al., 2019).

In Chouippe et al. (2019), we presented a framework for high-fidelity numerical simulations of heat and mass transfer around a falling ice particle which is not in thermal equilibrium with its surroundings. Even though the focus was predominantly set towards several methodological questions, it was shown qualitatively that the local supersaturation differs strongly from the simpler considerations of Fukuta and Lee (1986).

The present work aims to revisit the details of the supersaturation field around an idealized falling hydrometeor and to quantify the volume of air which is affected by the presence of hydrometeors. Furthermore, it is yet to be evaluated how significant meteor-induced enhancement of ice nucleation is in clouds (Korolev et al., 2020), and therefore, we attempt to link our results to heterogeneous nucleation of aerosol particles.

In order to assess the spatial structure of supersaturation around hydrometeors, knowledge on the flow around it is essential. This is not an easy task, since the features of the flow may strongly depend on parameters such as the size, shape and surface properties of the hydrometeor, which typically vary substantially, especially if ice particles are concerned.

In this work we utilize the numerical framework previously presented in Chouippe et al. (2019)
and model the hydrometeor as a sphere with constant diameter *D*, constant surface temperature *T*_{p} and
a constant vapor pressure *e*_{v,p} at its surface. The latter is determined by the assumption of local
equilibrium at the surface and, thus, corresponds to the saturation vapor pressure *e*_{sat,j} at *T*_{p}.
Here, *j* denotes the phase of water on the surface of the hydrometeor, where the subscript “*i*” (ice) will be used to
denote the solid phase and “*w*” (water) to denote the liquid phase.

The boundary conditions for the heat and mass transfer problem are motivated by the conditions expected
for riming ice particles in the wet-growth regime (Ludlam limit). A surface temperature of *T*_{p}=0 ^{∘}C
is assumed; i.e., the surface of the ice sphere is warmer than the local cloud temperature, since the latter is
generally below the freezing point in mixed-phase clouds. Please note that at this surface temperature
the equilibrium vapor pressure with respect to both the liquid and solid phase coincide, and hence, no assumptions on the
phase at the surface have to be made.
Since the actual riming process is not modeled explicitly
but only enters through the boundary conditions of the heat and mass transfer problem, the results of this work
are more general and can in principle be applied to any configuration which contains a mechanism leading to a warm hydrometeor,
such as a freezing droplet or rapid decreases in cloud temperature. Nonetheless, riming ice particles are likely the
most relevant application, since the riming mechanism allows for a temperature gap which is sustained and not transient.

One of the key parameters for the investigation of wake-induced ice nucleation is the distribution of the saturation ratio

where *e* denotes the local vapor pressure. A parametrization of *e*_{sat,j} with respect to
temperature for both phase equilibria in the temperature range of interest is adopted from Murphy and Koop (2005).
The spatial distribution of *S*_{j} can then be reconstructed from the temperature field *T*
and the vapor concentration *n*_{v}, which are governed by the transport equations

where 𝒟_{T} and ${\mathcal{D}}_{{n}_{\mathrm{v}}}$ denote the heat and vapor diffusivities, and $\mathit{u}=(u,v,w{)}^{T}$ is the
velocity field of the surrounding air. The partial pressure of water is linked to these transported quantities by

where *k*_{b} is Boltzmann's constant. Even though the diffusivities of both vapor and heat are similar in magnitude,
it was shown by Chouippe et al. (2019) that assuming equal diffusivities leads to an underestimation of the
saturation ratio, and hence, both equations will be treated separately.

Under the boundary conditions described earlier, Eq. (2) leads to an outward heat
flux because the surface temperature of the hydrometeor is higher than that of its surroundings. Furthermore,
at *T*_{p}=0 ^{∘}C the vapor concentration at the meteor's surface
is generally higher than in the ambient, since the latter is assumed to be in equilibrium with the liquid phase
at the corresponding cloud temperature *T*_{∞}<*T*_{p} due to the presence of
droplets in mixed-phase clouds; i.e., ${e}_{\mathrm{\infty}}={e}_{\mathrm{sat},w}\left({T}_{\mathrm{\infty}}\right)$. This results in an outward vapor flux, and hence,
the hydrometeor will also evaporate. It is therefore of primary importance to resolve the temporal and spatial variations
of both temperature and concentration fields.

In Chouippe et al. (2019) we showed that buoyant forces due to density variations within the fluid phase caused by variations in temperature and concentration are negligible at the parameter point of interest. The flow can then be approximated by the incompressible Navier–Stokes equations,

where *p* denotes the hydrodynamic pressure, *ρ*_{f} the density and *ν* the kinematic viscosity of the fluid.
The flow is driven by the falling motion of the hydrometer, which is assumed to be moving
with a constant velocity *v*_{p} equal to its terminal velocity *v*_{T}; i.e.,

where *e*_{z} is the unit vector in the *z* direction. The fluid ahead of the hydrometeor is assumed to be at rest.
This setup is equivalent to a system with a fixed sphere positioned in an
upcoming flow, which differs from that of a freely falling mobile sphere in the sense that fluctuations
in *v*_{p} are not permitted. The necessity to account for these variations was discussed in Chouippe et al. (2019),
and it was found that it has little influence in the present context due to the expected high value of the density ratio.

The Navier–Stokes and scalar transport equations are solved numerically in non-dimensional form.
In particular lengths are scaled by the particle diameter *D* and velocity components by the terminal velocity *v*_{T}, and the scalar
transport equations are formulated in terms of

which has the advantage that various
dimensional boundary conditions, e.g., various values of *T*_{∞}, can be studied in post-processing from a single simulation run.
Under the simplifications stated, the non-dimensional problem can be parametrized by the Reynolds number, $\mathit{Re}=\left|{\mathit{v}}_{p}\right|D/\mathit{\nu}$, the
Prandtl number, $\mathit{Pr}=\mathit{\nu}/{\mathcal{D}}_{T}$, and the Schmidt number, $\mathit{Sc}=\mathit{\nu}/{\mathcal{D}}_{{n}_{\mathrm{v}}}$. While the Reynolds number is varied
up to a value of *Re*=600, the Prandtl and Schmidt numbers
are set, respectively, to values of *Pr*=0.72 and *Sc*=0.63 in this study, which correspond to representative values expected for humid air
in the temperature range of interest.

The numerical method employed is based on the method of Jenny and Dušek (2004), which has been used by, e.g., Kotouč et al. (2008) and Chouippe et al. (2019) to simulate heat and mass transfer around spherical particles. It relies on a spectral/spectral-element discretization of the Navier–Stokes equations (Eqs. 5, 6) coupled to the transport of heat and mass (Eqs. 2, 3). The spherical particle is placed at the origin of a cylindrical domain Ω. On the surface of the particle, impermeability and no-slip boundary conditions are imposed, while the scalar fields are subject to the constant Dirichlet boundary conditions $\stackrel{\mathrm{\u0303}}{T}={\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{v}}=\mathrm{1}$. At the upstream boundary of the simulation domain, a uniform velocity condition and $\stackrel{\mathrm{\u0303}}{T}={\stackrel{\mathrm{\u0303}}{n}}_{\mathrm{v}}=\mathrm{0}$ are enforced, while the downstream boundary is subject to zero-gradient boundary conditions for velocity as well as the scalars. The lateral boundaries are stress-free, with zero-gradient boundary conditions for the scalars, and zero pressure is imposed. The axial–radial plane is decomposed using the spectral-element method of Patera (1984), while the homogeneous azimuthal direction is treated using Fourier decomposition. The numerical mesh used in this work is shown in Fig. 1. For the temporal integration, a time splitting method is used, which consists of an explicit third-order Adams–Bashforth discretization for the advective terms and a first-order fully implicit discretization of the diffusion terms (Rønquist, 1988).

In the vertical direction the domain has a total length of 62 *D*, where the inflow length is 12 *D* and the outflow length 50 *D*
measured from the center of the sphere. The diameter of the domain is 16 *D*.
The computational domain has therefore been extended in the rear of the sphere compared to our former simulations
presented in Chouippe et al. (2019).
In total, 463 two-dimensional elements, each containing [6×6] collocation points, have been distributed over the domain.
The azimuthal Fourier series is truncated at the seventh mode. This resolution is comparable to our previous work and has been shown to give good results
for the momentum as well as for heat and mass transfer. All relevant scales of the flow and the scalar fields are resolved.

For more details on the numerical framework, the reader is referred to the previous work of Chouippe et al. (2019).

## 3.1 Hydrometeor wake regimes in clouds

The dynamics of the flow around spherical objects are, under the assumptions stated in the previous section, fully
parametrized by the Reynolds number, which depends on the diameter and, indirectly, through modulation of the terminal velocity,
on the density of the hydrometeor. Depending on the Reynolds number value, various flow states emerge in the
wake (Johnson and Patel, 1999; Jenny et al., 2004; Kotouč et al., 2009).
At low values the wake is axisymmetric and steady. When the first critical value of the Reynolds number, ${\mathit{Re}}_{\mathrm{c},\mathrm{1}}\approx \mathrm{212}$, is exceeded,
the wake becomes oblique with respect to the falling direction (*steady oblique regime)*, and only planar symmetry may be observed (Ghidersa and Dušek, 2000).
At ${\mathit{Re}}_{\mathrm{c},\mathrm{2}}\approx \mathrm{273}$ a second instability of Hopf type occurs (Ghidersa and Dušek, 2000), which leads to unsteady periodic vortex shedding,
while planar symmetry is still maintained (*oscillating oblique regime*). Eventually the vortex shedding becomes chaotic at
${\mathit{Re}}_{\mathrm{c},\mathrm{3}}\approx \mathrm{360}$ (Ormières and Provansal, 1999), and all instantaneous symmetries are lost (*chaotic regime*).
The listed regimes differ significantly in their ability to transfer heat and mass (Chouippe et al., 2019), and thus, different characteristics
in producing local supersaturation are to be expected.

In order to estimate the distribution of ice particle wake regimes in clouds, we adopt the size distribution of Marshall and Palmer (1948),

where ${\stackrel{\mathrm{\u203e}}{N}}_{\mathrm{met}}$ is the number concentration density of hydrometeors per unit volume of air, and ${\stackrel{\mathrm{\u203e}}{N}}_{\mathrm{met},\mathrm{0}},\mathit{\lambda}$ are model parameters.
Equation (10) was originally developed for raindrops; however, its validity for sufficiently large
ice particles has been
demonstrated by Passarelli (1978), Houze et al. (1979), Gordon and Marwitz (1984), Herzegh and Hobbs (1985) and Patade et al. (2015)
for various types of natural clouds.
The values of the model parameters ${\stackrel{\mathrm{\u203e}}{N}}_{\mathrm{met},\mathrm{0}}$, *λ* are usually obtained by airborne measurements and given as
functions of the cloud temperature.
One such parametrization is provided by Houze et al. (1979) for frontal clouds
within a temperature range from −42 to +6 ^{∘}C, which will
be used in the following for the size distribution of primary ice and is given by
the constitutive equations

where *T* denotes the cloud temperature in ^{∘}C. For this cloud type, ice enhancement has been
previously reported in the presence of riming ice particles (Hogan et al., 2002).
Even though it is known that the number density of ice particles of sub-millimeter size might deviate significantly
from Eq. (10) in a sub-exponential or super-exponential manner (Passarelli, 1978),
depending on cloud conditions, we assume that Eq. (10) holds nonetheless for all ice particle sizes
for the sake of simplicity.

For reasons explained in Sect. 3.3, our quantity of interest is not the number concentration
but rather the volume fraction *ϕ* of hydrometeors in clouds, which can be obtained
by integrating the corresponding moment of the Marshall–Palmer distribution; i.e.,

In order to determine the wake regimes from the size distribution, the terminal velocity needs to be approximated as a function
of the ice particle diameter, and hence, further assumptions on the density of the ice particles have to be made.
For graupel particles, the density may take values which range between
0.05 and 0.89 g cm^{−3}, depending on growth conditions and history
(Pruppacher and Klett, 2010). For our estimation we choose a value of 0.6 g cm^{−3}, as we are primarily concerned with
ice particles in the wet-growth mode, which typically constitute the upper end of the density range.
Using this value and the empirical drag law of Schiller and Naumann (1933), the critical diameters for regime transition are calculated
from the critical Reynolds numbers assuming a fluid density of 1 kg m^{−3} and a
kinematic viscosity of $\mathrm{1.68}\times {\mathrm{10}}^{-\mathrm{5}}$ m^{2} s^{−1}, which corresponds
approximately to the values expected at 750 kPa air pressure and
−10 ^{∘}C ambient temperature (see Table 1).
The total solid volume fraction can then be subdivided
into the contributions by ice particles with a certain wake regime; i.e.,

where the index *j* denotes one of the four wake regimes (steady axisymmetric, steady oblique, oscillating oblique and chaotic) and ${D}_{j,min}$ and
${D}_{j,max}$ the corresponding lower and upper limit for the hydrometeor diameter to be in this regime.

Figure 2 shows the estimated volume fraction of frozen hydrometeors in clouds as a function
of cloud temperature. The values are typically smaller than 10^{−5} and decrease exponentially with decreasing
temperature. When comparing the various regimes, it can be observed that the largest contributions generally stem
from hydrometeors in the axisymmetric or chaotic regime. This becomes especially clear when looking at the relative
contribution of *ϕ*_{j} to the total volume fraction *ϕ* as shown in Fig. 3. As can be
seen, more than 80 % of ice particles (by volume) exhibit one of these two wake regimes. Chaotic wakes
are dominant at temperatures close to the freezing point, while axisymmetric wakes dominate when the temperature is very
low. We will therefore focus solely on the axisymmetric and chaotic regimes in the following.

## 3.2 Supersaturation in the wake of hydrometeors

We now examine the saturation profiles in the wake of hydrometeors in the two regimes of interest.
Figure 4 shows isocontours of supersaturation with respect to ice, defined as ${s}_{i}={S}_{i}-\mathrm{1}$, for two different
values of the Reynolds number, namely *Re*=75, which lies within the axisymmetric regime, and
*Re*=600 in the chaotic wake regime. The ambient temperature was set to a value of ${T}_{\mathrm{\infty}}=-\mathrm{15}$ ^{∘}C;
i.e., the temperature gap between the ambient air and the meteor surface is 15 K, which corresponds to
rather extreme cloud conditions with high liquid water contents (Greenan and List, 1995).
Since the ambient fluid is already supersaturated
with respect to ice, the threshold of the isocontours is given as an excess to the value in the ambient.

At *Re*=75, the flow is steady, and thus, so is the supersaturation field. In the chaotic regime,
the flow is characterized by time-dependent vortex shedding from the ice particle's boundary layer, and therefore
excess supersaturation appears intermittently. In both regimes, significant excess supersaturation with respect to the
ambient can be observed far downstream, and the volume of air which is affected by the wake is by far larger
than the volume of the meteor. This can be seen more clearly when averaging the fields in the azimuthal direction,
which is statistically homogeneous for both cases, as well as in time for the unsteady flow. The averaged
excess supersaturation with respect to ice is shown in Fig. 5 for both cases.
The saturation profiles differ substantially from the ones obtained by Fukuta and Lee (1986) for similar
boundary conditions and Reynolds numbers. This is presumably due to their strong simplification of potential flow,
which is incapable of reproducing the boundary layer and the flow in the wake of the ice particle correctly.
This indeed leads to strong modifications in the distribution
of supersaturated regions. Regions of high vapor content and relatively low temperatures, i.e., supersaturated regions,
which are created in the mixing layer close to the riming particle are transported further downstream by the detaching
vortices in the chaotic regime. This phenomenon does not occur in the simulations of Fukuta and Lee (1986), who observed
a complete decay of
supersaturation to the inflow value within few meteor diameters downstream (approximately within two diameters
at −20 ^{∘}C),
while Fig. 5 clearly indicates that regions with ${s}_{i}(\mathit{x},t)>{s}_{i,\mathrm{\infty}}$ may
be observed at downstream distances of the order of 50 *D* at −15 ^{∘}C.

In order to compare the induced supersaturation for different values of *T*_{∞} and between different regimes, we compute
the time-averaged volume of air which is supersaturated above a given threshold,

where *V*_{s} is the supersaturated volume, *s*^{*} the corresponding threshold and *H* the Heaviside step function. The integral is
to be understood as a volume integral over the simulation domain Ω, and 〈.〉_{t} denotes the time-averaging
operator. The supersaturated volume can be computed with respect to ice or the liquid phase, which will be indicated in the following by
a phase subscript for the threshold variable.

Figure 6 shows the normalized volume of air supersaturated with respect to ice (a) and liquid (b)
as a function of the thresholds ${s}_{i}^{*}$ and ${s}_{w}^{*}$, respectively.
It can be seen that in the axisymmetric regime, the supersaturated volume is generally larger than in the
chaotic regime. This is caused by the enhanced mixing properties of the vortical wake structures in the latter case,
which lead to a faster decay of the scalar field.
For both regimes, the range of observed values of supersaturation is similar.
In fact, differences in range may only occur due to the distinct diffusivities of the temperature and water vapor fields, which
indeed lead to subsaturations in the chaotic regime in some parts of the flow.
While the range can be estimated reasonably well using simple mixing arguments (Chouippe et al., 2019), we observe that
the highest values of supersaturation only occur in a very confined portion of the wake.
However, if the temperature difference
between the ambient and the ice particle is sufficiently high, a strong excess in supersaturation compared to the ambient
value may occur for both phases in a volume comparable to the size of the meteor, if no mechanism for depletion is
considered. The volume of air which is affected by the wake-induced supersaturation is of the order of 10^{3} particle
volumes for all temperature gaps of interest.

## 3.3 Ice enhancement due to meteor wakes

In the previous section it was shown that considerable supersaturations may be observed in the wake of hydrometeors and that the volume of air which is affected by it is far larger than what has been previously estimated in the literature. In order to analyze our spatially resolved data in terms of ice nucleation, a constitutive relation is required which links supersaturation to the concentration of activated ice nuclei (IN). A possible parametrization obtained from continuous-flow diffusion chamber (CFDC) measurements for natural aerosols is given by Meyers et al. (1992) and reads

Here, *N*_{IN} denotes the number concentration of IN activated by the deposition and condensation–freezing
mechanisms. Equation (16) is reported to be strictly valid in the parameter range

which corresponds to the range of the underlying CFDC data; albeit it has also been applied
outside of this range (Meyers et al., 1992).
The temperatures of interest in the current work fall well into the range of validity.
As the distribution of *s*_{i} in the wake (Fig. 6a) shows,
*s*_{i} does only slightly exceed 25 % at the lowest ambient temperature considered.
However, only a small portion of the wake exhibits such high values, and hence, imprecisions are
likely to be insignificant for integral quantities.
The volumetric distribution of *s*_{w} (Fig. 6b) indicates that
supersaturation with respect to liquid exceeds the CFDC data range in significant portions of the
wake when *T*_{∞}*⪅*−10 ^{∘}C. Therefore, the constitutive relation
may underestimate the contribution of the condensation–freezing mode, as this mode
exhibits increased activity under these conditions, as has been demonstrated by Schaller and Fukuta (1979) for
various substances. Having this in mind, we are not aware of any parametrization of
condensation–freezing nucleation for natural aerosols which can be directly applied under these conditions.

One way to quantify the enhancement of nucleation activity is the ice enhancement factor, i.e., the ratio between the observed number concentration of ice nuclei and that expected for a reference state. This concept has been used before by Baker (1991) to study meteor-induced ice enhancement. For our problem, the local ice enhancement factor can be expressed as

where *N*_{IN}(** x**,

*t*) is calculated from the instantaneous local supersaturation with respect to ice using Eq. (16), and

*N*

_{IN,∞}is the IN concentration evaluated for ambient conditions. Thus, the ice enhancement factor gives a quantitative measure on how many additional IN are activated in the wake of the hydrometeor compared to the ambient.

Figure 7 visualizes the local ice enhancement factor in the wake
of a hydrometeor at ${T}_{\mathrm{\infty}}=-\mathrm{15}$ ^{∘}C. At this rather favorable parameter point, we observe
that ice nucleation is considerably increased (up to a factor of 6) in the near vicinity of the
hydrometeor, and significant values may still be observed in the far wake. It is therefore conceivable that
a considerable amount of new ice may be formed in the wake, which is consistent with recent experimental results of
Prabhakaran et al. (2020), who observed nucleation of water droplets and ice particles in the wake of a hot drop
under cold conditions.

Using Eq. (16), the supersaturated volume can be expressed as a function of the ice enhancement factor
in order to gain insight into the volumetric distribution of *f*_{i}.
Figure 8 shows the distributions for various
temperatures in both wake regimes. The ice enhancement is generally higher in the axisymmetric regime than in the chaotic regime.
Please note, that for a more precise analysis of this behavior, the actual flow-driven distribution of APs in the wake should be considered.
Depending on the flow dynamics, ice-nucleating particles may accumulate in the wake or be expelled from it (Homann and Bec, 2015),
which might have a non-negligible impact on ice enhancement. However, the consideration of aerosol dynamics is outside the scope of
the present work.

The analysis of the local ice enhancement factor leads to the conclusion that the wake of a hydrometeor
may act as a site of increased nucleation activity
and that the volume in which this increased activity occurs is much larger than the volume of the
hydrometeor itself. Having identified wake-induced supersaturation as a possible ice multiplication mechanism
in the near surroundings of a hydrometeor, its significance in the larger context of clouds is yet to be
investigated.
In order to assess this, we suggest computing a coarse-grained ice enhancement factor, which takes into
account the spatial distribution of *f*_{i} obtained by the numerical simulations, as well as the concentration
of primary ice particles. Therefore, we define a control volume 𝒱 within
a cloud which is sufficiently large such that the hydrometeor size distribution follows Eq. (10)
but sufficiently small such that the saturation field is uniform (with a value of *s*_{i,∞}) if no hydrometeor is present.
The volume-averaged ice enhancement factor within this control volume is given by

where the integral has been decomposed into an excess contribution, which evaluates to zero outside of the wake, and a base contribution, which is equal to unity. At sufficiently small volume fractions, the wakes of individual hydrometeors can be assumed to not interact, and hence, the integral in the enumerator of the first term in Eq. (19) can be expressed as a sum of contributions from independent hydrometeors; i.e.,

where Ω^{(k)}⊆𝒱 indicates that the field *N*_{IN}(** x**,

*t*) is to be taken from a simulation with an appropriate value of the Reynolds number for meteor (

*k*). Assuming a continuous particle size distribution according to Eq. (10), the volume-averaged ice enhancement within 𝒱 can be expressed as

where Ω(*D*) denotes the simulation domain for a meteor with diameter *D*. An evaluation of this expression
requires knowledge on how the supersaturation field evolves with the Reynolds number. However, our results for
the two regimes investigated indicate that the values of the non-dimensional integral

are of the same order of magnitude for a wide range of Reynolds numbers, as they only differ approximately by a factor of 3 for the two Reynolds numbers investigated. For reasons of simplicity, we therefore assume that ${\stackrel{\mathrm{\u0303}}{I}}_{\mathrm{\Omega}\left(D\right)}$ can be approximated for a given regime by the value determined for a single Reynolds number within that regime, ${\stackrel{\mathrm{\u0303}}{I}}_{{\mathrm{\Omega}}_{j}}$, which allows us to rewrite Eq. (21) in the simplified form

Since ${\stackrel{\mathrm{\u0303}}{I}}_{{\mathrm{\Omega}}_{j}}$ is now a constant for each regime, this can be reexpressed using Eq. (14) to yield

which demonstrates the importance of the ice volume fraction on global ice enhancement. As has already been discussed and shown in Fig. 3, the axisymmetric and chaotic regimes contribute the most to the volume fraction. We therefore disregard the contribution of the two other regimes and adjust the threshold of regime transition accordingly to a value of 1.26 mm, which allows us to evaluate Eq. (24) from our flow simulation data.

Both *ϕ*_{j} and ${\stackrel{\mathrm{\u0303}}{I}}_{{\mathrm{\Omega}}_{j}}$ are functions of ambient temperature, and thus the global ice enhancement factor originating from
meteor wakes can be expressed as a function of cloud temperature. While the former decreases exponentially with decreasing temperature,
the latter exhibits a strong increase, leading to a counteracting effect. Furthermore, a large number of ice particles pertain to the axisymmetric
regime for low temperatures, which is slightly more favorable in terms of ice enhancement.

Please be aware that, even though our aim is to quantify the wake-induced ice formation, the volume fraction of ice will be regarded as constant. The reason for this is that the approach used in this study does not allow us to derive the timescales of growth for newly created ice, and thus, the time-dependent coupling with the size distribution is inaccessible. Therefore, the following considerations merely apply to the initial state of a possible rapid glaciation process.

Figure 9 shows the global ice enhancement factor as a function of ambient temperature. It can be seen that ice nucleation in a given cloud volume is only marginally increased due to the presence of meteor wakes. At the highest temperature gap considered, the induced supersaturation activates roughly 0.08 % more ice than the ambient, and hence, we conclude that under typical cloud conditions hydrometeor wakes do not considerably affect ice nucleation on a cloud scale, at least when the volume fraction of ice in the cloud is as low as predicted by the parametrization of Houze et al. (1979).

Since the volume fraction is the limiting factor for ice enhancement, we revisited our assumption on the size distribution of ice particles in order to assess the sensitivity to the parametrization. Using the parametrization of Heymsfield et al. (2002) for deep subtropical and tropical clouds, the ice volume fraction is of the same order of magnitude as for the presented results, such that the conclusions which have been drawn for frontal clouds also persist for this type of cloud.

## 3.4 Timescales of nucleation and history effects

The ice enhancement analysis presented in the previous subsection did not explicitly address two factors which may be of importance when investigating IN activation, namely the time of exposure required to actually activate an AP and history effects, i.e., the circumstance that IN which have been activated in the wake may stay activated once they are not exposed to the wake anymore. The latter leads to an accumulative effect, which suggests that it is more reasonable to consider the volume swept by the meteors rather than look at the volume within a cloud which is instantaneously exposed to high supersaturations. Indeed, Prabhakaran et al. (2020) used this argument to show that essentially all interstitial APs within a cloud volume of interest are exposed to high supersaturations within a few minutes. Following their arguments, we estimated the time required to expose a significant portion of the cloud to the wakes using

where ${\dot{V}}_{\mathrm{sweep}}=\mathit{\u03f5}{v}_{p}{D}^{\mathrm{2}}\mathit{\pi}/\mathrm{4}$ is the volumetric flow rate
of air which is swept by a hydrometeor with diameter *D* and velocity *v*_{p}, and *ϵ* is
an unknown factor assumed to be of the order of unity. Using the size distribution of
Houze et al. (1979) and the terminal velocity for smooth spheres, we obtain
*τ*_{sweep}≈110 s at −15 ^{∘}C, which fits the estimation
of Prabhakaran et al. (2020) very well. The use of an empirical law for the terminal velocity
of frozen hydrometeors of natural shape leads to larger values of *τ*_{sweep}, which are, however, found to
be of a similar order of magnitude.

In the following we apply the swept-volume argument to ice multiplication by quantifying
the excess ice nucleation rate induced by the meteor wakes, henceforth denoted as *j*_{met}.
Under the assumption that nucleation occurs sufficiently fast to achieve the IN concentrations predicted by
Eq. (16), the nucleation rate can be estimated from the number of IN activated
in the wake and the time it takes to replenish the volume of fluid affected by high supersaturations.
The former is obtained directly from our simulation data by computing the volume integral
${\int}_{\mathrm{\Omega}\left(D\right)}\left(\langle {N}_{\mathrm{IN}}(\mathit{x},t){\rangle}_{t}-{N}_{\mathrm{IN},\mathrm{\infty}}\right)\mathrm{d}\mathit{x}$, while the latter is
difficult to define objectively. We propose to estimate the timescale of wake renewal
by

where ${V}_{\mathrm{aff}}=\mathit{\gamma}{D}^{\mathrm{3}}\mathit{\pi}/\mathrm{6}$ is the volume affected by the wake of a hydrometeor of diameter *D*,
which should be proportional to the volume of the hydrometeor. The prefactor *γ* is currently unknown
but might be related to the concept of supersaturated volume. The timescale
*τ*_{expo} can be interpreted as the characteristic time a fluid volume is exposed to high supersaturations
and hence the subscript. From the definitions of *V*_{aff} and ${\dot{V}}_{\mathrm{sweep}}$, it follows that
${\mathit{\tau}}_{\mathrm{expo}}\propto D/{v}_{p}$,
with the constant of proportionality being referred to as *C*_{expo} hereafter. This new constant
contains both unknown coefficients *ϵ* and *γ* and might be understood as the non-dimensional
streamwise length of the wake. Again, this length is difficult to define rigorously
due to the asymptotic decay of supersaturation, but judging from Fig. 7,
it is likely of the order of 10.
The swept-volume-limited nucleation rate for an ensemble of meteors is then given by

under the assumptions that IN activate sufficiently fast during the exposure, that a sufficient number of
activatable APs are present and that those are homogeneously distributed within the wake. The resulting nucleation
rate is displayed in Fig. 10 for *C*_{expo}=10 (solid blue line) and the range
$\mathrm{1}<{C}_{\mathrm{expo}}<\mathrm{100}$ in order to pay regard to the uncertainties associated with
this tunable parameter (shaded blue area). Under the assumptions just stated, meteor wakes indeed
appear to be capable of activating a large number of IN in a short amount of time, as has been
concluded by Prabhakaran et al. (2020). In fact, at −10 ^{∘}C
it would only take around 40 s for the number concentration of wake-activated IN to match the
concentration of primary meteors.

This result strongly contradicts our conclusion that wake-induced supersaturation is of little significance
in clouds based on the ice enhancement factor. As will be demonstrated in the following, the resolution
to this apparent contradiction requires an examination
of the timescales of nucleation and exposure. For *C*_{expo}=10, the exposure times are of the order of
5 ms for all diameters of interest, which is substantially shorter than the timescales usually
relevant for cloud modeling (of the order of minutes). As Eq. (16)
has been developed for cloud modeling, the validity of the assumption that IN activation can
be regarded as instantaneous at the timescales considered should be brought into question.
Indeed, classical nucleation theory (Fletcher, 1958) suggests that nucleation is a
time-dependent process until the activated fraction of APs approaches unity. From concentrations of
IN obtained from continuous-flow diffusion chamber (CFDC) experiments, the nucleation rate may be estimated
by taking into account the residence time in the apparatus, *τ*_{nucl} (Hoose and Möhler, 2012).
Since Eq. (16) is based on CFDC data, we make the conjecture that the local nucleation rate
may be approximated by the relationship

In Hoose and Möhler (2012) residence times ranging from 1.6 to 120 s
are reported for various CFDC experiments. The primary data used to obtain Eq. (16) also
suggest that the peak concentration *N*_{IN} is achieved with residence times of approximately
10 s (Al-Naimi and Saunders, 1985) and that shorter residence times lead
to lower concentrations, which is in accordance with the arguments stated above.
The crucial assumption of Eq. (27) that IN concentrations predicted by
Eq. (16) can be achieved within the exposure time is disproved by acknowledging that
*τ*_{nucl}≫*τ*_{expo}. The large discrepancy in timescales suggests that the nucleation rate is not
limited by the rate at which interstitial APs are entrained into the wake (the swept volume) but rather by the timescale of the nucleation process
itself and thus by the instantaneously exposed volume. A more appropriate estimate of the nucleation rate is then given by

under the assumption that APs are entrained sufficiently fast into the wake, which seems reasonable
given the arguments by Prabhakaran et al. (2020), and that they are distributed homogeneously
within the wake. The new estimation of the nucleation rate is shown in Fig. 10 for
*τ*_{nucl}=10 s (solid red line) and the range 1 s $<{\mathit{\tau}}_{\mathrm{nucl}}<\mathrm{100}$ s
(shaded red area). As can be seen, the rates are several orders of magnitude smaller than what has been
previously estimated, and it would now take around 82 h for the number concentration of
wake-activated IN to match the concentration of primary meteors, which is in better agreement with the
conclusion obtained with the help of the ice enhancement factor. In fact, it is straightforward to demonstrate that
these two quantities are directly linked by the relationship

Furthermore, it can be shown that as soon as a significant cloud volume is swept quickly at low volume fractions,
the exposure time will automatically be short as *τ*_{expo}∝*τ*_{sweep} for a given meteor
concentration, and thus, ice enhancement is again found to be limited by the low volumetric concentration of ice in clouds.

In this study we have performed numerical simulations of momentum, heat and
mass transfer around a warm hydrometeor in order to assess the distribution
of supersaturation in its wake and, moreover, the implications on ice nucleation enhancement.
Our simulation method is based on a body-conforming spectral/spectral-element discretization,
and all relevant scales of the flow problem have been resolved. The hydrometeor is assumed to
be of spherical shape and to possess a uniform surface temperature of 0 ^{∘}C, while
the ambient temperature has been varied in a range between −15 and
0 ^{∘}C. The vapor concentration is kept at saturation value with respect to ice
at the particle surface, while it is saturated with respect to water in the ambient
(reflecting the presence of supercooled droplets).
Two different values of the Reynolds number have been simulated in order to
capture the characteristics of the most relevant wake regimes, namely *Re*=75, where
the wake is steady and axisymmetric, and *Re*=600, where the wake is chaotic.

We found that significant values of supersaturation can be attained in the wake of warm hydrometeors, which persist long enough to be observed at several tens of particle diameters downstream of the meteor for sufficiently high differences in temperature. The supersaturated volume of air exceeds the estimations by Fukuta and Lee (1986) by far, which is attributed to the more accurate representation of the flow in the current study. This is an important observation since one of the key arguments for disregarding wake-induced ice nucleation is the proclaimed small zone of influence (Baker, 1991). It should be borne in mind that the ambient air was considered to be quiescent in this study. Under atmospheric conditions, turbulent fluctuations with scales comparable in size with the hydrometeors are expected to be present, which will presumably lead to a faster decay of the vapor concentration and temperature in the wake (Bagchi and Kottam, 2008). It could be worthwhile to investigate the influence of turbulence on the supersaturated volume in the future, as well as the importance of fluctuations in the ambient saturation ratios in comparison to wake-induced fluctuations.

Using a constitutive relation provided by Meyers et al. (1992) for the ice nuclei concentration, we estimated that heterogeneous nucleation may be locally enhanced by a factor of up to 6 if the temperature gap between the ambient and the meteor's surface is sufficiently large. This enhancement increases strongly with decreasing cloud temperature and so does the affected volume of air, which is typically of the order of several hundred particle volumes. It is therefore conceivable that ice multiplication can be triggered by the meteor's wakes, which is in agreement with the recent experimental observations by Prabhakaran et al. (2020). However, as arguments on upscaling have shown, this local effect alone presumably has little significance in clouds due to the low volumetric concentration of primary ice. This conclusion appears contradictory to the arguments given by Prabhakaran et al. (2020), who assessed that hydrometeors are capable of sampling a large cloud volume with their wakes within a short amount of time. While this is indeed the case, we have shown that the time during which ice-nucleating particles are exposed to the wake is generally too short to activate all of them, even at high values of supersaturation. The nucleation rate is therefore not limited by the rate at which interstitial APs are encountered, and therefore, the instantaneously exposed volume is argued to be more relevant, supporting the results obtained with the help of the ice enhancement factor.

In order for wake-induced supersaturation to be a relevant SIP, the nucleation rate needs to be considerably higher than what has been estimated in this work. While it might be possible that the condensation–freezing mode has been underestimated in this work, it seems unlikely that this underestimation is significant enough to substantially affect the conclusion. As the relevance of this SIP is mainly determined by the volumetric concentration of ice particles, high contents of ice are necessary for it to be active. Furthermore, in order to achieve an adequate temperature difference between the ambient and the ice, high liquid water contents are also required. Hence, if this SIP occurs in natural clouds, convective clouds are the most likely candidates.

Nonetheless, it is conceivable that the present mechanism in conjunction with one (or several) secondary ice processes is of greater relevance to the problem of ice formation. After a rapid glaciation process has been triggered, wake-induced nucleation might become significant, as the ice concentration then increases to considerable values. Furthermore, the question of whether or not wake-induced nucleation alone might trigger such a process has not been fully resolved yet, since the feedback on the size distribution of ice particles has been disregarded in this study. In order to assess this dynamical behavior, further time-resolved information on the activation of ice nuclei in the wake appears to be necessary. By tracking the trajectories of APs, the estimation of the nucleation rate may be further improved, and a resulting parametrization may be added to existing cloud models with explicit microphysics in order to assess the relative importance of wake-induced nucleation in comparison to other secondary ice processes.

The datasets are available upon request to the corresponding author.

TL and MU conceptualized the idea. JD provided the simulation code and mesh. MK conducted the simulations. MK, AC and MU post-processed the data, and TL assisted in the interpretation of the results. MK prepared and revised the paper, and AC contributed to it. MU, AC, JD and TL proofread the paper.

The authors declare that they have no conflict of interest.

The authors thank Alexei Kiselev for helpful discussions. The constructive feedback by Alexei Korolev and an anonymous reviewer is greatly acknowledged as their suggestions helped to improve and clarify this paper. The simulations were partially performed at the Steinbuch Centre for Computing in Karlsruhe, and the computer resources, technical expertise and assistance provided by this center are thankfully acknowledged. We acknowledge support by the KIT-Publication Fund of the Karlsruhe Institute of Technology.

The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.

This paper was edited by Martina Krämer and reviewed by Alexei Korolev and one anonymous referee.

Al-Naimi, R. and Saunders, C. P. R.: Measurements of Natural Deposition and Condensation-Freezing Ice Nuclei with a Continuous Flow Chamber, Atmos. Environ., 19, 1871–1882, https://doi.org/10.1016/0004-6981(85)90012-5, 1985. a

Auer, A. H., Veal, D. L., and Marwitz, J. D.: Observations of Ice Crystal and Ice Nuclei Concentrations in Stable Cap Clouds, J. Atmos. Sci., 26, 1342–1343, https://doi.org/10.1175/1520-0469(1969)026<1342:OOICAI>2.0.CO;2, 1969. a

Bacon, N. J., Swanson, B. D., Baker, M. B., and Davis, E. J.: Breakup of Levitated Frost Particles, J. Geophys. Res.-Atmos., 103, 13763–13775, https://doi.org/10.1029/98JD01162, 1998. a

Bagchi, P. and Kottam, K.: Effect of Freestream Isotropic Turbulence on Heat Transfer from a Sphere, Phys. Fluids, 20, 073305, https://doi.org/10.1063/1.2963138, 2008. a

Bagchi, P., Ha, M. Y., and Balachandar, S.: Direct Numerical Simulation of Flow and Heat Transfer From a Sphere in a Uniform Cross-Flow, J. Fluids Eng., 123, 347–358, https://doi.org/10.1115/1.1358844, 2001. a

Baker, B. A.: On the Nucleation of Ice in Highly Supersaturated Regions of Clouds, J. Atmos. Sci., 48, 1904–1907, https://doi.org/10.1175/1520-0469(1991)048<1905:OTNOII>2.0.CO;2, 1991. a, b, c, d, e

Bouchet, G., Mebarek, M., and Dušek, J.: Hydrodynamic Forces Acting on a Rigid Fixed Sphere in Early Transitional Regimes, Eur. J. Mech.-B Fluid., 25, 321–336, https://doi.org/10.1016/j.euromechflu.2005.10.001, 2006. a

Chouippe, A., Krayer, M., Uhlmann, M., Dušek, J., Kiselev, A., and Leisner, T.: Heat and Water Vapor Transfer in the Wake of a Falling Ice Sphere and Its Implication for Secondary Ice Formation in Clouds, New J. Phys., 21, 043043, https://doi.org/10.1088/1367-2630/ab0a94, 2019. a, b, c, d, e, f, g, h, i, j, k

de Stadler, M. B., Rapaka, N. R., and Sarkar, S.: Large Eddy Simulation of the near to Intermediate Wake of a Heated Sphere at Re = 10,000, Int. J. Heat Fluid Fl., 49, 2–10, https://doi.org/10.1016/j.ijheatfluidflow.2014.05.013, 2014. a

Dye, J. E. and Hobbs, P. V.: The Influence of Environmental Parameters on the Freezing and Fragmentation of Suspended Water Drops, J. Atmos. Sci., 25, 82–96, https://doi.org/10.1175/1520-0469(1968)025<0082:TIOEPO>2.0.CO;2, 1968. a

Field, P. R., Lawson, R. P., Brown, P. R. A., Lloyd, G., Westbrook, C., Moisseev, D., Miltenberger, A., Nenes, A., Blyth, A., Choularton, T., Connolly, P., Buehl, J., Crosier, J., Cui, Z., Dearden, C., DeMott, P., Flossmann, A., Heymsfield, A., Huang, Y., Kalesse, H., Kanji, Z. A., Korolev, A., Kirchgaessner, A., Lasher-Trapp, S., Leisner, T., McFarquhar, G., Phillips, V., Stith, J., and Sullivan, S.: Secondary Ice Production: Current State of the Science and Recommendations for the Future, Meteorol. Monogr., 58, 7.1–7.20, https://doi.org/10.1175/AMSMONOGRAPHS-D-16-0014.1, 2017. a, b

Fletcher, N. H.: Size Effect in Heterogeneous Nucleation, J. Chem. Phys., 29, 572–576, https://doi.org/10.1063/1.1744540, 1958. a

Fukuta, N. and Lee, H. J.: A Numerical Study of the Supersaturation Field around Growing Graupel, J. Atmos. Sci., 43, 1833–1843, https://doi.org/10.1175/1520-0469(1986)043<1833:ANSOTS>2.0.CO;2, 1986. a, b, c, d, e, f, g

Gagin, A.: The Effect of Supersaturation on the Ice Crystal Production by Natural Aerosols, Journal de Recherches Atmosphériques, 6, 175–185, 1972. a, b

Ghidersa, B. and Dušek, J.: Breaking of Axisymmetry and Onset of Unsteadiness in the Wake of a Sphere, J. Fluid Mech., 423, 33–69, https://doi.org/10.1017/S0022112000001701, 2000. a, b

Gordon, G. L. and Marwitz, J. D.: An Airborne Comparison of Three PMS Probes, J. Atmos. Ocean. Tech., 1, 22–27, https://doi.org/10.1175/1520-0426(1984)001<0022:AACOTP>2.0.CO;2, 1984. a

Greenan, B. J. W. and List, R.: Experimental Closure of the Heat and Mass Transfer Theory of Spheroidal Hailstones, J. Atmos. Sci., 52, 3797–3815, https://doi.org/10.1175/1520-0469(1995)052<3797:ECOTHA>2.0.CO;2, 1995. a

Hallett, J. and Mossop, S. C.: Production of Secondary Ice Particles during the Riming Process, Nature, 249, 26–28, https://doi.org/10.1038/249026a0, 1974. a

Herzegh, P. H. and Hobbs, P. V.: Size Spectra of Ice Particles in Frontal Clouds: Correlations between Spectrum Shape and Cloud Conditions, Q. J. Roy. Meteor. Soc., 111, 463–477, https://doi.org/10.1002/qj.49711146810, 1985. a

Heymsfield, A. J., Bansemer, A., Field, P. R., Durden, S. L., Stith, J. L., Dye, J. E., Hall, W., and Grainger, C. A.: 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, https://doi.org/10.1175/1520-0469(2002)059<3457:OAPOPS>2.0.CO;2, 2002. a

Hobbs, P. V.: Ice Multiplication in Clouds, J. Atmos. Sci., 26, 315–318, https://doi.org/10.1175/1520-0469(1969)026<0315:IMIC>2.0.CO;2, 1969. a

Hobbs, P. V. and Alkezweeny, A. J.: The Fragmentation of Freezing Water Droplets in Free Fall, J. Atmos. Sci., 25, 881–888, https://doi.org/10.1175/1520-0469(1968)025<0881:TFOFWD>2.0.CO;2, 1968. a

Hobbs, P. V. and Rangno, A. L.: Ice Particle Concentrations in Clouds, J. Atmos. Sci., 42, 2523–2549, https://doi.org/10.1175/1520-0469(1985)042<2523:IPCIC>2.0.CO;2, 1985. a

Hogan, R. J., Field, P. R., Illingworth, A. J., Cotton, R. J., and Choularton, T. W.: Properties of Embedded Convection in Warm-Frontal Mixed-Phase Cloud from Aircraft and Polarimetric Radar, Q. J. Roy. Meteor. Soc., 128, 451–476, https://doi.org/10.1256/003590002321042054, 2002. a, b

Homann, H. and Bec, J.: Concentrations of Inertial Particles in the Turbulent Wake of an Immobile Sphere, Phys. Fluids, 27, 053301, https://doi.org/10.1063/1.4919723, 2015. a

Hoose, C. and Möhler, O.: Heterogeneous ice nucleation on atmospheric aerosols: a review of results from laboratory experiments, Atmos. Chem. Phys., 12, 9817–9854, https://doi.org/10.5194/acp-12-9817-2012, 2012. a, b

Houze, R. A., Hobbs, P. V., Herzegh, P. H., and Parsons, D. B.: Size Distributions of Precipitation Particles in Frontal Clouds, J. Atmos. Sci., 36, 156–162, https://doi.org/10.1175/1520-0469(1979)036<0156:SDOPPI>2.0.CO;2, 1979. a, b, c, d

Jenny, M. and Dušek, J.: Efficient Numerical Method for the Direct Numerical Simulation of the Flow Past a Single Light Moving Spherical Body in Transitional Regimes, J. Comput. Phys., 194, 215–232, https://doi.org/10.1016/j.jcp.2003.09.004, 2004. a

Jenny, M., Dušek, J., and Bouchet, G.: Instabilities and Transition of a Sphere Falling or Ascending Freely in a Newtonian Fluid, J. Fluid Mech., 508, 201–239, https://doi.org/10.1017/S0022112004009164, 2004. a

Johnson, T. A. and Patel, V. C.: Flow Past a Sphere up to a Reynolds Number of 300, J. Fluid Mech., 378, 19–70, https://doi.org/10.1017/S0022112098003206, 1999. a, b

Koenig, L. R.: The Glaciating Behavior of Small Cumulonimbus Clouds, J. Atmos. Sci., 20, 29–47, https://doi.org/10.1175/1520-0469(1963)020<0029:TGBOSC>2.0.CO;2, 1963. a

Korolev, A., Heckman, I., Wolde, M., Ackerman, A. S., Fridlind, A. M., Ladino, L. A., Lawson, R. P., Milbrandt, J., and Williams, E.: A new look at the environmental conditions favorable to secondary ice production, Atmos. Chem. Phys., 20, 1391–1429, https://doi.org/10.5194/acp-20-1391-2020, 2020. a

Kotouč, M., Bouchet, G., and Dušek, J.: Loss of Axisymmetry in the Mixed Convection, Assisting Flow Past a Heated Sphere, Int. J. Heat Mass Tran., 51, 2686–2700, https://doi.org/10.1016/j.ijheatmasstransfer.2007.10.005, 2008. a

Kotouč, M., Bouchet, G., and Dušek, J.: Transition to Turbulence in the Wake of a Fixed Sphere in Mixed Convection, J. Fluid Mech., 625, 205, https://doi.org/10.1017/S0022112008005557, 2009. a

Marshall, J. S. and Palmer, W. M. K.: The Distribution of Raindrops with Size, J. Meteorol., 5, 165–166, https://doi.org/10.1175/1520-0469(1948)005<0165:TDORWS>2.0.CO;2, 1948. a

Meyers, M. P., DeMott, P. J., and Cotton, W. R.: New Primary Ice-Nucleation Parameterizations in an Explicit Cloud Model, J. Appl. Meteorol., 31, 708–721, https://doi.org/10.1175/1520-0450(1992)031<0708:NPINPI>2.0.CO;2, 1992. a, b, c

Mossop, S. C.: The Origin and Concentration of Ice Crystals in Clouds, B. Am. Meteorol. Soc., 66, 264–273, https://doi.org/10.1175/1520-0477(1985)066<0264:TOACOI>2.0.CO;2, 1985. a

Murphy, D. M. and Koop, T.: Review of the Vapour Pressures of Ice and Supercooled Water for Atmospheric Applications, Q. J. Roy. Meteor. Soc., 131, 1539–1565, https://doi.org/10.1256/qj.04.94, 2005. a

Nix, N. and Fukuta, N.: Nonsteady-State Kinetics of Droplet Growth in Cloud Physics, J. Atmos. Sci., 31, 1334–1343, https://doi.org/10.1175/1520-0469(1974)031<1334:NSKODG>2.0.CO;2, 1974. a

Ormières, D. and Provansal, M.: Transition to Turbulence in the Wake of a Sphere, Phys. Rev. Lett., 83, 6–9, https://doi.org/10.1103/PhysRevLett.83.80, 1999. a

Passarelli, R. E.: Theoretical and Observational Study of Snow-Size Spectra and Snowflake Aggregation Efficiencies, J. Atmos. Sci., 35, 882–889, https://doi.org/10.1175/1520-0469(1978)035<0882:TAOSOS>2.0.CO;2, 1978. a, b

Patade, S., Prabha, T. V., Axisa, D., Gayatri, K., and Heymsfield, A.: Particle Size Distribution Properties in Mixed-Phase Monsoon Clouds from in Situ Measurements during CAIPEEX, J. Geophys. Res.-Atmos., 120, 10418–10440, https://doi.org/10.1002/2015JD023375, 2015. a

Patera, A. T.: A Spectral Element Method for Fluid Dynamics: Laminar Flow in a Channel Expansion, J. Comput. Phys., 54, 468–488, https://doi.org/10.1016/0021-9991(84)90128-1, 1984. a

Prabhakaran, P., Weiss, S., Krekhov, A., Pumir, A., and Bodenschatz, E.: Can Hail and Rain Nucleate Cloud Droplets?, Phys. Rev. Lett., 119, 128701, https://doi.org/10.1103/PhysRevLett.119.128701, 2017. a

Prabhakaran, P., Kinney, G., Cantrell, W., Shaw, R. A., and Bodenschatz, E.: High Supersaturation in the Wake of Falling Hydrometeors: Implications for Cloud Invigoration and Ice Nucleation, Geophys. Res. Lett., 47, e2020GL088055, https://doi.org/10.1029/2020GL088055, 2020. a, b, c, d, e, f, g, h, i, j

Pruppacher, H. R. and Klett, J. D.: Microphysics of Clouds and Precipitation, Springer, the Netherlands, https://doi.org/10.1007/978-0-306-48100-0, 2010. a, b

Rangno, A. L. and Hobbs, P. V.: Ice Particle Concentrations and Precipitation Development in Small Polar Maritime Cumuliform Clouds, Q. J. Roy. Meteor. Soc., 117, 207–241, https://doi.org/10.1002/qj.49711749710, 1991. a

Rønquist, E. M.: Optimal Spectral Element Methods for the Unsteady Three-Dimensional Incompressible Navier-Stokes Equations, PhD thesis, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA, 1988. a

Schaller, R. C. and Fukuta, N.: Ice Nucleation by Aerosol Particles: Experimental Studies Using a Wedge-Shaped Ice Thermal Diffusion Chamber, J. Atmos. Sci., 36, 1788–1802, https://doi.org/10.1175/1520-0469(1979)036<1788:INBAPE>2.0.CO;2, 1979. a

Schiller, L. and Naumann, A.: Über Die Grundlegenden Berechnungen Bei Der Schwerkraftaufbereitung, Z. Ver. Dtsch. Ing, 77, 318–320, 1933. a

Takahashi, T., Nagao, Y., and Kushiyama, Y.: Possible High Ice Particle Production during Graupel–Graupel Collisions, J. Atmos. Sci., 52, 4523–4527, https://doi.org/10.1175/1520-0469(1995)052<4523:PHIPPD>2.0.CO;2, 1995. a

Vardiman, L.: The Generation of Secondary Ice Particles in Clouds by Crystal–Crystal Collision, J. Atmos. Sci., 35, 2168–2180, https://doi.org/10.1175/1520-0469(1978)035<2168:TGOSIP>2.0.CO;2, 1978. a

Zhou, W. and Dušek, J.: Chaotic States and Order in the Chaos of the Paths of Freely Falling and Ascending Spheres, Int. J. Multiphas. Flow, 75, 205–223, https://doi.org/10.1016/j.ijmultiphaseflow.2015.05.010, 2015. a