8.3 Ocean Heat Content

The ocean has absorbed approximately 93% of the excess heat trapped by anthropogenic greenhouse gases since the 1970s. Ocean heat content (OHC) is the most reliable measure of Earth's energy imbalance and the key metric for quantifying global warming, thermosteric sea level rise, and committed future warming.

Ocean Heat Content: Definition and Measurement

Ocean heat content is defined as the vertically integrated thermal energy per unit area:

$$\text{OHC} = \int_{-H}^{0} \rho(z)\, c_p(z)\, T(z)\, dz$$

Units: J/m². Global OHC change often expressed in zettajoules (ZJ = 10²¹ J)

The change in OHC relative to a reference period is more meaningful than the absolute value:

$$\Delta\text{OHC} = \int_{-H}^{0} \rho\, c_p\, \Delta T(z)\, dz \approx \rho_0 c_p \int_{-H}^{0} \Delta T(z)\, dz$$

~93%

Of excess heat absorbed by the ocean

~0.5–1.0 W/m²

Current planetary energy imbalance

>400 ZJ

Added since 1970

The Argo Observing Network

The Argo array of ~4,000 profiling floats has revolutionized OHC monitoring since 2005. Each float cycles from 2000 m depth to the surface every 10 days, measuring temperature and salinity profiles:

Standard Argo

Profiles 0–2000 m. ~4,000 active floats. Near-global coverage (excluding ice-covered and marginal seas). Temperature accuracy ±0.002°C. 4–5 year battery life.

Deep Argo

Profiles 0–6000 m. Expanding network (>200 floats). Critical for measuring deep ocean warming below 2000 m, which accounts for ~10% of total OHC change but is poorly sampled.

Before Argo, OHC estimates relied on ship-based expendable bathythermographs (XBTs), which had systematic depth and temperature biases that required correction (Gouretski & Koltermann, 2007; Levitus et al., 2012). These XBT corrections are critical because uncorrected biases can produce spurious decadal variability in OHC time series.

Biogeochemical Argo (BGC-Argo)

An extension of the Argo program that adds sensors for dissolved oxygen, pH, nitrate, chlorophyll, backscatter, and irradiance. BGC-Argo provides simultaneous physical and biogeochemical profiles, enabling studies of ocean acidification, deoxygenation, and carbon uptake that are intrinsically linked to ocean heat content changes.

Heat Content Changes by Depth and Basin

Warming is not uniform in depth or space. The upper 700 m accounts for roughly 60% of total OHC change, while the 700–2000 m layer accounts for ~30%:

0–700 m: +8–10 ZJ/year

Best observed layer. Strong interannual variability from ENSO. Warming rate has increased since the 1990s. Atlantic accounts for disproportionate share due to deep overturning.

700–2000 m: +4–5 ZJ/year

Significant warming detected. Southern Ocean dominates heat uptake at these depths. Less affected by ENSO variability. Argo coverage since 2005.

>2000 m: +1–2 ZJ/year

Detectable warming in abyssal waters, particularly Antarctic Bottom Water. Measured by deep hydrographic repeat sections and emerging Deep Argo.

Thermosteric Sea Level Rise

Ocean warming causes thermal expansion, raising sea level even without adding water mass. The thermosteric sea level change is related to OHC through the thermal expansion coefficient:

$$\Delta \eta_{\text{thermo}} = \int_{-H}^{0} \frac{\alpha(T,S,p)}{\rho c_p}\, \Delta Q(z)\, dz \approx \frac{\alpha_0}{\rho_0 c_p}\, \Delta\text{OHC}$$

where $\alpha = -(1/\rho)\partial\rho/\partial T$ is the thermal expansion coefficient

Thermosteric expansion contributes approximately 1.3 mm/yr to current global mean sea level rise (~35–40% of the total). The expansion coefficient increases with temperature and pressure, so warming deep, warm waters causes more expansion per unit heat input.

