Lander Van Tricht of ETH Zürich and colleagues rebuilt the Global Glacier Evolution Model, known as GloGEM, so that it ingests daily temperature and precipitation rather than monthly averages. They ran both versions across more than 12,500 glaciers larger than five square kilometres, which together hold about 94 percent of the world’s glacier ice, and compared the results against more than 4,000 independent field measurements of annual mass balance. The daily version tracked the observations more closely overall, but its real advantage emerged in extreme years, when the gap between model and reality matters most.

During the record Alpine melt year of 2022, the daily model simulated a mass balance 0.31 metres of water equivalent more negative than the monthly version, landing substantially closer to what glaciologists actually measured. Similar improvements appeared in Western Canada and the United States in 2023 and in Svalbard in 2024, where the daily-monthly difference reached nearly a metre of water equivalent. For anomalies more than two standard deviations from the mean, the daily model cut the root mean squared error from 0.98 to 0.83 metres. In other words, the monthly model systematically underestimates the damage done by heatwaves and overestimates the relief provided by cold spells.

The physical reasons are straightforward. Monthly averaging smooths away the very days on which glaciers lose most of their ice. A single scorching week inside an otherwise mild month barely registers in a monthly mean, yet it can strip metres of snow and expose dark underlying ice, which absorbs more sunlight and melts faster in a runaway feedback. Daily forcing also captures the exact temperature at which precipitation flips between snow and rain, a threshold that determines whether a storm adds mass to a glacier or washes away. And because the daily model updates surface type, snow versus ice, every day rather than once a month, it reproduces the seasonal timing of these transitions far more faithfully.

Those differences cascade into projections. By 2100, the daily model projects global ice volume losses of 26 percent, plus or minus 10, under the low-emission SSP1-2.6 scenario and 47 percent, plus or minus 17, under the high-emission SSP5-8.5 scenario, relative to 2025. That is an additional 5 to 10 percent of mass loss compared with the monthly implementation, equivalent to 20 to 33 extra millimetres of global mean sea-level rise. The divergence grows steadily throughout the century and shows no sign of stabilising by 2100, which matters because glaciers respond on multi-centennial timescales.

The regional pattern is strikingly uneven. In low-relief, ice-cap-dominated regions such as Svalbard, the Russian Arctic and the Antarctic and Sub-Antarctic islands, the narrow elevation range amplifies the sensitivity of melt-albedo and surface-type feedbacks, and the daily model projects up to 25 percent greater ice loss by 2100 than the monthly version. Mid-latitude regions with limited ice, such as Central Europe and the Caucasus, show a transient divergence peaking around mid-century before both models converge on near-total deglaciation. Alaska, Scandinavia and the Southern Andes fall in between, with the daily model typically projecting 2 to 5 percent more loss.

Iceland is the lone exception, and it is instructive. There, heavy precipitation frequently arrives with warm south-westerly air masses from the North Atlantic, while cold northerly periods tend to be dry. Resolving this co-occurrence at daily timescales means more precipitation is classified as rain rather than snow, reducing accumulation. The monthly model smooths over those warm, wet days, overestimates snowfall, and consequently projects about 4 percent more ice loss than the daily version by 2100. A sensitivity experiment in which the monthly model received daily-based precipitation phase confirmed that this partitioning, not melt, drives Iceland’s reversed behaviour.

The study also reframes how glacier models talk to hydrologists. Runoff, defined as the sum of ice melt, snow melt and rainfall over glacierised terrain minus refreezing, is the primary channel through which glacier change reaches downstream water systems. The daily model resolves the seasonal peak in flow and its future shifts: in Western Canada and the United States, peak glacier runoff is projected to arrive roughly 46 days earlier by 2100, a shift the monthly model misrepresents by about 15 days. In High Mountain Asia, the daily model finds that ice melt supplies 70 to 80 percent of late-summer glacier runoff today, against 44 percent in the monthly configuration, a discrepancy with direct consequences for the nearly one-tenth of humanity that depends on glacier-fed water towers.

Daily resolution is not free. The authors report that the daily version takes roughly ten times longer to run than the monthly one, which helps explain why most state-of-the-art global models still use aggregated forcing. But the paper suggests cheaper fixes: correcting the assumption that sub-monthly temperature variability stays constant as the climate warms, and improving precipitation-phase partitioning in monthly models, could remove much of the structural bias without full daily computation. Remotely sensed snowline and snow-cover datasets offer a further route to constrain where models place the transient snowline, a quantity the daily model can be validated against directly.

The broader message is that temporal resolution is not a technical footnote but a modelling choice with measurable consequences for coastal flood exposure, erosion and saltwater intrusion. As heatwaves grow more frequent and more intense, models that average them away will keep underestimating how fast the world’s mountain ice is disappearing. The 20 to 33 extra millimetres of sea-level rise that emerge simply from reading the weather day by day are, the authors argue, a lower bound on what better-resolved physics will reveal as global glacier models continue to converge on the messy, day-to-day reality of a warming atmosphere.

Subject of Research: Effect of daily versus monthly meteorological forcing on global glacier mass loss and sea-level rise projections

Article Title: Resolving short-term meteorological variability amplifies projected glacier loss and sea-level rise

Article References: Van Tricht, L., Huss, M., Zekollari, H., Aguayo, R., Beer, J., van Tiel, M., von der Esch, A., & Farinotti, D. (2026). Resolving short-term meteorological variability amplifies projected glacier loss and sea-level rise. Nature Communications, 17(1), Article 10696. https://doi.org/10.1038/s41467-026-78260-y

Image Credits: AI Generated

DOI: 10.1038/s41467-026-78260-y

Keywords: glaciers, sea-level rise, climate modelling, GloGEM, mass balance, daily forcing, heatwaves, Svalbard, glacier runoff, CMIP6, cryosphere, Nature Communications

The stakes are high. The Greenland Ice Sheet has been responsible for roughly 20 percent of global sea level rise since 2002, and nearly half of that ice loss is moderated by the presence of ice mélange at the fronts of its tidewater glaciers. These floating aggregations of ice act as a physical brake on calving, the dramatic process by which chunks of glacier ice break off and become icebergs. When mélange is thick and rigid, it presses back against the glacier terminus, suppressing the release of large icebergs and slowing the glacier’s seaward flow. When mélange weakens or disappears, that restraint vanishes, and glaciers can accelerate their discharge of ice into the ocean.

Yet scientists have long struggled to explain why ice mélange behaves so differently from one Greenland fjord to another. Observations reveal a remarkable spread of behavior: some fjords see their mélange collapse seasonally each spring, disintegrating over just a few days, while others host mélange that persists year-round, seemingly indefinitely. Understanding what tips a mélange from durable to doomed has been a stubborn puzzle, and one with direct implications for forecasting how quickly Greenland’s glaciers will shed ice in a warming climate.

To crack the puzzle, the research team built a coupled model that brings together two worlds of physics that are usually simulated separately. The first component captures high-resolution ocean circulation within the fjord, including the turbulent plumes of meltwater that rise along glacier fronts and the transfer of heat from the ocean to the ice. The second component models the mélange itself as a dynamic granular material, a flowing collection of discrete ice blocks that jam, shear, and fracture much like sand or gravel. This granular treatment is essential, because ice mélange does not behave like a continuous sheet of ice; its strength emerges from the crowding and friction among thousands of individual bergs pressing against one another.

Running multi-year simulations with this coupled framework, the researchers were able to reproduce the observed seasonal variability of ice mélange in real fjords. In their simulations, the mélange follows an annual rhythm governed by the competition between ice supplied by calving at the glacier front and ice removed by melting and by export out of the fjord. During cold seasons, newly calved bergs consolidate into a rigid mass, and the mélange grows in extent and strength. As conditions warm, the balance shifts, and the fragile structure becomes vulnerable to rapid disintegration.

The crucial insight from the new work is the mechanism of collapse itself. Rather than eroding steadily in step with melt rates, the model reveals that mélange possesses a threshold behavior. Below a critical level of submarine melting, the mélange can persist indefinitely, repairing itself seasonally with fresh calved ice. But once melting from below pushes the system past that threshold, the internal structure of the mélange loses its rigidity. The granular fabric that gives the mélange its buttressing power unravels quickly, and the entire ice field can collapse in a matter of days, mirroring the abrupt springtime disintegrations that observers have documented in Greenlandic fjords.

