# Thermal history of the Earth

The thermal history of the Earth involves the study of the cooling history of Earth's interior. It is a sub-field of geophysics. (Thermal histories are also computed[by whom?] for the internal cooling of other planetary and stellar bodies.) The study of the thermal evolution of Earth's interior is uncertain and controversial in all aspects, from the interpretation of petrologic observations used to infer the temperature of the interior, to the fluid dynamics responsible for heat loss, to material properties that determine the efficiency of heat transport.

## Overview

Observations that can be used to infer the temperature of Earth's interior range from the oldest rocks on Earth to modern seismic images of the inner core size. Ancient volcanic rocks can be associated with a depth and temperature of melting through their geochemical composition. Using this technique and some geological inferences about the conditions under which the rock is preserved, the temperature of the mantle can be inferred. The mantle itself is fully convective, so that the temperature in the mantle is basically constant with depth outside the top and bottom thermal boundary layers. This is not quite true because the temperature in any convective body under pressure must increase along an adiabat, but the adiabatic temperature gradient is usually much smaller than the temperature jumps at the boundaries. Therefore, the mantle is usually associated with a single or potential temperature that refers to the mid-mantle temperature extrapolated along the adiabat to the surface. The potential temperature of the mantle is estimated to be about 1350 C today. There is an analogous potential temperature of the core but since there are no samples from the core its present-day temperature relies on extrapolating the temperature along an adiabat from the inner core boundary, where the iron solidus is somewhat constrained.

## Thermodynamics

The simplest mathematical formulation of the thermal history of Earth's interior involves the time evolution of the mid-mantle and mid-core temperatures. To derive these equations one must first write the energy balance for the mantle and the core separately. They are,

${\displaystyle Q_{\text{surf}}=Q_{\text{sec,man}}+Q_{\text{rad}}+Q_{\text{cmb}}}$

for the mantle, and

${\displaystyle Q_{\text{cmb}}=Q_{\text{sec,core}}+Q_{\text{L}}+Q_{\text{G}}}$

for the core. ${\displaystyle Q_{\text{surf}}}$ is the surface heat flow [W] at the surface of the Earth (and mantle), ${\displaystyle Q_{\text{sec,man}}=M_{\text{man}}c_{\text{man}}dT_{\text{man}}/dt}$ is the secular cooling heat from the mantle, and ${\displaystyle M_{\text{man}}}$, ${\displaystyle c_{\text{man}}}$, and ${\displaystyle T_{\text{man}}}$ are the mass, specific heat, and temperature of the mantle. ${\displaystyle Q_{\text{rad}}}$ is the radiogenic heat production in the mantle and ${\displaystyle Q_{\text{cmb}}}$ is the heat flow from the core mantle boundary. ${\displaystyle Q_{\text{sec,core}}=M_{\text{core}}c_{\text{core}}dT_{\text{core}}/dt}$ is the secular cooling heat from the core, and ${\displaystyle Q_{\text{L}}}$ and ${\displaystyle Q_{\text{G}}}$ are the latent and gravitational heat flow from the inner core boundary due to the solidification of iron.

Solving for ${\displaystyle dT_{\text{man}}/dt}$ and ${\displaystyle dT_{\text{core}}/dt}$ gives,

${\displaystyle {\frac {dT_{\text{man}}}{dt}}={\frac {3(-Q_{\text{surf}}-Q_{\text{cmb}})}{4\pi \rho _{\text{m}}c_{\text{m}}(R^{3}-R_{\text{c}}^{3})}}+{\frac {Q_{\text{rad}}}{V_{\text{m}}\rho _{\text{m}}c_{\text{m}}}}}$

and,

${\displaystyle {\frac {dT_{\text{core}}}{dt}}=Q_{\text{cmb}}\left[A_{\text{c}}(L+E_{G})\left({\frac {R_{i}}{R_{\text{c}}}}\right)^{2}\rho _{i}{\frac {dR_{i}}{dT_{\text{cmb}}\eta _{\text{c}}}}-{\frac {R_{\text{c}}^{3}-R_{i}^{3}}{3R_{\text{c}}^{3}}}\rho _{\text{c}}c_{\text{c}}\right]^{-1}}$

