4.2 Fronts

Fronts are boundaries between air masses of different temperature and humidity. They are zones of enhanced weather activity including precipitation, wind shifts, and temperature changes.

Types of Fronts

Cold Front

Cold air advances, lifts warm air. Steep slope (~1:50-1:100). Narrow band of intense weather.

Weather: towering cumulus, thunderstorms, heavy rain, wind shift

Warm Front

Warm air advances, rides over cold air. Gentle slope (~1:200-1:300). Wide band of weather.

Weather: cirrus → altostratus → nimbostratus, steady rain

Occluded Front

Cold front catches warm front, lifting warm sector aloft.

Can be cold or warm occlusion depending on relative temperatures

Stationary Front

Neither air mass advancing. Prolonged weather along front.

Frontal Slope & Margules Equation

The slope of a frontal surface is governed by the Margules equation, which balances the Coriolis force, pressure gradient, and buoyancy across the front:

$$\tan\alpha = \frac{f\,(T_2 v_2 - T_1 v_1)}{g\,(T_2 - T_1)}$$

where $\alpha$ is the frontal slope angle, $f$ is the Coriolis parameter,$T_1, T_2$ are the temperatures of the cold and warm air masses, and$v_1, v_2$ are the along-front wind components

The temperature gradient across a front intensifies sharply over a narrow zone. The magnitude of the horizontal temperature gradient can be expressed as:

$$|\nabla_H T| = \left|\frac{\partial T}{\partial n}\right| \approx \frac{\Delta T}{\Delta n}$$

where $n$ is the cross-front direction; typical values are 5--10 K per 100 km for active fronts

The rate of frontal lifting is determined by the frontal slope and the component of wind normal to the front:

$$w_{\text{front}} = u_n \tan\alpha$$

where $u_n$ is the cross-front wind component and $\alpha$ is the frontal slope angle

Frontogenesis

Frontogenesis is the process of intensifying a front by increasing the horizontal temperature gradient. The scalar frontogenesis function is defined as:

$$F = \frac{d}{dt}|\nabla_p \theta| = \frac{1}{|\nabla_p \theta|}\left(\nabla_p \theta \cdot \nabla_p \frac{d\theta}{dt} - \nabla_p \theta \cdot \nabla_p \mathbf{V} \cdot \nabla_p \theta\right)$$

$F > 0$ indicates frontogenesis (gradient tightening), $F < 0$ indicates frontolysis (gradient weakening)

The quasi-geostrophic (QG) frontogenesis drives a cross-frontal ageostrophic circulation. This thermally direct circulation is described by the Sawyer--Eliassen equation:

$$N^2 \frac{\partial^2 \psi}{\partial y^2} + f^2 \frac{\partial^2 \psi}{\partial z^2} = 2Q_{\text{front}}$$

where $\psi$ is the ageostrophic streamfunction, $N^2$ is the Brunt-Vaisala frequency squared, and $Q_{\text{front}}$ is the QG frontogenetic forcing

This circulation features ascending warm air on the warm side of the front and descending cold air on the cold side, reinforcing the frontal structure.

Interactive Simulation: Frontal Cross-Section Model

Python

Plots cross-sections through cold and warm fronts using the Margules slope equation. Shows temperature fields, frontal surface geometry, wind shear arrows, and cloud layer structure for both front types side by side.

frontal_cross_section.py106 lines

import numpy as np
import matplotlib.pyplot as plt

# Cross-sections through cold front and warm front using Margules slope
g = 9.81
f = 1.0e-4  # Coriolis parameter at ~45N

# Cold front parameters
T_cold = 275  # K (cold air)
T_warm = 285  # K (warm air)
v_cold = 5    # m/s (along-front wind, cold side)
v_warm = 15   # m/s (along-front wind, warm side)

# Margules slope: tan(alpha) = f*(T2*v2 - T1*v1) / (g*(T2 - T1))
tan_alpha_cold = f * (T_warm * v_warm - T_cold * v_cold) / (g * (T_warm - T_cold))
slope_cold = 1 / tan_alpha_cold  # 1:N ratio

# Warm front: gentler slope
v_cold_wf = 2; v_warm_wf = 10
tan_alpha_warm = f * (T_warm * v_warm_wf - T_cold * v_cold_wf) / (g * (T_warm - T_cold))
slope_warm = 1 / tan_alpha_warm

fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 5.5))
fig.patch.set_facecolor('#0f172a')
for ax in (ax1, ax2):
    ax.set_facecolor('#1e293b')
    ax.tick_params(colors='#cbd5e1')
    for spine in ax.spines.values():
        spine.set_color('#334155')

# --- Cold Front Cross-Section ---
x_cf = np.linspace(-200, 300, 500)  # km from front at surface
z_max = 10  # km
# Frontal surface: z = tan_alpha * x (for x > 0, front tilts over cold air)
z_front_cf = np.clip(tan_alpha_cold * (-x_cf) * 1000 / 1000, 0, z_max)  # km