To test whether this modeled threshold holds in the real world, the team parameterized the collapse condition and compared it against observations from five fjords across Greenland. The comparison showed that the model’s threshold captures the observed spread of mélange behavior, distinguishing fjords where mélange persists from those where it breaks down each year. This kind of validation matters, because it suggests the threshold could serve as a practical diagnostic, a way to assess which of Greenland’s glaciers are currently protected by stable mélange and which are hovering near the tipping point.

The findings carry an unsettling implication. Perennial ice mélange is currently present at the majority of the largest glaciers in Greenland, which means some of the ice sheet’s biggest contributors to sea level rise are being held in check by these floating braces. The new results indicate that many of these barriers may be capable of sudden collapse if subjected to increased submarine melt, the process by which relatively warm Atlantic-sourced waters, circulating deep in the fjords, erode ice from beneath the waterline. Because ocean heat content around Greenland has been rising, and because warm waters are increasingly able to reach the deep channels that feed glacial fjords, the reservoirs of melt potential pressing against these mélange fields are growing.

If a perennial mélange does collapse, the consequences would cascade. Freed from the back-pressure of the mélange, the glacier terminus becomes free to calve more vigorously, and the removal of resistive floating ice can allow the glacier to speed up its flow toward the sea. Faster flow means more ice delivered across the grounding line and into the ocean, accelerating glacier mass loss and, ultimately, sea level rise. The nonlinear nature of the collapse means this acceleration could arrive abruptly rather than gradually, complicating efforts to project Greenland’s contribution to future sea levels on decadal timescales.

The study, which was supported by the National Science Foundation through the GLACIOME project, represents a step forward in treating the ocean, the ice, and the fractured mélange between them as a single interacting system. It also highlights how much remains to be learned: fjord bathymetry, seasonal oceanography, and calving dynamics all shape where each glacier sits relative to its collapse threshold. But the central message is clear. Greenland’s glacial fjords are not destined to lose their ice mélange slowly and predictably. Some of these frozen barriers could hold for years and then vanish in days, and the warming ocean pressing against their undersides is steadily pushing more of them toward that edge.

Subject of Research: Threshold collapse of ice mélange at tidewater glaciers driven by increasing submarine melt

Article Title: Sudden collapse of ice mélange under increasing submarine melt at tidewater glaciers

Article References: Summers, P. T., Robel, A. A., Jackson, R. H., King, M. D., Amundson, J. M., & Burton, J. C. (2026). Sudden collapse of ice mélange under increasing submarine melt at tidewater glaciers. Nature Communications. https://doi.org/10.1038/s41467-026-78423-x

Image Credits: AI Generated

DOI: 10.1038/s41467-026-78423-x

Keywords: Greenland Ice Sheet, ice mélange, tidewater glaciers, submarine melting, calving, sea level rise, granular flow, fjord oceanography, glacier retreat, Nature Communications, climate change, cryosphere

The Himalayan cryosphere, often described as the water tower of Asia, supplies meltwater to billions of people downstream. But the same snowfields and glaciers that store this water also capture atmospheric deposition with remarkable efficiency. Snow scavenging, the process by which falling and lying snow removes aerosol particles and dissolved gases from the air, makes high-altitude snowpacks sensitive recorders of both natural mineral dust and anthropogenic emissions. The review collates studies from the Western, Central and Eastern Himalayas as well as the Tibetan Plateau, sometimes called the Third Pole, and reveals a chemical signature that is anything but pristine.

One of the clearest findings across the compiled datasets is the overwhelming dominance of crustal ions, particularly calcium (Ca2+) and bicarbonate (HCO3−). These species originate largely from the weathering of mineral dust, which is lifted from arid regions of Central and South Asia and deposited onto mountain snowfields. However, the review also documents notable enrichment of trace elements such as zinc (Zn), lead (Pb) and manganese (Mn), which cannot be explained by crustal sources alone. The combination points to a dual influence: natural geological weathering working alongside emissions from industry, fossil fuel combustion and biomass burning transported over long distances.

Spatial patterns emerging from the review are striking. Pollutant concentrations are consistently higher on the southern slopes of the Himalayas than on the northern side. This asymmetry reflects the geography of emission sources. The southern slopes face the densely populated Indo-Gangetic Plain, where industrial activity, vehicular traffic and agricultural burning generate vast quantities of aerosols and heavy metals. The northern slopes, facing the more sparsely populated Tibetan Plateau, receive comparatively less direct pollution, although long-range transport still delivers contaminants from Central Asia and beyond.

Seasonal variability adds another layer of complexity. The review indicates that accumulation of pollutants in snow is enhanced during non-monsoon periods, when dry conditions and reduced precipitation allow contaminants to build up in the snowpack without being flushed away. During the Indian summer monsoon, heavy rainfall and snowfall act as efficient scavenging agents, removing large amounts of material from the atmosphere but also depositing it onto the cryosphere in concentrated pulses. The interplay between monsoon dynamics and western disturbances, weather systems that bring moisture from the Mediterranean and Atlantic regions, plays a key role in transporting contaminants from South and Central Asia into the high mountains.

The implications of this contamination extend beyond the chemistry of the snow itself. Trace elements such as lead, mercury and cadmium are toxic even at low concentrations, and their presence in glacier snow raises concerns about what happens when that snow melts. Meltwater carries dissolved and particulate contaminants into downstream rivers, potentially affecting drinking water quality, aquatic ecosystems and agricultural soils. Studies cited in the review, including work on the Rongbuk Glacier on Mount Everest and the Zhadang Glacier on the southern Tibetan Plateau, have documented mercury export in glacier runoff, highlighting a direct pathway from atmospheric deposition to human exposure.

Black carbon and light-absorbing particles further complicate the picture. Although the review focuses primarily on ionic species and trace elements, it draws on related research showing that dark particles deposited on snow surfaces reduce albedo, the reflectivity of the snow, causing it to absorb more solar radiation and melt faster. This accelerates glacier retreat across the region, which has been documented extensively from the Dokriani Glacier in Garhwal to the Machoi Glacier in Kashmir and the glaciers of Sikkim. The combination of chemical contamination and physical darkening creates a feedback loop in which pollution both records and drives cryospheric change.

Melt-induced fractionation presents another technical challenge for interpreting the archive. As snowpacks warm and melt, ions and trace elements are not lost uniformly. Some species are eluted preferentially in the first meltwater fractions, while others remain concentrated in the remaining snow. This means that the chemical record preserved in a snowpit or ice core may not perfectly reflect the original atmospheric deposition. The review acknowledges this complication, drawing on studies from the Alps and other mountain regions that have quantified how meltwater redistributes solutes within the snowpack, and notes that such processes must be accounted for when reconstructing past atmospheric conditions from glacial archives.

Despite the growing body of research, the review identifies critical data gaps, particularly in the Eastern Himalayas. While the Western Himalayas, the Central Himalayas and the Tibetan Plateau have been relatively well studied, the Eastern Himalayan region, spanning parts of Arunachal Pradesh, Sikkim, Bhutan and northeastern India, remains underrepresented in the published record. This gap is significant because the Eastern Himalayas receive different atmospheric circulation patterns and may experience distinct pollution transport pathways. Without standardised monitoring across all sub-regions, it is difficult to build a complete picture of contaminant fluxes into the cryosphere or to predict how these will change under future emission scenarios and a warming climate.

The authors emphasise the need for expanded and standardised monitoring programmes to better assess future ecological and health risks associated with cryospheric contamination. As the climate warms and glaciers retreat, the pollutants stored in snow and ice will be released into meltwater over increasingly short timescales. Understanding the full inventory of contaminants locked in the Himalayan cryosphere, and the mechanisms by which they arrived there, is therefore not merely an academic exercise. It is a prerequisite for protecting the water security and public health of the more than a billion people who depend on the rivers born from these mountains. The review serves as both a synthesis of two decades of research and a call to action for the scientific community to fill the remaining gaps before the archive itself disappears.

Subject of Research: Major ionic species and trace element contamination in Himalayan snow and glacier ice

Article Title: A review of major ionic species and trace elements in the Himalayan cryosphere

Article References: Sharma, K. D., Dutta, J., & Jyethi, D. S. (2026). A review of major ionic species and trace elements in the Himalayan cryosphere. Environmental Monitoring and Assessment, 198(9), Article 991. https://doi.org/10.1007/s10661-026-15840-x

Image Credits: AI Generated

DOI: 10.1007/s10661-026-15840-x

Keywords: Himalayas, cryosphere, snow chemistry, trace elements, major ions, glacier meltwater, atmospheric deposition, Tibetan Plateau, black carbon, heavy metals, monsoon, pollution transport

