Hallett-Mossop Process & Secondary Ice Production: How Rime Splintering in Mixed-Phase Clouds Multiplies Ice Crystals and Triggers Rapid Glaciation

1. The Glaciation Cascade: An Afternoon by the Coast

Stand on a sun-drenched headland in late August, and the sky presents an illusion of stately permanence. The air at ground level is thick, maritime, and fragrant with the sharp geosmin of parched soil and the brine of tidal shallows. To the west, swollen thermals loft invisible parcels of moist air thousands of metres into the troposphere. At first, the emerging cumulus clouds appear as hard-edged sculptures: brilliant white domes of condensed liquid droplets with crisp, cauliflower-like contours, their boundaries carved sharply against the azure backdrop.

   [ CRISP CUMULUS TOP ]               [ GLACIATED FIBROUS ANVIL ]
   Liquid Supercooled Droplets          Secondary Ice Splinter Cascade
   (Sharply defined boundaries)         (Diffuse, fibrous cirriform veil)
            \                                      /
             \----->  Hallett-Mossop Zone  ------->
                      (-3°C to -8°C Layer)

Then, over the span of barely ten minutes, something extraordinary and violent occurs in the cloud’s upper interior. The rigid, boiled-cotton perimeter of the cloud dome softens. Its billowed ridges dissolve into a diffuse, fibrous veil, fraying at the edges like brushed silk. The barometric pressure drops by two millibars in a quiet, rhythmic flutter. The wind at your collar veers twenty degrees to the northwest, suddenly chilled by dozens of degrees Celsius, carrying the cold metallic scent of melting ice aloft.

High above, an unseen chain reaction has consumed the cloud. Millions of tons of liquid water, poised in an unstable liquid state well below zero degrees Celsius, have frozen almost instantaneously into a blinding shroud of crystalline ice. A torrential curtain of virga falls from the cloud base, swiftly reaching the ground as fat, heavy raindrops.

To classical thermodynamics, this abrupt metamorphosis ought to be an impossibility. The atmosphere does not contain enough freezing seeds to ignite such a rapid transformation. Yet here, playing out before the naked eye, is one of the most elegant, destructive, and perplexing feedback loops in atmospheric science: the secondary ice multiplication cascade.

2. What’s Actually Happening: The Ice Multiplication Paradox

To understand why this sudden freezing mystified meteorologists for half a century, one must first dismantle a common misconception: pure water does not automatically turn into ice at $0^\circ\text{C}$. In the clean, buoyant updrafts of maritime cumulus clouds, tiny water droplets can remain liquid down to temperatures as low as $-38^\circ\text{C}$—a metastable condition known as supercooling.

For a supercooled water drop to freeze at milder sub-zero temperatures (between $-2^\circ\text{C}$ and $-15^\circ\text{C}$), it requires an external template: a microscopic aerosol known as an Ice Nucleating Particle (INP). These rare particles—often grains of desert mineral dust, biological fragments, or soot—possess a lattice geometry that mimics the hexagonal symmetry of an ice crystal, providing a stable foundation upon which liquid water molecules can arrange themselves into a solid lattice.

Herein lies the profound historical puzzle known to meteorologists as the Ice Multiplication Paradox. When atmospheric research aircraft equipped with optical array probes flew through developing maritime convective clouds in the 1960s and 1970s, they measured primary ice nucleating particles in concentrations of roughly $1\text{ to }10\text{ particles per cubic metre}$ of air. Yet, within minutes of the cloud top reaching modest sub-zero altitudes, the observed concentration of pristine ice crystals surged to $10,000\text{ to }100,000\text{ crystals per cubic metre}$—a discrepancy of four to five orders of magnitude.

Where were these millions of unheralded ice crystals coming from?

The answer arrived in 1974 through a series of landmark laboratory experiments conducted by atmospheric physicists John Hallett and S. Colin Mossop at the Commonwealth Scientific and Industrial Research Organisation (CSIRO) and later detailed in their seminal research published through the World Meteorological Organization and leading journals. Hallett and Mossop demonstrated that the cloud was not relying solely on primary nucleation seeds brought up from the ground. Instead, under strict environmental conditions, existing ice particles were acting as microphysical factories, mechanically shattering water drops to spawn hundreds of daughter ice splinters in a self-sustaining chain reaction.

