The root of the decades-long scatter lies in an experimental trap that is easy to describe but fiendishly hard to escape. To measure how creep rate changes with temperature, researchers need to compare deformation at different temperatures. The traditional approach was to prepare several snow samples that were supposedly identical and test each one at a different temperature. The problem is that no two snow samples are ever truly identical. Snow is a two-phase material of sintered ice crystals and air, and its densification depends intimately on the spatial arrangement of the ice skeleton. Recent work by the same group showed that even a 5 percent difference in density between two samples, which is common in practice, can shift the strain rate by about 50 percent, equivalent to a temperature difference of roughly 4 degrees. Add to that the temperature fluctuations of typical cold rooms, which reach plus or minus 0.6 degrees Celsius due to defrost cycles and heating from X-ray sources, and the uncertainty swamps the signal being measured.

Védrine’s team solved the problem with an elegant inversion of the protocol: instead of comparing different samples at different temperatures, they tested the same sample at a sequence of temperatures. The centerpiece is a newly developed thermo-mechanical cell, a slim aluminium cylinder 20 millimeters in diameter with walls just half a millimeter thick so that X-rays can pass through. Two Peltier thermoelectric modules at the top and bottom of the cell control the sample temperature independently of the surrounding cold room, keeping temperature differences within the cell on the order of a thousandth of a kelvin, with absolute accuracy better than half a degree. A hollow piston applies a controlled load, a miniature load cell beneath the sample verifies that the force stays constant and that the piston moves without friction, and a laser sensor tracks the vertical displacement of the sample with a resolution of one micrometer. The entire assembly rotates inside an X-ray tomograph, which images the evolving microstructure every eight hours.

The experimental protocol was deliberately gentle. Natural snow collected in the field was stored for months at minus 20 degrees Celsius, where it matured into decomposing and fragmented precipitation particles, then sieved and sintered for at least two weeks. Five samples with initial densities between 278 and 320 kilograms per cubic meter were loaded with a stress of 1.25 kilopascals, comparable to the overburden of a 45-centimeter layer of snow of that density. Each experiment lasted five days and comprised five one-day temperature steps from minus 6 to minus 18 degrees Celsius in 3-degree increments. Because the applied stress was low, the samples deformed only slightly, with a mean final vertical strain of just 4 percent, keeping microstructural evolution to a minimum. Tomography confirmed that density changed by less than 4 percent and specific surface area by less than 7 percent over the whole test, meaning the ice framework remained essentially unchanged while the temperature varied.

One subtlety remained: even a small amount of densification causes geometric hardening, an increase in the reference stress of the ice matrix that slows creep even at constant temperature. To disentangle this effect from genuine temperature sensitivity, the team used a state-of-the-art snow viscoplastic model that describes strain rate as a function of stress, solid fraction, and a temperature-dependent prefactor. By dividing the measured strain rate by the mechanical component, they isolated the pure temperature response. The approach proved robust: neglecting geometric hardening entirely would have produced an activation energy about 10 percent too high, but within the physically meaningful range of model parameters the estimated temperature sensitivity varied by only 3 percent. Digital volume correlation applied to the tomographic images independently confirmed that the laser-based displacement measurements faithfully captured the bulk deformation of the samples.

The results are striking. When the corrected strain-rate prefactor is plotted against inverse temperature on a logarithmic scale, the data from all five samples collapse onto a single master curve, revealing a clear Arrhenius-type relationship, but with a kink. A single Arrhenius fit gives an apparent activation energy of 105 plus or minus 6 kilojoules per mole, yet it fails to capture a change in slope around minus 12 degrees. A two-regime fit is statistically preferred according to the Akaike Information Criterion and yields an activation energy of 126 plus or minus 6 kilojoules per mole above a transition temperature of minus 13 plus or minus 1 degree Celsius, and 51 plus or minus 18 kilojoules per mole below it. In other words, snow creeps far more readily per degree of warming near the melting point than in colder conditions, and the transition between the two behaviors is sharp.

Remarkably, this two-regime pattern mirrors what has long been known for polycrystalline ice. Classic creep studies of glacier ice dating back to Glen’s work in the 1950s reported activation energies of roughly 60 kilojoules per mole below minus 10 degrees and around 130 kilojoules per mole above, values that bracket the new snow measurements almost perfectly. This is intriguing because earlier research by Védrine and colleagues showed that ice inside snow does not deform like dense polycrystalline ice, since the pore network must accommodate the deformation and creep is dominated by basal glide, the sliding of crystal planes along their base. The authors caution that the apparent similarity may be partly deceptive. When they allowed the stress exponent to vary with temperature, following single-crystal measurements of basal slip reported by Jones and Brunet in 1978, the entire dataset could instead be described by a single effective activation energy of 78 plus or minus 6 kilojoules per mole, consistent with basal-slip creep. Part of the apparent bilinear behavior may therefore reflect a temperature-dependent stress exponent rather than a true change in the underlying deformation mechanism.

Whatever the microphysical explanation, the practical implications for snowpack modeling are immediate. The team compared their new parametrization with the settlement laws used in two widely deployed detailed snowpack models, Crocus and Snowpack. The Brun et al. 1992 formulation implemented in Crocus corresponds to an activation energy of about 62 kilojoules per mole, close to the new low-temperature value but far below the high-temperature one, and it exhibits the weakest temperature sensitivity of all. The Kojima 1975 law for fresh snow in Snowpack behaves similarly at low densities. Only the new parametrization captures two distinct temperature regimes. In simulations of a one-meter isothermal snowpack with an initial density of 100 kilograms per cubic meter, the differences matter: after seven days, a snowpack held at minus 20 degrees Celsius is predicted to be 20 to 34 percent taller than one at minus 6 degrees, depending on which law is used, and the new measurements generally imply faster compaction at low densities than either existing formulation predicts.

Accurate snow compaction is not a niche concern. Settlement rates shape snow stratigraphy, which in turn governs avalanche forecasting, estimates of snow water equivalent for hydrology and reservoir management, and even the interpretation of polar ice cores, where firn densification laws link trapped air to paleoclimate records. The study’s authors note that their experiments did not extend above minus 6 degrees Celsius, leaving open whether the creep rate accelerates nonlinearly as the melting point approaches, as often reported for polycrystalline ice, or whether a basal-slip-dominated material like snow behaves differently. They also emphasize that dedicated experiments are still needed to constrain how the stress exponent of basal slip varies with temperature. But the methodological lesson is already secure: by holding the microstructure fixed and letting only the temperature change, the Grenoble team has turned a seventy-year-old scatter plot into a coherent physical law, and given snowpack models their most rigorous temperature calibration to date.

Subject of Research: Temperature dependence of snow viscoplastic creep measured through controlled in-tomograph experiments

Article Title: Temperature sensitivity of snow viscoplasticity: Evidence from controlled creep experiments

Article References: Temperature sensitivity of snow viscoplasticity: Evidence from controlled creep experiments. (n.d.). https://doi.org/10.5194/tc-20-5813-2026

Image Credits: AI Generated

DOI: 10.5194/tc-20-5813-2026

Keywords: snow, viscoplasticity, creep, activation energy, Arrhenius law, snowpack modeling, X-ray tomography, ice mechanics, avalanche forecasting, The Cryosphere, Peltier temperature control, firn densification