The study, published in Nature Geoscience by C. Y. Fu of the University of Cambridge’s Scott Polar Research Institute, M. B. Osman of the University of Cambridge, and M. A. Aquino-López of the Centro de Investigación en Matemáticas in Mexico, tackles a problem that has plagued the latest generations of climate models. Although winter sea ice extent shows the least year-to-year variability in the observational record, the disagreement between models is widest precisely at this time of year, with the largest spread seen in March. The absolute bias in wintertime sea ice in the CMIP5 and CMIP6 model archives is substantially greater than the bias in September, meaning that the models collectively paint too much ice on the winter Arctic Ocean — and nobody knows exactly which models are to blame.

The consequences of this blind spot are not academic. Winter ice cover regulates the release of heat and moisture from the ocean to the atmosphere, and its reduction has been implicated in colder winters and heavier snowfall across Eurasia and North America, intensified Northern Hemisphere cyclones and shifting storm tracks, and extreme precipitation events in the Mediterranean. The 2025 annual maximum Arctic sea ice extent was the lowest on record, bringing the long-term winter decline into sharp focus. Yet efforts to refine projections have concentrated overwhelmingly on the summer minimum, potentially biasing model evaluation toward the season the models already handle relatively well.

Fu and colleagues’ insight was that the past can serve as an out-of-sample test of model performance. During the mid-Holocene, roughly six thousand years ago, orbitally driven changes enhanced summer insolation in the Northern Hemisphere, producing an increased annual mean radiative forcing in the Arctic comparable to the greenhouse gas forcing anticipated in the coming decades. Previous analyses of the Palaeoclimate Modeling Intercomparison Project have shown that Arctic warming in the mid-Holocene and in the future is driven by common physical processes from autumn through winter, with similar contributors to the multi-model spread across all seasons. Models that capture past sea ice dynamics, the reasoning goes, should also perform well in future projections — both periods being effectively unseen by the models during their development.

The obstacle has always been quantitative data. Palaeoclimate applications of sea ice proxies have remained largely qualitative, hampered by sparse records and large reconstruction uncertainties. To overcome this, the team compiled a new Arctic-wide database of mid-Holocene biomarker records based on the ln(PIP25) index, a fully quantitative sea ice proxy built on the relative concentrations of IP25 — an ice proxy with 25 carbon atoms produced by sea ice diatoms — and associated phytoplankton sterols such as brassicasterol and dinosterol, which serve as open-water indicators. The resulting compilation of 46 observations, spanning latitudes from roughly 54 to 83 degrees north, is to the researchers’ knowledge the largest assembly of palaeo-sea ice biomarker records ever assembled.

The technical heart of the study lies in a recently developed proxy system model called BaySIC, which translates simulated sea ice concentrations from 27 climate models of the PMIP3 and PMIP4 mid-Holocene ensembles into proxy units for direct, fully quantitative comparison with the sediment observations. Crucially, the forward modelling accounts for spatiotemporal variations in proxy seasonality: because the ln(PIP25) index carries a seasonal bias toward late winter and early spring, when sea ice extent reaches its annual maximum, the team used climatological sea ice concentrations from the months corresponding to each model’s simulated timing of local sea ice breakup rather than fixed calendar months. This matters enormously in high-latitude regions, where the timing of ice breakup can shift substantially between the modern climate and the mid-Holocene.

On top of this framework, the researchers built a Bayesian palaeoclimate model weighting scheme that adjusts the contribution of each model to the ensemble mean based on proxy evidence. The scheme maximizes the likelihood that the weighted multi-model mean of predicted proxy values reproduces the observations across all core sites, while implicitly accounting for model interdependence — related models that contribute similar information share weight, effectively penalizing redundancy without requiring explicit skill-dependence metrics. The team first derived prior weights from the models’ performance against historical observations from 1950 to 1999, then updated these with the mid-Holocene data to obtain posterior weights that maximize proxy-model agreement across both intervals.

The results are striking. In validation experiments, the weighted multi-model mean aligned with observations far more closely than the unweighted mean, and the uncertainty range was reduced by approximately 6.3 million square kilometres — an 85 percent narrowing — bringing the projection uncertainty down to the level of interannual variability. When historical and mid-Holocene constraints were combined, the weights effectively eliminated the mean positive bias in winter sea ice extent, leaving a residual of just 0.016 million square kilometres, a percentage error of 0.099 percent. The mid-Holocene period alone yielded more effective constraints than the historical record, supporting the notion that out-of-sample evaluation based on past warm periods provides information about model reliability beyond anything the short observational record can offer.

Applying the same weights to future projections of the same models under three emissions scenarios — SSP1, SSP2 and SSP5 in CMIP6, corresponding to RCP2.6, RCP4.5 and RCP8.5 in CMIP5 — the weighted mean winter sea ice extent was consistently lower than the unweighted estimate, by roughly 4 to 6 percent for the decade 2090 to 2099. For the high-emissions SSP5 scenario, the end-of-century projection was tightened by over 10 million square kilometres. The spatial pattern is equally consequential: the weighted projections show higher probabilities of sea ice loss in the Barents, Greenland and Labrador Seas, corresponding to northward ice-edge shifts of up to 5.5 degrees of latitude — regions already identified as vulnerable to Atlantic water intrusion and associated atmosphere-ocean-sea ice feedbacks.

The cascading implications reach far beyond the polar circle. Reduced winter sea ice in the Barents Sea has been linked to cold extremes in Eurasia through local atmospheric heating and weakening of the polar vortex, as well as altered cyclone tracks over the Nordic Seas and increased rainfall in the Mediterranean and East Asia. The team’s regression analyses suggest that late winter and early spring snowfall over Southeastern Greenland may be underestimated by more than 30 percent in current projections, with clear implications for Greenland Ice Sheet mass loss and future sea level rise rates, while end-of-century winter-spring precipitation in the Mediterranean may be twice what the CMIP ensemble anticipates — supporting elevated flood risks under fossil-fuelled development pathways. The palaeoclimate-weighted projections indicate that winter sea ice retreat, and the extreme mid-latitude weather it may fuel, are probably more severe than current models suggest — a finding that demands the same attention long given to the Arctic’s more famous summer decline.

Subject of Research: Cryospheric science

Article Title: Palaeoclimate constraints narrow uncertainty in future Arctic winter sea ice projections

Article References: Fu, C. Y., Osman, M. B., & Aquino-López, M. A. (2026). Palaeoclimate constraints narrow uncertainty in future Arctic winter sea ice projections. Nature Geoscience. https://doi.org/10.1038/s41561-026-02123-w

Image Credits: AI Generated

DOI: 10.1038/s41561-026-02123-w

Keywords: Arctic sea ice, palaeoclimate, climate models, mid-Holocene, Bayesian statistics, CMIP6, Nature Geoscience, winter sea ice, sea ice projections, cryosphere, constraints, narrow

The problem the index addresses is both scientific and political. In arid zones, water policy often proceeds without any regulatory framework for identifying which frozen landscapes matter most, and existing watershed classification schemes rely on morphometric analysis, general hydrology, or ecological indicators that rarely account for the cryosphere. Yet in the Copiapó basin, the cryosphere is everything. Long-term records from Chile’s water authority show that the river’s discharge, typically between roughly 1.49 and 1.82 cubic meters per second at the La Puerta gauging station, is strongly modulated by meltwater from its high-altitude tributaries. Centuries of tree-ring reconstructions reveal declining rainfall in northern Chile, and megadroughts have intensified since 2010, making the structural dependence of the river on snow and ice melt more consequential than ever.

Building the index required an enormous mapping effort. The researchers compiled a GIS-ready inventory of 1,862 cryospheric landforms across the watershed: 62 debris-free glaciers, 16 debris-covered glaciers, 173 rock glaciers, 848 gelifluction taluses, and 763 protalus lobes. The numerical dominance of periglacial features is itself a revelation, since gelifluction taluses and protalus lobes together account for 87 percent of inventoried landforms while glaciers make up less than 4 percent. Applying objective criteria, at least one active cryospheric landform plus measurable seasonal snow permanence, the team classified 12 of the basin’s 27 sub-watersheds as cryospheric. These 12 units, all situated at high elevation with a mean of 3,617 meters above sea level, contain 97 percent of the mapped landforms and overlap almost perfectly with high-probability permafrost zones predicted by an independent global permafrost model.

