34. Earth’s internal heat

Opening disclaimer: I’m going to be talking about “heat” a lot in this one. Formally, “heat” is defined as a process of energy flow, and not as an amount of thermal energy in a body. However to people who aren’t experts in thermodynamics (i.e. nearly everyone), “heat” is commonly understood as an “amount of hotness” or “amount of thermal energy”. To avoid the linguistic awkwardness of using the five-syllable phrase “thermal energy” in every single instance, I’m just going to use this colloquial meaning of “heat”. Even some of the papers I cite use “heat” in this colloquial sense. I’ve already done it in the title, which to be technically correct should be the more awkward and less pithy “Earth’s internal thermal energy”.

The interior of the Earth is hot. Miners know first hand that as you go deeper into the Earth, the temperature increases. The deepest mine on Earth is the TauTona gold mine in South Africa, reaching 3.9 kilometres below sea level. At this depth, the rock temperature is 60°C, and considerable cooling technology is required to bring the air temperature down to a level where the miners can survive. The Kola Superdeep Borehole in Russia reached a depth of 12.2 km, where it found the temperature to be 180°C.

Deeper in the Earth, the temperature gets hot enough to melt rock. The results are visible in the lava that emerges from volcanic eruptions. How did the interior of the Earth get that hot? And exactly how hot is it down there?

For many years, geologists have been measuring the amount of thermal energy flowing out of the Earth, at thousands of measuring stations across the planet. A 2013 paper analyses some 38,374 heat flow measurements across the globe to produce a map of the mean heat flow out of the Earth, shown below[1]:

From the map, you can see that most of Earth’s heat emerges at the mid-ocean ridges, deep underwater. This makes sense, as this is where rising plumes of magma from deep within the mantle are acting to bring new rock material to the crust. The coolest areas are generally geologically stable regions in the middle of tectonic plates.

Although the heat flow out of the Earth’s surface is of the order of milliwatts per square metre, the surface has a lot of square metres. The overall heat flow out of the Earth comes to a total of around 47 terawatts[2]. In contrast, the sun emits close to 4×10¹⁴ terawatts of energy in total, and the solar energy falling on the Earth’s surface is 1360 watts per square metre, over 10,000 times as much as the heat energy leaking out of the Earth itself. So the sun dominates Earth’s heating and weather systems by roughly that factor.

So the Earth generates some 47 TW of thermal power. Where does this huge amount of energy come from? To answer that, we need to go all the way back to when the Earth was formed, some 4.5 billion years ago.

Our sun formed from the gaseous and dusty material distributed throughout the Galaxy. This material is not distributed evenly, and where there is a denser concentration, gravity acts to draw in more material. As the material is pulled in, any small motions are amplified into an overall rotation. The result is an accretion disc, with matter spiralling into a growing mass at the centre. When the central concentration accumulates enough mass, the pressure ignites nuclear reactions and a star is born. Some of the leftover material continues to orbit the new star and forms smaller accretions that eventually become planets or smaller bodies.

The process of accreting matter generates thermal energy. Gravitational potential energy reduces as matter pulls closer together, and the resulting collisions between matter particles convert it into thermal energy, heating up the accumulating mass. Our Earth was born hot. As the matter settled into a solid body, the shrinking further heated the core through the Kelvin-Helmholtz mechanism. The total heat energy from the initial formation of the Earth dissipates only very slowly into space, and that process is still going on today, 4.5 billion years later.

It’s not known precisely how much of this primordial heat is left in Earth or how much flows out, but various different studies suggest it is somewhere in the range of 12-30 TW, roughly a quarter to two-thirds of Earth’s total measured heat flux[3]. So that’s not the only source of the heat energy flowing out of the Earth.

The other source of Earth’s internal heat is radioactive decay. Some of the matter in the primordial gas and dust cloud that formed the sun and planets was produced in the supernova explosions of previous generations of stars. These explosions produce atoms of radioactively unstable isotopes. Many of these decay relatively rapidly and are essentially gone by now. But some isotopes have very long half-lives, most importantly: potassium-40 (1.25 billion years), thorium-232 (14.05 billion years), uranium-235 (703.8 million years), and uranium-238 (4.47 billion years). These isotopes still exist in significant quantities inside the Earth, where they continue to decay, releasing energy.

We have a way of probing how much radioactive energy is released inside the Earth. The decay reactions produce neutrinos (which we’ve met before), and because they travel unhindered through the Earth these can be detected by neutrino observatories. These geoneutrinos have energy ranges that distinguish them from cosmic neutrino sources, and of course always emerge from underground. The observed decay rates from geoneutrinos correspond to a total radiothermal energy production of 10-30 TW, of the same order as the primordial heat flux. (The neutrinos themselves also carry away part of the energy from the radioactive decays, roughly 5 TW, but this is an additional component not deposited as thermal energy inside the Earth.)

