10.1 Mars Atmosphere

Mars possesses a thin, cold, CO₂-dominated atmosphere with a surface pressure only 0.6% that of Earth. Despite its tenuous nature, the Martian atmosphere drives global dust storms, sculpts polar ice caps through seasonal CO₂ sublimation cycles, and preserves evidence of a warmer, wetter past. Understanding Mars's atmosphere is central to planetary science and future human exploration.

Atmospheric Composition and Basic Properties

The Martian atmosphere is overwhelmingly carbon dioxide, with minor amounts of nitrogen, argon, and trace gases. The composition was precisely measured by the Mars Science Laboratory (Curiosity rover) and confirmed by earlier Viking lander mass spectrometry.

95.32%

CO₂ (carbon dioxide)

2.7%

N₂ (nitrogen)

1.6%

Ar (argon)

Trace Species

O₂: 0.174% (produced by UV photolysis of CO₂)
CO: 0.0747% (product of CO₂ photodissociation)
H₂O: 0.03% variable (seasonal, latitudinal dependence)
NO, O₃, H₂O₂: trace amounts critical for photochemistry
CH₄: debated detection at ~0.4 ppb (Curiosity TLS)

Key Physical Parameters

Mean surface pressure: ~610 Pa (6.1 mbar)
Mean surface temperature: ~210 K (-63°C)
Temperature range: 130 K to 308 K
Mean molecular weight: 43.34 g/mol
Surface gravity: 3.72 m/s²
Specific gas constant: R/M = 191.8 J/(kg·K)

The surface pressure on Mars varies seasonally by approximately 25% due to the condensation and sublimation of CO₂ at the polar caps. During southern hemisphere winter, the south polar cap absorbs a significant fraction of atmospheric CO₂, lowering global pressure from ~680 Pa to ~510 Pa. This seasonal "breathing" of the atmosphere has no analogue on Earth.

Scale Height and Hydrostatic Equilibrium

The vertical structure of any planetary atmosphere is governed by hydrostatic balance, where the upward pressure gradient force balances the downward gravitational force. On Mars, the lower gravity and higher mean molecular weight combine to produce a scale height somewhat larger than one might initially expect.

Hydrostatic equilibrium:

$$\frac{dp}{dz} = -\rho g = -\frac{p M g}{R T}$$

Integration yields the barometric formula with pressure scale height H:

Pressure scale height:

$$H = \frac{R T}{M g} = \frac{8.314 \times 210}{0.04334 \times 3.72} \approx 10.8 \text{ km}$$

Barometric formula:

$$p(z) = p_0 \exp\left(-\frac{z}{H}\right) = p_0 \exp\left(-\frac{M g z}{R T}\right)$$

Comparison with Earth

Earth: H ≈ 8.5 km (T=288K, M=0.029 kg/mol, g=9.81 m/s²). Mars has a larger scale height despite its colder temperature because the lower gravity (3.72 vs 9.81 m/s²) more than compensates for the heavier mean molecular weight (43.3 vs 29 g/mol).

Density Profile

For an isothermal atmosphere, the density also decreases exponentially: $\rho(z) = \rho_0 \exp(-z/H)$. Mars surface density is ~0.020 kg/m³, about 1.6% of Earth's. This extremely thin air limits aerodynamic braking during EDL (entry, descent, landing).

Greenhouse Effect on Mars

Despite its CO₂-rich atmosphere, Mars has a surprisingly weak greenhouse effect due to the low atmospheric mass and pressure. The effective temperature and actual surface temperature differ by only about 5 K.

Effective radiative temperature (no atmosphere):

$$T_{\text{eff}} = \left[\frac{S(1-A)}{4\sigma}\right]^{1/4} = \left[\frac{589 \times (1-0.25)}{4 \times 5.67 \times 10^{-8}}\right]^{1/4} \approx 210 \text{ K}$$

Greenhouse warming comparison:

$$\Delta T_{\text{greenhouse}} = T_{\text{surface}} - T_{\text{eff}}$$

Mars: $\Delta T \approx 5$ K | Earth: $\Delta T \approx 33$ K | Venus: $\Delta T \approx 510$ K

The weak Martian greenhouse arises because: (1) low total atmospheric mass means low total opacity despite high CO₂ mixing ratio; (2) pressure broadening of absorption lines is minimal at 6 mbar; (3) the near-surface atmosphere is optically thin in the thermal infrared window regions. Adding more CO₂ alone cannot raise Mars surface temperature enough for liquid water -- dust and ice cloud feedbacks would be needed.