Quantifying the ice hidden in these landforms demanded a unified physics-based framework rather than simple area-volume scaling. For debris-free glaciers, the team integrated the global consensus ice-thickness product of Farinotti and colleagues over Randolph Glacier Inventory outlines. For rock glaciers, protalus lobes, and gelifluction taluses, they applied a viscoplastic creep inversion that links surface velocity, basal shear stress, and ice rheology to estimate thickness, constraining the model with surface velocities derived from Sentinel-1 radar interferometry. The DInSAR processing, with interferometric pairs selected for short temporal and perpendicular baselines, achieved velocity uncertainties of roughly 3 to 4 millimeters per year, sufficient to detect the slow creep of frozen debris slopes. Ice-content fractions were treated as bounded priors, 0.5 plus or minus 0.10 for rock glaciers and 0.37 plus or minus 0.12 for the other periglacial landforms, and all errors were propagated quadratically through to ice mass.

Fieldwork in the Montosa and Del Potro sub-watersheds provided crucial ground truth. Electrical resistivity tomography surveys, using 600-meter profiles with electrodes spaced 5 meters apart, revealed high-resistivity zones exceeding 4 to 10 kilo-ohm meters, the signature of massive or ice-saturated permafrost. The geophysics indicated an ice-rich core roughly 15 meters thick in a gelifluction talus and 20 to 30 meters thick, with local maxima above 35 meters, in an active rock glacier. These thicknesses are an order of magnitude greater than values assumed for analogous landforms in parts of the existing literature, and they suggest that periglacial ice reserves in the arid Andes are being systematically undercounted. Notably, Chile’s national periglacial inventory currently systematizes only rock glaciers, omitting the gelifluction taluses and protalus lobes that dominate the Atacama cryosphere by extent.

The snow component of the index drew on 23 years of daily MODIS fractional snow cover imagery from 2000 to 2022. Rather than relying on isolated satellite scenes, the researchers computed snow permanence, the fraction of days each 500-meter pixel remained snow-covered, and integrated it into mean snow-covered area per sub-watershed. Critically, they validated this variable hydrologically: September snow extent correlated robustly with October-to-March discharge at the Pastillo gauging station, with a coefficient of determination of 0.58 across 22 years. The regression intercept, 0.55 cubic meters per second, represents melt-season flow sustained independently of the current year’s snowpack, presumably by longer-residence ice storage. The ratio of this intercept to each year’s observed discharge yielded the empirical weighting factor that balances ice versus snow in the index, averaging 0.43 but shifting to 0.27 in wet years and 0.62 in dry years.

The resulting rankings reveal a watershed split into two hydrological personalities. The Del Potro River sub-watershed holds the largest ice mass at 0.901 gigatons, followed closely by Montosa with 0.873 gigatons and Manflas with 0.218 gigatons. Under average and dry conditions, Montosa tops the NCI at 0.72 plus or minus 0.12, followed by Estero Come Caballos at 0.62 plus or minus 0.23, Del Potro at 0.61 plus or minus 0.18, Manflas at 0.59 plus or minus 0.17, and Figueroa at 0.50 plus or minus 0.21. In dry-year simulations, Montosa and Del Potro occupy the top two positions in practically every Monte Carlo run, and together the top-ranked basins concentrate roughly 78 percent of the inventoried ice mass. In wet years, however, the picture inverts: northern, snow-rich basins such as Estero Come Caballos surge up the ranking on the strength of their persistent snow cover, despite hosting relatively few cryospheric landforms.

This dual behavior is the central management insight of the index. Southern, ice-rich sub-watersheds function as drought insurance, sustaining baseflow when snowfall fails, while northern, snow-driven basins generate the bulk of runoff in high-precipitation years and may require infrastructure to capture the increased meltwater that warming will deliver. The authors also traced the geography behind the pattern to the South American Arid Diagonal, the belt of extreme aridity that fragments the cryosphere in the north into small, isolated bodies. In the central sector between Cadillal and Caserones hills, orographic blocking by the higher Argentine Andes starves the Chilean slopes of snowfall, producing a conspicuous gap with no cryospheric landforms at all. The very slow DInSAR velocities of Copiapó rock glaciers, millimeters per year, one to two orders of magnitude slower than counterparts near 30 degrees south, are physically consistent with colder, drier permafrost and the cryogenic, rather than glacigenic, origin of these landforms.

The implications extend well beyond one Chilean basin. Projections for the Southern Andes indicate that glacier meltwater contributions may decline markedly during future megadroughts, which raises the urgency of protecting all cryospheric water reserves, including the periglacial stores that global inventories routinely ignore. The NCI framework, with its explicit Monte Carlo propagation of uncertainty and hydrologically grounded weighting, offers a reproducible template for any arid mountain region where frozen landscapes sustain downstream populations. The authors caution that the index is calibrated for the Copiapó watershed and requires complete resource data to transfer elsewhere, and that under the ongoing megadrought the ice term should gain weight in future applications. Still, the message is unambiguous: five sub-watersheds, Montosa, Del Potro, Manflas, Estero Come Caballos, and Figueroa, host 83.8 percent of the basin’s cryospheric ice volume and constitute the foundation of continuous water production for the Copiapó River. Identifying them, the researchers argue, provides the technical basis for integrating glacial and periglacial reserves into Chile’s water management strategy, and potentially into a legal instrument to protect the country’s cryosphere before the frozen reserves that keep a desert river flowing quietly disappear.

Subject of Research: A normalized cryospheric index for ranking sub-watersheds by ice storage and snow permanence in arid Andean environments

Article Title: Introducing a new normalized cryospheric index (NCI) to categorize sub-watersheds on arid environments

Article References: Ulloa, C., García, A., & Beniest, A. (2026). Introducing a new normalized cryospheric index (NCI) to categorize sub-watersheds on arid environments. The Cryosphere, 20(10), 5829-5866. https://doi.org/10.5194/tc-20-5829-2026

Image Credits: AI Generated

DOI: 10.5194/tc-20-5829-2026

Keywords: cryosphere, Normalized Cryospheric Index, Copiapó watershed, Atacama Desert, rock glaciers, permafrost, snow permanence, water resources, DInSAR, Andes, drought, watershed management

Permafrost is conventionally defined as soil or rock that has remained at or below 0 degrees Celsius for at least two consecutive years, and it blankets more than 20 percent of the northern hemisphere’s land area. Yet the definition conceals a subtlety that has fascinated cryospheric scientists for decades: permafrost is not necessarily completely frozen. Salts dissolved in pore water depress the freezing point, and capillary forces plus adsorption onto fine-grained mineral surfaces allow thin films of liquid water to persist at temperatures well below zero. Those unfrozen films can become highways for ions and fluids, permitting redistribution and migration of water within the cryotic pore space. Near-surface saline permafrost, in which soluble salts are abundant enough to significantly lower the freezing point, may cover as much as 35 percent of the continuous permafrost zone of the circumpolar north.

The research team drilled three cores at the ADE site in Adventdalen, a U-shaped fjord valley carved by glaciers and later flooded by seawater during the Early Holocene. Two boreholes spaced half a meter apart were drilled in March 2017 and treated as a single core, ADE-17, while two more, ADE-1 and ADE-2, were drilled in March 2022 within roughly 30 meters of the first. The valley fill at the site consists of fine aeolian deposits over fluviatile gravel and roughly twelve meters of deltaic sediments laid down when the sea reached this far up the valley, all resting on Lower Cretaceous shales. After the sea retreated following rapid glacial rebound, the exposed marine sediments froze from the top down, while later fluviatile and aeolian layers froze in place as they accumulated, creating a shallow syngenetic cap above the deeper epigenetic permafrost.

What the cores contained stunned the researchers. Ground ice chemistry differed dramatically between boreholes separated by no more than a few tens of meters. In ADE-17 and ADE-1, the shallow syngenetic ice was essentially fresh water, while the deeper epigenetic sections were brackish to saline, with chloride concentrations between 450 and 3230 milligrams per liter. ADE-2 was another story entirely: salinity reached 7960 milligrams of chloride per liter, with one sample exceeding 23,000 milligrams per liter and total dissolved solids above 40,000 milligrams per liter, saltier than seawater. Even the syngenetic section of ADE-2, deposited only two to four thousand years ago, carried total dissolved solids above 6000 milligrams per liter, a clue that something had pushed saline fluid upward into sediments that had never known the sea.

