28.01.03 · astronomy / solar-system

Planetary interiors and surfaces: comparative planetology, magnetic fields

stub3 tiersLean: nonepending prereqs

Anchor (Master): Stevenson, D. J. — Planetary magnetic fields (2003)

Intuition Beginner

Every planet has a layered structure: a dense metallic core at its centre, a rocky mantle around it, and a thin crust on the outside. This layering formed when the young planets were hot and molten. Heavy iron sank toward the middle, while lighter rock floated up toward the surface. This separation by density is called differentiation.

Earth is unique among the rocky planets in having plate tectonics. Its crust is broken into moving plates that recycle material, build mountains, and cause earthquakes. Mars is too small to keep this activity going, so its ancient surface is preserved. Venus is close to Earth in size but has no moving plates, possibly because its surface is too dry and stiff.

Earth also has a strong magnetic field. It is generated by churning liquid metal deep in the outer core. This invisible shield deflects harmful charged particles streaming outward from the Sun, protecting the atmosphere and life at the surface.

Visual Beginner

The table below compares the interior structure, present magnetic field, and a defining surface feature for the major bodies of the solar system.

Body	Interior structure	Magnetic field today	Defining surface feature
Earth	Iron core, silicate mantle, thin crust	Strong, active	Moving tectonic plates
Moon	Small iron core, thick rocky mantle	None	Dark basalt plains (the maria)
Mars	Iron-rich core, rocky mantle	None today	Olympus Mons, Valles Marineris
Venus	Iron core, thick rocky shell	None	Volcanic plains, no plate motion
Mercury	Very large iron core	Weak, active	Curved cliffs (lobate scarps)
Jupiter	Layer of metallic hydrogen	Strongest of any planet	Great Red Spot storm
Europa	Rocky interior, global water ocean	None (induced field)	Cracked icy shell

One Earth radius is about 6,371 kilometres. Mercury's iron core makes up roughly 75 percent of that planet's radius, whereas Earth's core fills only about 55 percent.

Worked example Beginner

Consider why small planets like Mars cooled much faster than Earth. A planet stores heat throughout its entire volume, but that heat can escape only through its outer surface. When a planet's radius doubles, its surface area grows by a factor of four (because $2 \times 2 = 4$ ), while its volume grows by a factor of eight (because $2 \times 2 \times 2 = 8$ ). The larger planet therefore has only half as much cooling surface for each unit of stored heat.

Mars has about half of Earth's radius. By this scaling, it loses its internal heat roughly twice as fast per unit of mass. Over billions of years, its metal core cooled enough that the churning stopped, and its magnetic field died out. Earth, larger and still hot inside, keeps its outer core churning and its protective field active. This single idea, that cooling rate depends on size, explains why Earth and Mars look so different today.

Check your understanding Beginner

Formal definition Intermediate+

Planetary differentiation and the iron catastrophe

A planet is said to be differentiated when its interior is stratified by density: a dense metallic core overlain by a silicate mantle and a low-density crust. Differentiation is driven by the iron catastrophe, the epoch in which heat from accretion and short-lived radioactivity (chiefly $^{26} Al$ during the first million years, and $^{40} K$ , $^{232} Th$ , $^{235} U$ , $^{238} U$ over billions of years thereafter) raised the interior temperature past the melting point of iron-rich alloys. Once melting began, dense iron-nickel sank toward the centre and lighter silicates rose, releasing further gravitational potential energy as heat in a runaway feedback ^{[Ch. 1 Planetary science]}.

The four principal thermal energy sources shaping a planetary interior are:

Accretion energy, the gravitational binding energy released as planetesimals collide and merge. For Earth this is of order $G M^{2} / R \sim 4 \times 1 0^{32} J$ .
Differentiation energy, the additional energy released as iron sinks through the mantle to form the core.
Radioactive decay, dominated by the long-lived isotopes $^{40} K$ , $^{232} Th$ , $^{235} U$ , $^{238} U$ in the present era, and by $^{26} Al$ in the first few million years.
Tidal heating, generated by time-varying distortion of a body by a massive neighbour, important for Io, Europa, and Enceladus.

The structure equations for a rocky planet

Under the assumption of hydrostatic equilibrium and spherical symmetry, the interior of a planet satisfies three coupled first-order equations. Conservation of mass gives