3. The Science of the Hallett-Mossop Process

The Hallett-Mossop (H-M) process, often termed rime splintering, is a secondary ice production mechanism that operates within a remarkably narrow microphysical envelope. It does not occur across the entire cloud, but is instead confined to a specific thermal and kinematic sweet spot:

1. The Strict Thermal Window: Secondary ice production is active exclusively between $-3^\circ\text{C}$ and $-8^\circ\text{C}$, exhibiting a sharp, Gaussian-like peak in efficiency at precisely $-5^\circ\text{C}$. Warmer than $-3^\circ\text{C}$, accreted water freezes too slowly without building sufficient internal stress; colder than $-8^\circ\text{C}$, the freezing shell is too ductile and resilient to shatter.
2. The Bimodal Droplet Spectra: The cloud must possess a broad droplet size distribution featuring the simultaneous coexistence of: - Large supercooled cloud droplets with diameters $d > 24\,\mu\text{m}$ (which deliver the mass and liquid volume required for explosive freezing), and - Small cloud droplets with diameters $d < 13\,\mu\text{m}$ (which pack into the interstices of the accreted riming structure, creating a rigid mechanical anchor).
3. Accretor Kinematics: An existing ice particle—typically a millimeter-sized graupel pellet (soft hail)—must be falling through the updraft at a relative impact velocity between $0.5\text{ m s}^{-1}\text{ and }3.0\text{ m s}^{-1}$. If the velocity is too low, the droplet coalesces without fragmentation; if the velocity is too high, kinetic impact energy dissipates the internal pressure prematurely.

The Microphysics of Hydrostatic Rupture

When a supercooled water droplet ($d > 24\,\mu\text{m}$) collides with the frozen surface of a falling graupel particle within the $-5^\circ\text{C}$ zone, it undergoes asymmetric outward-inward freezing. The base of the droplet in direct contact with the cold ice substrate freezes first. Concurrently, heat transfer to the surrounding sub-zero air causes a rigid, hemispherical outer ice shell to freeze across the droplet's exterior surface within milliseconds, encapsulating a core of remaining liquid water.

As the trapped interior water continues to freeze, it undergoes the well-known anomalous density transition of water: expanding by approximately $9\%$ in volume as its molecules rearrange from a disordered liquid into an open hexagonal crystalline structure ($\rho_{\text{liquid}} \approx 1.000\text{ g cm}^{-3} \to \rho_{\text{ice}} \approx 0.917\text{ g cm}^{-3}$).

Because the outer ice shell is rigid, unyielding, and mechanically pinned to the graupel surface, this $9\%$ volumetric expansion generates immense internal hydrostatic pressures exceeding tens of megapascals. When the tangential hoop stress surpasses the critical tensile strength of the thin outer ice shell, the shell ruptures catastrophically. The explosive release ejects tiny fragments of the freezing shell into the ambient air stream—microscopic acicular ice needles measuring between $5\,\mu\text{m}\text{ and }50\,\mu\text{m}$ in length.

For every milligram of rime accreted under optimal $-5^\circ\text{C}$ conditions, the Hallett-Mossop process liberates approximately $350\text{ secondary ice splinters}$.

Mathematical Formulation of Secondary Splinter Generation

To quantify this runaway microphysical process within numerical weather models, atmospheric dynamicists express the time-dependent production rate of secondary ice crystals ($\frac{dN_s}{dt}$) per unit volume of air through the following kinematic accretion integral:

$$\frac{dN_s}{dt} = \Psi(T) \cdot \pi r_g^2 \cdot |v_g - v_d| \cdot E_{\text{coll}} \cdot n_L(d > 24\,\mu\text{m})$$

