Correctly representing the snow on sea-ice in coupled numerical weather prediction models has great potential to improve weather forecast and climate monitoring applications, such as climate reanalyses, which are usually produced using such systems. In this study two different methodologies to account for the effect of the snowpack accumulating over the sea-ice in 5-day global forecasts are compared. This dataset contains the high-resolution (dx=9km) model data for the different simulations, with and without the snow over sea-ice, co-located with observations at the locations of the SHEBA and N-ICE2015 field campaigns. The dataset also contains the daily averaged skin temperature, for the different simulations (interpolated on 0.25x0.25) ...