Ocean Thermal Inertia and Committed Warming

The ocean's enormous heat capacity creates thermal inertia that delays the full temperature response to radiative forcing. The effective climate sensitivity can be estimated from the energy budget:

$$\Delta F = \lambda \Delta T + \frac{d(\text{OHC})}{dt \cdot A_{\text{ocean}}}$$

$\Delta F$ = radiative forcing, $\lambda$ = feedback parameter, $\Delta T$ = surface warming

Even if atmospheric greenhouse gas concentrations were stabilized today, the ocean would continue to warm for centuries as it approaches equilibrium. This "committed warming" is estimated at 0.4–0.6°C beyond current temperatures, arising from the ~0.7 W/m² planetary energy imbalance stored primarily in the ocean.

Climate Sensitivity from OHC

The effective climate sensitivity (ECS) can be estimated from the energy budget: $\text{ECS} = F_{2\times CO_2} / \lambda$ where $\lambda = (\Delta F - \dot{Q}_{ocean})/\Delta T$. Using observed OHC changes yields ECS estimates of 2.5–4.0°C, consistent with model-based estimates.

Ocean Heat Uptake Efficiency

The ocean heat uptake efficiency $\kappa = Q_{ocean}/\Delta T$ (typically 0.5–1.5 W/m²/K) determines how much of the radiative forcing goes into warming the ocean vs the atmosphere. Higher κ means slower surface warming but greater committed future warming.

Python: OHC Computation, Time Series & Thermosteric Sea Level

Python

!/usr/bin/env python3

script.py84 lines

#!/usr/bin/env python3
"""ocean_heat_content.py - OHC from profiles, time series, steric SLR"""
import numpy as np
import matplotlib.pyplot as plt

rho0 = 1025.0   # reference density (kg/m^3)
cp = 3990.0      # specific heat (J/kg/K)
alpha_T = 2.0e-4 # thermal expansion coefficient (1/K)
A_ocean = 3.61e14  # ocean area (m^2)

# --- 1. Compute OHC from temperature profiles ---
z = np.linspace(0, 2000, 200)  # depth (m)
dz = z[1] - z[0]

# Reference profile (climatology)
T_ref = 18 * np.exp(-z / 500) + 2

# Warming scenarios (anomaly profiles)
dT_2000 = 0.3 * np.exp(-z / 300)
dT_2010 = 0.5 * np.exp(-z / 350)
dT_2020 = 0.8 * np.exp(-z / 400)

fig, axes = plt.subplots(2, 2, figsize=(14, 10))

ax1 = axes[0, 0]
ax1.plot(T_ref, z, 'b-', lw=2, label='Reference')
ax1.plot(T_ref + dT_2000, z, 'g--', lw=1.5, label='Year 2000')
ax1.plot(T_ref + dT_2010, z, 'orange', lw=1.5, label='Year 2010')
ax1.plot(T_ref + dT_2020, z, 'r-', lw=2, label='Year 2020')
ax1.invert_yaxis()
ax1.set_xlabel('Temperature (C)'); ax1.set_ylabel('Depth (m)')
ax1.set_title('Temperature Profiles with Warming')
ax1.legend(); ax1.grid(True, alpha=0.3)

# --- 2. OHC time series ---
ax2 = axes[0, 1]
years = np.arange(1960, 2025)
# Approximate OHC trend (ZJ, relative to 1960)
ohc_trend = np.where(years < 1970, 0,
            5 * (years - 1970) + 0.08 * (years - 1970)**2)
# Add ENSO variability
enso_var = 8 * np.sin(2*np.pi*years/3.7 + 1.2) + \
           5 * np.sin(2*np.pi*years/5.2)
np.random.seed(42)
noise = np.cumsum(np.random.randn(len(years)) * 2)
ohc_series = ohc_trend + enso_var + noise
ohc_series -= ohc_series[0]