Where the physical variables represent: - $\Psi(T)$: The empirical splinter production efficiency function (splinters generated per unit mass of accreted rime, peaking sharply at $\sim 3.5 \times 10^5\text{ splinters g}^{-1}\text{ at } T = -5^\circ\text{C}$, dropping linearly to zero at $-3^\circ\text{C}$ and $-8^\circ\text{C}$). - $r_g$: The effective cross-sectional radius of the falling graupel accretor ($\text{m}$). - $|v_g - v_d|$: The differential terminal fall velocity between the heavy graupel particle and the suspended supercooled liquid droplets ($\text{m s}^{-1}$). - $E_{\text{coll}}$: The hydrodynamic collision-coalescence efficiency (the probability that a droplet in the geometric sweep path actually impacts the graupel rather than being deflected along aerodynamic streamlines). - $n_L(d > 24\,\mu\text{m})$: The number concentration of large supercooled droplets per unit volume ($\text{m}^{-3}$).

Worked Numerical Example: The 10-Minute Glaciation Engine

Consider a vigorous maritime convective updraft over the North Atlantic, sampled by airborne instruments operating under protocols defined by the UK Met Office Cloud Physics Research Branch.

Let us isolate the microphysical kinetics occurring around a single graupel particle and trace its multiplier effect across a cubic metre of cloud air:

1. Input Cloud Parameters at $T = -5^\circ\text{C}$:

• Graupel radius ($r_g$): $1.0\text{ mm} = 1.0 \times 10^{-3}\text{ m}$
• Graupel cross-sectional area ($A = \pi r_g^2$): $3.14 \times 10^{-6}\text{ m}^2$
• Differential fall velocity ($|v_g - v_d|$): $2.0\text{ m s}^{-1}$
• Collision efficiency ($E_{\text{coll}}$): $0.80$
• Large supercooled droplet concentration ($n_L$, $d > 24\,\mu\text{m}$): $5.0 \times 10^7\text{ droplets m}^{-3}$ ($50\text{ cm}^{-3}$)
• Splinter generation yield ($\Psi$): $1\text{ secondary splinter per } 250\text{ large droplet impacts}$ ($4.0 \times 10^{-3}\text{ splinters/droplet}$)

2. Swept Volume Rate:

The geometric volume of cloud air swept out by the falling graupel grain per second is:

$$\dot{V}_{\text{sweep}} = A \cdot |v_g - v_d| = (3.14 \times 10^{-6}\text{ m}^2) \times (2.0\text{ m s}^{-1}) = 6.28 \times 10^{-6}\text{ m}^3\text{ s}^{-1}$$

3. Droplet Accretion Rate:

Accounting for hydrodynamic collision efficiency ($E_{\text{coll}} = 0.80$), the number of large droplets accreted per second ($\dot{N}_{\text{drop}}$) is:

$$\dot{N}{\text{drop}} = \dot{V}{\text{sweep}} \cdot E_{\text{coll}} \cdot n_L = (6.28 \times 10^{-6}) \times 0.80 \times (5.0 \times 10^7) \approx 251.2\text{ large droplets s}^{-1}$$

4. Splinter Production Rate:

Multiplying the accreted droplets by the secondary splinter yield:

$$\frac{dN_s}{dt} = 251.2\text{ droplets s}^{-1} \times (4.0 \times 10^{-3}\text{ splinters/droplet}) \approx 1.005\text{ splinters per second}$$

A single millimeter-sized graupel grain produces $\approx 1\text{ fresh ice crystal every second}$.

While one splinter per second may sound modest, consider the convective feedback loop. Within a modest convective cloud updraft of $3.0\text{ m s}^{-1}$, these newly ejected secondary needles are far lighter than the parent graupel. Their terminal velocity is less than $0.1\text{ m s}^{-1}$, meaning they are instantly lofted back up through the $-5^\circ\text{C}$ zone.

As they ascend, they feed on ambient water vapor via the Bergeron-Findeisen process—growing rapidly into full-sized stellar dendrites and rimed pellets within three to five minutes.

Each daughter crystal becomes a new graupel accretor. The population of secondary ice ($N_s$) evolves exponentially according to the classic branching cascade:

$$N_s(t) = N_0 \cdot e^{\lambda t}$$

Where the multiplication growth constant $\lambda \approx 0.008\text{ to }0.015\text{ s}^{-1}$.

At $t = 600\text{ seconds}$ ($10\text{ minutes}$), a starting seed concentration of $N_0 = 1\text{ primary crystal m}^{-3}$ compounds to:

$$N_s(600) = 1 \cdot e^{0.012 \times 600} = e^{7.2} \approx 1,340\text{ crystals m}^{-3}$$