## Thermal Catastrophe

In 1862, Lord Kelvin calculated the age of the Earth by assuming that Earth had formed as a completely molten object, and determined the amount of time it would take for the near-surface to cool to its present temperature. Since uniformitarianism required a much older Earth, there was a contradiction. Eventually, the additional heat sources within the Earth were discovered, allowing for a much older age. This section is about a similar paradox in current geology, called the thermal catastrophe.

The thermal catastrophe of the Earth can be demonstrated by solving the above equations for the evolution of the mantle with ${\displaystyle Q_{\text{cmb}}=0}$. The catastrophe is defined as when the mean mantle temperature ${\displaystyle T_{\text{man}}}$ exceeds the mantle solidus so that the entire mantle melts. Using the geochemically preferred Urey ratio of ${\displaystyle Ur=1/3}$ and the geodynamically preferred cooling exponent of ${\displaystyle {\text{beta}}=1/3}$ the mantle temperature reaches the mantle solidus (i.e. a catastrophe) in 1-2 Ga. This result is clearly unacceptable because geologic evidence for a solid mantle exists as far back as 4 Ga (and possibly further). Hence, the thermal catastrophe problem is the foremost paradox in the thermal history of the Earth.

The "New Core Paradox" [1] posits that the new upward revisions to the empirically measured thermal conductivity of iron [2][3][4] at the pressure and temperature conditions of Earth's core imply that the dynamo is thermally stratified at present, driven solely by compositional convection associated with the solidification of the inner core. However, wide spread paleomagnetic evidence for a geodynamo [5] older than the likely age of the inner core (~1 Gyr) creates a paradox as to what powered the geodynamo prior to inner core nucleation. Recently it has been proposed that a higher core cooling rate and lower mantle cooling rate can resolve the paradox in part.[6][7][8] However, the paradox remains unresolved.

Two additional constraints have been recently proposed. Numerical simulations of the material properties of high pressure-temperature iron [9] claim an upper limit of 105 W/m/K to the thermal conductivity. This downward revision to the conductivity partially alleviates the issues of the New Core Paradox by lowering the adiabatic core heat flow required to keep the core thermally convective. Also, recent geochemical experiments [10] have led to the proposal that radiogenic heat in the core is larger than previously thought. This revision, if true, would also alleviate issues with the core heat budget by providing an additional energy source back in time.

• Boehler, Reinhard (1996). "MELTING TEMPERATURE OF THE EARTH'S MANTLE AND CORE: Earth's Thermal Structure". Annual Review of Earth and Planetary Sciences. 24 (1): 15–40. Bibcode:1996AREPS..24...15B. doi:10.1146/annurev.earth.24.1.15.
• Davies, Geoffrey F. (2001). Dynamic earth : plates, plumes and mantle convection (Repr. ed.). Cambridge: Cambridge University Press. ISBN 9780521599337.
• Fowler, C.M.R. (2006). "7. Heat". The solid earth : an introduction to global geophysics (2nd ed.). Cambridge, UK: Cambridge University Press. pp. 269–325. ISBN 9780521893077.
• Jacobs, J.A. (1992). "4. The thermal history of the Earth". Deep interior of the earth (1st ed.). London: Chapman & Hall. ISBN 9780412365706.
• McKenzie, Dan; Weiss, Nigel (1975). "Speculations on the Thermal and Tectonic History of the Earth". Geophysical Journal of the Royal Astronomical Society. 42 (1): 131–174. doi:10.1111/j.1365-246X.1975.tb05855.x.
• Pollack, Henry N.; Hurter, Suzanne J.; Johnson, Jeffrey R. (1993). "Heat flow from the Earth's interior: Analysis of the global data set". Reviews of Geophysics. 31 (3): 267. Bibcode:1993RvGeo..31..267P. doi:10.1029/93RG01249.
• Sharpe, H. N.; Peltier, W. R. (1978). "Parameterized mantle convection and the Earth's thermal history". Geophysical Research Letters. 5 (9): 737–740. Bibcode:1978GeoRL...5..737S. doi:10.1029/GL005i009p00737.
• Williams, Quentin (6 October 1997). "Why is the earth's core so hot? And how do scientists measure its temperature?". Ask the experts. Scientific American. Retrieved 6 April 2013.

