Letters from the SkySnow crystals, grown here from water vapour
A snow crystal starts as a speck of ice in a cloud and grows from the vapour around it for half an hour or so as it falls. Its shape records the temperature and humidity it met along the way. On this page a crystal grows from the published equations of a lattice model, one water-vapour step at a time.
Snowflakes are letters sent from heaven.Ukichiro Nakaya
This page needs WebGL2 with float render targets.
Part one · the crystal
Ice, written by the weather
This page needs WebGL2 with float render targets.
View
01 · 1611
Why six?
In 1611 Johannes Kepler gave his patron a small book for the new year, On the Six-Cornered Snowflake, asking why every snowflake has six points. He thought about stacked spheres and could go no further. By 1929 X-ray studies had shown how water molecules lock together in ice: in puckered hexagonal rings, layer on layer.
Grown slowly, ice takes the shape of its own lattice, a hexagonal prism. Here a single cell of ice sits in air at −15 °C that holds a little more vapour than ice can bear, and grows into a plate.
ring spacing a = 4.52 Å, layer c = 7.36 Å
02 · 1885
The snowflake man
Wilson Bentley, a farm boy from Jericho, Vermont, fixed a microscope to a bellows camera and in 1885, aged 19, made the first photograph of a single snow crystal. He kept at it for 46 years and more than 5,000 crystals; the book he made with William Humphreys in 1931 holds over 2,000 of them, white on black.
This is the same model drawn in his style. The fine lines inside the arms are ridges where ice piled up a little thicker as the edge swept past.
03 · instability
Why it branches
A crystal feeds on vapour that diffuses in from the air around it, and it eats faster than the air can refill. The vapour thins out around the crystal. Any corner that sticks out reaches into richer air, grows faster and sticks out further. W. W. Mullins and R. F. Sekerka worked out this runaway in 1964; on a snow crystal it turns the six corners of a plate into six arms, and every bump on an arm into a side branch.
This view shows the vapour: dark where the crystal has drawn it down, bright where the cloud is still full.
vapour at the tips … of the cloud’s
04 · kinetics
Why it has flat faces
Branching has an opponent. A molecule landing on a flat face of ice has nowhere to hold on, so flat faces grow slowly and stay flat. When vapour is scarce the facets win and you get a hexagonal plate; with more vapour the corners outrun them, first as ribs and sectors, then as branches, then as ferns.
Move the humidity while the crystal grows. Each change is written into the edge it reaches.
form …
05 · 1936
Nakaya’s diagram
On 12 March 1936, in a cold room at Hokkaido University, Ukichiro Nakaya grew the first snow crystal in a laboratory, on a single rabbit hair. Varying temperature and humidity, he found that the shape depends sensitively on temperature. His diagram, later extended by Kobayashi (1961), runs: plates near −2 °C, columns and needles near −5, plates and stars again from −10 to −22, columns below that.
Tap the diagram to grow a crystal there. In the column bands the crystal grows along its axis; switch to 3D to see it from the side. Nobody has fully explained why the habit flips with temperature like this; it is still studied today.
here …water-saturated air at −15 °C is … supersaturated over ice
06 · the journey
Letters from the sky
A falling crystal passes through air of different temperature and humidity, and each stretch adds a different kind of growth to its edge. Nakaya saw that you can read the journey back from the shape. This view colours each part of the crystal by the kind of growth the air called for when that part froze, with a faint ring every equal stretch of time.
The six arms come out alike for the same reason: all six made the same journey at the same moment.
07 · chance
Twins and strangers
In the model so far every arm is identical to the last digit, because nothing distinguishes one from another. Real air is never that even. Add a little turbulence, small random flickers in the vapour from cell to cell, and the arms begin to differ in their details while keeping the same overall plan.
Two crystals that fell side by side through the same air can look much alike. Two that took different paths almost never will, because the number of possible journeys is enormous.
largest arm-to-arm difference …
08 · the cloud chamber
Your cloud
Tap or drag on the cloud chart under the crystal to change the air it grows in, at any moment. Or draw a path and let the crystal fall along it.
Model parameters
The six published sets the chart chooses from. Picking one here fixes the parameters until you touch the chart again.
form …diameter …supersaturation over ice …water saturation here …model steps …column length …
Part two · working
How the crystal is grown
The model. Janko Gravner and David Griffeath (2008) treat a flat snow crystal as a hexagonal lattice of cells; this page takes a cell to be 6 µm across. A cell holds ice, a thin quasi-liquid film on the surface, and vapour. Every step the vapour diffuses; at the crystal’s edge part of it freezes and the rest joins the surface film; edge cells join the crystal when enough has gathered, and a little melts back. Whether an edge cell joins depends on how many ice neighbours it has: tips and flat faces (one or two neighbours) need more mass than nooks (three), which is how the model separates facets from branches.
On the GPU. The lattice is stored in axial coordinates in a float texture, 512 cells across on a computer and 384 on a phone, inside a hexagonal domain whose rim is held at the cloud’s vapour density. One step is three shader passes. Neighbour sums are added in an order that does not depend on direction, so in single precision the six arms stay identical to the last bit.
The cloud. Saturation vapour pressures over ice and supercooled water follow Murphy and Koop (2005). The chart’s vertical axis is the excess vapour density over ice saturation, as in Nakaya’s and Kobayashi’s diagrams. Each point of the chart picks one of six parameter sets that Gravner and Griffeath published, each of which grows a form seen in nature: hexagonal plate, plate with ribs, sectored plate, stellar plate, stellar dendrite and fernlike dendrite. That choice follows the observed diagram; the model itself knows nothing about temperature.
Columns. The lattice is two-dimensional, so it only grows the flat habits. In the column bands the lattice runs at half the vapour and a column length is added along the axis at a rate set by the humidity, peaking near −5 °C. The 3D view draws that column with the lattice crystal at each end, which is how capped columns form.
Drawing. Ice is clear, so the photo view renders it as a microscope sees it: the crystal’s outline and the ridges of thicker ice bend light from two coloured lamps. The record view stores, for every cell, the moment it froze and the temperature and humidity at that moment.
One step, at an edge cell x (Gravner–Griffeath)
d′(x) = (1/7) Σy∈N(x) d(y)
b ← b + (1−κ)d, c ← c + κd, d ← 0
join if n ≤ 2 and b ≥ β; n = 3 and (b ≥ 1 or Σd < θ, b ≥ α); n ≥ 4
b ← (1−μ)b, c ← (1−γ)c
The cloud (Murphy & Koop)
ln ei = 9.5504 − 5723.3/T + 3.5307 ln T − 0.0072833 T
σ = e/ei − 1
Checks
Computed when the page loads, except where marked.
Quantity
This page
Reference
Part three · people
Four centuries of looking closely
Sources
J. Kepler, Strena seu de nive sexangula (On the Six-Cornered Snowflake), Frankfurt (1611).
W. H. Barnes, “The crystal structure of ice between 0 °C and −183 °C”, Proc. R. Soc. A 125, 670 (1929).
W. A. Bentley and W. J. Humphreys, Snow Crystals, McGraw-Hill (1931).
U. Nakaya, Snow Crystals: Natural and Artificial, Harvard University Press (1954).
T. Kobayashi, “The growth of snow crystals at low supersaturations”, Phil. Mag. 6, 1363 (1961).
W. W. Mullins and R. F. Sekerka, “Stability of a planar interface during solidification of a dilute binary alloy”, J. Appl. Phys. 35, 444 (1964).
C. Magono and C. W. Lee, “Meteorological classification of natural snow crystals”, J. Fac. Sci. Hokkaido Univ. Ser. VII, 2, 321 (1966).
D. M. Murphy and T. Koop, “Review of the vapour pressures of ice and supercooled water for atmospheric applications”, Q. J. R. Meteorol. Soc. 131, 1539 (2005).
K. G. Libbrecht, “The physics of snow crystals”, Rep. Prog. Phys. 68, 855 (2005).
J. Gravner and D. Griffeath, “Modeling snow crystal growth II: A mesoscopic lattice map with plausible dynamics”, Physica D 237, 385 (2008).
K. Kikuchi, T. Kameda, K. Higuchi and A. Yamashita, “A global classification of snow crystals, ice crystals, and solid precipitation based on observations from middle latitudes to polar regions”, Atmos. Res. 132–133, 460 (2013).
K. G. Libbrecht, Snow Crystals: A Case Study in Spontaneous Structure Formation, Princeton University Press (2022).