Dust Storms and Dust Cycle

Dust is the dominant aerosol on Mars and plays a central role in atmospheric dynamics, radiative transfer, and climate variability. Dust storms range from small local events to planet-encircling global events that can persist for months.

Local Storms

Diameter < 2000 km. Frequent, especially near polar caps during spring. Driven by surface heating and topographic channeling. Optical depth $\tau_{\text{dust}} \sim 0.5-2$.

Regional Storms

Cover significant latitude bands. Can merge from local storms. Duration: weeks. Optical depth $\tau_{\text{dust}} \sim 2-5$. Often triggered near perihelion (L_s ~ 250°).

Global Storms

Planet-encircling, obscure surface for months. Occur ~1 per 3 Mars years. Notable: 1971, 2001, 2007, 2018. Optical depth $\tau_{\text{dust}} > 5$, can reach 10+.

Dust lifting threshold -- friction velocity:

$$u_* > u_{*t} = A\sqrt{\frac{\rho_p g d_p}{\rho_a}\left(1 + \frac{6 \times 10^{-7}}{\rho_p g d_p^{2.5}}\right)}$$

where A ≈ 0.1, $\rho_p$ is particle density (~3000 kg/m³), $d_p$ is particle diameter (~few $\mu$m), $\rho_a$ is air density

Dust Devil Mechanics

Dust devils are convective vortices driven by superadiabatic surface heating. On Mars, dust devils can reach heights of 8-10 km (vs ~1 km on Earth) due to the deeper convective boundary layer. The pressure deficit at the vortex core is: $\Delta p \sim \rho v_\theta^2$ where $v_\theta$ is the tangential velocity (up to 30 m/s on Mars). Dust devils are responsible for significant dust lifting and were observed cleaning solar panels on the Spirit rover, extending its mission.

CO₂ Ice Caps and Seasonal Cycle

Mars has permanent water ice caps at both poles, overlaid by seasonal CO₂ ice (dry ice) that condenses directly from the atmosphere during winter. This process sequesters approximately 25% of the total atmospheric mass each year, making it the largest known seasonal mass transfer in the solar system.

CO₂ condensation temperature at Mars surface pressure:

$$T_{\text{cond}} = \frac{L_{\text{sub}}}{R \ln(p_{\text{sat}}/p)} \approx 148 \text{ K at } 610 \text{ Pa}$$

North Polar Cap

Permanent residual cap: water ice, ~1000 km diameter, 2 km thick. Seasonal CO₂ frost extends to ~65°N in winter. CO₂ ice thickness: up to 1-2 m. Sublimes completely by late spring.

South Polar Cap

Permanent residual cap: thin CO₂ ice layer over deep water ice (up to 3.7 km thick). Seasonal frost extends to ~50°S. Asymmetry due to Mars orbital eccentricity: southern winter is longer and colder (aphelion).

Seasonal mass budget (condensation rate):

$$\frac{dm}{dt} = \frac{F_{\text{IR}} - F_{\text{solar}}(1-A_{\text{ice}})}{L_{\text{sub,CO}_2}}$$

where $L_{\text{sub,CO}_2} = 5.9 \times 10^5$ J/kg and $A_{\text{ice}} \approx 0.6-0.8$ for fresh CO₂ frost

Atmospheric Escape and Evolution

Mars has lost the vast majority of its early atmosphere over ~4 billion years. Multiple escape mechanisms operate, with their relative importance varying over solar system history. The MAVEN (Mars Atmosphere and Volatile EvolutioN) mission has provided the most detailed measurements of present-day escape rates.

Escape velocity from Mars:

$$v_{\text{esc}} = \sqrt{\frac{2GM}{R}} = \sqrt{\frac{2 \times 6.674 \times 10^{-11} \times 6.417 \times 10^{23}}{3.390 \times 10^6}} \approx 5.03 \text{ km/s}$$

Jeans (thermal) escape flux for species i:

$$\Phi_J = \frac{n_c v_{\text{th}}}{2\sqrt{\pi}} (1 + \lambda_c) e^{-\lambda_c}, \quad \lambda_c = \frac{v_{\text{esc}}^2}{v_{\text{th}}^2} = \frac{GM m_i}{k_B T_c r_c}$$

where $v_{\text{th}} = \sqrt{2k_BT/m_i}$ is the most probable thermal speed, $n_c, T_c, r_c$ are evaluated at the exobase

Jeans (Thermal) Escape

