Greenland ice sheet mass loss continues to accelerate as global temperatures increase. The surface albedo of the ice sheet determines the amount of absorbed solar energy, which is a key factor in driving surface snow and ice melting. Satellite-retrieved snow albedo allows us to compare and optimise modelled albedo over the entirety of the ice sheet. We optimise the parameters of the albedo scheme in the ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystems) land surface model for 3 random years taken over the 2000–2017 period and validate over the remaining years. In particular, we want to improve the albedo at the edges of the ice sheet, since they correspond to ablation areas and show the greatest variations in runoff and surface mass balance. By giving a larger weight to points at the ice sheet's edge, we improve the model–data fit by reducing the root-mean-square deviation by over 25 % for the whole ice sheet for the summer months. This improvement is consistent for all years, even those not used in the calibration step. We also show the optimisation successfully improves the model–data fit at 87.5 % of in situ sites from the PROMICE (Programme for Monitoring of the Greenland Ice Sheet) network. We conclude by showing which additional model outputs are impacted by changes to the albedo parameters, encouraging future work using multiple data streams when optimising these parameters.