The ionic ratios told a story of two different waters. In ADE-17 and ADE-1, sodium-to-chloride and sulfate-to-chloride ratios far exceeded those of seawater, indicating that porewater had undergone substantial water-rock interaction, dissolving minerals or exchanging ions during long residence in the valley-fill sediments. In ADE-2, by contrast, the ratios closely matched seawater, with sodium-to-chloride values between 0.84 and 1.02, suggesting a marine-derived saline component that had experienced little chemical modification on its journey. Because freezing and salt expulsion would have produced sodium-to-chloride ratios lower than seawater rather than higher, the team ruled out brine formation by freezing, pointing instead to an intrusion of saline fluid that mixed with the original, fresher ground ice.

The decisive evidence came from the decay chains of uranium and thorium. Radium has four naturally occurring isotopes with half-lives spanning an extraordinary range: radium-223 and radium-224 decay within days, radium-228 within 5.75 years, and radium-226 within 1601 years. All are produced by the alpha decay of thorium isotopes, which bind tightly to mineral grains, while radium is comparatively mobile and tends to enter the dissolved phase. Alpha recoil during decay ejects radium from the mineral lattice into pore space. Because short-lived isotopes reach steady-state activities within weeks, the ratio of a long-lived isotope such as radium-226 to a short-lived one such as radium-223 depends almost entirely on how long the ice or fluid has sat in place, growing toward equilibrium with its parent over five to six half-lives, or roughly eight to ten thousand years for radium-226.

In ADE-17 and ADE-1, the long-to-short-lived radium ratios were high, mostly between 20 and 70, close to or above the secular equilibrium value of 21.7 and approaching the parent ratios measured in the bulk sediment. This implies that the ground ice in these cores is ancient, potentially as old as the Early Holocene permafrost that formed roughly 9000 to 9500 years ago after the sediments emerged from the sea. The team attributes the elevated long-lived radium activities to slow diffusion out of the mineral grains themselves, likely through partially liquidized nano-pores, where thin films of unfrozen water along grain boundaries allow ions to escape over years to decades. Simulations using an effective diffusion coefficient of 10 to the minus 13 square meters per second showed that radium could escape silt-sized grains within years by this pathway, enriching the pore space in radium-226 and radium-228 without requiring any contribution from the short-lived isotopes.

ADE-2 defied this pattern. Its long-to-short-lived radium ratios were strikingly low, between 4 and 14, and matched the parent isotope ratios found on grain surfaces in the exchangeable fraction of the sediment rather than in the bulk grains. In other words, the radium clock in this core had been reset recently, before the longer-lived isotopes had any chance to diffuse out from inside the grains or accumulate significantly from recoil. Modeling the buildup of the radium-226 to radium-223 ratio showed that with no adsorption, reaching equilibrium would take several thousand years, but with adsorption coefficients typical of seawater salinity, the timescale collapses to centuries or decades. The team concludes that the saline ground ice in ADE-2 is younger than 100 years, and possibly only years to decades old, the product of a Late Holocene, and likely very recent, intrusion of saline fluids that mixed with the original ground ice in a fingering style.

Where did those fluids come from? A direct connection to Adventfjorden roughly ten kilometers away seems implausible given the low hydraulic gradient and shale-dominated geology, which would produce residence times of thousands of years. The researchers instead point to a local sandstone aquifer beneath the permafrost, most likely the Helvetiafjellet Formation more than 70 meters below the surface, or possibly the deeper fractured Geerdalen Formation. Fractured sandstone typically carries a high water-to-rock ratio and low thorium content, which would explain both the seawater-like chemistry of the intruding fluid and its capacity to reset the radium clock through enhanced desorption. Artesian discharge through pingo structures nearby demonstrates that water does ascend from the sub-permafrost zone in this region, driven by hydraulic heads from surrounding mountains or by basal permafrost aggradation, and the same pressure may have forced saline fluids upward into the shallow permafrost at the ADE site.

The findings carry weight well beyond one Arctic valley. They demonstrate that saline permafrost is heterogeneous, dynamic, and permeable to fluid migration on human timescales, not the inert deep-frozen archive it is often assumed to be. As Arctic air temperatures continue to climb, thawing of saline permafrost may mobilize subsurface brines, threatening infrastructure built on frozen ground and potentially accelerating the release of stored carbon to the atmosphere. The same radioactive clocks that revealed a decades-old intrusion in Adventdalen now offer scientists a tool to detect such hidden flow elsewhere in the circumpolar north, turning the quiet chemistry of decaying atoms into an early-warning system for a warming Arctic.

Subject of Research: Ground ice residence times and intra-permafrost fluid flow in Svalbard permafrost traced with radium and thorium isotopes

Article Title: The history of ground ice formation and intra-permafrost fluid flow as documented by Ra and Th isotopes

Article References: Rotem, D., Weinstein, Y., Harlavan, Y., Torfstein, A., & Christiansen, H. H. (2026). The history of ground ice formation and intra-permafrost fluid flow as documented by Ra and Th isotopes. The Cryosphere, 20(9), 5303-5326. https://doi.org/10.5194/tc-20-5303-2026

Image Credits: AI Generated

DOI: 10.5194/tc-20-5303-2026

Keywords: permafrost, ground ice, radium isotopes, thorium isotopes, Svalbard, Adventdalen, saline permafrost, fluid migration, geochronology, cryosphere, Arctic warming, hydrogeology

The tool is the Snow Microwave Radiative Transfer model, known as SMRT, which has been developed and refined by the snow remote-sensing community since 2015. Originally built to simulate passive microwave emission from snow, SMRT evolved into a general-purpose framework capable of handling a wide range of cryospheric environments and sensors. It already included a module for conventional low-resolution-mode altimetry, the older style of radar altimetry in which the satellite illuminates a broad, pulse-limited footprint. What it lacked was support for the delay-Doppler, or synthetic aperture radar, processing mode that most modern sensors, including CryoSat-2 and Sentinel-3, actually use. The new version, SMRT 1.7, fills that gap with a dedicated SAR altimetry module described in the journal Geoscientific Model Development.

The significance of the delay-Doppler technique lies in how it reshapes the radar echo. Rather than letting the pulse spread across a wide oval footprint, the satellite exploits the Doppler shift of the returning signal to synthetically focus the beam along its flight direction, shrinking the illuminated strip to roughly 300 metres in the along-track direction. This improves the precision of elevation measurements and reduces the blurring effect of sloped or rugged terrain. But the focusing comes at a price: the resulting waveform has a distinctly different shape from the classic ocean-altimetry curve described by Brown’s famous 1977 model, and interpreting it over snow and ice demands purpose-built physics.

Picard and his colleagues, working with partners in Canada, Norway, Italy, and France, took a deliberately modular approach. Their key insight is that the echo can be decomposed into two nearly independent pieces that are later combined. The first piece describes the vertical journey of the wave into the snowpack, capturing the echoes produced by the surface, by internal interfaces between layers, and by volume scattering within the snow itself. The second piece describes the horizontal spreading of the spherical radar wave over the curved Earth and the filtering imposed by the SAR processing. The first component already existed in SMRT; the second required the team to survey the scientific literature and implement not one but eight different delay-Doppler formulations, published between 2004 and 2019 by groups including Wingham and colleagues, Halimi and colleagues, Ray and colleagues, Boy and colleagues, Buchhaupt and colleagues, Dinardo and colleagues, and Landy and colleagues.

Why eight? Because, as the team’s careful review makes clear, there is no single best model. The formulations differ in how they represent surface topography, whether they treat the terrain statistically as a probability distribution of heights or deterministically as a digital elevation model, how they approximate the antenna pattern and the compressed pulse shape, and how they balance analytical elegance against computational speed. Some models, such as those of Dinardo and Buchhaupt, are extremely fast and well suited to physical retracking algorithms that fit observed waveforms in near real time. Others, such as the facet-based model of Landy, can ingest a real digital elevation model and capture the electromagnetic bias caused by skewed slope distributions, but at a far greater computational cost. By bundling all eight into a single plug-and-play framework, SMRT lets researchers compare them directly and pick the right tool for the job.

The technical machinery works in three steps. First, existing SMRT modules compute the backscatter of every layer and interface, drawing on theories such as the Improved Born Approximation for volume scattering and the Geometrical Optics approximation for rough surfaces. A new addition handles coherent, quasi-specular backscatter, which matters when the surface is smooth, as it often is over sea ice, and which depends in a subtle way on the altitude of the sensor. Second, the user-selected delay-Doppler model computes a delay-Doppler map for each interface and for the volume, applying a slant-range correction that migrates echoes arriving at off-nadir angles back to their correct range. Third, the pieces are combined into the final delay-Doppler map, which is summed along the Doppler dimension to yield the waveform a satellite would record.

