A machine learning model accurately predicts the spatial distribution of geothermal gradients in Colombia using global-scale geological and geophysical datasets.