We developed a method to estimate the aerodynamic properties of the Greenland Ice Sheet surface using either UAV or ICESat-2 elevation data. We show that this new method is able to reproduce the important spatiotemporal variability in surface aerodynamic roughness, measured by the field observations. The new maps of surface roughness can be used in atmospheric models to improve simulations of surface turbulent heat fluxes and therefore surface energy and mass balance over rough ice worldwide.