Modelling temperature in the Earth’s crust is accomplished by populating a geological model with thermal properties, such as thermal conductivity and rates of heat production, and solving a numerical model of thermal diffusion with assigned boundary conditions.