Molecules in the high-velocity tail of the Maxwell-Boltzmann distribution exceed escape velocity. Significant for H and D atoms at the exobase (~200 km altitude, T ~ 200-350 K). The Jeans parameter $\lambda_c$ for atomic hydrogen on Mars is ~3-5, making thermal escape efficient. For heavier species (O, C, N), $\lambda_c > 50$, so Jeans escape is negligible.

Solar Wind Stripping

Mars lacks a global magnetic field (lost ~4 Gyr ago as the core cooled). Solar wind ions directly interact with the upper atmosphere, picking up and accelerating ionospheric ions (O⁺, CO₂⁺, O₂⁺). MAVEN measured current loss rate: ~100 g/s of atmosphere. During solar storms, rates can increase by 10-20x. Over 4 Gyr, estimated total loss: 0.5-1 bar equivalent.

Sputtering

Energetic ions (from solar wind or magnetospheric pickup) impact the upper atmosphere and eject neutral atoms with enough energy to escape. Particularly important for heavier species (O, C, N) that cannot escape thermally. Sputtering yield depends on incident ion energy and atmospheric composition.

Photochemical (Dissociative Recombination) Escape

Dissociative recombination of O₂⁺ produces two O atoms, each carrying ~2.5 eV of kinetic energy: $\text{O}_2^+ + e^- \rightarrow \text{O}^* + \text{O}^*$. If the O atoms are produced above the exobase with sufficient upward velocity, they escape. This is the dominant oxygen escape mechanism on Mars today (~50% of total O loss).

Evidence for a Warmer, Wetter Past

Multiple lines of geological, mineralogical, and isotopic evidence point to Mars having had a thicker atmosphere and liquid surface water in its first ~1 billion years (Noachian period, >3.7 Gyr ago).

Geomorphological Evidence

Dendritic valley networks (Noachian-age terrain)
Outflow channels from catastrophic floods (Hesperian)
Deltas in crater basins (Jezero crater, Perseverance landing site)
Possible shorelines of a northern ocean

Mineralogical Evidence

Phyllosilicates (clays): require sustained liquid water
Sulfates (jarosite, gypsum): evaporite deposits
Hematite "blueberries" at Meridiani Planum
Carbonates in ancient rock units

Isotopic Evidence (D/H ratio)

The deuterium-to-hydrogen ratio (D/H) in the Martian atmosphere is ~5-6x the terrestrial SMOW value. This enrichment indicates preferential loss of lighter H atoms relative to D through atmospheric escape. From mass balance: $\frac{(D/H)_{\text{now}}}{(D/H)_{\text{initial}}} = f^{(1-1/\alpha)}$ where f is the fraction remaining and $\alpha$ is the fractionation factor. The D/H enrichment implies Mars has lost >80% of its original water inventory, corresponding to a global equivalent layer of ~20-40 m.

MAVEN Mission Results

NASA's MAVEN spacecraft (arrived September 2014) has been dedicated to understanding Mars atmospheric loss. Key findings include:

Present-Day Loss Rates

Total atmospheric loss: ~100 g/s under quiet solar conditions. During solar storms and CME impacts, loss rates increase by factors of 10-20. Ion escape (O⁺, CO₂⁺, O₂⁺) dominates. Extrapolated over 4 Gyr: loss equivalent to ~0.5-1 bar of CO₂, consistent with a thicker early atmosphere.

Upper Atmosphere Structure

MAVEN mapped the thermosphere (100-250 km) and ionosphere in detail. Exobase temperature: 200-350 K (varying with solar cycle). Discovered highly variable "breathing" of the upper atmosphere in response to solar EUV variability. Identified discrete aurora patches linked to crustal magnetic field remnants.

Terraforming Considerations

Proposals to make Mars habitable require dramatically increasing atmospheric pressure and temperature. The fundamental challenge is the mass budget.

Required atmospheric mass for 1 bar surface pressure:

$$m_{\text{atm}} = \frac{p_{\text{target}} \times 4\pi R^2}{g} = \frac{10^5 \times 4\pi \times (3.39 \times 10^6)^2}{3.72} \approx 3.9 \times 10^{18} \text{ kg}$$

The current Mars atmosphere has mass ~2.5 x 10¹⁶ kg. Known CO₂ reservoirs on Mars (polar caps, regolith adsorption, carbonate rocks) could contribute at most ~0.1-0.2 bar even if fully released. Studies by Jakosky and Edwards (2018) concluded that terraforming Mars using currently accessible CO₂ is not feasible with present-day technology. Importing volatiles (e.g., redirecting comets) would require ~10¹⁷ kg of material -- equivalent to thousands of 10-km comets.

