10.3 Giant Planet Atmospheres

The gas giants (Jupiter, Saturn) and ice giants (Uranus, Neptune) possess deep hydrogen-helium atmospheres with no solid surface. These worlds exhibit spectacular banded cloud systems, powerful zonal jets, long-lived vortex storms, and significant internal heat sources that drive atmospheric dynamics in fundamentally different regimes from terrestrial planets.

Jupiter's Atmosphere: Composition & Cloud Layers

Jupiter is composed primarily of H₂ (~89.8%) and He (~10.2%) by volume, with trace amounts of CH₄, NH₃, H₂O, and disequilibrium species (PH₃, GeH₄, AsH₃). The visible atmosphere is organized into alternating bright zones and dark belts, bounded by zonal jets.

Layer 1 (Uppermost): Ammonia Ice (NH₃)

Cloud base at ~0.5-0.7 bar, T ~ 130-150 K. White ammonia ice crystals form the visible cloud tops. These define the bright zones (upwelling, NH₃-rich) and dark belts (subsidence, NH₃-depleted). Optical depth $\tau \sim 5$-20 depending on latitude.

Layer 2 (Middle): Ammonium Hydrosulfide (NH₄SH)

Cloud base at ~2-3 bar, T ~ 200-230 K. Formed by the reaction $\text{NH}_3 + \text{H}_2\text{S} \to \text{NH}_4\text{SH}$. This reddish-brown cloud layer may contribute to Jupiter's characteristic coloration. Less directly observed than the NH₃ layer.

Layer 3 (Deepest): Water (H₂O)

Cloud base at ~5-7 bar, T ~ 270-300 K. The most massive cloud layer. Drives convective thunderstorms with lightning detected by Juno. The Galileo probe entered an anomalously dry hot spot, but Juno MWR observations show water abundance increasing with depth to at least 3x solar below 20 bar.

Cloud condensation follows the Clausius-Clapeyron equation:

$$\frac{d\ln p_{\text{sat}}}{dT} = \frac{L}{RT^2} \quad \Rightarrow \quad p_{\text{sat}}(T) = p_0 \exp\left[-\frac{L}{R}\left(\frac{1}{T} - \frac{1}{T_0}\right)\right]$$

Cloud base forms where partial pressure of the condensable equals saturation vapor pressure

Jovian Circulation: Geostrophic Turbulence & Alternating Jets

Jupiter's ~30 alternating zonal jets (speeds up to 180 m/s) arise from geostrophic turbulence on a rapidly rotating planet. In 2D turbulence, energy cascades to larger scales (inverse cascade), while the $\beta$-effect organizes this energy into zonal bands at the Rhines scale:

$$L_{\beta} = \pi\sqrt{\frac{2U}{\beta}}, \quad \beta = \frac{2\Omega\cos\phi}{R_p}$$

For Jupiter: $\beta \approx 4.5 \times 10^{-12}$ m⁻¹s⁻¹, U ~ 50 m/s gives Lβ ~ 15,000 km

Thermal wind equation for giant planets:

$$\frac{\partial u}{\partial z} = -\frac{g}{fT}\frac{\partial T}{\partial y}$$

Meridional temperature gradients couple to vertical wind shear; warm zones are anticyclonic, cool belts cyclonic

Zones and Belts

Zones (light bands): high cloud tops, upwelling, anticyclonic vorticity. Belts (dark bands): lower cloud tops, subsidence, cyclonic vorticity. The equatorial zone is strongly prograde (~100 m/s). Juno gravity data revealed jets extend ~3000 km deep before rigid-body rotation takes over.

Shallow vs. Deep Models

Two paradigms: (1) "Shallow" -- jets from inverse cascade in a thin weather layer at the Rhines scale; (2) "Deep" -- convective columns parallel to the rotation axis (Taylor-Proudman theorem). Juno confirmed jets reach ~3000 km, suggesting both mechanisms play a role.

Great Red Spot: Anticyclonic Vortex Dynamics

The GRS is the largest storm in the solar system -- an anticyclonic vortex at 22 deg S observed continuously for over 350 years. Its persistence is explained by 2D vortex dynamics:

~16,000 km

East-west extent (shrinking)

~12,000 km

North-south extent

~190 m/s

Peak wind speed at edges

2D barotropic vorticity equation (Euler equation on sphere):

$$\frac{\partial \zeta}{\partial t} + J(\psi, \zeta + f) = F - D, \quad \zeta = \nabla^2\psi$$

J = Jacobian advection, F = forcing, D = dissipation, f = Coriolis parameter

