ASTR 320 Lecture 04 Notes: Planetary Heat Budgets, Radiogenic Heating, and Convection

Learning Objectives

Identify the sources of a planet's internal heat: accretional heat, differentiation heat, radiogenic decay, and (for some giant planets) ongoing gravitational contraction.
Derive the exponential decay law for radiogenic heat production and apply it to Earth's major heat-producing isotopes.
Explain the Rayleigh number criterion for the onset of solid-state mantle convection.
Compare surface heat flow across terrestrial planets and connect it to their present-day geological activity.

Terms Developed in Context

These terms organize the lecture's progression from measured planetary properties to physical interpretation. Each entry expands the matching explorer panel and connects it to the lecture's evidence, model, or calculation.

radiogenic heating

Tidal heating is flagged here, ahead of its full quantitative treatment in Lecture 13, specifically to prevent an overgeneralization of this lecture's surface-to-volume cooling argument: Io is smaller than the Moon yet is the most volcanically active body in the solar system, precisely because it receives an enormous non-radiogenic heat input from tidal flexing that the simple size-scaling argument does not capture.

Use and relationship: The Rayleigh Number and Radiogenic Heating Combined The worked example and the lecture's comparison cases show how the term changes an observable or model prediction.

Boundary or confusion: The Rayleigh number is the same dimensionless number used in fluid dynamics and atmospheric science generally (it also governs, for example, convection in a pot of heated water or in a star's envelope); its application here to a nominally "solid" mantle is possible only because geological timescales (millions of years) are long enough for silicate rock to behave as an extremely viscous fluid, an idea first proposed to explain continental drift and now confirmed by both seismic tomography (imaging actual mantle plumes) and geodetic measurement of present-day plate motions.

half-life

Radiogenic heat, decaying exponentially with each isotope's own half-life, is Earth's dominant present-day internal heat source, though secular cooling contributes comparably.

Use and relationship: Because ⁴⁴K has the shortest half-life of the four major heat producers, Earth's radiogenic heat production was substantially higher in the past and has been declining ever since formation. The worked example and the lecture's comparison cases show how the term changes an observable or model prediction.

Rayleigh number

Earth's mantle Rayleigh number is estimated at 10⁶-10⁷, vastly super-critical, meaning vigorous convection rather than marginal stability is expected -- consistent with plate tectonics.

Use and relationship: The Rayleigh number is the same dimensionless number used in fluid dynamics and atmospheric science generally (it also governs, for example, convection in a pot of heated water or in a star's envelope); its application here to a nominally "solid" mantle is possible only because geological timescales (millions of years) are long enough for silicate rock to behave as an extremely viscous fluid, an idea first proposed to explain continental drift and now confirmed by both seismic tomography (imaging actual mantle plumes) and geodetic measurement of present-day plate motions. The worked example and the lecture's comparison cases show how the term changes an observable or model prediction.

mantle convection

Distinguishing these heat sources matters because they have very different timescales: accretional and differentiation heat were front-loaded within the first ∼100 Myr and have long since been radiated away from small bodies, while radiogenic heat decays on the multi-billion-year half-lives of ²³⁶U, ²³&sup8;U, ²³²Th, and ⁴⁴K, meaning it is still the dominant heat source powering Earth's mantle convection and dynamo today, 4.5 Gyr after formation.

Use and relationship: The exponential decay law here is identical in form to the radioactive decay students have seen in introductory nuclear physics, but its planetary-science consequence is important: Earth's mantle was heated more vigorously by radioactivity in the Archean era (roughly 2.5-4 Gyr ago) than today, which is part of the explanation for why plate tectonics, mantle convection vigor, and volcanic/magnetic activity generally were more intense earlier in Earth's history and have been gradually declining since. The worked example and the lecture's comparison cases show how the term changes an observable or model prediction.

stagnant lid

Stagnant-lid convection (Venus, Mars, Mercury, the Moon, most other rocky bodies): the cold surface layer is too strong (or too dry) to break into mobile plates and instead sits as a single immobile shell atop a convecting interior.

Use and relationship: Mantle Convection Regimes: Stagnant Lid vs. The worked example and the lecture's comparison cases show how the term changes an observable or model prediction.

surface-to-volume scaling

This surface-to-volume scaling argument is one of the most powerful and general results in all of planetary science: it requires no detailed knowledge of a specific body's composition or history to predict, correctly, that smaller bodies cool and become geologically inactive faster than larger ones, and it is the same scaling argument (applied to stars rather than planets) that explains why low-mass stars live far longer than high-mass stars in ASTR 330.

Use and relationship: Tidal heating is flagged here, ahead of its full quantitative treatment in Lecture 13, specifically to prevent an overgeneralization of this lecture's surface-to-volume cooling argument: Io is smaller than the Moon yet is the most volcanically active body in the solar system, precisely because it receives an enormous non-radiogenic heat input from tidal flexing that the simple size-scaling argument does not capture. The worked example and the lecture's comparison cases show how the term changes an observable or model prediction.

Sources of Planetary Interior Heat

Accretional heat: kinetic energy of infalling planetesimals converted to heat during formation, largely radiated away or retained depending on burial depth and accretion rate.
Differentiation heat: gravitational potential energy released as dense material (iron) sinks to form a core (Lecture 03), a one-time but substantial heat source.
Radiogenic heat: ongoing heat production from the decay of long-lived radioactive isotopes (⁴⁴⁴K, ²³⁵U, ²³⁾U, ²³⁆Th) distributed through the mantle and crust.
Giant planets add a further source: slow gravitational (Kelvin-Helmholtz) contraction, and for Saturn, possible helium rain releasing gravitational energy (Lecture 11).

Radiogenic Heat Production: The Decay Law

Each radioactive isotope's heat production rate decays exponentially: H(t) = H₀ e⁻λt, with decay constant λ = ln2 / t₁⁄₂.
Earth's present-day radiogenic heat budget is dominated by four isotopes with half-lives ranging from 1.25 Gyr (⁴⁴K) to 14.0 Gyr (²³²Th).
Because ⁴⁴K has the shortest half-life of the four major heat producers, Earth's radiogenic heat production was substantially higher in the past and has been declining ever since formation.

The exponential decay law here is identical in form to the radioactive decay students have seen in introductory nuclear physics, but its planetary-science consequence is important: Earth's mantle was heated more vigorously by radioactivity in the Archean era (roughly 2.5-4 Gyr ago) than today, which is part of the explanation for why plate tectonics, mantle convection vigor, and volcanic/magnetic activity generally were more intense earlier in Earth's history and have been gradually declining since.

Present-Day Radiogenic Heat Production Rates

⁴⁴K (t₁⁄₂ = 1.25 Gyr): contributes roughly 40% of Earth's present bulk silicate Earth radiogenic heat production, despite modest abundance, because of its shorter half-life.
²³⁶U (t₁⁄₂ = 4.47 Gyr) and ²³²Th (t₁⁄₂ = 14.0 Gyr): together contribute most of the remainder, with Th typically contributing somewhat more heat than U owing to greater cosmochemical abundance.
Total present-day radiogenic heat production for the bulk silicate Earth is estimated at roughly 20 TW, compared to a total observed surface heat flow of roughly 47 TW -- implying a substantial non-radiogenic (secular cooling) contribution as well.

The roughly 20 TW versus 47 TW gap between present radiogenic heat production and total observed surface heat flow is itself an important and actively debated result: it implies that Earth is still cooling from its formation (secular cooling of the core and mantle) at a rate comparable to its radiogenic heat production, a conclusion drawn from combining bulk silicate Earth geochemical models with direct heat-flow measurements and, since 2005, direct measurement of Earth's geoneutrino flux (from KamLAND and Borexino), which independently confirms the radiogenic heat production estimate.

The Rayleigh Number and the Onset of Convection

The Rayleigh number Ra = (ρ g α ΔT d³) / (κ μ) compares the destabilizing effect of thermal buoyancy to the stabilizing effects of viscous dissipation and thermal diffusion.
Convection begins once Ra exceeds a critical value Ra_c ∼ 1000-2000 (depending on boundary conditions); solid-state mantle convection occurs because Earth's mantle, though nominally solid, deforms plastically over geological timescales at very low effective viscosity (∼10²¹ Pa·s).
Earth's mantle Rayleigh number is estimated at 10⁶-10⁷, vastly super-critical, meaning vigorous convection rather than marginal stability is expected -- consistent with plate tectonics.

The Rayleigh number is the same dimensionless number used in fluid dynamics and atmospheric science generally (it also governs, for example, convection in a pot of heated water or in a star's envelope); its application here to a nominally "solid" mantle is possible only because geological timescales (millions of years) are long enough for silicate rock to behave as an extremely viscous fluid, an idea first proposed to explain continental drift and now confirmed by both seismic tomography (imaging actual mantle plumes) and geodetic measurement of present-day plate motions.

Mantle Convection Regimes: Stagnant Lid vs. Mobile Lid

Mobile-lid convection (Earth's plate tectonics): the cold surface layer participates directly in convection, recycling into the mantle at subduction zones.
Stagnant-lid convection (Venus, Mars, Mercury, the Moon, most other rocky bodies): the cold surface layer is too strong (or too dry) to break into mobile plates and instead sits as a single immobile shell atop a convecting interior.
Earth appears to be the only terrestrial planet with active plate tectonics today, making it observationally anomalous rather than the default expectation.

Understanding why Earth alone sustains mobile-lid convection is one of the biggest open questions in comparative planetology (raised in Lecture 01) and is thought to depend on a combination of factors including surface water (which weakens rock rheology, promoting plate boundary failure), a particular range of mantle temperature and viscosity, and possibly even the presence of pre-existing weaknesses from a specific accretional history -- no single factor has been shown to be both necessary and sufficient, which is why Lecture 05 treats plate tectonics as a genuinely open comparative problem rather than a solved one.

Heat Flow Measurements Across the Terrestrial Planets

Earth: mean surface heat flow ∼91-92 mW/m², with strong contrast between young oceanic ridges (>200 mW/m²) and old continental cratons (∼40-60 mW/m²).
Moon: Apollo heat-flow probes measured 16-21 mW/m², an order of magnitude below Earth, consistent with its small size, thick stagnant lid, and largely extinct interior activity.
Mars: InSight's heat-flow probe (HP³) failed to achieve its planned burial depth, but independent geophysical estimates suggest ∼20-24 mW/m², intermediate between Earth and the Moon.

The strong contrast between Earth's young oceanic ridges and old continental interiors is itself a heat-flow signature of plate tectonics: new oceanic lithosphere is thin and hot immediately after formation at a mid-ocean ridge and then cools and thickens as it ages and moves away from the ridge, so heat flow measured at the seafloor is a direct, continuously renewed record of the plate-recycling process that Lecture 05 develops in full.

Why Small Bodies Cool Faster: Surface-to-Volume Scaling

A body's total heat generation scales with its volume (∝ R³), while its capacity to lose heat scales with its surface area (∝ R²).
The ratio of heat loss to heat generation therefore scales as R²/R³ = 1/R: smaller bodies lose heat proportionally faster and cool on shorter timescales.
This single scaling argument explains why the Moon and Mercury are geologically quiescent today while Earth (and, differently, Venus) remain active.

The Rayleigh Number and Radiogenic Heating Combined

A planet's convective vigor depends jointly on its radiogenic heat production (this lecture's exponential decay law) and its size (this lecture's surface-to-volume scaling).
Mars, intermediate in size between the Moon and Earth, is a natural test case: it appears to have lost most active mantle convection and any global magnetic dynamo (Lecture 10) roughly 4 Gyr ago, consistent with faster cooling than Earth but slower than the Moon.
This combined argument previews Lab 02's quantitative heat-budget calculation for Mars versus Earth.

Mars sits at an instructive middle ground in almost every heat-budget comparison in this course: too small to sustain Earth-like plate tectonics or an active dynamo today, but large enough to have retained substantial early geological and magnetic activity (recorded in ancient crustal magnetization, discussed further in Lecture 10), making it the clearest natural test of the surface-to-volume cooling argument developed in this section.

Tidal Heating as a Non-Radiogenic Heat Source

For satellites in eccentric or resonant orbits (Io, Europa, Enceladus), periodic tidal flexing dissipates orbital and rotational energy as internal heat.
Tidal heating can dominate over radiogenic heating in small icy or rocky satellites, as this course shows quantitatively for Io in Lecture 13.
This heat source has no analog in isolated planets and depends on orbital dynamics (eccentricity, orbital resonance) rather than only on a body's size and radioactive inventory.

Reading and Scope Note

OpenStax Astronomy 2e Section 9.1-9.2 (Earth's interior and plate tectonics) introduces heat flow and mantle convection at an introductory level.
The Rayleigh-number formalism and the quantitative radiogenic decay-law treatment go beyond OpenStax and follow Turcotte & Schubert, Geodynamics, Chapters 4 and 6.

Students should treat this lecture's Rayleigh-number and decay-law derivations as the primary quantitative content, with OpenStax serving only as descriptive background, consistent with this course's stated policy of exceeding the introductory textbook's depth wherever a genuinely upper-division derivation is available.

Worked Example: Earth's Present-Day Radiogenic Heat Production

Model each isotope's heat production as H(t) = H₀ e⁻λt with λ = ln2/t₁⁄₂, using standard bulk-silicate-Earth present-day heat production rates (per kilogram of isotope) and Earth mantle abundances.

Half-lives used: ⁴⁴K = 1.25 Gyr (λ = 0.555 Gyr⁻&sup9;), ²³⁶U = 4.47 Gyr (λ = 0.155 Gyr⁻&sup9;), ²³²Th = 14.0 Gyr (λ = 0.050 Gyr⁻&sup9;).

At t = 4.5 Gyr after formation, the surviving fraction of each isotope's initial inventory is e⁻λt: ⁴⁴K retains 8.25%, ²³⁶U retains 49.8%, and ²³²Th retains 80.0% of its original abundance.

Result: ⁴⁴K has decayed by more than a factor of 8 over Earth history, while ²³²Th has barely decayed at all -- exactly why ⁴⁴K, despite being a minor heat contributor today, dominated Earth's radiogenic heat budget much more strongly in the distant past, and why the present-day total (∼20 TW) is well below the ∼60-80 TW radiogenic heat production Earth likely had shortly after formation.

This calculation uses only the decay law H(t)=H0 e^(-lambda t) applied to three of Earth's four major heat-producing isotopes; Lab 02 extends it to a full quantitative comparison of Earth versus Mars total heat budgets using each planet's estimated bulk composition.

Summary

Radiogenic heat, decaying exponentially with each isotope's own half-life, is Earth's dominant present-day internal heat source, though secular cooling contributes comparably.
The Rayleigh number formalizes when a planetary interior convects; Earth's mantle is vastly super-critical, consistent with observed plate tectonics.
Surface-to-volume scaling (heat loss ∝ 1/R relative to heat generation) explains why small bodies cool and become geologically quiescent faster than large ones, though tidal heating is an important exception developed later.

Looking ahead: Lecture 05 examines the surface expression of this convecting interior: plate tectonics on Earth and its absence, so far, everywhere else in the solar system.