Verification proceeded on several fronts. Under simplified conditions, all eight models produced waveforms that peaked within the same range gate, with only small differences in the leading and trailing edges attributable to numerical and theoretical approximations. The team’s implementation of the Landy model matched the original Matlab code to within one percent. And in a particularly elegant test, the researchers disabled the slant-range correction to produce pseudo-low-resolution-mode waveforms and compared them against SMRT’s independent LRM module, finding near-perfect agreement with the Wingham18 formulation. These cross-checks matter because a modelling error of even a few percent in the peak position could be misread as a real change in surface elevation.

The real-world test came on the Antarctic ice sheet. The team drove the model with in-situ measurements collected during four field campaigns, including density and specific surface area profiles from snow cores, temperature readings, and photogrammetric estimates of centimetre-scale surface roughness, supplemented by the Reference Elevation Model of Antarctica for large-scale terrain statistics. Comparing the simulations with Sentinel-3 Ku-band waveforms at eighteen sites, they found that the model reproduced the observed latitudinal gradient in waveform amplitude, with stronger echoes over the smooth, wind-scoured plateau near Dome C and weaker ones in the rougher coastal margins. When the poorly constrained radar-scale roughness was optimized at each site, the agreement improved markedly, confirming that surface roughness exerts first-order control on how much of the echo comes from the snow surface versus the snow volume beneath it.

The study is candid about its limits. The current implementation assumes the medium is homogeneous across a footprint spanning many kilometres, neglects multiple scattering between layers, and covers only unfocused SAR processing, leaving fully focused and interferometric modes for future work. Multiple scattering is likely to become significant at the Ka-band frequency of the upcoming Copernicus CRISTAL mission, where the team’s calculations suggest single scattering accounts for less than half the backscatter in the top metre of Antarctic snow. Even so, the payoff is substantial. With a verified, open-source framework that links detailed snow physics to realistic SAR waveforms, scientists can now attack long-standing problems such as the penetration bias that offsets ice-sheet elevation measurements, the validity of assuming that radar echoes from sea ice originate at the snow-ice interface, and the design of retrieval algorithms for snow depth on sea ice and thickness of lake ice. As new missions like CRISTAL and Sentinel-6 come online, the ability to simulate exactly what their radars should see, given what we know about the snow and ice below, may prove decisive for turning polar radar echoes into climate-grade measurements.

Subject of Research: Simulation of SAR radar altimeter echoes from snow and ice surfaces using the SMRT model

Article Title: Simulating SAR altimeter echoes from cryospheric surfaces with the Snow Microwave Radiative Transfer (SMRT) model version 1.7

Article References: Picard, G., Murfitt, J., Zakharova, E., Zeiger, P., Arnaud, L., Aublanc, J., Landy, J. C., Scagliola, M., & Duguay, C. (2026). Simulating SAR altimeter echoes from cryospheric surfaces with the Snow Microwave Radiative Transfer (SMRT) model version 1.7. Geoscientific Model Development, 19(18), 9103-9129. https://doi.org/10.5194/gmd-19-9103-2026

Image Credits: AI Generated

DOI: 10.5194/gmd-19-9103-2026

Keywords: radar altimetry, SMRT model, cryosphere, Antarctica, Sentinel-3, delay-Doppler, sea ice thickness, ice sheet elevation, snow microwave radiative transfer, waveform retracking, Ku-band radar, CRISTAL

The product, called GMM-SWE v1.0, covers the Northern Hemisphere between 35° and 85° north latitude for every month from 2003 to 2022. Rather than relying on a single measurement technique, the framework fuses two very different streams of information: satellite measurements of tiny changes in Earth’s gravity field from the GRACE and GRACE-FO missions, and high-resolution meteorological data from the ERA5-Land reanalysis. The two are combined within a physics-constrained data assimilation scheme built around an adaptive Kalman filter, a mathematical engine that continuously reconciles model predictions with observations by weighting each according to its uncertainty. The result is a spatially continuous, physically consistent record of three fundamental snow properties: snow water equivalent, snow depth, and snow density.

The gravitational component is the conceptual heart of the approach. The GRACE satellites, flying in tandem, detect minuscule variations in the distance between one another as they pass over regions of changing mass on the surface below. When snow accumulates, the ground beneath it becomes fractionally heavier, and the satellites register the shift. The signal, however, is not clean. Changes in total terrestrial water storage also include groundwater fluctuations, soil moisture, surface water, and other hydrological components. To isolate the snow contribution, the researchers explicitly subtracted these non-snow terms, derived from an ensemble of land surface and hydrological models. This mass-balance approach means the snow estimates are anchored to actual physical mass rather than to statistical relationships that may break down outside the conditions in which they were calibrated.

Gravity alone, however, cannot tell the whole story. The GRACE measurements are coarse in space and time, and they say nothing directly about how deep the snowpack is or how dense it has become. GMM-SWE therefore incorporates high-resolution meteorological forcing and global topographic information to capture elevation-dependent melt dynamics and temperature thresholds, a refinement that matters enormously in complex mountain terrain where snow behavior can change dramatically over a few hundred meters of altitude. Once the assimilated snow mass has been estimated, it is further refined with a Rauch–Tung–Striebel smoother, a technique that reduces temporal noise and ensures that the seasonal evolution of the snowpack remains physically coherent rather than jittering from month to month.

Perhaps the most innovative element of the dataset is its dynamic decoupling of snow mass and snow volume. Traditional large-scale products often assume a static relationship between depth and density, an assumption that ignores the fact that a snowpack is a living, evolving material. Fresh powder may sit at densities below 100 kilograms per cubic meter, while a spring snowpack that has survived repeated melt–freeze cycles can exceed 500. GMM-SWE retrieves snow depth and density by coupling the smoothed snow water equivalent estimates with a physics-based SNOW-17 compaction scheme resolved at sub-daily time steps. This allows the dataset to explicitly represent densification and metamorphic processes within the snowpack, capturing how grains round, settle, and transform under their own weight and through temperature fluctuations, without leaning on fixed density assumptions that vary wildly across climates and seasons.

The validation effort behind the dataset is substantial. The researchers compared their estimates against an extensive network of 13,135 in-situ ground stations spanning a wide range of climatic and topographic regimes, as well as against established products including ERA5-Land, GLDAS, WGHM, and GlobSnow. The results show consistently strong performance, with particular skill in the high-latitude continental interiors of Siberia and North America, where seasonal snow dominates the annual water cycle. Crucially, the framework substantially improves the representation of snow water equivalent in complex mountain regions, such as high mountain Asia, where passive microwave satellite retrievals, the backbone of many conventional snow products, are severely limited by deep snowpacks, wet snow conditions, and rugged topography that scramble the microwave signal.

Sensitivity analyses add an important layer of credibility. The researchers found that the large-scale patterns of snow water equivalent and their year-to-year variability are governed primarily by the physical mass-balance constraints at the core of the method, rather than by finely tuned parameter choices. In other words, the dataset’s performance does not hinge on arbitrary calibration decisions that might not generalize; it emerges from the underlying physics of mass conservation. This robustness is exactly what is needed for a product intended to serve hydrological applications and climate-scale assessments of cryospheric change, where users must be able to trust that the numbers reflect reality rather than modeler preference.

The practical implications reach far beyond the scientific community. Snow water equivalent is the single most important variable for seasonal water forecasting in snow-dominated basins: it determines how much meltwater will arrive in rivers and reservoirs months later. Water managers across the western United States, Central Asia, and the Himalayan headwaters make decisions worth billions of dollars based on snow estimates, and errors in those estimates propagate directly into flood risk, drought planning, and hydropower generation. A unified, physically informed dataset that spans two decades and the entire hemisphere gives these communities a consistent baseline, and gives climate scientists the long, homogeneous record they need to detect trends in seasonal snow storage and to evaluate how well climate models capture them.

The dataset also arrives at a moment of rapid change. Warming temperatures are altering the timing of snowfall, shortening the snow season, and shifting precipitation from snow to rain in many mid-latitude regions, with pronounced regional differences in how seasonal snow storage is evolving. Disentangling those regional trajectories requires observations that are consistent across borders and biomes, something no network of ground stations alone can provide. By combining the global reach of satellite gravimetry with the local realism of meteorological forcing and snow physics, GMM-SWE offers a template for how the next generation of cryospheric monitoring products might be built, and the authors have released both the dataset and the MATLAB production code openly, inviting scrutiny, reuse, and improvement by the wider community.