To within the uncertainties, the sum of the estimated primordial and measured radiothermal energy fluxes is equal to the total measured 47 TW flux. So that’s good.

Once you know how much heat is being generated inside the Earth, you can start to apply heat transfer equations, knowing the thermodynamic properties of rock and iron, how much conduction and convection can be expected, and cross-referencing it with our knowledge of the physical state of these materials under different temperature and pressure conditions. There’s also additional information about the internal structure of the Earth that we get from seismology, but that’s a story for a future article. Putting it all together, you end up with a linked series of equations which you can solve to determine the temperature profile of the Earth as a function of depth.

The results are all self-consistent, with observations such as the temperature of the rock in deep mine shafts and the rate of detection of geoneutrinos, with structural constraints provided by seismology, and with the temperature constraints and known modes of heat flow from the core to the surface of the Earth.

That is, they’re all consistent assuming the Earth is a spherical body of rock and iron. If the Earth were flat, the thermal transport equations would need to be changed to reflect the different geometry. As a first approximation, assume the flat Earth is relatively thin (i.e. a cylinder with the radius larger than the height). We still measure the same amount of heat flux emerging from the Earth’s surface, so the same amount of heat has to be either (a) generated inside it, or (b) being input from some external energy source underneath the flat Earth. However geoneutrino energy ranges indicate that they come from radioactive decay of Earthly minerals, so it makes sense to conclude that radiothermal heating is significant.

If radioactive decay is producing heat within the bulk of the flat Earth, then half of the produced neutrinos will emerge from the underside, and thus be undetectable. So the total heat production should be double that deduced from neutrino observations, or somewhere in the range 20-60 TW. To produce twice the energy, you need twice the mass of the Earth. If the flat Earth is a disc with radius 20,000 km (the distance from the North Pole to the South Pole), then to have the same volume as the spherical Earth it would need to be 859 km thick. But we need twice as much mass to produce the observed thermal energy flux, so it should be approximately 1720 km thick. Some fraction of the geoneutrinos will escape from the sides of the cylinder of this thickness, which means we need to add more rock to produce a bit more energy to compensate, so the final result will be a bit thicker.

There’s no obvious reason to suppose that a flat Earth can’t be a bit over 1700 km thick, as opposed to any other thickness. With over twice as much mass as our spherical Earth, the surface gravity of this thermodynamically correct flat Earth would be over 2 Gs (i.e. twice the gravity we experience), which is obviously wrong, but then many flat Earth models deny Newton’s law of gravity anyway (because it causes so many problems for the model).

But just as in the spherical Earth model the observed geoneutrino flux only accounts for roughly half the observed surface heat flux. The other half could potentially come from primordial heat left over from the flat Earth’s formation – although as we’ve already seen, what we know about planetary formation precludes the formation of a flat Earth in the first place. The other option is (b) that the missing half of the energy is coming from some source underneath the flat Earth, heating it like a hotplate. What this source of extra energy is is mysterious. No flat Earth model that I’ve seen addresses this problem, let alone proposes a solution.

What’s more, if such a source of energy under the flat Earth existed, then it would most likely also radiate into space around the edges of the flat Earth, and have observable effects on the objects in the sky above us. What we’re left with, if we trust the sciences of radioactive decay and thermal energy transfer, is a strong constraint on the thickness of the flat Earth, plus a mysterious unspecified energy source underneath – neither of which are mentioned in standard flat Earth models.

References:

[1] Davies, J. H. “Global map of solid Earth surface heat flow”. Geochemistry, Geophysics, Geosystems, 14(10), p.4608-4622, 2013. https://doi.org/10.1002/ggge.20271

[2] Davies, J.H., Davies, D.R. “Earth’s surface heat flux”. Solid Earth, 1(1), p.5-24, 2010. https://doi.org/10.5194/se-1-5-2010

[3] Dye, S.T. “Geoneutrinos and the radioactive power of the Earth”. Reviews of Geophysics, 50(3), 2012. https://doi.org/10.1029/2012RG000400

[4] Arevalo Jr, R., McDonough, W.F., Luong, M. “The K/U ratio of the silicate Earth: Insights into mantle composition, structure and thermal evolution”. Earth and Planetary Science Letters, 278(3-4), p.361-369, 2009. https://doi.org/10.1016/j.epsl.2008.12.023

[5] Boehler, R. “Melting temperature of the Earth’s mantle and core: Earth’s thermal structure”. Annual Review of Earth and Planetary Sciences, 24(1), p.15-40, 1996. https://doi.org/10.1146/annurev.earth.24.1.15