\frac{d m}{d r} = 4 π r^{2} ρ (r),

and balance between pressure and self-gravity (hydrostatic equilibrium) gives

\frac{d P}{d r} = - \frac{G m ( r ) ρ ( r )}{r ^{2}} .

Together with an equation of state $P = P (ρ, T)$ and the requirement that the temperature profile respect energy transport by radiation, conduction, or convection, these equations determine the radial profiles of density, pressure, and temperature. A terrestrial planet is conventionally divided into a solid inner core (iron-nickel alloy), a liquid outer core (iron alloyed with light elements such as sulfur, oxygen, and silicon), a silicate mantle (olivine and pyroxene at depth, transforming to higher-pressure phases such as perovskite near the core boundary), and a thin crust ^{[Ch. 1-6 Interiors and magnetic fields]}.

Comparative surfaces of the terrestrial bodies

Earth is the only terrestrial planet with plate tectonics. Its lithosphere is broken into about a dozen rigid plates that move over the ductile asthenosphere. New crust forms at mid-ocean ridges and descends back into the mantle at subduction zones, driven primarily by the pull of cold, dense descending slabs (slab pull) together with a smaller contribution from gravitational sliding off ridges (ridge pull). The cycle in which supercontinents assemble, rift apart, and reassemble is the Wilson cycle, with a characteristic period of roughly 400 to 500 million years ^{[Ch. 1-4 Mantle convection and plate tectonics]}.

