To improve surface mass balance (SMB) estimates for the Greenland Ice Sheet (GrIS), we developed a 5ĝ€̄km resolution regional climate model combining the Japan Meteorological Agency Non-Hydrostatic atmospheric Model and the Snow Metamorphism and Albedo Process model (NHM-SMAP) with an output interval of 1ĝ€̄h, forced by the Japanese 55-year reanalysis (JRA-55). We used in situ data to evaluate NHM-SMAP in the GrIS during the 2011-2014 mass balance years. We investigated two options for the lower boundary conditions of the atmosphere: an offline configuration using snow, firn, and ice albedo, surface temperature data from JRA-55, and an online configuration using values from SMAP. The online configuration improved model performance in simulating 2ĝ€̄m air temperature, suggesting that the surface analysis provided by JRA-55 is inadequate for the GrIS and that SMAP results can better simulate physical conditions of snow/firn/ice. It also reproduced the measured features of the GrIS climate, diurnal variations, and even a strong mesoscale wind event. In particular, it successfully reproduced the temporal evolution of the GrIS surface melt area extent as well as the record melt event around 12 July 2012, at which time the simulated melt area extent reached 92.4ĝ€̄%. Sensitivity tests showed that the choice of calculation schemes for vertical water movement in snow and firn has an effect as great as 200ĝ€̄Gtĝ€̄year-1 in the GrIS-wide accumulated SMB estimates
a scheme based on the Richards equation provided the best performance.