Abstract. The polar version of the regional climate model RACMO2, version 2.3p1, is used to study the effect of model resolution on the simulated climate and surface mass balance (SMB) of South Greenland. The model data on resolutions of 60, 20, 6.6 and 2.2 km are intercompared and compared to SMB observations using three different data refinement methods: nearest neighbour, bilinear interpolation and a statistical downscaling method utilising the local dependency of fields on elevation. Furthermore, it is estimated how the errors induced by model resolution compare to errors induced by the model physics and initialisation. The results affirm earlier studies that SMB components which are tightly linked to elevation, like runoff, can be refined successfully, as soon as the ablation zone is reasonably well resembled in the source dataset. Precipitation fields are also highly elevation dependent, but precipitation has no systematic correlation with elevation, which inhibits statistical downscaling to work well. If refined component-wise, 20 km resolution model simulations can reproduce the SMB ablation observations almost as good as the finer resolution model simulations. Nonetheless, statistical downscaling and regional climate modelling are complementary, the best results are obtained when high resolution RACMO2 data are statistically refined. Model estimates in the accumulation zone do not benefit from statistical downscaling; hence, a resolution of about 25 km is sufficient to resolve the accumulation zone of the Greenland Ice Sheet. Finally, additional tests show that model resolution is as important as properly resolving spatial albedo patterns, correctly initialising the firn column and uncertainties in the modelled precipitation and turbulent exchange.