Geophysics & Planetary Science

Geophysics and planetary science investigate the physical processes governing solid planets, their interiors, surfaces, atmospheres, and magnetic fields — from the iron core of the Earth to the hydrogen envelopes of gas giants and the exotic atmospheres of exoplanets orbiting distant stars. This discipline sits at the intersection of classical mechanics, thermodynamics, fluid dynamics, and electromagnetism, applied to some of the most complex natural systems known. It is also one of the most consequential branches of physics: understanding plate tectonics, climate dynamics, and planetary habitability bears directly on questions of practical urgency for civilization and of profound philosophical weight about life in the universe.

Earth’s Interior and Seismology

Our knowledge of Earth’s interior comes almost entirely from seismology — the study of elastic waves generated by earthquakes and artificial sources as they propagate through the planet. The fundamental wave equation for an elastic, isotropic medium yields two types of body waves: P-waves (primary, or compressional waves) with velocity $v_P = \sqrt{(\lambda + 2\mu)/\rho}$, and S-waves (secondary, or shear waves) with velocity $v_S = \sqrt{\mu/\rho}$, where $\lambda$ and $\mu$ are the Lame parameters and $\rho$ is the density. Since S-waves require a medium with shear rigidity ($\mu > 0$), they cannot propagate through liquids — a fact that proved decisive in establishing Earth’s internal structure.

In 1906, Richard Dixon Oldham used the absence of S-waves in a shadow zone to demonstrate the existence of a liquid outer core. In 1936, Inge Lehmann discovered the solid inner core by analyzing P-wave arrivals that had been refracted through a deeper boundary. These discoveries, together with systematic travel-time analysis, revealed Earth’s layered structure: a thin silicate crust (5-70 km thick), a viscous mantle extending to 2900 km depth, a liquid iron-nickel outer core, and a solid inner core with a radius of approximately 1220 km. The pressure at Earth’s center reaches roughly $360\;\text{GPa}$ and the temperature is estimated at $5000\text{-}6000\;\text{K}$.

Modern seismic tomography uses arrival-time data from thousands of earthquakes recorded at global seismic networks (such as the IRIS Global Seismographic Network) to construct three-dimensional maps of seismic velocity anomalies within the mantle. Regions of faster-than-average velocity correspond to cold, dense material (such as subducted oceanic slabs), while slower regions indicate hot upwelling material (such as mantle plumes). Full-waveform inversion techniques, which fit the entire seismogram rather than just arrival times, are revealing increasingly fine-grained structure, including the ultra-low-velocity zones at the core-mantle boundary that may represent partial melt or compositionally distinct material.

Plate Tectonics and Geodynamics

The theory of plate tectonics, synthesized in the 1960s from decades of observations by Alfred Wegener (continental drift, 1912), Harry Hess (seafloor spreading, 1962), and many others, describes the outer rigid layer of the Earth — the lithosphere — as a mosaic of approximately fifteen major tectonic plates that move relative to one another over the more ductile asthenosphere beneath. Plates are created at mid-ocean ridges, where upwelling mantle material generates new oceanic crust, and destroyed at subduction zones, where one plate dives beneath another into the deep mantle.

The driving forces of plate motion include ridge push (gravitational sliding of plates away from elevated ridges), slab pull (the negative buoyancy of cold, dense subducting slabs), and basal drag from underlying mantle convection. The relative importance of these forces remains debated, but slab pull is generally regarded as the dominant mechanism. Plate velocities, measured today by GPS at centimeter-per-year precision, range from about $1$ to $15\;\text{cm/yr}$.

The engine underlying plate tectonics is mantle convection — the slow, thermally driven circulation of the solid but ductile mantle. Whether convection occurs as a single layer spanning the entire mantle or in two distinct layers separated at the 660-km seismic discontinuity has been a central question in geodynamics. Seismic tomography images of subducted slabs penetrating into the lower mantle favor whole-mantle convection, though some slabs appear to stagnate at intermediate depths. The vigor of convection is characterized by the Rayleigh number:

$\text{Ra} = \frac{\alpha g \Delta T\,d^3}{\kappa\nu}$

where $\alpha$ is the thermal expansion coefficient, $g$ is gravitational acceleration, $\Delta T$ is the temperature contrast across the mantle, $d$ is the mantle depth, $\kappa$ is thermal diffusivity, and $\nu$ is kinematic viscosity. Earth’s mantle Rayleigh number ($\sim 10^7$) is well above the critical value for convection ($\sim 10^3$), ensuring vigorous flow.

