Geothermal heat flux (GHF) is a crucial boundary condition for making accurate predictions of ice sheet mass loss, yet it is poorly known in Greenland due to inaccessibility of the bedrock. Here we use a machine learning algorithm on a large collection of relevant geologic features and global GHF measurements and produce a GHF map of Greenland that we argue is within ∼15% accuracy. The main features of our predicted GHF map include a large region with high GHF in central-north Greenland surrounding the NorthGRIP ice core site, and hot spots in the Jakobshavn Isbrae catchment, upstream of Petermann Gletscher, and near the terminus of Nioghalvfjerdsfjorden glacier. Our model also captures the trajectory of Greenland movement over the Icelandic plume by predicting a stripe of elevated GHF in central-east Greenland. Finally, we show that our model can produce substantially more accurate predictions if additional measurements of GHF in Greenland are provided. Plain Language Summary The heat generated at the interior regions of Earth (geothermal heat flux, GHF) can be high enough to melt the bottom layers of ice sheets, decrease friction between ice and bedrock, and increase ice discharge to the ocean. This heat, however, cannot be directly measured in ice sheets because the bedrock is inaccessible. Here we present a novel approach to estimate this heat. We combine all the available geologic, tectonic, and GHF data that are available on all continents. We then establish a complex relationship between GHF and all the geologic-tectonic features using machine learning techniques and then predict the GHF for the Greenland Ice Sheet. We utilize all information from available ice cores and bedrock boreholes to improve the GHF prediction in Greenland. Thus, the new GHF map honors tectonic settings, regional geology, and measurements from ice cores and can be used as an important input parameter to numerical ice sheet models that aim at lowering the uncertainties of future sea level rise predictions. This study derives a new map of GHF for the GrIS using statistical relationships between global heat flux observations and the combined influence of local geology and regional tectonic setting. Compilations of global RESEARCH LETTER 10.1002/2017GL075661 Key Points: • A new geothermal heat flux map of Greenland is obtained within ∼15% accuracy using machine learning techniques • The new map honors regional geology, tectonic settings, and ice core measurements • Pockets of high heat flux are predicted in central-north Greenland and upstream of several fast-flowing outlet glaciers