## References

1. ^ Olson, P. (24 October 2013). "The New Core Paradox". Science. 342 (6157): 431–432. Bibcode:2013Sci...342..431O. doi:10.1126/science.1243477.
2. ^ de Koker, N.; Steinle-Neumann, G.; Vlcek, V. (28 February 2012). "Electrical resistivity and thermal conductivity of liquid Fe alloys at high P and T, and heat flux in Earth's core". Proceedings of the National Academy of Sciences. 109 (11): 4070–4073. Bibcode:2012PNAS..109.4070D. doi:10.1073/pnas.1111841109. PMC 3306690. PMID 22375035.
3. ^ Pozzo, Monica; Davies, Chris; Gubbins, David; Alfè, Dario (11 April 2012). "Thermal and electrical conductivity of iron at Earth's core conditions". Nature. 485 (7398): 355–358. arXiv:1203.4970. Bibcode:2012Natur.485..355P. doi:10.1038/nature11031. PMID 22495307.
4. ^ Gomi, Hitoshi; Ohta, Kenji; Hirose, Kei; Labrosse, Stéphane; Caracas, Razvan; Verstraete, Matthieu J.; Hernlund, John W. (November 2013). "The high conductivity of iron and thermal evolution of the Earth's core". Physics of the Earth and Planetary Interiors. 224: 88–103. Bibcode:2013PEPI..224...88G. doi:10.1016/j.pepi.2013.07.010.
5. ^ Tarduno, J. A.; Cottrell, R. D.; Watkeys, M. K.; Hofmann, A.; Doubrovine, P. V.; Mamajek, E. E.; Liu, D.; Sibeck, D. G.; Neukirch, L. P.; Usui, Y. (4 March 2010). "Geodynamo, Solar Wind, and Magnetopause 3.4 to 3.45 Billion Years Ago". Science. 327 (5970): 1238–1240. Bibcode:2010Sci...327.1238T. doi:10.1126/science.1183445. PMID 20203044.
6. ^ Driscoll, P.; Bercovici, D. (November 2014). "On the thermal and magnetic histories of Earth and Venus: Influences of melting, radioactivity, and conductivity". Physics of the Earth and Planetary Interiors. 236: 36–51. Bibcode:2014PEPI..236...36D. doi:10.1016/j.pepi.2014.08.004.
7. ^ Labrosse, Stéphane (February 2015). "Thermal evolution of the core with a high thermal conductivity". Physics of the Earth and Planetary Interiors. 247: 36–55. Bibcode:2015PEPI..247...36L. doi:10.1016/j.pepi.2015.02.002.
8. ^ Davies, Christopher J. (April 2015). "Cooling history of Earth's core with high thermal conductivity". Physics of the Earth and Planetary Interiors. 247: 65–79. Bibcode:2015PEPI..247...65D. doi:10.1016/j.pepi.2015.03.007.
9. ^ Zhang, Peng; Cohen, R. E.; Haule, K. (28 January 2015). "Effects of electron correlations on transport properties of iron at Earth's core conditions". Nature. 517 (7536): 605–607. Bibcode:2015Natur.517..605Z. doi:10.1038/nature14090.
10. ^ Wohlers, Anke; Wood, Bernard J. (15 April 2015). "A Mercury-like component of early Earth yields uranium in the core and high mantle 142Nd". Nature. 520 (7547): 337–340. Bibcode:2015Natur.520..337W. doi:10.1038/nature14350. PMC 4413371.