As with any ambitious synthesis, the work is a preprint currently under peer review, and its assumptions, particularly the model-based subtraction of non-snow water storage components, will face scrutiny from specialists. But the underlying idea is undeniably elegant: the same planet-watching satellites originally designed to track melting ice sheets and draining aquifers can, with the right physics bolted on, be repurposed to weigh the winter snowpack of an entire hemisphere. In doing so, they transform one of hydrology’s most stubborn uncertainties into a monthly, measurable quantity, one that may prove indispensable as the world’s snowlines climb and the freshwater stored in snow becomes an ever more precious and contested resource.

Subject of Research: A GRACE and meteorology data assimilation dataset of Northern Hemisphere snow water equivalent, snow depth, and snow density for 2003–2022

Article Title: GMM-SWE v1.0: A GRACE–Meteorology-Constrained Dataset of Snow Water Equivalent, Snow Depth, and Snow Density over the Northern Hemisphere, 2003–2022

Article References: Mohasseb, H., & Yi, S. (2026). GMM-SWE v1.0: A GRACE–Meteorology-Constrained Dataset of Snow Water Equivalent, Snow Depth, and Snow Density over the Northern Hemisphere, 2003–2022. https://doi.org/10.5194/essd-2026-675

Image Credits: AI Generated

DOI: 10.5194/essd-2026-675

Keywords: snow water equivalent, GRACE, GRACE-FO, data assimilation, snow density, snow depth, ERA5-Land, Northern Hemisphere, cryosphere, hydrology, satellite gravimetry, climate change

Snow algae are found on every continent, predominantly in polar and alpine settings. They are keystone primary producers, acting as terrestrial carbon sinks, facilitating nutrient cycling, shaping microbial communities and providing nutrition for higher organisms. But they also have a darker climatic significance: by darkening snow surfaces, they reduce albedo, the fraction of incoming solar radiation that snow reflects, and thereby accelerate melting. Snow algae have been shown to reduce snow albedo by up to 44 percent. In one Antarctic Peninsula study, albedo fell from 0.85 at algae-free control sites to 0.44 over green-dominated blooms and 0.65 over red ones, corresponding to roughly 40 percent and 20 percent reductions respectively. Despite this ecological weight, the mechanisms governing how snow algae migrate within snow remain poorly understood, which is precisely the gap the review sets out to fill.

The authors frame snow algae as biologically active particles within the framework of active matter physics, a class of non-equilibrium soft matter whose motion cannot be explained by equilibrium physics alone. Unlike passive particles, microalgae harvest energy from their environment and convert it into directed motion. Their propulsion comes from flagellar beating, a process in which adenosine triphosphate is converted to adenosine diphosphate, releasing energy that powers microtubule sliding within the flagella. The motile species examined in the review are typically biflagellated, and the symmetry and phase relationships of their two flagella produce distinctive swimming gaits, including the helical swimming pattern well characterised in the model freshwater alga Chlamydomonas reinhardtii, which the authors use as a proxy where snow algae data are lacking.

Scale matters enormously in this world. The review walks through the dimensionless Reynolds number, the ratio of inertial to viscous forces, to illustrate just how alien the algal fluid environment is. For Chlamydomonas reinhardtii, with a swimming speed of roughly 130 micrometres per second and a cell length of about 10 micrometres, the Reynolds number is around 1.3 times ten to the minus three. For a human swimmer it is about 1.7 times ten to the sixth, roughly a billion times larger. At such low Reynolds numbers, inertia is negligible and reciprocal strokes produce no net motion, which is why microalgae have evolved helical flagella and whip-like techniques to navigate. Any understanding of how algae move through snow must therefore begin with the physics of viscous-dominated swimming.

The snowpack itself is a dynamic, porous and evolving medium. Snow crystals, typically between half a millimetre and three millimetres across, are up to three orders of magnitude larger than the algae, which measure roughly 5 to 40 micrometres. Within the pack, crystals continuously metamorphose under temperature and pressure gradients, forming rounded grains in the upper layers and large depth hoar crystals at the base. In melting snow, where algae flourish, water films bind crystals into clusters. At the macroscopic scale, the algae follow a seasonal cycle: they overwinter as dormant cysts, germinate into green motile flagellated cells with the spring melt, migrate upward toward the surface in response to meltwater, light and released nutrients, and then transform back into non-motile, pigmented cysts. Field experiments on an Alaskan ice field suggest that at the peak of the growing season, actively resurfacing cells account for about 65 percent of surface algal abundance, with passive dispersal by wind, water and birds making up the remaining 35 percent.

At the microscopic scale, the story becomes subtler. Snow algae were long thought to inhabit the quasi-liquid layer, a thin film of liquid water that coats snow crystals even at sub-zero temperatures. But that layer is astonishingly thin, ranging from a few molecular layers, roughly 0.37 nanometres each, up to about 10 nanometres as temperature rises. Motile snow algae cells are six to twenty micrometres across, several orders of magnitude too large to swim within it. In non-melting snow, moreover, about 80 percent of the limited liquid water is held in menisci at crystal contact points. This size disparity suggests that active swimming within the quasi-liquid layer is unlikely, and helps explain why algal migration and bloom development peak during snowmelt, when larger interconnected water channels form. X-ray tomography of the red snow alga Sanguina nivaloides has shown that its dormant cysts occur only in the liquid water fraction of the snowpack, never inside ice grain cores.

Fluid flow within the pack adds another layer of physics. Preferential flow paths carry meltwater at speeds between 12 and 30 millimetres per second, far exceeding the reported mean swimming speed of the snow alga Chlamydomonas nivalis, a mere 0.061 millimetres per second. Whether algae can swim against such flows depends on the viscous diffusion timescale, which scales with the square of channel width divided by kinematic viscosity. In narrow channels or saturated capillary zones, where flow is effectively stagnant, algae can move independently of the bulk water. In wider channels they are likely transported passively by advection. The review also highlights interfacial pre-melting and thermal regelation as possible migration aids: when an algal cell is embedded in ice near its melting point, surface forces induce a thin melted film, and a temperature gradient can drive the cell toward warmer regions through cycles of melting and refreezing. Some Chlamydomonas species even produce ice-binding proteins that hinder ice crystal growth, and exopolymeric substances and antifreeze glycoproteins may further enhance survival and motility in icy conditions.

Much of the review is devoted to tactic behaviour, the biased swimming of cells toward or away from stimuli. Phototaxis is the best documented. Motile snow algae swim toward light at low to moderate intensities but switch to negative phototaxis above a critical threshold; in a Japanese alpine snowpack, motile cells ascended nearly to the surface during low-light hours and descended 10 to 20 centimetres when solar radiation peaked at up to 755 watts per square metre. Cell density at the surface was negatively correlated with solar radiation and air temperature, and because solute gradients showed no day-night variation, light rather than nutrients drove the migration. Intriguingly, some snow algal species lack eyespots yet still exhibit phototaxis, while close relatives with eyespots do not, raising questions about how light sensing works in snow, where reflections off countless crystals could produce confusing signals. Ultraviolet-B radiation, meanwhile, appears to act mainly as a physiological stressor that suppresses motility rather than as a directional cue.

Chemotaxis, gravitaxis, gyrotaxis and thermotaxis complete the picture, though each is far less understood in snow algae specifically. Snowpacks are oligotrophic, and laboratory experiments suggest phosphorus, rather than nitrogen, is the limiting nutrient for the snow alga Chloromonas typhlos, making phosphorus gradients a plausible but untested driver of chemotactic movement. Gravitaxis, the orientation of cells relative to gravity, remains debated between passive mechanisms such as bottom-heavy mass distribution and differential sedimentation, and active physiological sensing; mutant studies in Chlamydomonas point to active, membrane-excitability-dependent signalling. Gyrotaxis, the combination of gravitational and viscous torques in flowing fluid, has been shown experimentally to generate bioconvective patterns in snow algal suspensions, and may focus swimming algae into downward meltwater flows, accelerating their transport to the base of the pack. Thermotaxis has never been demonstrated in snow algae, though psychrophilic species show optimal swimming speeds below 10 degrees Celsius, in sharp contrast to mesophilic species that peak above 20 degrees.