Earthquakes occur when accumulated elastic strain along a fault is released suddenly. The Mohr-Coulomb failure criterion describes the conditions for fracture: failure occurs when the shear stress $\tau$ on a plane exceeds the frictional resistance, $\tau = c + \mu_f\,\sigma_n$, where $c$ is the cohesive strength, $\mu_f$ is the friction coefficient, and $\sigma_n$ is the normal stress. The seismic moment $M_0 = \mu A \bar{D}$ (where $\mu$ is the shear modulus, $A$ is the rupture area, and $\bar{D}$ is the average slip) quantifies the size of an earthquake and defines the moment magnitude scale $M_w = \frac{2}{3}\log_{10}M_0 - 10.7$, which has replaced the older Richter scale for large events.

Geomagnetism and the Dynamo

Earth possesses a global magnetic field that is approximately dipolar at the surface, with a current strength of about $3 \times 10^{-5}\;\text{T}$ at the equator. This field is generated by the geodynamo — convective motions of electrically conducting liquid iron in the outer core that amplify and sustain magnetic field lines through electromagnetic induction. The governing equation is the magnetic induction equation:

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

where $\mathbf{B}$ is the magnetic field, $\mathbf{v}$ is the fluid velocity, and $\eta = 1/(\mu_0\sigma)$ is the magnetic diffusivity ($\sigma$ being electrical conductivity). The first term on the right represents advection of field lines by the flow (the dynamo effect); the second represents ohmic diffusion. For the dynamo to operate, the magnetic Reynolds number $\text{Rm} = vL/\eta$ must exceed a critical value of order $10\text{-}100$.

One of the most striking features of Earth’s magnetic field is its history of polarity reversals — episodes in which the north and south magnetic poles swap, recorded in the magnetization of ocean-floor basalts and continental volcanic rocks. The paleomagnetic record reveals that reversals occur at irregular intervals averaging a few hundred thousand years, with the last major reversal (the Brunhes-Matuyama reversal) occurring approximately $780{,}000$ years ago. The magnetic polarity timescale, established by systematic paleomagnetic dating, provided key evidence for seafloor spreading and was instrumental in the acceptance of plate tectonics in the 1960s. Modern numerical dynamo simulations, solving the coupled magnetohydrodynamic equations on supercomputers, can reproduce many features of the observed field, including secular variation, dipole tilt, and spontaneous polarity reversals.

Atmospheric Physics and Climate

A planet’s atmosphere is the thin gaseous envelope held by gravity, and its structure is determined by the interplay of radiative transfer, thermodynamics, and fluid dynamics. Earth’s atmosphere is divided into layers defined by the temperature profile: the troposphere (surface to $\sim 12\;\text{km}$, where temperature decreases with altitude), the stratosphere (up to $\sim 50\;\text{km}$, warmed by ozone absorption of UV radiation), the mesosphere, and the thermosphere. The pressure at any altitude $z$ in hydrostatic equilibrium is:

$P(z) = P_0\,\exp\left(-\frac{z}{H}\right), \qquad H = \frac{k_B T}{mg}$

where $H$ is the scale height (about $8.5\;\text{km}$ for Earth), $m$ is the mean molecular mass, and $g$ is the gravitational acceleration.

The global energy balance is set by the competition between incoming solar radiation and outgoing infrared emission. Earth absorbs solar radiation at a rate $\pi R^2 (1 - A) S$, where $S \approx 1361\;\text{W/m}^2$ is the solar constant and $A \approx 0.30$ is the planetary albedo. In radiative equilibrium, this is balanced by thermal emission $4\pi R^2 \sigma T_e^4$, yielding an effective temperature $T_e \approx 255\;\text{K}$ — well below the actual surface temperature of $\sim 288\;\text{K}$. The $33\;\text{K}$ difference is the greenhouse effect, caused by atmospheric gases (primarily $\text{CO}_2$, $\text{H}_2\text{O}$, $\text{CH}_4$) that are transparent to incoming shortwave radiation but absorb and re-emit outgoing longwave radiation.

