10.2 Venus Atmosphere
Venus is Earth's "evil twin" -- nearly identical in size and mass, yet possessing a hellish atmosphere with surface temperatures of 737 K (464 C), pressures of 92 bar, and perpetual sulfuric acid cloud cover. Venus represents the most dramatic example of a runaway greenhouse effect in our solar system and serves as a critical cautionary tale for understanding climate tipping points on terrestrial planets.
Atmospheric Composition and Physical Properties
The Venusian atmosphere is overwhelmingly carbon dioxide with a small nitrogen component. What distinguishes it from Mars (also CO₂-dominated) is the sheer mass: Venus has roughly 90 times more atmosphere than Earth by mass, and nearly all of it is a potent greenhouse gas.
96.5%
CO₂ (carbon dioxide)
3.5%
N₂ (nitrogen)
92 bar
Surface pressure
Trace Species
- SO₂: 150 ppm (drives sulfuric acid cloud formation)
- H₂O: 20-30 ppm (extremely desiccated)
- CO: 17 ppm (product of CO₂ photodissociation)
- HCl: 0.1-0.6 ppm
- HF: 1-5 ppb
- OCS (carbonyl sulfide): 4 ppm below clouds
- Noble gases: ³⁶Ar, ²⁰Ne enriched relative to Earth
Key Physical Parameters
- Surface temperature: 737 K (464°C)
- Surface pressure: 9.2 MPa (92 bar)
- Surface density: ~65 kg/m³ (supercritical CO₂)
- Mean molecular weight: 43.45 g/mol
- Surface gravity: 8.87 m/s²
- Scale height (surface): H ≈ 15.9 km
- Total atmospheric mass: 4.8 x 10²⁰ kg (~93x Earth)
At the surface, the CO₂ exists in a supercritical state -- above both the critical temperature (304 K) and critical pressure (73.8 bar). This means the distinction between liquid and gas phases vanishes, and the atmosphere near the surface behaves more like a dense fluid than what we normally consider a "gas."
The Runaway Greenhouse Effect
Venus exemplifies the catastrophic endpoint of a positive feedback loop between surface temperature, atmospheric water vapor, and infrared opacity. The runaway greenhouse mechanism explains how a planet with initially Earth-like conditions can evolve to a hothouse state.
Effective temperature without greenhouse (using Bond albedo A = 0.77):
$$T_{\text{eff}} = \left[\frac{S_0(1-A)}{4\sigma}\right]^{1/4} = \left[\frac{2601 \times (1-0.77)}{4 \times 5.67 \times 10^{-8}}\right]^{1/4} \approx 227 \text{ K}$$
Actual greenhouse warming:
$$\Delta T_{\text{GH}} = T_{\text{surface}} - T_{\text{eff}} = 737 - 227 = 510 \text{ K}$$
The Positive Feedback Loop
1. Early Venus received ~40% more solar flux than Earth (closer to Sun). 2. Higher insolation increased surface temperature, evaporating more ocean water. 3. Water vapor is a powerful greenhouse gas, trapping more infrared radiation. 4. Higher temperature evaporated more water -- positive feedback. 5. Eventually all ocean water entered the atmosphere. 6. UV photodissociation split H₂O; hydrogen escaped to space, oxygen was consumed by surface oxidation. 7. Without liquid water, carbonate-silicate weathering cycle shut down. 8. Volcanic CO₂ accumulated with no removal mechanism, locking in extreme greenhouse.
Komabayashi-Ingersoll Limit
There exists a maximum outgoing longwave radiation (OLR) for a planet with a water-vapor-saturated atmosphere. This is known as the Komabayashi-Ingersoll or Simpson-Nakajima limit:
$$F_{\text{OLR,max}} \approx 282 \text{ W/m}^2$$
If absorbed solar radiation exceeds this limit, the planet cannot radiate energy fast enough to maintain equilibrium. Temperature rises without bound until all water is in the vapor phase. For Venus, $F_{\text{abs}} = S_0(1-A)/4 \approx 150$ W/m² (below the limit due to high albedo), but early Venus with lower albedo would have exceeded it.
Atmospheric Super-Rotation
One of the most puzzling features of Venus's atmosphere is its super-rotation: the atmosphere at cloud level rotates 60 times faster than the solid planet beneath it. This is fundamentally different from Earth, where atmospheric angular momentum is less than the solid body's.
Cloud-top zonal wind speed:
~100 m/s
4-day circumnavigation at cloud top
Surface rotation speed:
~1.6 m/s
243-day retrograde sidereal period
Atmospheric angular momentum excess:
$$L_{\text{atm}} = \int \rho u R \cos\phi \, dV \gg \Omega R^2 \int \rho \cos^2\phi \, dV$$
The atmosphere carries ~10³ times more angular momentum per unit mass than the surface rotation would supply
Proposed Mechanisms
The Gierasch-Rossow-Williams mechanism proposes that the super-rotation is maintained by: (1) meridional Hadley circulation transporting angular momentum equatorward and upward; (2) planetary-scale waves (thermal tides, Kelvin waves, Rossby waves) transporting angular momentum from mid-latitudes toward the equator. The thermal tide component arises because solar heating of the cloud layer creates a pressure asymmetry that drives retrograde flow. Recent GCM simulations have partially reproduced super-rotation, but the precise balance of mechanisms remains debated.
Vertical Wind Profile
Near the surface, winds are very slow (~0.5-2 m/s), measured by Venera landers. Winds increase roughly linearly with altitude through the cloud layer, reaching ~100 m/s at 65-70 km. Above the clouds, winds decrease. The vertical shear is: $du/dz \approx 100/70000 \approx 1.4 \times 10^{-3}$ s⁻¹. Despite slow surface winds, the dense atmosphere exerts significant dynamic pressure: $q = \frac{1}{2}\rho v^2 \approx \frac{1}{2}(65)(1)^2 \approx 32$ Pa, comparable to a gentle Earth breeze of ~7 m/s.
Cloud Layers and Sulfuric Acid Cycle
Venus is permanently shrouded in thick clouds that extend from roughly 48 to 70 km altitude. These clouds are composed primarily of concentrated sulfuric acid (H₂SO₄) droplets, not water, and play a dominant role in the planet's radiative and chemical balance.
Upper Cloud (58-70 km)
H₂SO₄ droplets (~75% concentration by weight), particle radii ~1 μm. Temperature: 200-260 K. This is the cloud top visible from space. UV absorber (possibly S₂O or FeCl₃) creates dark markings.
Middle Cloud (50-58 km)
Largest H₂SO₄ droplets, radii up to ~3.5 μm. Highest optical depth layer. Temperature: 260-340 K. Concentrated sulfuric acid (~85% H₂SO₄ by weight).
Lower Cloud (48-50 km)
Larger, fewer particles. Temperature: 340-380 K. Below ~48 km, sulfuric acid droplets evaporate due to heat. "Rain" never reaches the surface (virga).
Sulfuric acid formation cycle:
$$\text{SO}_2 + \text{O} \rightarrow \text{SO}_3 \quad ; \quad \text{SO}_3 + \text{H}_2\text{O} \rightarrow \text{H}_2\text{SO}_4$$
Below clouds, thermal decomposition reverses the process: $\text{H}_2\text{SO}_4 \rightarrow \text{SO}_3 + \text{H}_2\text{O} \rightarrow \text{SO}_2 + \text{O} + \text{H}_2\text{O}$
The cloud layer has a total visible optical depth of $\tau \approx 25-40$, reflecting ~77% of incident sunlight (giving Venus its high Bond albedo). Only ~2.5% of solar radiation reaches the surface. However, the thermal opacity of the sub-cloud CO₂ is so enormous that even this small fraction maintains surface temperatures of 737 K.
Venus vs Earth: Divergent Evolution
Venus and Earth formed from the same region of the proto-solar disk, with similar compositions, sizes, and initial volatile inventories. Their dramatically different present states illustrate how small differences in initial conditions can lead to divergent atmospheric evolution.
Total CO₂ Inventory
Venus atmosphere: ~90 bar CO₂. Earth: ~90 bar equivalent of CO₂ is locked in carbonate rocks (limestone, dolomite) and dissolved in oceans. The total carbon inventories are remarkably similar -- the difference is phase, not quantity. Earth's liquid water enables the carbonate-silicate weathering cycle that sequesters CO₂ from the atmosphere into rocks over geological time.
Water Loss
Venus D/H ratio: ~150x terrestrial. This extreme enrichment indicates Venus once had substantial water (possibly an ocean) that was lost through: (1) UV photodissociation of H₂O in the upper atmosphere; (2) hydrogen escape to space (hydrodynamic and Jeans escape); (3) oxygen consumed by surface mineral oxidation. Estimated original water: 5-500 m global equivalent layer (GEL).
Carbonate-silicate cycle (operative on Earth, absent on Venus):
$$\text{CaSiO}_3 + \text{CO}_2 \rightleftharpoons \text{CaCO}_3 + \text{SiO}_2$$
This weathering thermostat removes CO₂ when warm and releases it when cold, stabilizing Earth's climate over geological time
Hadley Circulation on a Slowly Rotating Planet
Venus's extremely slow rotation (243 Earth days, retrograde) means the Coriolis parameter is negligibly small. This fundamentally changes the large-scale atmospheric circulation compared to Earth.
Coriolis parameter comparison:
$$f = 2\Omega \sin\phi \quad ; \quad \frac{f_{\text{Venus}}}{f_{\text{Earth}}} = \frac{\Omega_V}{\Omega_E} \approx \frac{1}{243} \approx 0.004$$
Single Hadley Cell
On Earth, the Hadley cell extends to ~30° latitude before Coriolis deflection creates the subtropical jet. On Venus, the negligible Coriolis force allows a single thermally-direct Hadley cell to extend from equator to pole in each hemisphere. Warm air rises at the equator, flows poleward at altitude, descends at the poles, and returns near the surface.
Polar Vortices
The descending branch of the Hadley cell creates persistent polar vortices at both poles. These are warm-core structures (due to adiabatic heating of descending air) with distinctive S-shaped or dipole morphology observed by Venus Express. The polar vortex temperature at 65 km is ~40 K warmer than surrounding regions.
Atmospheric and Thermal Tides
On Venus, thermal tides (pressure oscillations driven by solar heating) are a dominant dynamical feature, unlike Earth where they are secondary to baroclinic instability.
Solar thermal tide forcing:
$$Q(\lambda, \phi, z) = Q_0(z) \sum_{s} J_s(\phi, z) \cos\left(s\lambda - s\Omega_{\text{sun}} t\right)$$
where s is the zonal wavenumber, $\Omega_{\text{sun}}$ is the apparent solar angular velocity, and $J_s$ is the heating function
The diurnal (s=1) and semidiurnal (s=2) thermal tides play a crucial role in angular momentum transport. On slowly rotating planets, thermal tides can transfer angular momentum from the mean flow to the atmosphere, potentially contributing to super-rotation. The solar semidiurnal tide on Venus is particularly important because its phase speed is prograde relative to the surface, allowing it to accelerate the super-rotating flow.
Phosphine Detection and Exploration Missions
Phosphine (PH₃) Controversy
In September 2020, Greaves et al. reported detection of PH₃ at ~20 ppb in Venus's cloud layer via JCMT and ALMA submillimeter observations. On Earth, PH₃ is associated with anaerobic biological processes or industrial chemistry. No known abiotic mechanism can produce PH₃ in Venus's oxidizing atmosphere. However, subsequent reanalyses disputed the detection, finding possible calibration artifacts. The ALMA data, when recalibrated, showed the signal reduced to ~1-5 ppb or absent entirely. The debate remains unresolved and has motivated new Venus missions.
Cloud Habitability Hypothesis
The Venus cloud layer at 48-60 km altitude has temperatures (300-350 K) and pressures (0.4-2 bar) compatible with liquid water chemistry. Some astrobiologists hypothesize extremophile microorganisms could exist as aeroplankton within sulfuric acid cloud droplets, using UV-absorbing pigments as energy sources. The "unknown UV absorber" in Venus's clouds has never been chemically identified and remains a potential (though speculative) biosignature.
Past Missions
- Venera 7-14 (1970-1982): Surface landers, survived 23-127 min
- Pioneer Venus (1978): Orbiter + 4 atmospheric probes
- Vega 1-2 (1985): Balloons in cloud layer (~54 km)
- Magellan (1990-94): Radar mapping of entire surface
- Venus Express (2006-2014): Comprehensive atmospheric study
- Akatsuki (2015-present): Cloud dynamics, thermal structure
Upcoming Missions
- DAVINCI (NASA, ~2031): Atmospheric descent probe
- VERITAS (NASA): Orbital radar and emissivity mapper
- EnVision (ESA, ~2031): Orbital geological study
- Shukrayaan (ISRO): Orbiter with SAR
Radiative Transfer in the Dense Atmosphere
The extreme optical thickness of Venus's atmosphere presents unique challenges for radiative transfer. In the thermal infrared, the CO₂ absorption bands are so pressure-broadened that virtually no window regions remain.
Pressure-broadened line half-width (Lorentz profile):
$$\gamma_L = \gamma_0 \left(\frac{p}{p_0}\right)\left(\frac{T_0}{T}\right)^n$$
At 92 bar, absorption lines are ~1000x broader than at 1 bar, causing massive overlap and near-continuous opacity
Surface energy balance:
$$F_{\text{solar,surface}} \approx 17 \text{ W/m}^2 = \sigma T_s^4 (1 - \epsilon_{\text{atm}}) \implies \epsilon_{\text{atm}} \approx 0.9999$$
The atmosphere has an effective emissivity of 99.99% -- nearly a perfect blackbody blanket
Fortran: Pressure-Temperature Profile with Real Gas EOS
This Fortran program computes the Venus atmospheric pressure-temperature profile from the surface to 100 km altitude, using the van der Waals equation of state to account for real gas behavior at the extreme pressures near the surface.
program venus_pt_profile
implicit none
! === Constants ===
double precision, parameter :: R_gas = 8.314d0 ! J/(mol·K)
double precision, parameter :: g = 8.87d0 ! m/s^2 Venus gravity
double precision, parameter :: M_co2 = 0.04401d0 ! kg/mol (CO2)
double precision, parameter :: cp_co2 = 850.0d0 ! J/(kg·K) at ~400K
double precision, parameter :: pi = 3.14159265d0
! Van der Waals constants for CO2
double precision, parameter :: a_vdw = 0.3658d0 ! Pa·m^6/mol^2
double precision, parameter :: b_vdw = 4.286d-5 ! m^3/mol
! === Surface conditions ===
double precision, parameter :: T_surf = 737.0d0 ! K
double precision, parameter :: p_surf = 9.2d6 ! Pa (92 bar)
! === Profile computation ===
integer, parameter :: nlevels = 1001
double precision :: z(nlevels), T(nlevels), p(nlevels), rho(nlevels)
double precision :: z_max, dz, T_lapse, Vm, rho_ideal, rho_vdw
double precision :: f_vdw, df_vdw, p_vdw
integer :: i, iter
double precision :: z_km
z_max = 100.0d3 ! 100 km
dz = z_max / dble(nlevels - 1)
! Initialize surface values
z(1) = 0.0d0
T(1) = T_surf
p(1) = p_surf
! Compute density at surface using van der Waals EOS
! (p + a/Vm^2)(Vm - b) = RT -> solve for Vm (molar volume)
! Newton iteration
Vm = R_gas * T(1) / p(1) ! ideal gas initial guess
do iter = 1, 50
f_vdw = (p(1) + a_vdw / Vm**2) * (Vm - b_vdw) - R_gas * T(1)
df_vdw = (p(1) + a_vdw / Vm**2) - 2.0d0 * a_vdw * (Vm - b_vdw) / Vm**3
Vm = Vm - f_vdw / df_vdw
if (abs(f_vdw) < 1.0d-6) exit
end do
rho(1) = M_co2 / Vm ! density in kg/m^3
! Dry adiabatic lapse rate for CO2
T_lapse = g / cp_co2 ! ~0.0104 K/m = 10.4 K/km
write(*,'(A)') '======================================================='
write(*,'(A)') ' Venus Atmospheric Profile (van der Waals EOS)'
write(*,'(A)') '======================================================='
write(*,'(A8, A12, A14, A14, A14)') &
'z(km)', 'T(K)', 'p(Pa)', 'p(bar)', 'rho(kg/m3)'
write(*,'(F8.1, F12.1, ES14.4, F14.3, F14.3)') &
z(1)/1.0d3, T(1), p(1), p(1)/1.0d5, rho(1)
! === Integrate upward ===
do i = 2, nlevels
z(i) = dble(i-1) * dz
! Temperature: follow adiabat in lower atmosphere,
! then transition to observed profile
if (z(i) < 60.0d3) then
! Troposphere: ~adiabatic with reduced lapse rate
! Observed lapse rate is sub-adiabatic above 30 km
if (z(i) < 30.0d3) then
T(i) = T(i-1) - T_lapse * dz ! near-adiabatic
else
T(i) = T(i-1) - 0.5d0 * T_lapse * dz ! reduced lapse
end if
else
! Above 60 km: roughly isothermal to mesosphere
T(i) = max(T(i-1) - 2.0d-3 * dz, 170.0d0)
end if
! Pressure: hydrostatic equilibrium
! dp = -rho * g * dz = -(p * M / (R * T)) * g * dz [ideal gas approx]
p(i) = p(i-1) * exp(-M_co2 * g * dz / (R_gas * T(i-1)))
! Density from van der Waals (Newton iteration)
Vm = R_gas * T(i) / max(p(i), 1.0d0) ! initial guess
if (Vm < b_vdw * 1.1d0) Vm = b_vdw * 1.1d0
do iter = 1, 50
f_vdw = (p(i) + a_vdw / Vm**2) * (Vm - b_vdw) - R_gas * T(i)
df_vdw = (p(i) + a_vdw / Vm**2) - 2.0d0 * a_vdw * (Vm - b_vdw) / Vm**3
if (abs(df_vdw) < 1.0d-30) exit
Vm = Vm - f_vdw / df_vdw
if (Vm < b_vdw * 1.01d0) Vm = b_vdw * 1.01d0
if (abs(f_vdw) < 1.0d-6) exit
end do
rho(i) = M_co2 / Vm
! Print every 10 km
z_km = z(i) / 1.0d3
if (mod(nint(z_km), 10) == 0 .and. abs(z_km - nint(z_km)) < 0.2d0) then
write(*,'(F8.1, F12.1, ES14.4, F14.3, F14.3)') &
z_km, T(i), p(i), p(i)/1.0d5, rho(i)
end if
end do
! === Compare ideal vs van der Waals at surface ===
rho_ideal = p_surf * M_co2 / (R_gas * T_surf)
write(*,'(A)') ''
write(*,'(A,F10.3,A)') ' Surface density (ideal gas): ', rho_ideal, ' kg/m³'
write(*,'(A,F10.3,A)') ' Surface density (van der Waals): ', rho(1), ' kg/m³'
write(*,'(A,F8.2,A)') ' Real gas correction: ', &
(rho(1) - rho_ideal)/rho_ideal * 100.0d0, ' %'
end program venus_pt_profileInteractive Simulation: Runaway Greenhouse Effect
PythonModels surface temperature vs optical depth for Venus, Earth, and Mars using T_surface = T_eff * (1 + 3*tau/4)^(1/4). Shows how Venus's thick CO2 atmosphere creates extreme greenhouse warming of ~510 K, and compares temperature profiles for all three planets.
Click Run to execute the Python code
Code will be executed with Python 3 on the server