Simulated surface conditions of the Goddard Earth Observing System model, version 5 (GEOS-5), atmospheric general circulation model (AGCM) are examined for the contemporary Greenland Ice Sheet (GrIS). A surface parameterization that explicitly models surface processes including snow compaction, meltwater percolation and refreezing, and surface albedo is found to remedy an erroneous deficit in the annual net surface energy flux and provide an adequate representation of surface mass balance (SMB) in an evaluation using simulations at two spatial resolutions. The simulated 1980–2008 GrIS SMB average is 24.7 ± 4.5 cm yr−1 water-equivalent (w.e.) at ½° model grid spacing, and 18.2 ± 3.3 cm yr−1 w.e. for 2° grid spacing. The spatial variability and seasonal cycle of the ½° simulation compare favorably to recent studies using regional climate models, while results from 2° integrations reproduce the primary features of the SMB field. In comparison to historical glaciological observations, the coarser-resolution model overestimates accumulation in the southern areas of the GrIS, while the overall SMB is underestimated. These changes relate to the sensitivity of accumulation and melt to the resolution of topography. The GEOS-5 SMB fields contrast with available corresponding atmospheric models simulations from phase 5 of the Coupled Model Intercomparison Project (CMIP5). It is found that only a few of the CMIP5 AGCMs examined provide significant summertime runoff, a dominant feature of the GrIS seasonal cycle. This is a condition that will need to be remedied if potential contributions to future eustatic change from polar ice sheets are to be examined with GCMs.