Atmospheric dynamics are governed by the Navier-Stokes equations on a rotating sphere. The Coriolis force deflects large-scale flows, establishing the pattern of Hadley, Ferrel, and polar cells that dominate the general circulation. Mid-latitude weather is driven by baroclinic instability — the growth of wave-like disturbances that extract energy from the meridional temperature gradient. Climate sensitivity — the equilibrium warming per doubling of $\text{CO}_2$ — is amplified by feedbacks including water vapor, ice-albedo, cloud, and lapse-rate effects; current estimates place the equilibrium climate sensitivity (ECS) at $2.5\text{-}4.0\;\text{K}$, with paleoclimate evidence from ice cores and ocean sediments constraining this range.

Planetary Formation and Solar System Bodies

The solar system formed approximately $4.57$ billion years ago from the gravitational collapse of a molecular cloud, producing a rotating protoplanetary disk of gas and dust around the young Sun. Close to the Sun, where temperatures exceeded the condensation points of volatile compounds, only refractory silicates and metals could condense, giving rise to the terrestrial planets (Mercury, Venus, Earth, Mars). Beyond the snow line (roughly $2.7\;\text{AU}$), water ice condensed, greatly increasing the available solid mass and enabling the formation of massive planetary cores that accreted thick hydrogen-helium envelopes to become the giant planets (Jupiter, Saturn, Uranus, Neptune).

Planet formation proceeds through a hierarchy of stages: microscopic dust grains coagulate into millimeter-sized aggregates, which grow through the streaming instability (a collective aerodynamic effect) into kilometer-sized planetesimals, which then merge through gravitational encounters into planetary embryos and ultimately planets. The formation of Jupiter’s core must have occurred within a few million years, before the gas disk dissipated. Planetary migration — the gravitational interaction of forming planets with the disk — reshuffles orbital architectures. The Nice model proposes that a late instability in the orbits of the giant planets triggered the Late Heavy Bombardment ($\sim 3.9$ billion years ago), scattering planetesimals throughout the inner solar system.

Each solar system body tells a distinct physical story. Venus exhibits a runaway greenhouse effect, with surface temperatures reaching $735\;\text{K}$ and a dense $\text{CO}_2$ atmosphere at $90\;\text{atm}$. Mars once had liquid water on its surface but lost most of its atmosphere due to the absence of a global magnetic field and the consequent stripping by the solar wind. The moons of the outer planets harbor some of the most intriguing environments in the solar system: Europa and Enceladus possess subsurface liquid water oceans maintained by tidal heating, making them prime targets in the search for extraterrestrial life.

Exoplanets and Habitability

The discovery of the first exoplanet orbiting a Sun-like star — 51 Pegasi b, detected by Michel Mayor and Didier Queloz in 1995 via the radial velocity method — inaugurated a new era in planetary science. As of 2026, more than $5{,}700$ exoplanets have been confirmed, primarily through the transit method (pioneered by the Kepler space telescope) and the radial velocity method. The transit of a planet across its host star produces a dip in the light curve with depth:

$\frac{\Delta F}{F} \approx \left(\frac{R_p}{R_\star}\right)^2$

where $R_p$ and $R_\star$ are the planetary and stellar radii. Combined with radial velocity measurements, which yield the planetary mass through the semi-amplitude $K = \frac{M_p\sin i}{(M_\star + M_p)^{2/3}}\left(\frac{2\pi G}{P}\right)^{1/3}$, one obtains the bulk density and hence constraints on composition.

The diversity of exoplanets far exceeds what the solar system suggested: hot Jupiters orbiting their stars in days, super-Earths with no solar system analogue, compact multi-planet systems in orbital resonance, and free-floating planets unbound to any star. The concept of the habitable zone — the range of orbital distances where a planet could maintain liquid water on its surface — provides a first-order framework for assessing potential habitability, though the actual habitability of any world depends on atmospheric composition, magnetic field strength, tidal locking, and many other factors.

Transmission spectroscopy — measuring the wavelength-dependent absorption of starlight passing through a transiting planet’s atmosphere — has detected molecular species including water vapor, carbon dioxide, methane, and sulfur dioxide in exoplanet atmospheres. The James Webb Space Telescope, launched in 2021, has dramatically advanced atmospheric characterization of both giant and rocky exoplanets. Detecting unambiguous biosignatures — atmospheric signatures that can only be explained by biological activity, such as simultaneous oxygen and methane in thermodynamic disequilibrium — remains one of the grand challenges of twenty-first-century science, connecting geophysics and planetary science to the most fundamental question of all: whether life exists beyond Earth.