By $t = 900\text{ seconds}$ ($15\text{ minutes}$):

$$N_s(900) = 1 \cdot e^{0.012 \times 900} = e^{10.8} \approx 49,000\text{ crystals m}^{-3}$$

In less than a quarter of an hour, the microphysical regime shifts completely. The liquid droplet reservoir is exhausted, and the cloud transitions into total glaciation.

4. Observer Signatures, Polarimetric Radar & Meteorological Impact

The secondary ice multiplication cascade is not merely a theoretical curiosity of microphysics; it dictates the structural evolution of severe weather, shapes aviation safety protocols, and governs global climate models.

Dual-Polarization Radar Signatures

Modern meteorological networks, such as the NOAA National Severe Storms Laboratory Dual-Pol Radar Network, detect the signature of the Hallett-Mossop process in real time using dual-polarization Doppler radar. Polarimetric radars emit radio pulses with both horizontal and vertical wave orientations, measuring how hydrometeors distort the reflected signal.

When the Hallett-Mossop cascade ignites within the $-3^\circ\text{C}$ to $-8^\circ\text{C}$ isotherm, it produces two distinct polarimetric signatures:

1. Differential Reflectivity ($Z_{dr}$) Peaks: Because the ejected secondary ice fragments grow preferentially as elongated, horizontally-oriented needle crystals (acicular prisms), they reflect stronger horizontal pulses than vertical pulses. This generates a localized column of elevated $Z_{dr}$ ($+2.0\text{ to }+4.0\text{ dB}$) right above the $-5^\circ\text{C}$ altitude.
2. Specific Differential Phase ($K_{dp}$) Columns: As millions of microscopic needles align under aerodynamic forces, they induce a phase shift in the radar beam, creating a narrow, high-density $K_{dp}$ column ($> 1.5^\circ\text{ km}^{-1}$). Meteorologists tracking severe storms treat these polarimetric spikes as unambiguous warnings of rapid cloud electrification, imminent hail growth, and severe downdraft generation.

Electrification and Severe Weather

Secondary ice production is the primary catalyst for lightning generation in convective storms. When rising secondary ice splinters collide with descending graupel pellets in the presence of supercooled liquid water, non-inductive charge separation occurs: the heavier graupel strips electrons from the ice crystals, acquiring a negative charge, while the lighter ice splinters carry positive charge to the top of the cloud. Without the rapid million-fold multiplication of ice splinters provided by the Hallett-Mossop mechanism, thunderstorm electrification would proceed too slowly to generate frequent cloud-to-ground lightning.

Hazards to Modern Aviation

For commercial and military aviation, the Hallett-Mossop thermal zone represents one of the most hazardous flight regimes in atmospheric physics. As detailed by safety advisories from the Federal Aviation Administration (FAA) and the International Civil Aviation Organization (ICAO), aircraft flying through the $-3^\circ\text{C}$ to $-8^\circ\text{C}$ layer encounter a volatile mixture of high supercooled liquid water content and dense ice crystal concentrations.

                                  [ AIRFRAME ICING HAZARD ]
                                  Mixed-Phase Accretion Zone
                                     (-3°C to -8°C Layer)
                                      /                 \
                                     v                   v
                        [ Structural Freezing ]     [ High Ice Water Content ]
                         Rapid glaze ice forms       Ice crystals melt on warm
                         on leading wing edges;      engine turbines, refreeze
                         degrades aerodynamic lift   in compressor stages (core stall)

Supercooled drops strike the airframe, rapidly forming heavy structural rime and clear ice on leading wing edges, while microscopic ice needles bypass engine inertial separators. These ingested crystals melt upon hitting warm engine compressor surfaces, evaporatively cooling the metal until subsequent crystals refreeze into solid ice sheets inside the engine core, causing uncommanded power loss and turbine stalls.

Parameterization in Numerical Weather Prediction

Because global numerical weather prediction models—such as the Integrated Forecasting System (IFS) at the European Centre for Medium-Range Weather Forecasts (ECMWF) and the NOAA/NCAR Weather Research and Forecasting (WRF) model—operate on grid scales of several kilometres, they cannot resolve individual $24\,\mu\text{m}$ droplet collisions directly.