ax2.plot(years, ohc_series, 'r-', lw=1.5)
ax2.fill_between(years, 0, ohc_series, alpha=0.3, color='red')
ax2.set_xlabel('Year'); ax2.set_ylabel('OHC anomaly (ZJ)')
ax2.set_title('Global Ocean Heat Content (0-2000m)')
ax2.grid(True, alpha=0.3)

# --- 3. Thermosteric sea level contribution ---
ax3 = axes[1, 0]
# Compute steric height change from temperature anomaly
steric_2000 = alpha_T * np.sum(dT_2000) * dz * 1000  # mm
steric_2010 = alpha_T * np.sum(dT_2010) * dz * 1000
steric_2020 = alpha_T * np.sum(dT_2020) * dz * 1000

years_steric = np.arange(1970, 2025)
steric_ts = alpha_T * (years_steric - 1970)**1.3 * 0.3  # mm
ax3.plot(years_steric, steric_ts, 'b-', lw=2)
ax3.set_xlabel('Year'); ax3.set_ylabel('Thermosteric SLR (mm)')
ax3.set_title('Thermosteric Sea Level Rise')
ax3.grid(True, alpha=0.3)

# --- 4. Depth distribution of warming ---
ax4 = axes[1, 1]
depth_layers = ['0-300m', '300-700m', '700-2000m', '>2000m']
ohc_frac = [35, 25, 30, 10]  # approximate % of total OHC change
colors = ['#ff4444', '#ff8844', '#ffbb44', '#ffdd88']

bars = ax4.barh(depth_layers, ohc_frac, color=colors, edgecolor='white')
ax4.set_xlabel('Fraction of Total OHC Change (%)')
ax4.set_title('Distribution of Ocean Warming by Depth')
for bar, val in zip(bars, ohc_frac):
    ax4.text(bar.get_width() + 1, bar.get_y() + bar.get_height()/2,
             f'{val}%', va='center', color='white')

plt.tight_layout()
plt.savefig('ocean_heat_content.png', dpi=150)
plt.show()

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Fortran: 1D Upwelling-Diffusion Ocean Heat Uptake Model

Fortran

1D upwelling-diffusion model for ocean heat uptake

program.f9095 lines

program ocean_heat_uptake
  ! 1D upwelling-diffusion model for ocean heat uptake
  ! Used for estimating climate response and OHC evolution
  implicit none
  integer, parameter :: dp = selected_real_kind(15)
  integer, parameter :: nz = 100
  real(dp), parameter :: pi = 3.141592653589793_dp

real(dp) :: z(nz), T(nz), T_new(nz), dz, dt
  real(dp) :: kappa_v, w_up, rho0, cp0
  real(dp) :: Q_surface, lambda_fb, dT_surf
  real(dp) :: forcing, ohc, ohc_old, dOHC_dt
  real(dp) :: steric_slr, alpha_T
  integer :: i, step, nsteps, year
  real(dp) :: total_depth, t_years

! Model parameters
  total_depth = 4000.0_dp     ! ocean depth (m)
  dz = total_depth / real(nz-1, dp)
  kappa_v = 1.0e-4_dp         ! vertical diffusivity (m^2/s)
  w_up = 4.0e-7_dp            ! upwelling velocity (m/s, ~12 m/yr)
  rho0 = 1025.0_dp; cp0 = 3990.0_dp
  lambda_fb = 1.2_dp          ! climate feedback (W/m^2/K)
  alpha_T = 2.0e-4_dp         ! thermal expansion (1/K)

! Time stepping
  dt = 0.5_dp * dz**2 / kappa_v  ! stability condition
  dt = min(dt, 86400.0_dp * 30.0_dp)  ! max 30 days
  nsteps = nint(300.0_dp * 365.25_dp * 86400.0_dp / dt)  ! 300 years

! Depth grid
  do i = 1, nz
    z(i) = real(i-1, dp) * dz
  end do

! Initial condition: uniform temperature
  T = 4.0_dp  ! deep ocean reference (C)