Fortran: Jeans Escape Flux Calculation

This Fortran program computes the Jeans thermal escape flux for several atmospheric species (H, D, He, O, C, N₂, CO₂) at the Martian exobase, demonstrating why only light species escape thermally.

program jeans_escape_mars
  implicit none

  ! === Constants ===
  double precision, parameter :: kB    = 1.381d-23   ! Boltzmann constant (J/K)
  double precision, parameter :: G_grav = 6.674d-11  ! gravitational constant
  double precision, parameter :: pi    = 3.14159265358979d0
  double precision, parameter :: amu   = 1.6605d-27  ! atomic mass unit (kg)

  ! === Mars parameters ===
  double precision, parameter :: M_mars = 6.417d23   ! Mars mass (kg)
  double precision, parameter :: R_mars = 3.390d6    ! Mars radius (m)
  double precision, parameter :: z_exo = 200.0d3     ! exobase altitude (m)
  double precision, parameter :: T_exo = 250.0d0     ! exobase temperature (K)
  double precision, parameter :: r_exo = R_mars + z_exo ! exobase radius

  ! === Species data (name, mass in amu, exobase density in m^-3) ===
  integer, parameter :: nspecies = 7
  character(len=4) :: names(nspecies)
  double precision :: mass_amu(nspecies), n_exo(nspecies)
  double precision :: m_kg, v_esc, v_th, lambda_c, phi_jeans
  integer :: i

  ! Initialize species
  names    = (/ 'H   ', 'D   ', 'He  ', 'O   ', 'C   ', 'N2  ', 'CO2 ' /)
  mass_amu = (/ 1.0d0, 2.0d0, 4.0d0, 16.0d0, 12.0d0, 28.0d0, 44.0d0 /)
  n_exo    = (/ 1.0d10, 1.0d7, 1.0d8, 1.0d12, 1.0d9, 1.0d10, 1.0d11 /)

  ! === Escape velocity at exobase ===
  v_esc = sqrt(2.0d0 * G_grav * M_mars / r_exo)

  write(*,'(A)') '================================================='
  write(*,'(A)') ' Jeans Escape from Mars Exobase'
  write(*,'(A,F6.0,A)') ' Exobase altitude: ', z_exo/1.0d3, ' km'
  write(*,'(A,F6.1,A)') ' Exobase temperature: ', T_exo, ' K'
  write(*,'(A,F8.2,A)') ' Escape velocity: ', v_esc, ' m/s'
  write(*,'(A)') '================================================='
  write(*,'(A)') ''
  write(*,'(A6, A12, A12, A10, A15)') &
       'Spec', 'v_th(m/s)', 'lambda_c', 'exp(-lam)', 'Phi(m-2s-1)'

  do i = 1, nspecies
    m_kg = mass_amu(i) * amu

    ! Most probable thermal speed
    v_th = sqrt(2.0d0 * kB * T_exo / m_kg)

    ! Jeans parameter (escape parameter)
    lambda_c = (v_esc / v_th)**2

    ! Jeans escape flux (particles per m^2 per s)
    ! Phi = n * v_th / (2*sqrt(pi)) * (1 + lambda_c) * exp(-lambda_c)
    if (lambda_c < 500.0d0) then
      phi_jeans = n_exo(i) * v_th / (2.0d0 * sqrt(pi)) &
                  * (1.0d0 + lambda_c) * exp(-lambda_c)
    else
      phi_jeans = 0.0d0  ! negligibly small
    end if

    write(*,'(A6, F12.1, F12.2, ES10.2, ES15.4)') &
         names(i), v_th, lambda_c, exp(-min(lambda_c,500.0d0)), phi_jeans
  end do

  write(*,'(A)') ''
  write(*,'(A)') 'Note: Only H and D have significant Jeans escape.'
  write(*,'(A)') 'Heavier species require non-thermal escape mechanisms.'

end program jeans_escape_mars

Interactive Simulation: Mars Atmospheric Profile

Python

Plots temperature, pressure, and density profiles for Mars atmosphere (0-120 km) compared with Earth on the same scale. Calculates scale height H = kT/(mg) for both planets and shows the CO2 sublimation curve marking polar cap regions.

mars_atmospheric_profile.py87 lines

import numpy as np
import matplotlib.pyplot as plt

# Mars Atmospheric Profile compared with Earth
R_gas = 8.314  # J/(mol K)

