The temperature distribution in the crust, and associated uncertainty, was simulated from the ensemble of Curie depth realisations assigned to a lower thermal boundary condition of a crustal model (inc. sedimentary thickness, Moho depth, heat production, thermal conductivity), constructed from various geophysical and geochemical data sets.