# Temperature field
X_cf, Z_cf = np.meshgrid(x_cf, np.linspace(0, z_max, 100))
# Transition zone across front
front_z_grid = np.clip(tan_alpha_cold * (-X_cf) * 1000 / 1000, 0, z_max)
sigma_cf = 0.5  # transition width (km)
blend = 0.5 * (1 + np.tanh((Z_cf - front_z_grid) / sigma_cf))
T_field_cf = T_cold + blend * (T_warm - T_cold) - 6.5 * Z_cf

cf1 = ax1.contourf(X_cf, Z_cf, T_field_cf, levels=20, cmap='coolwarm', alpha=0.8)
cbar1 = fig.colorbar(cf1, ax=ax1, shrink=0.85)
cbar1.set_label('Temperature (K)', color='#cbd5e1')
cbar1.ax.tick_params(colors='#cbd5e1')
ax1.plot(x_cf, z_front_cf, color='#38bdf8', linewidth=3)
ax1.fill_between(x_cf, 0, z_front_cf, color='#38bdf8', alpha=0.15)
# Wind shear arrows
for xp in [-100, -50, 50, 150]:
    zp = 3
    if xp < 0:
        ax1.annotate('', xy=(xp + 20, zp), xytext=(xp, zp),
                     arrowprops=dict(arrowstyle='->', color='#60a5fa', lw=2))
    else:
        ax1.annotate('', xy=(xp + 30, zp + 0.3), xytext=(xp, zp),
                     arrowprops=dict(arrowstyle='->', color='#f87171', lw=2))

ax1.text(-150, 1.5, 'COLD AIR', color='#60a5fa', fontsize=13, fontweight='bold')
ax1.text(150, 6, 'WARM AIR', color='#f87171', fontsize=13, fontweight='bold')
ax1.set_xlabel('Distance from surface front (km)', color='#cbd5e1')
ax1.set_ylabel('Height (km)', color='#cbd5e1')
ax1.set_title(f'Cold Front (slope ~1:{abs(slope_cold):.0f})', color='#cbd5e1', fontsize=13)

# --- Warm Front Cross-Section ---
x_wf = np.linspace(-400, 200, 500)
z_front_wf = np.clip(tan_alpha_warm * x_wf * 1000 / 1000, 0, z_max)

X_wf, Z_wf = np.meshgrid(x_wf, np.linspace(0, z_max, 100))
front_z_wf = np.clip(tan_alpha_warm * X_wf * 1000 / 1000, 0, z_max)
blend_wf = 0.5 * (1 + np.tanh((Z_wf - front_z_wf) / sigma_cf))
T_field_wf = T_cold + blend_wf * (T_warm - T_cold) - 6.5 * Z_wf

cf2 = ax2.contourf(X_wf, Z_wf, T_field_wf, levels=20, cmap='coolwarm', alpha=0.8)
cbar2 = fig.colorbar(cf2, ax=ax2, shrink=0.85)
cbar2.set_label('Temperature (K)', color='#cbd5e1')
cbar2.ax.tick_params(colors='#cbd5e1')
ax2.plot(x_wf, z_front_wf, color='#f87171', linewidth=3)
ax2.fill_between(x_wf, 0, z_front_wf, color='#38bdf8', alpha=0.15)

# Cloud layers
cloud_labels = ['Ci', 'Cs', 'As', 'Ns']
cloud_x = [-300, -200, -100, -20]
cloud_z = [8, 6.5, 4.5, 2]
for cx, cz, cl in zip(cloud_x, cloud_z, cloud_labels):
    ax2.scatter(cx, cz, s=1500, c='#94a3b8', alpha=0.4, marker='o')
    ax2.text(cx, cz - 0.5, cl, ha='center', fontsize=9, color='#e2e8f0')

ax2.text(-300, 1.5, 'COLD AIR', color='#60a5fa', fontsize=13, fontweight='bold')
ax2.text(100, 6, 'WARM AIR', color='#f87171', fontsize=13, fontweight='bold')
ax2.set_xlabel('Distance from surface front (km)', color='#cbd5e1')
ax2.set_ylabel('Height (km)', color='#cbd5e1')
ax2.set_title(f'Warm Front (slope ~1:{abs(slope_warm):.0f})', color='#cbd5e1', fontsize=13)

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

print("=== Frontal Cross-Section Model ===")
print(f"Margules equation: tan(a) = f(T2*v2 - T1*v1) / [g(T2-T1)]")
print(f"Cold front: slope = 1:{abs(slope_cold):.0f}, tan(a) = {tan_alpha_cold:.5f}")
print(f"Warm front: slope = 1:{abs(slope_warm):.0f}, tan(a) = {tan_alpha_warm:.5f}")
print(f"Temperature contrast: {T_warm - T_cold:.0f} K")
print(f"Cold front is {abs(slope_warm/slope_cold):.1f}x steeper than warm front")

Click Run to execute the Python code

Code will be executed with Python 3 on the server

← 4.1 Air Masses 4.3 Cyclones →