Climate Model Fundamentals

What Are Climate Models?

Climate models are mathematical representations of the climate system based on physical, chemical, and biological principles. They simulate interactions between the atmosphere, ocean, land surface, and ice.

History of Climate Modeling

The first climate models were developed in the 1960s by Syukuro Manabe and Joseph Smagorinsky. Early models were simple energy balance models with just a few equations.

Where is incoming solar radiation and is outgoing longwave radiation.

Energy Balance Model

The simplest climate model is the zero-dimensional energy balance model:

Where:

• is the heat capacity of the climate system ()
• is the global mean temperature ()
• is the solar constant ()
• is the planetary albedo ()
• is the effective emissivity ()
• is the Stefan-Boltzmann constant ()

Python Implementation

import numpy as np
import matplotlib.pyplot as plt

def energy_balance(dt, T, S=1361, alpha=0.30, epsilon=0.61, C=4.0e8):
    sigma = 5.67e-8
    Q_net = S * (1 - alpha) / 4 - epsilon * sigma * T**4
    return Q_net / C

# Time stepping
T = 288.0  # Initial temperature in Kelvin
dt = 86400 * 365  # One year in seconds
years = 1000
temps = []
for _ in range(years):
    T += energy_balance(dt, T) * dt
    temps.append(T - 273.15)

plt.plot(range(years), temps)
plt.xlabel('Year')
plt.ylabel('Temperature (°C)')
plt.title('Energy Balance Model Convergence')
plt.show()

Model Components

A modern climate model includes:

1. Atmosphere module — radiation, convection, clouds
2. Ocean module — currents, mixing, sea ice
3. Land surface module — hydrology, vegetation, soil
4. Sea ice module — dynamics and thermodynamics

Radiative Transfer Equation

The fundamental equation governing radiative transfer in the atmosphere:

Where is the spectral radiance, is the absorption coefficient, is the density, and is the emission coefficient.

The Role of CO2

Carbon dioxide absorbs outgoing longwave radiation. The relationship between CO2 concentration and radiative forcing:

Where is the change in radiative forcing (), is the CO2 concentration, and is the pre-industrial concentration (280 ppm).

Model Grid Structure

Climate models divide the Earth into a three-dimensional grid. A typical GCM has:

• Horizontal resolution: 1 x 1 degree (~100 km)
• Vertical layers: 20-60 levels
• Time step: 15-30 minutes

The number of grid cells:

Climate Sensitivity

Equilibrium Climate Sensitivity (ECS) is defined as:

Where for doubling CO2.

Validation Metrics

Climate models are validated using:

• Global mean temperature
• Spatial patterns of temperature and precipitation
• Seasonal cycles
• Interannual variability (e.g., ENSO)
• Response to volcanic eruptions

Key Equations Summary

Exercises

1. Modify the Python code to include a linear ice-albedo feedback:
2. Calculate the equilibrium temperature for 560 ppm CO2
3. Plot the forcing curve for CO2 from 280 to 1000 ppm