print *, '=== 1D Upwelling-Diffusion Ocean Heat Uptake Model ==='
  print '(A,E10.3,A)', ' kappa_v = ', kappa_v, ' m^2/s'
  print '(A,F8.4,A)', ' w_up = ', w_up*3.15e7_dp, ' m/yr'
  print '(A,F6.1,A)', ' dt = ', dt/86400.0_dp, ' days'

print *, ''
  print *, ' Year    dT_surf(K)   OHC(ZJ)    Steric(mm)   Q_ocean(W/m2)'

ohc_old = 0.0_dp
  do step = 1, nsteps
    t_years = real(step, dp) * dt / (365.25_dp * 86400.0_dp)

! Radiative forcing: 1% CO2 increase per year (simplified)
    forcing = 3.7_dp * log(1.0_dp + 0.01_dp * min(t_years, 200.0_dp)) / &
              log(2.0_dp)
    forcing = min(forcing, 8.0_dp)  ! cap at ~4xCO2

! Surface boundary: Q = F - lambda * dT_surf
    dT_surf = T(1) - 4.0_dp
    Q_surface = forcing - lambda_fb * dT_surf

! Diffusion + upwelling (interior)
    do i = 2, nz-1
      T_new(i) = T(i) + dt * ( &
        kappa_v * (T(i+1) - 2.0_dp*T(i) + T(i-1)) / dz**2 - &
        w_up * (T(i) - T(i-1)) / dz )
    end do

! Surface: apply heat flux
    T_new(1) = T(1) + dt * Q_surface / (rho0 * cp0 * dz) + &
               dt * kappa_v * (T(2) - T(1)) / dz**2
    T_new(nz) = T(nz)  ! fixed bottom temperature

T = T_new

! Compute OHC and steric SLR
    ohc = 0.0_dp
    steric_slr = 0.0_dp
    do i = 1, nz
      ohc = ohc + rho0 * cp0 * (T(i) - 4.0_dp) * dz
      steric_slr = steric_slr + alpha_T * (T(i) - 4.0_dp) * dz
    end do
    ! Convert to global ZJ
    ohc = ohc * 3.61e14_dp / 1.0e21_dp

! Output every 10 years
    year = nint(t_years)
    if (mod(step, nint(10.0_dp*365.25_dp*86400.0_dp/dt)) == 0) then
      dOHC_dt = (ohc - ohc_old) / (10.0_dp * 365.25_dp * 86400.0_dp) * &
                1.0e21_dp / 3.61e14_dp
      print '(I5, F12.3, F12.2, F12.1, F12.3)', &
            year, dT_surf, ohc, steric_slr*1000.0_dp, dOHC_dt
      ohc_old = ohc
    end if
  end do

end program ocean_heat_uptake

Click Run to execute the Fortran code

Code will be compiled with gfortran and executed on the server

Major OHC Analysis Products

Several research groups produce global OHC estimates using different statistical methods, quality control procedures, and gap-filling techniques:

Levitus / NOAA NCEI

Objective analysis on a 1° grid. Quarterly updates. Extended back to 1955 using ship-based observations. Conservative gap-filling: unmeasured regions assumed zero anomaly.

Ishii / JMA

Japanese Meteorological Agency analysis. Independent XBT bias corrections. Extends from 1945. Good coverage in the Pacific basin.

Cheng / IAP-CAS

Institute of Atmospheric Physics, Chinese Academy of Sciences. Advanced mapping method fills gaps using ocean model covariances. Generally shows larger OHC trends than NOAA, likely more accurate.

Domingues / CSIRO

Australian analysis with careful attention to instrumental biases. Extended time coverage using historical data rescue. Important for understanding decadal variability.

All analyses agree on the robust long-term warming trend, but differ in magnitude and interannual variability due to different treatments of data-sparse regions and XBT fall-rate corrections.

Marine Heatwaves and Ecosystem Impacts

