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 and
indicated that it might be irrelevant, the more recent
study of 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.
Introduction
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 , which has been observed
for various cloud types
.
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 .
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 , fragmentation of ice
and freezing drops .
Most of these mechanisms have been implemented into cloud models with explicit ice microphysics
(see , 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, 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
. This hypothesis was corroborated by , 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 .
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 , 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 that the overall ice enhancement is likely to be negligible.
While the previously mentioned studies focused on the quantification of the supersaturation field,
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
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, 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 and
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 .
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
and thus strongly affect the distribution of supersaturation around hydrometeors .
In , 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 .
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
, and therefore, we attempt to link our results to heterogeneous nucleation of aerosol particles.
Methodology
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
and model the hydrometeor as a sphere with constant diameter D, constant surface temperature Tp and
a constant vapor pressure ev,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 esat,j at Tp.
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 Tp=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
Sj(e,T)=e/esat,j(T),
where e denotes the local vapor pressure. A parametrization of esat,j with respect to
temperature for both phase equilibria in the temperature range of interest is adopted from .
The spatial distribution of Sj can then be reconstructed from the temperature field T
and the vapor concentration nv, which are governed by the transport equations
2∂T∂t+u⋅∇T=DT∇2T,3∂nv∂t+u⋅∇nv=Dnv∇2nv,
where DT and Dnv denote the heat and vapor diffusivities, and 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
e=nvkbT,
where kb is Boltzmann's constant. Even though the diffusivities of both vapor and heat are similar in magnitude,
it was shown by 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. () leads to an outward heat
flux because the surface temperature of the hydrometeor is higher than that of its surroundings. Furthermore,
at Tp=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∞<Tp due to the presence of
droplets in mixed-phase clouds; i.e., e∞=esat,w(T∞). 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 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,
5∇⋅u=0,6∂u∂t+(u⋅∇)u=-1ρf∇p+ν∇2u,
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 vp equal to its terminal velocity vT; i.e.,
vp=-vTez,
where ez 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 vp are not permitted. The necessity to account for these variations was discussed in ,
and it was found that it has little influence in the present context due to the expected high value of the density ratio.
Spectral-element mesh in the axial–radial plane, where the lines depict the boundaries of the 463
elements. Each element contains [6×6] collocation points. The three-dimensionality is introduced
by Fourier expansion in the azimuthal direction, which is truncated at the seventh Fourier mode,
resulting in a cylindrical domain. A uniform velocity profile with constant temperature and vapor content is
imposed at the upstream boundary, while the downstream boundary is subject to zero-gradient boundary conditions (BCs) for both velocity and scalars.
The lateral boundaries are stress-free at zero pressure, with zero-gradient BCs for the scalars.
The surface of the spherical particle, whose center is located at the origin of the coordinate system, is impermeable, and a no-slip BC
for velocity is imposed, as well as constant Dirichlet BCs for the scalar fields.
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 vT, and the scalar
transport equations are formulated in terms of
8T̃=(T-T∞)/(Tp-T∞),9ñv=(nv-nv,∞)/(nv,p-nv,∞),
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, Re=vpD/ν, the
Prandtl number, Pr=ν/DT, and the Schmidt number, Sc=ν/Dnv. 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 , which has been used by, e.g.,
and to simulate heat and mass transfer around spherical particles.
It relies on a spectral/spectral-element discretization of the Navier–Stokes equations (Eqs. , )
coupled to the transport of heat and mass (Eqs. , ).
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 T̃=ñv=1. At the upstream boundary of the simulation domain,
a uniform velocity condition and T̃=ñv=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 , while the homogeneous
azimuthal direction is treated using Fourier decomposition. The numerical mesh used in this work is shown in Fig. .
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 .
In the vertical direction the domain has a total length of 62D, where the inflow length is 12D and the outflow length 50D
measured from the center of the sphere. The diameter of the domain is 16D.
The computational domain has therefore been extended in the rear of the sphere compared to our former simulations
presented in .
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 .
ResultsHydrometeor 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 .
At low values the wake is axisymmetric and steady. When the first critical value of the Reynolds number, Rec,1≈212, is exceeded,
the wake becomes oblique with respect to the falling direction (steady oblique regime), and only planar symmetry may be observed .
At Rec,2≈273 a second instability of Hopf type occurs , which leads to unsteady periodic vortex shedding,
while planar symmetry is still maintained (oscillating oblique regime). Eventually the vortex shedding becomes chaotic at
Rec,3≈360, and all instantaneous symmetries are lost (chaotic regime).
The listed regimes differ significantly in their ability to transfer heat and mass , 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 ,
N‾met=N‾met,0exp(-λD),
where N‾met is the number concentration density of hydrometeors per unit volume of air, and N‾met,0,λ are model parameters.
Equation () was originally developed for raindrops; however, its validity for sufficiently large
ice particles has been
demonstrated by , , , and
for various types of natural clouds.
The values of the model parameters N‾met,0, λ are usually obtained by airborne measurements and given as
functions of the cloud temperature.
One such parametrization is provided by 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
11N‾met,0(T)=5.5×106exp-0.09Tm-4,12λ(T)=9.6×10-1exp-0.056Tmm-1,
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 .
Even though it is known that the number density of ice particles of sub-millimeter size might deviate significantly
from Eq. () in a sub-exponential or super-exponential manner ,
depending on cloud conditions, we assume that Eq. () holds nonetheless for all ice particle sizes
for the sake of simplicity.
For reasons explained in Sect. , 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.,
ϕ=∫0∞πD36N‾met(D)dD.
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 gcm-3, depending on growth conditions and history
. For our estimation we choose a value of 0.6 gcm-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 , the critical diameters for regime transition are calculated
from the critical Reynolds numbers assuming a fluid density of 1 kgm-3 and a
kinematic viscosity of 1.68×10-5m2s-1, which corresponds
approximately to the values expected at 750 kPa air pressure and
-10∘C ambient temperature (see Table ).
The total solid volume fraction can then be subdivided
into the contributions by ice particles with a certain wake regime; i.e.,
ϕ≡∑jϕj=∑j∫Dj,minDj,maxπD36N‾met(D)dD,
where the index j denotes one of the four wake regimes (steady axisymmetric, steady oblique, oscillating oblique and chaotic) and Dj,min and
Dj,max the corresponding lower and upper limit for the hydrometeor diameter to be in this regime.
Summary of the critical Reynolds numbers for regime transition and the threshold diameters used to determine
the partial volume fractions.
Volume fraction of ice particles in clouds as a function of ambient temperature. The total volume fraction
ϕ is displayed by a solid black line, while the contribution by the regimes is
given by the colored lines.
Linestyles: axisymmetric regime (blue), steady oblique regime (red),
oscillating oblique regime (yellow) and chaotic regime (purple).
Figure 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. . 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.
Relative contribution to the total volume fraction by hydrometeors in a certain regime.
Linestyles: axisymmetric regime (blue), steady oblique regime (red),
oscillating oblique regime (yellow) and chaotic regime (purple).
Supersaturation in the wake of hydrometeors
We now examine the saturation profiles in the wake of hydrometeors in the two regimes of interest.
Figure shows isocontours of supersaturation with respect to ice, defined as si=Si-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∞=-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 .
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.
Isosurfaces of supersaturation in the wake at T∞=-15∘C. The
value of the isocontour is si,∞+0.02, i.e., two percentage points higher than the
ambient supersaturation. Two different wake regimes are depicted, which correspond to two different
hydrometeor sizes in our framework. (a) Axisymmetric regime at Re=75. (b) Chaotic regime at Re=600.
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. for both cases.
The saturation profiles differ substantially from the ones obtained by 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 , 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. clearly indicates that regions with si(x,t)>si,∞ may
be observed at downstream distances of the order of 50D at -15∘C.
Contours of excess supersaturation in the wake, averaged over time and azimuthal direction at T∞=-15∘C.
(a) Axisymmetric regime at Re=75. (b) Chaotic regime at Re=600.
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,
Vs(s*)≡∫Ω〈H(s(x,t)-s*)〉tdx,
where Vs 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 shows the normalized volume of air supersaturated with respect to ice (a) and liquid (b)
as a function of the thresholds si* and sw*, 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 , 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 103 particle
volumes for all temperature gaps of interest.
Volume of air where supersaturation exceeds a given threshold as a function of the
threshold. The volume supersaturated with respect to the solid phase is shown in (a), whereas (b)
shows the volume supersaturated with respect to the liquid phase.
The volume of the ice particle is used for normalization and four different
ambient temperatures are shown:
T∞=-6∘C (purple),
T∞=-9∘C (yellow),
T∞=-12∘C (red) and
T∞=-15∘C (blue).
Solid lines correspond to Re=600 (chaotic regime, time average), while dashed lines show the data obtained for
Re=75 (axisymmetric regime).
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 and reads
NIN=0.528exp12.96sim-3.
Here, NIN denotes the number concentration of IN activated by the deposition and condensation–freezing
mechanisms. Equation () is reported to be strictly valid in the parameter range
-20∘C<T<-7∘C,0.02<si<0.25,-0.05<sw<0.045,
which corresponds to the range of the underlying CFDC data; albeit it has also been applied
outside of this range .
The temperatures of interest in the current work fall well into the range of validity.
As the distribution of si in the wake (Fig. a) shows,
si 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 sw (Fig. b) 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 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 to study meteor-induced ice enhancement.
For our problem, the local ice enhancement factor can be expressed as
fi(x)=〈NIN(x,t)〉tNIN,∞,
where NIN(x,t) is calculated from the instantaneous local supersaturation with respect to ice using
Eq. (), and NIN,∞ 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 visualizes the local ice enhancement factor in the wake
of a hydrometeor at T∞=-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
, who observed nucleation of water droplets and ice particles in the wake of a hot drop
under cold conditions.
Contours of local ice enhancement factor in the wake, averaged over time and azimuthal direction at T∞=-15∘C.
(a) Axisymmetric regime at Re=75. (b) Chaotic regime at Re=600.
Using Eq. (), the supersaturated volume can be expressed as a function of the ice enhancement factor
in order to gain insight into the volumetric distribution of fi.
Figure 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 ,
which might have a non-negligible impact on ice enhancement. However, the consideration of aerosol dynamics is outside the scope of
the present work.
Volume of air with supersaturation above a given threshold as a function of the
ice enhancement factor. The volume is normalized by the volume of the ice particle, and four different temperatures
are shown:
T∞=-6∘C (purple),
T∞=-9∘C (yellow),
T∞=-12∘C (red) and
T∞=-15∘C (blue).
Solid lines correspond to Re=600 (chaotic regime, time average), while dashed lines show the data obtained for
Re=75 (axisymmetric regime).
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 fi obtained by the numerical simulations, as well as the concentration
of primary ice particles. Therefore, we define a control volume V within
a cloud which is sufficiently large such that the hydrometeor size distribution follows Eq. ()
but sufficiently small such that the saturation field is uniform (with a value of si,∞) if no hydrometeor is present.
The volume-averaged ice enhancement factor within this control volume is given by
〈fi〉V=∫V〈NIN(x,t)〉t-NIN,∞dx∫VNIN,∞dx+1,
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. () can be expressed
as a sum of contributions from independent hydrometeors; i.e.,
∫V〈NIN(x,t)〉t-NIN,∞dx=∑(k)∫Ω(k)〈NIN(x,t)〉t-NIN,∞dx,
where Ω(k)⊆V indicates that the field NIN(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. (), the volume-averaged
ice enhancement within V can be expressed as
〈fi〉V=1+1NIN,∞∫0∞N‾met(D)×∫Ω(D)〈NIN(x,t)〉t-NIN,∞dxdD,
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
ĨΩ(D)≡6π∫Ω(D)〈NIN(x̃,t)〉t-NIN,∞dx̃
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 ĨΩ(D)
can be approximated for a given regime by the value determined for a single Reynolds number
within that regime, ĨΩj, which allows us to rewrite Eq. () in the simplified form
〈fi〉V=1+1NIN,∞∑j∫Dj,minDj,maxπD36N‾met(D)ĨΩjdD.
Since ĨΩj is now a constant for each regime, this can be reexpressed using Eq. () to yield
〈fi〉V=1+∑jϕjĨΩj/NIN,∞,
which demonstrates the importance of the ice volume fraction on global ice enhancement.
As has already been discussed and shown in Fig. ,
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. ()
from our flow simulation data.
Both ϕj and ĨΩ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.
Global ice enhancement factor as a function of cloud temperature.
The inset shows the same data but on a semi-logarithmic scale.
Figure 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 .
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 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.
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, 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
τsweep=∫0∞N‾met(D)V˙sweep(D)dD-1,
where V˙sweep=ϵvpD2π/4 is the volumetric flow rate
of air which is swept by a hydrometeor with diameter D and velocity vp, and ϵ is
an unknown factor assumed to be of the order of unity. Using the size distribution of
and the terminal velocity for smooth spheres, we obtain
τsweep≈110s at -15∘C, which fits the estimation
of 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 jmet.
Under the assumption that nucleation occurs sufficiently fast to achieve the IN concentrations predicted by
Eq. (), 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
∫Ω(D)〈NIN(x,t)〉t-NIN,∞dx, while the latter is
difficult to define objectively. We propose to estimate the timescale of wake renewal
by
τexpo≡Vaff/V˙sweep,
where Vaff=γD3π/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 Vaff and V˙sweep, it follows that
τexpo∝D/vp,
with the constant of proportionality being referred to as Cexpo 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. ,
it is likely of the order of 10.
The swept-volume-limited nucleation rate for an ensemble of meteors is then given by
jmetexpo=∫0∞N‾met(D)1τexpo(D)×∫Ω(D)〈NIN(x,t)〉t-NIN,∞dxdD
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. for Cexpo=10 (solid blue line) and the range
1<Cexpo<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 . 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.
Limiting cases for the nucleation rate jmet. The swept-volume-limited estimate based on the arguments
of and evaluated from the present data according to Eq. () is shown
for Cexpo=10 (solid blue line), with the shaded area depicting the values obtained for 1<Cexpo<100.
The exposure-time-limited estimate defined in Eq. (), which is directly linked to the ice enhancement factor,
is shown for τnucl=10s (solid red line) and the
range 1s<τnucl<100s (shaded area).
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 Cexpo=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. ()
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 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.
Since Eq. () is based on CFDC data, we make the conjecture that the local nucleation rate
may be approximated by the relationship
jIN(x,t)≈NIN(x,t)/τnucl.
In residence times ranging from 1.6 to 120 s
are reported for various CFDC experiments. The primary data used to obtain Eq. () also
suggest that the peak concentration NIN is achieved with residence times of approximately
10 s and that shorter residence times lead
to lower concentrations, which is in accordance with the arguments stated above.
The crucial assumption of Eq. () that IN concentrations predicted by
Eq. () 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
jmetnucl=∫0∞N‾met(D)1τnucl×∫Ω(D)〈NIN(x,t)〉t-NIN,∞dxdD
under the assumption that APs are entrained sufficiently fast into the wake, which seems reasonable
given the arguments by , and that they are distributed homogeneously
within the wake. The new estimation of the nucleation rate is shown in Fig. for
τnucl=10s (solid red line) and the range 1 s<τnucl<100s
(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
〈fi〉V=jmetnuclNIN,∞/τnucl+1.
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.
Conclusion
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 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 . 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 .
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 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
.
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 , 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.
Data availability
The datasets are available upon request to the corresponding author.
Author contributions
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.
Competing interests
The authors declare that they have no conflict of interest.
Acknowledgements
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.
Financial support
The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
Review statement
This paper was edited by Martina Krämer and reviewed by Alexei Korolev and one anonymous referee.
ReferencesAl-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, 10.1016/0004-6981(85)90012-5, 1985.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,
10.1175/1520-0469(1969)026<1342:OOICAI>2.0.CO;2, 1969.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, 10.1029/98JD01162, 1998.Bagchi, P. and Kottam, K.: Effect of Freestream Isotropic Turbulence on Heat
Transfer from a Sphere, Phys. Fluids, 20, 073305,
10.1063/1.2963138, 2008.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,
10.1115/1.1358844, 2001.Baker, B. A.: On the Nucleation of Ice in Highly Supersaturated
Regions of Clouds, J. Atmos. Sci., 48, 1904–1907,
10.1175/1520-0469(1991)048<1905:OTNOII>2.0.CO;2, 1991.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, 10.1016/j.euromechflu.2005.10.001,
2006.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, 10.1088/1367-2630/ab0a94, 2019.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,
10.1016/j.ijheatfluidflow.2014.05.013, 2014.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,
10.1175/1520-0469(1968)025<0082:TIOEPO>2.0.CO;2, 1968.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,
10.1175/AMSMONOGRAPHS-D-16-0014.1, 2017.Fletcher, N. H.: Size Effect in Heterogeneous Nucleation, J. Chem. Phys., 29, 572–576, 10.1063/1.1744540, 1958.Fukuta, N. and Lee, H. J.: A Numerical Study of the Supersaturation
Field around Growing Graupel, J. Atmos. Sci., 43,
1833–1843, 10.1175/1520-0469(1986)043<1833:ANSOTS>2.0.CO;2, 1986.
Gagin, A.: The Effect of Supersaturation on the Ice Crystal Production by
Natural Aerosols, Journal de Recherches Atmosphériques, 6, 175–185, 1972.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, 10.1017/S0022112000001701, 2000.Gordon, G. L. and Marwitz, J. D.: An Airborne Comparison of Three PMS
Probes, J. Atmos. Ocean. Tech., 1, 22–27,
10.1175/1520-0426(1984)001<0022:AACOTP>2.0.CO;2, 1984.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,
10.1175/1520-0469(1995)052<3797:ECOTHA>2.0.CO;2, 1995.Hallett, J. and Mossop, S. C.: Production of Secondary Ice Particles during the
Riming Process, Nature, 249, 26–28, 10.1038/249026a0, 1974.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,
10.1002/qj.49711146810, 1985.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,
10.1175/1520-0469(2002)059<3457:OAPOPS>2.0.CO;2, 2002.Hobbs, P. V.: Ice Multiplication in Clouds, J. Atmos.
Sci., 26, 315–318, 10.1175/1520-0469(1969)026<0315:IMIC>2.0.CO;2,
1969.Hobbs, P. V. and Alkezweeny, A. J.: The Fragmentation of Freezing Water
Droplets in Free Fall, J. Atmos. Sci., 25,
881–888, 10.1175/1520-0469(1968)025<0881:TFOFWD>2.0.CO;2, 1968.Hobbs, P. V. and Rangno, A. L.: Ice Particle Concentrations in Clouds,
J. Atmos. Sci., 42, 2523–2549,
10.1175/1520-0469(1985)042<2523:IPCIC>2.0.CO;2, 1985.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, 10.1256/003590002321042054,
2002.Homann, H. and Bec, J.: Concentrations of Inertial Particles in the Turbulent
Wake of an Immobile Sphere, Phys. Fluids, 27, 053301,
10.1063/1.4919723, 2015.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, 10.5194/acp-12-9817-2012, 2012.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,
10.1175/1520-0469(1979)036<0156:SDOPPI>2.0.CO;2, 1979.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,
10.1016/j.jcp.2003.09.004, 2004.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, 10.1017/S0022112004009164, 2004.Johnson, T. A. and Patel, V. C.: Flow Past a Sphere up to a Reynolds Number
of 300, J. Fluid Mech., 378, 19–70,
10.1017/S0022112098003206, 1999.Koenig, L. R.: The Glaciating Behavior of Small Cumulonimbus Clouds,
J. Atmos. Sci., 20, 29–47,
10.1175/1520-0469(1963)020<0029:TGBOSC>2.0.CO;2, 1963.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, 10.5194/acp-20-1391-2020, 2020.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,
10.1016/j.ijheatmasstransfer.2007.10.005, 2008.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, 10.1017/S0022112008005557, 2009.Marshall, J. S. and Palmer, W. M. K.: The Distribution of Raindrops with Size,
J. Meteorol., 5, 165–166,
10.1175/1520-0469(1948)005<0165:TDORWS>2.0.CO;2, 1948.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,
10.1175/1520-0450(1992)031<0708:NPINPI>2.0.CO;2, 1992.Mossop, S. C.: The Origin and Concentration of Ice Crystals in
Clouds, B. Am. Meteorol. Soc., 66, 264–273,
10.1175/1520-0477(1985)066<0264:TOACOI>2.0.CO;2, 1985.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, 10.1256/qj.04.94, 2005.Nix, N. and Fukuta, N.: Nonsteady-State Kinetics of Droplet Growth in
Cloud Physics, J. Atmos. Sci., 31, 1334–1343,
10.1175/1520-0469(1974)031<1334:NSKODG>2.0.CO;2, 1974.Ormières, D. and Provansal, M.: Transition to Turbulence in the
Wake of a Sphere, Phys. Rev. Lett., 83, 6–9,
10.1103/PhysRevLett.83.80, 1999.Passarelli, R. E.: Theoretical and Observational Study of Snow-Size
Spectra and Snowflake Aggregation Efficiencies, J.
Atmos. Sci., 35, 882–889,
10.1175/1520-0469(1978)035<0882:TAOSOS>2.0.CO;2, 1978.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, 10.1002/2015JD023375, 2015.Patera, A. T.: A Spectral Element Method for Fluid Dynamics: Laminar Flow
in a Channel Expansion, J. Comput. Phys., 54, 468–488,
10.1016/0021-9991(84)90128-1, 1984.Prabhakaran, P., Weiss, S., Krekhov, A., Pumir, A., and Bodenschatz, E.: Can
Hail and Rain Nucleate Cloud Droplets?, Phys. Rev. Lett., 119,
128701, 10.1103/PhysRevLett.119.128701, 2017.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,
10.1029/2020GL088055, 2020.Pruppacher, H. R. and Klett, J. D.: Microphysics of Clouds and
Precipitation, Springer, the Netherlands, 10.1007/978-0-306-48100-0,
2010.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,
10.1002/qj.49711749710, 1991.
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.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,
10.1175/1520-0469(1979)036<1788:INBAPE>2.0.CO;2, 1979.
Schiller, L. and Naumann, A.: Über Die Grundlegenden Berechnungen Bei Der
Schwerkraftaufbereitung, Z. Ver. Dtsch. Ing, 77, 318–320, 1933.Takahashi, T., Nagao, Y., and Kushiyama, Y.: Possible High Ice Particle
Production during Graupel–Graupel Collisions, J.
Atmos. Sci., 52, 4523–4527,
10.1175/1520-0469(1995)052<4523:PHIPPD>2.0.CO;2, 1995.
Vardiman, L.: The Generation of Secondary Ice Particles in Clouds
by Crystal–Crystal Collision, J. Atmos.
Sci., 35, 2168–2180,
10.1175/1520-0469(1978)035<2168:TGOSIP>2.0.CO;2, 1978.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, 10.1016/j.ijmultiphaseflow.2015.05.010, 2015.