Why Anticyclones Persist

In 2D turbulence, energy cascades upward (inverse cascade) while enstrophy cascades downward. Large anticyclones are preferentially stable because: (1) they sit in anticyclonic shear zones between jets; (2) they merge with smaller same-sign vortices (vortex merger); (3) anticyclones on a $\beta$-plane radiate Rossby waves less efficiently than cyclones, making them more robust. The GRS continuously feeds on smaller vortices drifting into it.

Jupiter's Internal Heat: Kelvin-Helmholtz Contraction

Jupiter emits 1.67 times as much thermal radiation as it absorbs from the Sun. This excess comes from slow gravitational contraction (Kelvin-Helmholtz mechanism) -- a remnant of formation heat:

$$F_{\text{int}} = \sigma T_{\text{eff}}^4 - \frac{S_0(1-A)}{4} = 13.6 - 8.1 = 5.5 \text{ W/m}^2$$

T_eff = 124.4 K, T_eq = 110.0 K; total internal luminosity ~ 3.5 x 10¹⁷ W

Kelvin-Helmholtz cooling timescale:

$$\tau_{\text{KH}} = \frac{GM^2}{RL} \approx 10^{10} \text{ yr}$$

Jupiter is still contracting at ~2 cm/year, converting gravitational PE to thermal radiation

Adiabatic Interior Profile

The deep interior follows a dry adiabat with lapse rate $\Gamma = g/c_p \approx 24.8/14300 \approx 1.7$ K/km for the H₂/He mixture. Ortho-para hydrogen conversion modifies this in the 50-150 K range, releasing latent heat (1062 J/mol) that reduces the effective lapse rate by up to 40%.

Saturn: Helium Rain & Hexagonal Polar Vortex

Saturn has a similar H₂/He composition but with severely depleted helium (He ~ 3.25% vs. solar ~13.6%). This depletion is explained by helium immiscibility in metallic hydrogen at ~1-2 Mbar, causing "helium rain" that releases gravitational energy and accounts for Saturn's excess heat (emitted/absorbed ratio = 1.78).

Helium Rain Process

At pressures of 1-2 Mbar and T < 10,000 K, He becomes immiscible in metallic H. Helium-rich droplets form and sink, releasing gravitational PE as heat. This contributes ~50% of Saturn's excess luminosity beyond simple K-H contraction. Mean density is only 687 kg/m³ (less than water).

Hexagonal Polar Vortex

Saturn's north pole features a persistent hexagonal jet pattern (wavenumber-6) at ~78 deg N, discovered by Voyager (1981) and still present in Cassini images. This is a standing Rossby wave on the polar jet, with wavelength matching the Rossby deformation radius: $L_d = NH/f$. The hexagon spans ~30,000 km and extends at least 100 km deep.

Saturn's equatorial jet reaches ~450 m/s (far broader and faster than Jupiter's ~100 m/s equatorial jet). Great White Spots -- massive convective outbursts -- erupt roughly once per Saturn year (29.5 Earth years).

Uranus: Extreme Obliquity & Minimal Internal Heat

Uranus (H₂ ~82.5%, He ~15.2%, CH₄ ~2.3%) has an extreme obliquity of 97.8 deg, so each pole alternately faces the Sun for ~42 years. Despite this extreme seasonal asymmetry, the atmosphere maintains zonal banding.

Blue Color: Methane Absorption

Uranus's blue-green color arises from methane (CH₄) in the upper atmosphere, which absorbs red wavelengths ($\lambda > 600$ nm) while scattering blue light. The ~2.3% methane abundance (compared to Jupiter's 0.18%) creates much stronger absorption. The interior likely contains an "icy" mantle of H₂O, NH₃, and CH₄ under extreme pressure.

Anomalous Low Internal Heat

Uranus has $F_{\text{int}} < 0.042$ W/m² (emitted/absorbed ratio ~ 1.06) -- anomalously low compared to all other giants. Hypotheses: (1) a compositional gradient inhibits convection, trapping primordial heat; (2) a giant impact during formation caused rapid early cooling; (3) the symmetric circulation (pole-to-pole) efficiently redistributes heat.

Neptune: Strongest Winds & Great Dark Spot

Neptune (H₂ ~80%, He ~19%, CH₄ ~1.5%) hosts the fastest winds in the solar system despite receiving only 1/900th of Earth's solar flux. Its equatorial winds reach ~580 m/s retrograde (2100 km/h), which is marginally supersonic:

Sound speed at cloud level:

