Physics

HydroFirn: A numerical model for large-scale multidimensional firn hydrology

How the science connects

Computational mode…GlaciologyHydrology

AI Insight

Researchers have developed HydroFirn, a new computational model that simulates meltwater movement through firn (compacted snow) in multiple dimensions on the Greenland Ice Sheet. Unlike existing models that only calculate water flow in one dimension, this model captures the complex lateral and vertical dynamics of meltwater percolation and ice layer formation. The model demonstrates that lateral variations in firn properties significantly affect how deep meltwater penetrates and where impermeable ice layers form, which is critical for understanding surface mass balance.


Accurate modeling of firn hydrology is essential for converting satellite measurements of ice sheet elevation changes into mass loss estimates, which directly impacts sea-level rise projections. Better understanding of multidimensional meltwater processes will improve predictions of freshwater discharge into the ocean under climate warming scenarios.


⚠️ Preprint – Noch nicht peer-reviewed

Dieser Artikel wurde noch nicht von unabhängigen Experten begutachtet. Die Ergebnisse sind vorläufig und sollten mit Vorsicht interpretiert werden.

Abstract: Observations show the multidimensional dynamics of meltwater and distribution of ice layers in the firn on the Greenland Ice Sheet. However, state-of-the-art large-scale models for firn hydrology are essentially one-dimensional, limiting their ability to explain observed datasets and failing to reduce uncertainty in surface mass balance and sea-level rise estimates. Here, we present a large-scale, multidimensional, multiphase, and thermomechanical model to simulate firn hydrology. The model is highly efficient due to a novel algorithm in which an extra equation for pressure is solved only in saturated regions. Furthermore, the model can apply spatially heterogeneous boundary conditions to the unsaturated-saturated domain and allows for the dynamic formation of fully impermeable ice layers. The numerical results show excellent comparisons against analytic solutions to one- and two-dimensional problems that involve coupled unsaturated-saturated flows, thermodynamics, and phase change. We further apply the model to investigate field data from southwest Greenland and find that lateral heterogeneities strongly influence the depth of melt percolation and ice layer formation. Improved understanding of these local, multidimensional processes will provide physics-based constraints on firn densification, reduce uncertainty in converting altimetric elevation change to mass change, and improve estimates of freshwater fluxes to the ocean under a warming climate.

Source: HydroFirn: A numerical model for large-scale multidimensional firn hydrology