The firn meltwater Retention Model Intercomparison Project (RetMIP): Evaluation of nine firn models at four weather station sites on the Greenland ice sheet

Abstract. Perennial snow, or firn, covers 80 % of the Greenland ice sheet and has the capacity to retain part of the surface meltwater, buffering the ice sheet’s contribution to sea level. Multi-layer firn models are traditionally used to simulate the firn processes and estimate meltwater retention. We present the output from nine firn models, forced by weather-station-derived mass and energy fluxes at four sites representative of the dry snow, percolation, ice slab and firn aquifer areas. We compare the model outputs and evaluate them against in situ observations. Models that explicitly account for deep meltwater percolation overestimate percolation depth and consequently firn temperature at the percolation and ice slab sites although they accurately simulate the recharge of the firn aquifer. Models using Darcy's law and a bucket scheme compare favourably to observations at the percolation site but only the Darcy models accurately simulate firn temperature and thus meltwater percolation at the ice slab site. We find that Eulerian models, that transfer firn through fixed layers, smooth sharp gradients in firn temperature and density over time. From the model spread, we find that simulated densities (respectively temperature) have an uncertainty envelope of ±60 kg m−3 (resp. ±14 °C) in the dry snow area and up to ±280 kg m−3 (resp. ±15–18 °C) at warmer sites.