$$c_s = \sqrt{\frac{\gamma k_B T}{\mu m_H}} = \sqrt{\frac{1.4 \times 1.381 \times 10^{-23} \times 70}{2 \times 1.67 \times 10^{-27}}} \approx 570 \text{ m/s}$$

Significant Internal Heat

Neptune emits 2.6x absorbed solar radiation, providing a strong internal energy source that drives vigorous dynamics despite the planet's great distance from the Sun. This internal heat is consistent with ongoing K-H contraction and possibly ongoing differentiation of interior ices.

Great Dark Spot

Voyager 2 (1989) observed an anticyclonic vortex similar to Jupiter's GRS. However, it was transient -- gone by 1994 when Hubble looked. Companion bright clouds of methane ice form at vortex edges where air is forced upward. New dark spots have appeared since, indicating dynamic vortex generation.

Metallic Hydrogen & Magnetic Field Generation

At pressures above ~2 Mbar and temperatures above ~2000 K, molecular H₂ dissociates and electrons become delocalized, forming metallic liquid hydrogen. The transition is now thought to be gradual rather than a sharp phase boundary:

$$P_{\text{transition}} \sim 1\text{-}2 \text{ Mbar} \quad (10^{11}\text{-}2\times10^{11} \text{ Pa})$$

Occurs at ~15,000 km depth in Jupiter, ~30,000 km in Saturn

This metallic hydrogen layer is the largest electrical conductor in the solar system and generates Jupiter's powerful magnetic field (~14x Earth's equatorial field) through dynamo action. Juno revealed the magnetic field is far more complex than a simple dipole, with localized anomalies near the surface. The ice giants have very different, off-center tilted dipole fields generated in their ionic water layers.

Fortran: 1D Radiative-Convective Model for H₂/He Atmosphere

A 1D radiative-convective equilibrium model for a hydrogen-helium atmosphere with CH₄ and NH₃ absorption:

! giant_planet_rc.f90
! 1D Radiative-Convective model for H2/He atmosphere
! with methane and ammonia gray absorption
!
! Compile: gfortran -O2 -o giant_rc giant_planet_rc.f90
! Run:     ./giant_rc

