A new method of sea ice model evaluation is demonstrated. Data from the network of Arctic ice mass balance buoys (IMBs) are used to estimate distributions of vertical energy fluxes over sea ice in two densely sampled regions – the North Pole and Beaufort Sea. The resulting dataset captures seasonal variability in sea ice energy fluxes well, and it captures spatial variability to a lesser extent. The dataset is used to evaluate a coupled climate model, HadGEM2-ES (Hadley Centre Global Environment Model, version 2, Earth System), in the two regions. The evaluation shows HadGEM2-ES to simulate too much top melting in summer and too much basal conduction in winter. These results are consistent with a previous study of sea ice state and surface ...