The review closes with a call to arms. The seasonal cyst-to-flagellate cycle is broadly accepted, but the quantitative biophysics of how algae interact with quasi-liquid layers, navigate meltwater channels and respond to evolving snowpack structure remains undocumented. Filling these gaps, the authors argue, would improve predictions of snow algal blooms and their climatic consequences, refine hydrological and cryospheric models, advance our understanding of microswimmer behaviour in active matter physics, and even inform emerging biotechnologies, from low-temperature algal cultivation for astaxanthin production to medical microrobotics inspired by algal motility. The pink patches staining the world’s snowfields, it turns out, are not just a curiosity but a living laboratory for physics, ecology and climate science all at once.

Subject of Research: Biophysical mechanisms of snow algae motility and migration within snowpacks

Article Title: Reviews and syntheses: Snow algae on the move – biased motility and snowpack interaction from a biophysics perspective

Article References: de Vries, C. S., Sandells, M. J., Davey, M. P., Caldwell, G. S., & Croze, O. A. (2026). Reviews and syntheses: Snow algae on the move – biased motility and snowpack interaction from a biophysics perspective. Biogeosciences, 23(19), 6835-6855. https://doi.org/10.5194/bg-23-6835-2026

Image Credits: AI Generated

DOI: 10.5194/bg-23-6835-2026

Keywords: snow algae, biophysics, active matter, phototaxis, chemotaxis, gravitaxis, gyrotaxis, snowpack, albedo, cryosphere, microswimmers, flagella

The problem LADDIE addresses is one of scale. The most realistic simulations of ice–ocean interaction come from fully coupled, three-dimensional ocean models, but these demand enormous computing resources. Studies cited by the team indicate that ocean models need a resolution of roughly two kilometres or finer to accurately represent melting near the grounding line, the critical zone where ice first begins to float, and even higher resolution to capture the narrow basal channels that channelise melting on many ice shelves. Continental-scale coupled simulations today typically run at resolutions of four to fourteen kilometres below ice shelves, and only regional configurations achieve the two-kilometre threshold. Running such models over the multiple centuries relevant to ice sheet dynamics remains unfeasible, which is why most Antarctic sea-level projections still rely on stand-alone ice sheet models with simplified melt parameterisations.

Those parameterisations, which translate offshore ocean temperatures into fields of melt rates beneath the ice, come in several flavours, including quadratic temperature dependencies, the box model PICO, and one-dimensional plume models. Yet detailed comparisons against three-dimensional ocean models and observations have shown that they struggle to reproduce the horizontal pattern of melting, and earlier work by the same group demonstrated that ice sheet dynamics are highly sensitive to precisely that spatial pattern. LADDIE was conceived as a middle way: a two-dimensional model of the upper ocean mixed layer beneath ice shelves that resolves the influence of topographic steering and Coriolis deflection on meltwater flow, without the cost of simulating the full three-dimensional cavity.

Version 2.0 represents a deep technical overhaul. The original LADDIE was written in interpreted Python on a square grid with finite difference numerics, which limited it to coarse pan-Antarctic resolutions. The new version has been translated into compiled Fortran and rebuilt on an unstructured mesh of triangles and Voronoi cells using finite volume numerics, drawing on the Utrecht Polar SYstem toolbox developed for the UFEMISM ice sheet model. Ocean velocities are solved on the triangles, while melt rates and other variables are resolved on the cells surrounding the triangle vertices, an arrangement equivalent to the Arakawa-B discretisation used in the three-dimensional ocean model FESOM. This design gives the velocity field an effective resolution roughly twice as fine as that of the melt field, sharpening the simulated melt plumes.

The time stepping scheme has also been replaced. Where version 1.0 used a modified LeapFrog scheme with a Robert–Asselin filter, the new version employs a third-order Forward–Backward Runge–Kutta scheme, a generalised form of the scheme used in NEMO 5. With carefully chosen forward–backward weights, this scheme maintains numerical stability with time steps approximately twice as large as standard Runge–Kutta approaches. Combined with parallelisation across up to 128 processor cores, these changes cut computation time by an order of magnitude compared with the previous version. In a benchmark simulation of the Amundsen Sea region at roughly one-kilometre resolution, the team found that doubling the core count reduced computing time by about a third up to 64 cores, with 16 cores arguably the optimal trade-off for that particular domain.

The evaluation against three-dimensional ocean models is striking. In the idealised ISOMIP+ Ocean1 experiment, which prescribes an 80-kilometre-wide ice shelf under warm Amundsen Sea-like forcing, LADDIE 2.0 was compared with an ensemble of twelve ocean models. Its melt pattern closely resembled the multi-model mean, with a mean difference of just 0.3 metres per year and a root mean squared difference of 5.5 metres per year, placing it among the models closest to the ensemble average. The model also reproduced the melt patterns of the four models that cluster tightly around the mean, namely MITgcm-BAS, MITgcm-JPL, MPAS-Ocean and NEMO-CNRS, suggesting it captures the essential physics that establish those patterns.

At continental scale, the team configured LADDIE at approximately two-kilometre resolution, a mesh of more than 450,000 vertices, and compared it with the RISE multi-model ensemble of nine pan-Antarctic ocean models. Bulk metrics such as melt sensitivities, which describe how strongly melting responds to ocean warming, fell close to the ensemble median and mean. A comparison against four satellite-derived estimates of sub-shelf melting showed good overall agreement in integrated melt rates per ice shelf, notably without any regional tuning. The largest discrepancies, along the Bellingshausen Sea and at the Borchgrevink ice shelf, appear to stem partly from a warm bias in the ISMIP6 ocean forcing dataset, and preliminary tests with updated forcing significantly reduce the mismatch. The model does underestimate refreezing in cold cavities such as the Filchner–Ronne and Ross ice shelves, likely because it underestimates barotropic currents and tides there.

Perhaps the most visually compelling result comes from the Pine Island ice shelf, one of the fastest-melting glaciers in West Antarctica. Run at a nominal resolution of 120 metres, with a mesh of over 430,000 vertices, LADDIE reproduced the fine-grained network of basal channels observed in high-resolution satellite studies, including the well-known Y-shaped channel and the branching of wide along-flow channels into narrower ones deflected by Coriolis forces and topographic steering. This capability matters because channelised melting is underrepresented in large-scale remote sensing products, and because satellite altimetry cannot measure melt rates close to the grounding line, where bridging stresses break the assumption of hydrostatic balance and where melting may matter most for ice sheet stability.

The consequences for sea-level projections could be substantial. In an idealised coupled simulation integrating LADDIE with the UFEMISM ice sheet model on the same adaptive mesh, the team compared the dynamic melt model against a conventional quadratic melt parameterisation under identical high-melt ocean forcing. Although initial melt rates were tuned to be nearly identical, the coupled simulation produced a threefold increase in grounding line retreat and volume above floatation loss. The reason lies in feedbacks: with the parameterisation, ice shelf thinning lifts the ice out of the warm water layer and suppresses melting, a negative feedback, whereas LADDIE’s simulated meltwater flow sustains strong melting along steep basal slopes and shear margins, a partial positive feedback that amplifies the loss of buttressing.

The authors are careful to stress that two-dimensional melt models cannot replace three-dimensional ocean models, and that LADDIE’s rapid equilibration and simplified physics carry known limitations, including the absence of seasonal intrusions of warm water into shallow cavity regions. But by making a physically detailed, open-source melt model fast enough for centennial-scale, continent-wide simulations, and by fully integrating it with an ice sheet model on a shared mesh, the team hopes to push future intercomparison efforts such as ISMIP7 toward more realistic representations of ocean-driven melting. As the ocean continues to gnaw at the underside of Antarctica’s floating ice, tools like LADDIE may prove essential for turning scattered observations into credible forecasts of how much the sea will rise.

Subject of Research: A two-dimensional ocean mixed-layer model for simulating ocean-driven melting beneath Antarctic ice shelves and its coupling to ice sheet dynamics

Article Title: The one-Layer Antarctic model for Dynamical Downscaling of Ice–ocean Exchanges (LADDIE) version 2.0

Article References: Lambert, E., Jesse, F., & Berends, C. J. (2026). The one-Layer Antarctic model for Dynamical Downscaling of Ice–ocean Exchanges (LADDIE) version 2.0. Geoscientific Model Development, 19(19), 9441-9461. https://doi.org/10.5194/gmd-19-9441-2026

Image Credits: AI Generated

DOI: 10.5194/gmd-19-9441-2026

Keywords: Antarctica, ice shelves, sub-shelf melting, sea-level rise, ocean modelling, LADDIE, ice sheet dynamics, grounding line, basal channels, UFEMISM, model intercomparison, cryosphere