program giant_planet_rc
  implicit none

  integer, parameter :: dp = selected_real_kind(15)
  integer, parameter :: nlev = 80           ! number of levels
  integer, parameter :: niter = 50000       ! iteration steps

  ! Physical constants
  real(dp), parameter :: sigma = 5.67d-8    ! Stefan-Boltzmann
  real(dp), parameter :: R_gas = 8.314d0    ! J/(mol K)
  real(dp), parameter :: g_jup = 24.79d0    ! Jupiter gravity (m/s^2)
  real(dp), parameter :: cp = 14300.0d0     ! J/(kg K) for H2/He
  real(dp), parameter :: mu = 2.22d-3       ! mean molecular weight (kg/mol)

  ! Model parameters
  real(dp), parameter :: p_top = 1.0d3      ! top pressure (Pa = 0.01 bar)
  real(dp), parameter :: p_bot = 1.0d7      ! bottom pressure (Pa = 100 bar)
  real(dp), parameter :: F_int = 5.5d0      ! internal heat flux (W/m^2)
  real(dp), parameter :: F_sol = 8.1d0      ! absorbed solar flux (W/m^2)
  real(dp), parameter :: tau_ir0 = 5.0d0    ! IR optical depth at p_bot
  real(dp), parameter :: kappa_sw = 0.002d0 ! SW absorption coefficient

  ! Arrays
  real(dp) :: p(nlev), T(nlev), tau_ir(nlev)
  real(dp) :: F_up(nlev), F_dn(nlev), F_net(nlev)
  real(dp) :: dT_rad(nlev), lapse_rate
  real(dp) :: dp_lev, T_skin, dt_relax
  integer  :: k, iter

  ! Set up pressure grid (log-spaced)
  do k = 1, nlev
    p(k) = p_top * (p_bot / p_top) ** (dble(k-1) / dble(nlev-1))
  end do

  ! IR optical depth (proportional to pressure for well-mixed absorber)
  do k = 1, nlev
    tau_ir(k) = tau_ir0 * (p(k) / p_bot)
  end do

  ! Initial temperature: adiabat from 165 K at 1 bar
  do k = 1, nlev
    T(k) = 165.0d0 * (p(k) / 1.0d5) ** (R_gas / (mu * cp))
    T(k) = max(T(k), 50.0d0)
  end do

  ! Relaxation parameter
  dt_relax = 0.005d0

  ! === Iterative radiative-convective adjustment ===
  do iter = 1, niter

    ! --- Radiative fluxes (gray two-stream) ---
    ! Upward flux: F_up = sigma*T^4 * (1 - exp(-tau))
    ! Simplified: use Eddington approximation
    do k = 1, nlev
      F_up(k) = sigma * T(k)**4
      F_dn(k) = 0.5d0 * sigma * T(k)**4 * (1.0d0 - exp(-tau_ir(k)))
    end do

    ! Add internal heat flux at bottom
    F_up(nlev) = F_up(nlev) + F_int

    ! Add absorbed solar (deposited with depth)
    do k = 1, nlev
      F_dn(k) = F_dn(k) + F_sol * exp(-kappa_sw * p(k) / g_jup)
    end do

    ! Net flux and heating rate
    do k = 2, nlev - 1
      dp_lev = p(k+1) - p(k-1)
      F_net(k) = F_up(k) - F_dn(k)
      dT_rad(k) = -g_jup / cp * (F_net(k+1) - F_net(k-1)) / dp_lev
      T(k) = T(k) + dt_relax * dT_rad(k)
      T(k) = max(T(k), 40.0d0)
    end do

    ! --- Convective adjustment ---
    ! Enforce adiabatic lapse rate where radiative profile is unstable
    lapse_rate = g_jup / cp  ! K/m -> need K per pressure level
    do k = 2, nlev
      ! dT/dp for adiabat: (R T)/(mu cp p)
      if (T(k) - T(k-1) > lapse_rate * R_gas * T(k) / &
          (mu * cp * g_jup) * (p(k) - p(k-1))) then
        ! Unstable: adjust to adiabat
        T(k) = T(k-1) + R_gas * T(k-1) / (mu * cp) * &
               log(p(k) / p(k-1))
      end if
    end do

    ! Print progress
    if (mod(iter, 10000) == 0) then
      write(*,'(A,I7,A,F8.2,A,F8.2,A,F8.2)') &
        ' iter=', iter, &
        '  T(top)=', T(1), &
        '  T(1bar)=', T(nlev/2), &
        '  T(bot)=', T(nlev)
    end if
  end do

  ! === Output final profile ===
  write(*,'(A)') ''
  write(*,'(A)') '=== Giant Planet Radiative-Convective Profile ==='
  write(*,'(A)') '  P(bar)      T(K)      tau_IR'
  write(*,'(A)') '  ------      ----      ------'
  do k = 1, nlev, 4
    write(*,'(F10.4, F10.2, F10.4)') p(k)/1.0d5, T(k), tau_ir(k)
  end do

  write(*,'(A)') ''
  write(*,'(A,F8.2,A)') ' Skin temperature: ', T(1), ' K'
  write(*,'(A,F8.2,A)') ' 1-bar temperature: ', T(nlev/2), ' K'
  write(*,'(A,F8.2,A)') ' Deep temperature: ', T(nlev), ' K'

end program giant_planet_rc

Interactive Simulation: Jupiter & Saturn Atmospheric Structure

Python

Plots temperature-pressure profiles for Jupiter, Saturn, Uranus, and Neptune. Shows cloud condensation levels for NH3 and H2O, compares scale heights and tropopause locations, and uses adiabatic lapse rates with different compositions.

giant_planet_atmosphere.py104 lines

import numpy as np
import matplotlib.pyplot as plt

# Jupiter & Saturn Atmospheric Structure: T-P profiles and cloud layers
R = 8.314  # J/(mol K)

planets = {
    'Jupiter': {'T1bar': 165, 'g': 24.79, 'mu': 2.22e-3, 'cp': 14300,
                'color': '#f97316', 'x_NH3': 3e-4, 'x_H2O': 1e-3},
    'Saturn':  {'T1bar': 134, 'g': 10.44, 'mu': 2.14e-3, 'cp': 14500,
                'color': '#eab308', 'x_NH3': 2e-4, 'x_H2O': 8e-4},
    'Uranus':  {'T1bar': 76,  'g': 8.87,  'mu': 2.64e-3, 'cp': 13000,
                'color': '#06b6d4', 'x_NH3': 1e-5, 'x_H2O': 5e-4},
    'Neptune': {'T1bar': 72,  'g': 11.15, 'mu': 2.53e-3, 'cp': 13200,
                'color': '#3b82f6', 'x_NH3': 1e-5, 'x_H2O': 5e-4},
}

def p_sat_NH3(T):
    return np.exp(27.92 - 3754.0 / np.clip(T, 50, 500))