# Mars parameters
g_mars, M_mars, T0_mars, p0_mars = 3.72, 0.04334, 210.0, 610.0
H_mars = R_gas * T0_mars / (M_mars * g_mars)  # scale height

# Earth parameters
g_earth, M_earth, T0_earth, p0_earth = 9.81, 0.029, 288.15, 101325.0
H_earth = R_gas * T0_earth / (M_earth * g_earth)

z = np.linspace(0, 120000, 1000)  # 0-120 km
z_km = z / 1000

# Mars: piece-wise temperature profile
T_mars = np.where(z_km < 40, T0_mars - 1.5 * z_km,
         np.where(z_km < 80, T0_mars - 1.5*40 - 0.5*(z_km-40), 130.0))
T_mars = np.maximum(T_mars, 120)
p_mars = p0_mars * np.exp(-z / H_mars)
rho_mars = p_mars * M_mars / (R_gas * T_mars)

# Earth: troposphere lapse rate then isothermal stratosphere
T_earth = np.where(z_km < 11, T0_earth - 6.5*z_km,
          np.where(z_km < 50, 216.65, 216.65 + 1.5*(z_km-50)))
T_earth = np.where(z_km > 80, 216.65 + 1.5*30 - 3*(z_km-80), T_earth)
p_earth = p0_earth * np.exp(-z / H_earth)
rho_earth = p_earth * M_earth / (R_gas * T_earth)

# CO2 sublimation curve (approximate)
T_co2_sub = np.linspace(140, 195, 100)
p_co2_sub = 101325 * np.exp(-25000/R_gas * (1/T_co2_sub - 1/194.65))

fig, axes = plt.subplots(1, 3, figsize=(16, 7))
fig.patch.set_facecolor('#0f172a')
for ax in axes:
    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)

# Temperature profiles
axes[0].plot(T_mars, z_km, '#ef4444', lw=2.5, label='Mars')
axes[0].plot(T_earth, z_km, '#38bdf8', lw=2.5, label='Earth')
axes[0].axhspan(0, 5, alpha=0.1, color='#c084fc', label='Mars polar cap zone')
axes[0].axvline(148, color='#c084fc', ls='--', alpha=0.5, label='CO2 frost point')
axes[0].set_xlabel('Temperature (K)')
axes[0].set_ylabel('Altitude (km)')
axes[0].set_title('Temperature Profile')
axes[0].legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#cbd5e1', fontsize=8)

# Pressure profiles (log scale)
axes[1].semilogx(p_mars, z_km, '#ef4444', lw=2.5, label=f'Mars (H={H_mars/1000:.1f} km)')
axes[1].semilogx(p_earth, z_km, '#38bdf8', lw=2.5, label=f'Earth (H={H_earth/1000:.1f} km)')
axes[1].set_xlabel('Pressure (Pa)')
axes[1].set_title('Pressure Profile')
axes[1].legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#cbd5e1', fontsize=8)

# Density profiles (log scale)
axes[2].semilogx(rho_mars, z_km, '#ef4444', lw=2.5, label='Mars')
axes[2].semilogx(rho_earth, z_km, '#38bdf8', lw=2.5, label='Earth')
axes[2].set_xlabel('Density (kg/m^3)')
axes[2].set_title('Density Profile')
axes[2].legend(facecolor='#1e293b', edgecolor='#334155', labelcolor='#cbd5e1', fontsize=8)

for ax in axes:
    ax.set_ylim(0, 120)

plt.suptitle('Mars vs Earth Atmospheric Profiles', color='#cbd5e1', fontsize=14, y=1.02)
plt.tight_layout()
plt.savefig('output.png', dpi=150, bbox_inches='tight', facecolor=fig.get_facecolor())

rho0_mars = p0_mars * M_mars / (R_gas * T0_mars)
rho0_earth = p0_earth * M_earth / (R_gas * T0_earth)
print("=== Mars Atmospheric Profile ===")
print(f"Mars scale height: H = {H_mars/1000:.2f} km")
print(f"Earth scale height: H = {H_earth/1000:.2f} km")
print(f"Mars surface: T={T0_mars} K, P={p0_mars} Pa, rho={rho0_mars:.4f} kg/m^3")
print(f"Earth surface: T={T0_earth} K, P={p0_earth} Pa, rho={rho0_earth:.4f} kg/m^3")
print(f"Density ratio (Earth/Mars): {rho0_earth/rho0_mars:.0f}x")
print(f"CO2 frost temperature at Mars: ~148 K")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← Part 10 10.2 Venus Atmosphere →