Questions: Parameterized Thermal Models of Planetary Interiors

Question 1

In a parameterized thermal model, a planet's mantle becomes hotter due to an initial abundance of radioactive elements. What happens to the cooling rate?

• A. Cooling slows — hotter materials have higher thermal conductivity, retaining heat
• B. Cooling accelerates — a hotter mantle drives more vigorous convection, increasing heat transport
• C. Cooling rate is unchanged — it depends only on surface temperature, not interior temperature
• D. The planet heats up further in a runaway feedback loop

Answer: (B)

The Nusselt-Rayleigh parameterization captures a key negative feedback: a hotter mantle has lower viscosity, which increases the Rayleigh number (convective vigor), which increases the Nusselt number (ratio of actual to conductive heat transport). More vigorous convection transports heat more efficiently to the surface, where it is radiated to space — cooling the mantle faster. This self-regulating feedback prevents runaway heating and governs the long-term thermal evolution. A planet with more radiogenic heating initially cools faster, not slower, because the thermal feedback is stronger.

Question 2

Two planets are identical except that Planet A is larger. All else equal, which is more likely to retain active plate tectonics after 4 billion years?

• A. Planet B (smaller) — less mass means faster cooling, maintaining vigorous convection
• B. Planet A (larger) — more radioactive elements in a larger volume means more sustained heating
• C. They cool at identical rates because surface-area-to-volume ratio doesn't affect thermal models
• D. Planet B — smaller planets have thinner lithospheres that are easier to subduct

Answer: (B)

A larger planet has more radioactive elements (total inventory scales with volume), sustaining internal heating for longer. More critically, a larger planet has a lower surface-area-to-volume ratio, meaning heat escapes more slowly relative to the volume that must cool. Both effects keep the mantle hotter and convection more vigorous for longer, maintaining the conditions for plate tectonics. Mars (smaller than Earth) is thought to have lost plate tectonics early partly for this reason. Planet size is a primary predictor of thermal longevity.

Question 3

In planetary thermal models, a hotter mantle loses heat more slowly because high-temperature materials conduct heat less efficiently.

• T. True
• F. False

Answer: False

The opposite is true for the convective component, which dominates heat transport in planetary mantles. A hotter mantle has lower viscosity, which increases the Rayleigh number and drives more vigorous convection. The Nusselt number increases with Rayleigh number, meaning total heat transport (dominated by convection) increases with mantle temperature. The negative feedback between temperature and cooling rate is what stabilizes planetary thermal evolution. Conduction alone (which does not strongly depend on temperature) is only dominant in the rigid lithosphere.

Question 4

Parameterized thermal models predict that as a planet cools over billions of years, its lithosphere tends to thicken.

• T. True
• F. False

Answer: True

The lithosphere is the rigid outer shell where conduction dominates over convection. As the mantle cools, viscosity increases and the temperature at which material transitions from ductile (convecting) to rigid (conductive) deepens. This progressive thickening has major geological consequences: it suppresses volcanic activity, can shut down plate tectonics (transition to stagnant-lid regime), and reduces the likelihood of generating a magnetic dynamo. Mars is the archetypal example — its early lithosphere thickened rapidly, freezing the planet into a single-plate state.

Question 5

Explain the negative feedback loop that regulates planetary cooling in parameterized thermal models, and why it matters for long-term thermal evolution.

Model answer: A hotter mantle has lower viscosity, increasing the Rayleigh number (ratio of convective to diffusive heat transport). A higher Rayleigh number increases the Nusselt number via the power-law parameterization, meaning the mantle convects more vigorously and transports heat to the surface more efficiently. Faster heat loss cools the mantle, raising viscosity and reducing the Rayleigh number — a classic negative feedback. This self-regulation prevents runaway heating and stabilizes the cooling trajectory, but it also means a planet cannot indefinitely maintain vigorous convection: as it cools, convection weakens, the lithosphere thickens, and geological activity fades.

This feedback loop is encoded in the coupled ODEs — one for mantle temperature, one for core temperature — that parameterized models solve forward in time. The power-law Nusselt-Rayleigh relationship is the key parameterization that replaces the full 3D convection calculation with a tractable 1D description. Its accuracy is validated by comparing model predictions (volcanism timing, lithospheric thickness, magnetic field lifetime) against observational constraints from crater counts and remote sensing.