EGUsphere EGUsphere EGUsphere EGUsphere Copernicus Publications Göttingen, Germany 10.5194/egusphere-2026-2341 Evaluation of MARv3.14 over the Greenland Ice Sheet Timmermans Guillaume https://orcid.org/0009-0007-0725-1725 1 Noël Brice https://orcid.org/0000-0002-7159-5369 1 Kittel Christoph https://orcid.org/0000-0001-6586-9784 1 2 Dethinne Thomas https://orcid.org/0000-0003-3114-5248 1 3 Ghilain Nicolas 1 4 Fettweis Xavier https://orcid.org/0000-0002-4140-3813 1 Laboratory of Climatology, Department of Geography, SPHERES research unit, University of Liège, Liège, Belgium Physical Geography research group, Department of Geography, Vrije Universiteit Brussel, Brussels, Belgium Geospatial Data Science and City Information Modeling, Department of Geography, SPHERES research unit, University of Liège, Liège, Belgium Royal Meteorological Institute of Belgium, Uccle, Belgium 02 07 2026 2026 1 30 Copyright: © 2026 Guillaume Timmermans et al. 2026 This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/ This article is available from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2341/ The full text article is available as a PDF file from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2341/egusphere-2026-2341.pdf

Accurately estimating the surface mass balance (SMB) of the Greenland Ice Sheet (GrIS) is essential to quantify its contribution to sea-level rise. The polar regional atmospheric climate model MAR is widely used to simulate GrIS SMB and to force ice sheet dynamics models, highlighting the need for a thorough evaluation. Here, we evaluate the latest MAR version (MARv3.14) over Greenland at 5 km spatial resolution and examine the impact of coarser resolutions (10–30 km) on the simulated SMB and its components. MAR outputs are compared to a range of independent observations, including in situ SMB measurements, automatic weather station (AWS) records of near-surface meteorological variables, satellite-derived melt extent, and albedo products. At 5 km, MAR reproduces the observed SMB with a root-mean-square error (RMSE) of 0.51 m and a correlation of 0.93. For near-surface meteorological variables and surface energy budget fluxes, the model RMSE is smaller than the corresponding observed natural variability (i.e., standard deviation), indicating non-significant model error. Prescribing bare-ice albedo improves the model performance in the ablation zone, while biases remain in the accumulation zone, suggesting that further improvements are required in the snow albedo scheme. In addition, simulated melt timing is consistent with satellite-based melt extent products. Sensitivity experiments reveal that discrepancies between simulations at different spatial resolutions are mostly limited to the ice sheet margins where strong SMB and topographic gradients occur, notably in the southeast of Greenland where precipitation peaks. Differences in integrated SMB generally remain small and mostly non-significant, while individual components, i.e., precipitation and runoff, exhibit larger resolution-dependent variations. We find that a resolution of at least 10 km is required to accurately capture the GrIS climate and SMB. Reducing computation time about 10-fold relative to 5 km simulation, a 10 km grid makes a good compromise for long-term climate projections.