The other terrestrial bodies exhibit stagnant-lid tectonics: a single, thick, immobile lithosphere through which heat escapes by conduction and mantle convection proceeds beneath a rigid lid. Mars preserves this style, and its small size (about half Earth's radius) led to rapid cooling. Its surface records the consequences: the immense volcanic pile of the Tharsis bulge, the 21-kilometre-tall shield volcano Olympus Mons, and the 4,000-kilometre canyon system Valles Marineris. Mars has no global magnetic field today, but the Mars Global Surveyor mission discovered strong remanent magnetisation in the oldest crust, proof that a dynamo operated before roughly 4 billion years ago.

Venus is nearly Earth's twin in size and bulk composition yet lacks plate tectonics and a magnetic field. Its lithosphere is thick and dry, plausibly because surface water, which on Earth weakens rock and enables subduction, was lost to the runaway greenhouse effect. Venus is resurfaced episodically: a catastrophic resurfacing event roughly 300 to 500 million years ago buried earlier terrain under vast volcanic plains. Mercury possesses an unusually large iron core (about 75 percent of its radius) and a weak but active magnetic field; its surface bears curved cliffs called lobate scarps, formed as the planet cooled and its interior contracted. The BepiColombo mission (ESA and JAXA) is presently en route to refine these measurements ^{[Ch. 1 Planetary science]}.

The Moon formed from debris ejected by a giant impact early in Earth's history. Its small core and rapid cooling produced a geologically dead body whose dark, basalt-flooded plains (the maria) record volcanic activity that ceased over three billion years ago. Permanently shadowed craters at the lunar poles trap water ice, a resource of considerable interest for future exploration.

Giant planets and icy moons

The gas giants Jupiter and Saturn are composed predominantly of hydrogen and helium. Deep within them, pressure exceeds roughly 1.4 million atmospheres and hydrogen enters a metallic phase, conducting electricity as electrons move freely through the compressed fluid. This vast conducting region, combined with rapid rotation and internal convection, generates the strongest magnetic fields in the solar system. Jupiter's field, some twenty thousand times stronger than Earth's, channels charged particles into intense radiation belts and is associated with long-lived storms such as the Great Red Spot.

The ice giants Uranus and Neptune contain substantial mantles of water, ammonia, and methane in hot, dense, conductive states beneath their hydrogen-helium envelopes; their magnetic fields are offset and highly tilted. Among the icy satellites, comparative planetology reaches its richest expression. Europa, Jupiter's fourth-largest moon, harbours a global ocean of liquid water beneath an ice shell. Enceladus, a small moon of Saturn, jets water vapour and ice from fractures at its south pole. Ganymede is the only moon known to generate its own magnetic field, sustained by a liquid iron core kept warm by tidal dissipation ^{[Ch. 1-6 Interiors and magnetic fields]}.

Key result: the dynamo condition and the moment-of-inertia constraint Intermediate+

Conditions for a self-sustaining planetary dynamo

A planetary magnetic field is not a frozen relic but the product of a hydromagnetic dynamo: convective motion of an electrically conducting fluid twists and amplifies any seed magnetic field through magnetic induction. The evolution of the field $B$ in a conducting fluid moving with velocity $v$ obeys the induction equation

\frac{\partial B}{\partial t} = \nabla \times (v \times B) + η \nabla^{2} B,

where $η = (μ_{0} σ)^{- 1}$ is the magnetic diffusivity, with $μ_{0}$ the permeability of free space and $σ$ the electrical conductivity. The first term on the right amplifies and transports the field through fluid motion; the second diffuses it away. A dynamo is possible only when induction dominates diffusion, that is, when the magnetic Reynolds number $R m = v L / η$ is large ^{[Earth Planet. Sci. Lett. 208 (2003) 1-11]}.

Three conditions must hold simultaneously for an active dynamo:

A sufficiently large volume of electrically conducting fluid. In terrestrial planets this is the liquid outer core; in gas giants it is the metallic hydrogen layer.
Convection vigorous enough to move the fluid and stretch the field lines. Convection requires the core to be cooled from above so that buoyancy drives circulation.
Rotation fast enough that the Coriolis force organises the flow into helical columns aligned with the rotation axis. Without rotation the induction is too disordered to sustain a coherent field.

The natural dimensionless measure of the balance between Lorentz and Coriolis forces is the Elsasser number

Λ = \frac{σ B ^{2}}{ρ Ω},

where $B$ is the field strength, $ρ$ the fluid density, and $Ω$ the angular velocity of rotation. Numerical and laboratory dynamos tend to saturate at $Λ$ of order unity, the so-called magnetostrophic balance, which predicts field strengths broadly consistent with those of Earth, Jupiter, and Saturn.

This framework explains the striking diversity of planetary magnetism. Earth satisfies all three conditions and sustains a strong dipole. Venus lacks a dynamo despite its size, most likely because its core, deprived of light-element buoyancy by a dry interior and possibly lacking convection, cannot drive the required motion; its slow rotation (243 days) further weakens the Coriolis effect. Mars once had a dynamo, recorded in remanent crustal magnetisation, but its core cooled and solidified enough to shut convection down roughly four billion years ago. Mercury sustains a weak field, about one percent of Earth's, a puzzle because its small core should have cooled: the leading explanations are an unusually thin liquid outer-core shell driven by strong inner-core boundary convection, or a thermoelectric dynamo. Ganymede is the only satellite with an internally generated field, kept active by tidal heating that maintains a liquid iron core ^{[Earth Planet. Sci. Lett. 208 (2003) 1-11]}.

The moment of inertia as a probe of internal structure

Direct sampling of a planetary interior is impossible beyond a few kilometres' depth, so internal structure is inferred from bulk constraints. The most powerful of these is the dimensionless polar moment of inertia $C / (M R^{2})$ , where $C$ is the moment of inertia about the rotation axis, $M$ the mass, and $R$ the radius. For a uniform sphere $C / (M R^{2}) = 0.4$ ; any value below this signals that mass is concentrated toward the centre. Measured values are Earth $0.3307$ , Mars $0.365$ , the Moon $0.393$ , and Mercury $0.346$ . The low value for Earth immediately demands a dense core, while Mercury's value, combined with its large iron core, fixes the core radius at roughly three-quarters of the planetary radius. Spacecraft tracking of the gravity field and rotation state (libration) supplies $C / (M R^{2})$ for bodies where it can be measured ^{[Ch. 1-6 Interiors and magnetic fields]}.

Exercises Intermediate+

Exercise 1 (easy, multiple choice).

Which combination of three ingredients is required to sustain an active planetary dynamo?

A. A solid crust, a thick atmosphere, and slow rotation

B. A conducting fluid, convection, and sufficiently rapid rotation

C. A metallic hydrogen layer, a ring system, and a molten mantle

D. Abundant radioactive isotopes, a frozen core, and tidal locking

Hint

A dynamo is an electromagnetic induction process. Identify what carries the current, what moves it, and what organises the flow.

Answer

Option B. A self-sustaining dynamo requires (1) an electrically conducting fluid to carry currents, (2) convection to move that fluid and stretch the magnetic field lines, and (3) rotation rapid enough that the Coriolis force organises the convective flow into coherent helical structures. Without any one of these, as in the cooled core of Mars or the dry, slowly rotating Venus, no dynamo operates.

Exercise 2 (easy, multiple choice).

Earth has plate tectonics but Venus does not, even though the two planets are similar in size. What is the most widely accepted explanation?

A. Venus rotates too fast for its lithosphere to break into plates

B. Venus lost its surface water to a runaway greenhouse, drying and stiffening its lithosphere so that subduction cannot operate

C. Venus has no iron core and therefore lacks the convection that drives plate motion

D. Venus is too small to retain the internal heat needed for mantle convection

Hint

On Earth, water plays a key role in weakening rock and enabling the plates to bend and descend.

Answer

Option B. Water trapped in minerals dramatically lowers rock strength and allows the lithosphere to bend and break into subducting plates. Venus underwent a runaway greenhouse effect that stripped its surface water. The resulting dry, stiff lithosphere cannot subduct, so Venus is trapped in a stagnant-lid regime and resurfaces episodically through catastrophic volcanism rather than steady plate recycling.

Exercise 4 (medium, multiple choice).

Mars Global Surveyor detected strong patches of magnetisation locked into the oldest Martian crust, but Mars has no global magnetic field today. What conclusion follows?

A. Mars never had a dynamo; the magnetisation was implanted by the solar wind

B. Mars once had an active dynamo that switched off after its core cooled, leaving only fossil remanent magnetisation in ancient crust

C. Mars's dynamo is still active but hidden beneath a non-magnetic surface layer

D. The crustal fields were produced by lightning strikes during volcanic eruptions

Hint

Crustal magnetisation is a permanent record of the field that existed when the rocks cooled through their Curie temperature.

Answer

Option B. When basaltic crust cools through its Curie temperature, it locks in a record of the ambient magnetic field, producing remanent magnetisation. The existence of strongly magnetised ancient crust on Mars, together with the absence of a present global field, shows that a dynamo once operated but ceased. Because Mars is small it cooled quickly, its core stopped convecting, and the dynamo died approximately four billion years ago.

Exercise 5 (medium, short answer).

Mercury is far smaller than Earth, so its core should have cooled substantially, yet Mercury sustains a weak magnetic field about one percent as strong as Earth's. Outline two distinct mechanisms that have been proposed to explain how Mercury maintains a dynamo despite its small size.

Hint

Consider the geometry of the core (how thick is the liquid outer layer?) and alternative sources of electrical currents in a stratified fluid.

Answer

Two leading explanations are (1) an unusually thin liquid outer-core shell: if only a narrow layer of liquid iron alloy remains above a large solid inner core, the steep temperature gradients at the inner-core boundary can drive compositional convection vigorous enough to sustain a weak dynamo; and (2) a thermoelectric dynamo, in which temperature differences between adjacent regions of a stratified, electrically conducting core generate currents through the thermoelectric effect, producing a field without the same convective requirements as a conventional dynamo. MESSENGER and BepiColombo data are used to distinguish between these scenarios.

Exercise 6 (hard, short answer).

The magnetic Reynolds number is $R m = v L / η$ , where $v$ is a characteristic velocity, $L$ a length scale, and $η$ the magnetic diffusivity. Explain qualitatively what $R m$ measures, why a large $R m$ is necessary for a dynamo, and what would happen to a seed field in a planet whose core had $R m ≪ 1$ .

Hint

The induction equation has two right-hand terms: one from fluid motion and one from diffusion. $R m$ is their ratio.

Answer

The magnetic Reynolds number compares the rate at which fluid motion stretches and amplifies the magnetic field (the induction term $\nabla \times (v \times B)$ , of size $v B / L$ ) to the rate at which ohmic diffusion erases it (the term $η \nabla^{2} B$ , of size $η B / L^{2}$ ); their ratio is $v L / η$ . When $R m ≫ 1$ , induction dominates and any seed field can be stretched, twisted, and amplified faster than it diffuses away, allowing a self-sustaining dynamo. When $R m ≪ 1$ , diffusion overwhelms induction: fluid motion is too slow or the region too small to organise the field, and any seed magnetic field decays exponentially on the diffusive timescale $L^{2} / η$ . Such a planet would have no internally generated field.

Advanced results Master

Dynamo theory: induction, Cowling's theorem, and the alpha-omega mechanism

A rigorous theory of planetary magnetism is built on the magnetohydrodynamic (MHD) equations: the induction equation coupling the velocity field $v$ to the magnetic field $B$ , the Navier-Stokes equation including the Lorentz force $J \times B$ and the Coriolis force $2 ρ Ω \times v$ , and the equations of thermal convection. A central negative result is Cowling's theorem (1934): no axisymmetric (purely two-dimensional) velocity field can sustain an axisymmetric magnetic field against ohmic decay. A working dynamo must therefore involve genuinely three-dimensional, non-axisymmetric motion. This theorem explains why dynamos require turbulent, helical convection rather than simple circulation patterns.

The alpha-omega dynamo is the paradigmatic mean-field model that resolves Cowling's obstruction. Differential rotation (the omega-effect) winds a poloidal field into a strong toroidal field, stretching field lines in the azimuthal direction. Helical convective turbulence (the alpha-effect) twists toroidal field back into a poloidal component, closing the cycle. The competition between these two effects, modulated by the magnetic back-reaction on the flow, produces the field reversals seen in the solar cycle and, on far longer timescales, in the geodynamo. The geometry of the resulting large-scale field is broadly dipolar for the rapidly rotating planets (Earth, Jupiter, Saturn) but can be dominated by higher multipoles when the dynamo operates outside the magnetostrophic regime or the convective region is unusually deep; this is invoked to explain the strongly tilted and offset fields of Uranus and Neptune ^{[Earth Planet. Sci. Lett. 208 (2003) 1-11]}.

Equation of state and the Adams-Williamson equation

The radial density profile of a planet is constrained by seismology (for Earth) and by gravity-field measurements (for other bodies). Under the assumption of an adiabatic, compositionally uniform interior, the density satisfies the Adams-Williamson equation

\frac{d ρ}{d r} = - \frac{G m ( r ) ρ ( r ) ^{2}}{K _{S} ( r ) r ^{2}},

where $K_{S}$ is the adiabatic bulk modulus, itself determined from seismic wave velocities through $K_{S} = ρ (v_{p}^{2} - \frac{4}{3} v_{s}^{2})$ . For Earth this was assembled into the Preliminary Reference Earth Model (PREM) (Dziewonski and Anderson 1981), a self-consistent table of density, pressure, and seismic velocities as functions of radius, anchored by the moment of inertia and total mass. PREM reveals a sharp jump in density at the core-mantle boundary (from about 5.5 to 9.9 g/cm $^{3}$ ) and the solid-liquid transition at the inner-core boundary, confirmed by the absence of shear-wave transmission through the outer core.

Tidal heating and the Laplace resonance

For the icy satellites, the dominant heat source is tidal dissipation. The canonical expression for the time-averaged tidal heat generated in a synchronously rotating satellite of radius $R$ orbiting a planet of mass $M_{p}$ at mean distance $a$ , with eccentricity $e$ and mean motion $n$ , is

\dot{E}_{tidal} = \frac{21}{2} \frac{k _{2}}{Q} \frac{G M _{p}^{2} R ^{5} n e ^{2}}{a ^{6}},

where $k_{2}$ is the Love number measuring the body's elastic response and $Q$ the specific dissipation factor. Io, Europa, and Ganymede are locked in the Laplace resonance, a three-body configuration that maintains non-zero eccentricities and thereby sustains tidal heat. In Io this powers the most active volcanism in the solar system; in Europa it maintains a liquid water ocean beneath an ice shell; in Enceladus (held in resonance with Dione) it drives the south-polar geysers detected by Cassini.

Comparative tectonics and thermal evolution

The style of heat loss from a terrestrial planet is governed by the viscosity contrast across its mantle. On Earth the contrast is moderate, allowing the cold surface layer to subduct and recycle; this is active-lid (plate) tectonics. When the viscosity contrast is extreme, the surface forms a rigid, immobile stagnant lid through which heat escapes only by conduction, as on Venus, Mars, and Mercury. Numerical convection models identify a third, intermittent regime in which the lid episodically founders and resurfaces the planet, plausibly explaining the catastrophic resurfacing of Venus.

The long-term energy budget is summarised by the Urey ratio

U r = \frac{Q _{conv}}{Q _{rad}},

the ratio of convective heat flow out of the mantle to the heat produced by radioactive decay within it. A value $U r < 1$ indicates that the planet is cooling secularly, drawing on primordial heat; Earth's Urey ratio is estimated at roughly 0.2 to 0.5, implying significant secular cooling over four and a half billion years. Mars and Mercury, with much smaller mantles, cooled far faster and now lie deep in the stagnant-lid regime.

Crustal and ice-shell thickness from gravity and libration

Two complementary techniques probe the thickness of a planet's outer layer. For rocky bodies, combining topography with a gravity map yields crustal thickness through a topographic-isostatic inversion; the GRAIL mission applied this to the Moon and found a mean crustal thickness of 34 to 43 kilometres, far thicker than Earth's 6 to 7 kilometre oceanic crust. For icy satellites, the shell thickness is constrained by measuring the body's forced libration, the small oscillation in rotation rate driven by an eccentric orbit. A decoupled shell floating on a global ocean librates with a measurable amplitude; the detection of such libration for Europa and Enceladus confirmed the presence of subsurface oceans and bounded the ice-shell thickness at roughly 15 to 25 kilometres for Europa and a few tens of kilometres for Enceladus ^{[Ch. 1-6 Interiors and magnetic fields]}.

Connections Master

28.01.01 — The solar system: planets, moons, and small bodies. The inventory and bulk properties (mass, radius, density) catalogued there are the raw data from which interior models are constructed; the dense iron cores and the magnetic field inventory described here follow directly from those bulk constraints.
28.01.02 pending — Formation of the solar system. Accretion energy and the $^{26} Al$ heating of early planetesimals, treated there as drivers of melting, are precisely the energy sources that triggered the iron catastrophe and differentiation analysed here. The Moon-forming giant impact is the boundary condition for Earth's early thermal state.
28.02.01 — Stars and stellar evolution. The convective, rotating, electrically conducting interiors of planets obey the same MHD equations used in stellar dynamo theory; the solar dynamo and the geodynamo are sibling problems differing mainly in the conducting fluid (plasma versus liquid metal) and the geometry.
[27.01] — Earth science: geodynamics and plate tectonics. The Wilson cycle, slab pull, mantle convection, and the Urey ratio treated here are the dynamical core of terrestrial geology; this unit supplies the planetary-mechanics derivation, while the earth-science sequence develops the observational and historical record.
28.05.01 — Exoplanets: detection methods and habitability. Interior structure and magnetic field generation determine whether a rocky exoplanet can retain an atmosphere and surface water over geological time, and therefore whether it can be habitable; the dynamo and stagnant-lid criteria developed here are applied to characterise detected exoplanets.

Historical and philosophical context Master

The recognition that planets possess layered interiors emerged from two converging threads. The first was seismology. After Andrija Mohorovicic identified the crust-mantle boundary in 1909, Beno Gutenberg in 1914 determined the radius of the Earth's core from the shadow cast by seismic waves, and Inge Lehmann in 1936 discovered the solid inner core by analysing the arrivals of waves that should not have existed. These seismic results, synthesised into global reference models culminating in PREM, transformed the interior of the Earth from speculation into a measured structure.

The second thread was the theory of continental drift. Alfred Wegener proposed in 1912 that the continents move, but his mechanism lacked a convincing driving force. The breakthrough came half a century later. Fred Vine and Drummond Matthews, and independently Lawrence Morley, showed in 1963 that the symmetric pattern of magnetic stripes flanking mid-ocean ridges recorded the reversals of the geomagnetic field as new crust formed and spread apart. This single result married palaeomagnetism to sea-floor spreading and established plate tectonics as the unifying framework of geology within a few years. The recognition that Earth is unique among the terrestrial planets in sustaining this regime transformed comparative planetology into a science of contrasts: the same equations, applied to bodies of different size and composition, yield plate tectonics, episodic resurfacing, or a frozen stagnant lid.

The theory of planetary magnetism began with Joseph Larmor's 1919 proposal that the Sun's magnetic field is generated by a self-exciting dynamo. Thomas Cowling's 1934 theorem, proving that no axisymmetric dynamo can exist, appeared to doom the idea until Walter Elsasser and Edward Bullard in the 1940s and 1950s showed that three-dimensional, turbulent convection in the Earth's liquid outer core could evade the theorem and sustain the field. The mean-field alpha-omega formalism developed by Eugene Parker and by Max Steenbeck and Friedrich Krause in the 1950s and 1960s supplied the theoretical machinery still used today. The space age converted planetary magnetism from an Earth-bound subject into a comparative one: Mariner 10 discovered Mercury's field in 1974, Pioneer and Voyager mapped the enormous magnetospheres of the gas giants, Galileo detected Ganymede's field in 1996, and Mars Global Surveyor found the fossil crustal magnetisation of Mars in 1999. David Stevenson's 2003 synthesis consolidated these data into a unified account of why some planets are magnetised and others are not, framing the open questions that the BepiColombo, JUICE, and Europa Clipper missions now address.

Bibliography Master

Dziewonski, A. M. and Anderson, D. L. (1981). "Preliminary reference Earth model." Physics of the Earth and Planetary Interiors, 25, 297-356. The canonical seismically constrained radial model of Earth's density, pressure, and wave velocities, still the reference standard for interior structure.
Stevenson, D. J. (2003). "Planetary magnetic fields." Earth and Planetary Science Letters, 208, 1-11. The definitive comparative synthesis of why the planets do or do not sustain dynamos, framing the roles of core size, convection, and rotation.
de Pater, I. and Lissauer, J. J. (2015). Planetary Sciences, 2nd ed. Cambridge University Press. Comprehensive graduate text covering interiors, magnetic fields, atmospheres, and small bodies; chapters 1-6 develop the structure and dynamo material of this unit.
Turcotte, D. L. and Schubert, G. (2014). Geodynamics, 3rd ed. Cambridge University Press. The standard reference for mantle convection, plate tectonics, and the thermal evolution of terrestrial planets.
Vine, F. J. and Matthews, D. H. (1963). "Magnetic anomalies over oceanic ridges." Nature, 199, 947-949. The paper that linked marine magnetic stripes to sea-floor spreading and geomagnetic reversals, founding modern plate tectonics.
Lehmann, I. (1936). "'P'." Publications du Bureau Central Seismologique International, Serie A, 14, 87-115. Discovery of the Earth's solid inner core from the analysis of seismic wave arrivals.
Cowling, T. G. (1934). "The magnetic field of sunspots." Monthly Notices of the Royal Astronomical Society, 94, 39-48. The anti-dynamo theorem proving that no axisymmetric flow can maintain an axisymmetric field, which set the agenda for three-dimensional dynamo theory.
Parker, E. N. (1955). "Hydromagnetic dynamo models." Astrophysical Journal, 122, 293-314. Introduces the alpha-effect and the mean-field dynamo framework that underpins modern stellar and planetary dynamo theory.
Acuna, M. H. et al. (1999). "Global distribution of crustal magnetization discovered by the Mars Global Surveyor MAG/ER experiment." Science, 284, 790-793. Discovery of intense remanent crustal magnetisation on Mars, evidence for an ancient dynamo that has since ceased.
Ness, N. F., Behannon, K. W., Lepping, R. P., Whang, Y. C., and Schatten, K. H. (1974). "Magnetic field observations near Mercury: preliminary results from Mariner 10." Science, 185, 151-160. First detection of Mercury's weak planetary magnetic field.
Gurnett, D. A., Kurth, W. S., Roux, A., Bolton, S. J., and Kennel, C. F. (1996). "Galileo plasma wave observations during Ganymede encounters." Nature, 384, 535-537. Evidence from the Galileo mission for an internally generated magnetic field at Ganymede, the only moon known to possess one.
Porco, C. C. et al. (2006). "Cassini observes the active south pole of Enceladus." Science, 311, 1393-1401. Discovery of cryovolcanic plumes from the south-polar fractures of Enceladus, driven by tidal dissipation.
Wieczorek, M. A. et al. (2013). "The crust of the Moon as seen by GRAIL." Science, 339, 671-675. GRAIL gravity data yielding the global crustal thickness map of the Moon and its bulk composition.
Nimmo, F. and Stevenson, D. J. (2000). "Influence of early plate tectonics on the thermal evolution and magnetic field of Earth." Earth and Planetary Science Letters, 183, 367-377. Analysis linking the presence or absence of plate tectonics to core heat flow and the sustainability of a planetary dynamo.