def p_sat_H2O(T):
    return 101325 * np.exp(-51058/R * (1/np.clip(T, 100, 700) - 1/373.15))

def p_sat_NH4SH(T):
    return np.exp(34.0 - 10000.0 / np.clip(T, 100, 500))

p_bar = np.logspace(-2, 2, 2000)
p_Pa = p_bar * 1e5

fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 7))
fig.patch.set_facecolor('#0f172a')
for ax in [ax1, ax2]:
    ax.set_facecolor('#1e293b')
    ax.tick_params(colors='#cbd5e1')
    ax.xaxis.label.set_color('#cbd5e1')
    ax.yaxis.label.set_color('#cbd5e1')
    ax.title.set_color('#cbd5e1')
    for spine in ax.spines.values():
        spine.set_color('#334155')
    ax.grid(True, color='#334155', alpha=0.4)

# Panel 1: T-P profiles with cloud condensation
cloud_info = []
for name, par in planets.items():
    kappa = R / (par['mu'] * par['cp'] / par['mu'])
    # Use correct adiabatic exponent: R/(mu*cp) but mu is already in kg/mol
    kappa = R / (par['cp'] * par['mu'])
    T_prof = par['T1bar'] * (p_bar / 1.0) ** kappa
    ax1.semilogy(T_prof, p_bar, color=par['color'], lw=2.5, label=name)

# NH3 cloud base
    p_nh3_partial = par['x_NH3'] * p_Pa
    p_nh3_sat = p_sat_NH3(T_prof)
    idx = np.where(p_nh3_partial >= p_nh3_sat)[0]
    if len(idx) > 0:
        i = idx[0]
        ax1.plot(T_prof[i], p_bar[i], 's', color=par['color'], ms=10, zorder=5)
        cloud_info.append(f"{name} NH3 cloud: {p_bar[i]:.2f} bar, {T_prof[i]:.0f} K")

# H2O cloud base
    p_h2o_partial = par['x_H2O'] * p_Pa
    p_h2o_sat = p_sat_H2O(T_prof)
    idx = np.where(p_h2o_partial >= p_h2o_sat)[0]
    if len(idx) > 0:
        i = idx[0]
        ax1.plot(T_prof[i], p_bar[i], 'o', color=par['color'], ms=10, zorder=5)
        cloud_info.append(f"{name} H2O cloud: {p_bar[i]:.1f} bar, {T_prof[i]:.0f} K")

ax1.set_ylim(100, 0.01)
ax1.set_xlabel('Temperature (K)')
ax1.set_ylabel('Pressure (bar)')
ax1.set_title('T-P Profiles & Cloud Condensation')
ax1.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#cbd5e1')
ax1.text(80, 0.02, 'squares = NH3 clouds\ncircles = H2O clouds', color='#94a3b8', fontsize=8)

# Panel 2: Scale heights and energy balance
names = list(planets.keys())
T_eff = [124.4, 95.0, 59.1, 59.3]
T_eq  = [110.0, 81.3, 58.2, 46.6]
ratio = [1.67, 1.78, 1.06, 2.61]
colors = [planets[n]['color'] for n in names]
H_vals = [R * planets[n]['T1bar'] / (planets[n]['mu'] * planets[n]['g']) / 1000 for n in names]

x_pos = np.arange(4)
bars1 = ax2.bar(x_pos - 0.15, T_eff, 0.3, color=colors, alpha=0.9, label='T_eff')
bars2 = ax2.bar(x_pos + 0.15, T_eq, 0.3, color=colors, alpha=0.4, label='T_eq')
for i, (r, h) in enumerate(zip(ratio, H_vals)):
    ax2.text(i, max(T_eff[i], T_eq[i]) + 4, f'{r:.2f}x\nH={h:.0f}km',
             ha='center', fontsize=9, color=colors[i])
ax2.set_xticks(x_pos)
ax2.set_xticklabels(names, color='#cbd5e1')
ax2.set_ylabel('Temperature (K)')
ax2.set_title('Effective vs Equilibrium Temperature\n(ratio = emitted/absorbed, H = scale height)')
ax2.legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#cbd5e1')

plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor=fig.get_facecolor())

print("=== Giant Planet Atmospheric Structure ===")
for info in cloud_info:
    print(info)
for n in names:
    H = R * planets[n]['T1bar'] / (planets[n]['mu'] * planets[n]['g']) / 1000
    print(f"{n}: T(1bar)={planets[n]['T1bar']} K, Scale height H={H:.0f} km")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← 10.2 Venus 10.4 Exoplanets →