As ocean heat content increases, marine heatwaves (MHWs) — prolonged periods of anomalously warm ocean temperatures — have become more frequent and intense. A marine heatwave is defined when SST exceeds the local 90th percentile for at least 5 consecutive days:

$$\text{MHW intensity} = \overline{T}_{\text{event}} - T_{\text{climatology}}$$

Ecological Impacts

Mass coral bleaching (Great Barrier Reef 2016, 2017, 2020, 2024), kelp forest die-offs, harmful algal blooms, fishery collapses, and seabird mass mortality events.

Trend

MHW frequency has doubled since the 1980s. Duration increased by ~17%. Projected to further intensify under continued warming. Category IV (extreme) events becoming more common.

OHC in Climate Models and Detection-Attribution

OHC changes provide one of the strongest detection-attribution signals for anthropogenic climate change. Climate models reproduce observed OHC trends only when anthropogenic forcings (greenhouse gases, aerosols) are included:

$$\text{OHC}_{\text{obs}} = \beta_{\text{GHG}} \cdot \text{OHC}_{\text{GHG}} + \beta_{\text{nat}} \cdot \text{OHC}_{\text{nat}} + \epsilon$$

Optimal fingerprinting attributes observed OHC changes to anthropogenic and natural forcing patterns

Signal Detection

The anthropogenic OHC signal is detectable above natural variability at the 5-sigma level. Natural forcing (volcanic, solar) alone cannot explain the observed warming. The spatial pattern of warming (enhanced in the Atlantic, Southern Ocean) matches model predictions.

Model Spread

CMIP6 models show a factor-of-2 spread in OHC trends, largely driven by differences in climate sensitivity and ocean heat uptake efficiency. The ocean uptake efficiency $\kappa = \Delta OHC / (\Delta F \cdot A_{ocean} \cdot t)$ ranges from 0.5–1.5 W/m²/K across models.

Meridional Ocean Heat Transport

The ocean transports heat from the tropics toward the poles through both wind-driven (Ekman/gyre) and thermohaline circulation. The total meridional heat transport can be computed from ocean state estimates:

$$H(\phi) = \int_{-H}^{0}\int_{x_w}^{x_e} \rho c_p v T \, dx\, dz$$

Maximum poleward heat transport: ~1.5 PW at ~15°N in the Atlantic, ~0.7 PW at ~15°S in the Indo-Pacific

Atlantic Overturning

The AMOC transports ~1.3 PW northward, accounting for most of the Atlantic's heat transport. Unlike the Pacific, heat transport is northward at all latitudes due to the deep overturning cell. A weakening AMOC could reduce European warming by 0.5–1°C.

Indo-Pacific Gyres

Wind-driven subtropical gyres carry warm water poleward in western boundary currents (Kuroshio, East Australian Current). The shallow overturning circulation (subtropical cells) connects subtropical subduction to tropical upwelling zones.

The total ocean-atmosphere system transports approximately 5.5 PW poleward at mid-latitudes. The ocean contribution is roughly one-third of the total but dominates at tropical latitudes. Changes in ocean heat transport redistribute warming patterns and affect regional climate projections.

Key Concepts Summary

OHC Definition

Vertically integrated thermal energy $\int \rho c_p T\, dz$. Global OHC increase >400 ZJ since 1970. Rate ~14 ZJ/yr and accelerating.

Argo Revolution

~4,000 profiling floats sampling 0–2000 m globally every 10 days. Deep Argo extending to 6000 m. Transformed OHC monitoring since 2005.

Thermosteric SLR

Thermal expansion contributes ~1.3 mm/yr to sea level rise (~35% of total). Expansion coefficient increases with temperature, amplifying warm-water warming effects.

Committed Warming

Ocean thermal inertia means 0.4–0.6°C additional surface warming even with stabilized forcing. Earth's energy imbalance (~0.7 W/m²) stored in the ocean.

← 8.2 El Niño & La Niña 8.4 Sea Level Rise →