Instead, modelers implement bulk microphysics parameterizations (such as the Morrison, Thompson, or Milbrandt-Yau schemes) that mathematically represent the Hallett-Mossop process as a piecewise empirical threshold function. Accurately tuning this single non-linear rate equation is what allows supercomputers to correctly forecast whether a storm system will produce light stratiform drizzle or explosive convective hail.

5. Practical Outdoor Guidance: Reading the Glaciation Boundary

You do not need a multi-million-pound dual-polarization Doppler radar array or a high-performance research aircraft to witness secondary ice multiplication in real time. Any attentive observer, hiker, gardener, or sailor can track the signature of the Hallett-Mossop process using sensory observations and basic field instruments.

What to Look For in the Sky

1. The Dissolving Boundary (Liquid to Solid Phase Transition): Focus your gaze on the rapidly ascending towers of a Cumulus congestus cloud on a warm afternoon. While the cloud is composed of supercooled liquid droplets, its upper margins will maintain crisp, bulging, hard-edged boundaries that look like freshly sculpted marble. The moment the cloud top punches through the $-3^\circ\text{C}\text{ to }-8^\circ\text{C}$ isotherm (typically between $3,500\text{ and }5,000\text{ metres}$ altitude in mid-latitudes), the hard margins will suddenly blur, turn fibrous, and take on a soft, feathered appearance. This optical transition—from sharp liquid Mie scattering to diffuse crystalline scattering—marks the ignition of the secondary ice cascade.
2. The Emergence of Halo Phenomena: As secondary ice splinters grow into pristine hexagonal plates and columnar prisms, they begin to refract sunlight at exact geometric angles. If the glaciating top of the cloud drifts across the path of the sun, watch for the immediate appearance of Sun Dogs (Parhelia) or a faint $22^\circ\text{ halo}$. These optical phenomena cannot form in liquid clouds; their appearance is definitive proof that crystalline ice has conquered the cloud volume.
3. Virga and Precipitation Cascades: Watch the cloud base directly beneath a glaciating tower. Within five to eight minutes of the cloud top turning fibrous, dark vertical streaks of falling precipitation (virga) will emerge from the base. In maritime air masses, this is the direct surface consequence of the Hallett-Mossop process: riming graupel has grown large enough to overcome the updraft and fall toward the earth.

Instrument Readings to Monitor

• The Barometer: Keep a close watch on a high-resolution digital or aneroid barometer. When an isolated cumulus cloud undergoes explosive glaciation nearby, the latent heat of freezing released aloft ($\approx 3.34 \times 10^5\text{ J kg}^{-1}$ of frozen water) causes a pulse of rapid thermal expansion, driving a localized pressure perturbation: a momentary 1 to 2 millibar flutter, followed by a rapid, steady decline indicating the approach of a convective precipitation downdraft.
• The Thermometer and Anemometer: If you are downwind of a developing cell, watch for a sudden drop in ambient dry-bulb temperature ($3^\circ\text{C}\text{ to }8^\circ\text{C}$) accompanied by a crisp, gusting wind shift. This is the gust front, propelled by air that has been cooled by melting secondary ice aloft and dragged to the surface by the descending precipitation core.

6. Today’s Meteorological Rule of Thumb

When a towering cumulus loses its crisp, carved edges and frays into a soft, fibrous veil, the cloud has ignited its secondary ice engine: expect rain or lightning at the surface within fifteen minutes.

Key Takeaways for the Curious Observer

• Primary ice nucleating particles are astonishingly rare in the lower atmosphere, often numbering fewer than ten per cubic metre.
• The Hallett-Mossop process bridges this deficit by mechanically multiplying ice crystals up to 100,000-fold via the explosive hydrostatic rupture of supercooled droplets freezing onto falling graupel.
• The reaction is exquisitely temperature-dependent, operating solely between $-3^\circ\text{C}$ and $-8^\circ\text{C}$, and reaching maximum splintering efficiency at precisely $-5^\circ\text{C}$.
• Next time you watch a storm cloud fray at its summit, you are not merely watching passive water evaporate—you are witnessing a violent, self-accelerating chain reaction of microscopic explosions that transforms the physics of the sky.