Antipodal Vortex Model for Atmospheric Blocking

• The paper introduces an analytical point-vortex model that characterizes atmospheric blockings as interacting antipodal vortex pairs on a rotating sphere.
• It employs a Hamiltonian framework to derive exact criteria for stationarity, stability, and typical oscillation periods consistent with observed blocking lifetimes.
• By incorporating polar vortices, the model elucidates mechanisms for abrupt onset and breakdown of blockings, matching key climatological patterns.

A point-vortex model for atmospheric blockings is an analytically tractable representation of the large-scale, quasi-stationary high- and low-pressure dipoles that characterize atmospheric blocking events. In recent developments, such models on a rotating sphere focus on antipodal point vortex pairs (APV), which rigorously satisfy the governing hydrodynamic equations for thin atmospheric shells. These approaches provide a finite-dimensional, Hamiltonian dynamical system for vortex interaction, yielding exact criteria for blockade stationarity, stability, and lifetimes in realistic, planetary-scale flows (Mokhov et al., 2020, Mokhov et al., 2016).

1. Antipodal Point Vortex Dynamics on the Rotating Sphere

The foundational formulation is based on the absolute vorticity conservation law for a thin, rotating spherical shell of radius $R$ at angular velocity $\Omega$: $$\frac{\partial\Omega_a}{\partial t} + \frac{1}{R\sin\theta}\frac{\partial}{\partial\varphi}(\psi_\theta \Omega_a) - \frac{1}{R}\frac{\partial}{\partial\theta}(\psi_\varphi \Omega_a) = 0$$ where $\psi(\theta,\varphi,t)$ is the streamfunction, the velocity components are derived from its derivatives, and $\Omega_a = \zeta + f$ with $\zeta = -\Delta_{S}\psi$ and $f=2\Omega\cos\theta$.

A critical construction is the background “zero-absolute-vorticity” streamfunction $\psi_0(\theta) = -\Omega R^2\cos\theta$, which exactly cancels the Coriolis term $f$; this reduces the nonlinear coupling, allowing the total streamfunction to be written as

$$\psi(\theta, \varphi, t) = \psi_0(\theta) + \psi_\text{PV}(\theta, \varphi, t),$$

with $\psi_\text{PV}$ a superposition of APV singularities,

$$\psi_\text{APV}(\theta, \varphi) = \frac{\Gamma}{4\pi}\ln\left(\frac{1+\cos\gamma_i}{1-\cos\gamma_i}\right),$$

where each APV consists of a cyclone–anticyclone pair at antipodal locations $(\theta_i, \varphi_i)$ and $(\pi-\theta_i, \varphi_i+\pi)$.

For $N$ APVs, the generalized streamfunction is

$$\psi(\theta, \varphi, t) = -\Omega R^2\cos\theta + \sum_{i=1}^N \frac{\Gamma_i}{4\pi} \ln\frac{1+\cos\gamma_i(t)}{1-\cos\gamma_i(t)},$$

where $\Gamma_i$ denotes the circulation and $\cos\gamma_i$ is the chordal cosine between field and vortex positions (Mokhov et al., 2020, Mokhov et al., 2016).

2. Hamiltonian Structure and Equations of Motion

The APV system is exactly Hamiltonian, with vortex coordinates $(\theta_i, \varphi_i)$ evolving according to

$$\dot\theta_i = \frac{1}{R^2\sin\theta_i}\frac{\partial H}{\partial\varphi_i}, \qquad \dot\varphi_i = -\frac{1}{R^2\sin\theta_i}\frac{\partial H}{\partial\theta_i},$$

with the Hamiltonian

$$H = -\frac{1}{8\pi}\sum_{i\ne j}\Gamma_i\Gamma_j\ln[1-\cos\gamma_{ij}] - \Omega\sum_{i=1}^N \Gamma_i R^2\cos\theta_i,$$

$\cos\gamma_{ij}$ being the mutual chord angle.

This system conserves total kinetic energy, the $z$-component of angular momentum $M_z = \sum_{i=1}^N \Gamma_i R^2\cos\theta_i$, and vector impulse (Mokhov et al., 2020). The APV model is distinct in requiring vortex pairs of equal and opposite circulation at antipodal sites, ensuring that the streamfunction solves the spherical hydrodynamic equations exactly for all times (Mokhov et al., 2016).

3. Equilibrium, Stationarity, and Criteria for Blocking

Specializing to $N=2$, a steady (stationary) vortex pair satisfies

$$\dot\theta_1 = \dot\theta_2 = 0, \qquad \dot\varphi_1 = \dot\varphi_2 = 0.$$

Two key conditions emerge:

• Relative equilibrium requires the vortex strengths and co-latitudes to satisfy

$$\frac{\Gamma_1}{\Gamma_2} = -\frac{\sin\theta_2}{\sin\theta_1}$$

(enforcing opposite sign circulations on the same meridian).

• Zero drift (absolute stationarity) imposes

$$\Gamma_2 = -\Omega \, 4\pi R^2 \frac{\sin\theta_1\sin\theta_2}{\sin(\theta_2-\theta_1)}, \quad \Gamma_1 / \Gamma_2 < 0.$$

These link APV positions and strengths to planetary rotation, producing blocking configurations with one stronger midlatitude cyclone ($\Gamma_1 > 0$), balanced by a weaker high-latitude anticyclone ($\Gamma_2 < 0$). Observed vortex ratios for the Icelandic Low and Azores High ($\Gamma_1/\Gamma_2 \approx -1.9$) select $\theta_1 \approx 25^\circ$, $\theta_2 \approx 55^\circ$, consistent with climatological blocking centers (Mokhov et al., 2020).

Incorporating a polar APV modifies the stationarity and yields regimes where the two-cell equilibrium becomes exponentially unstable, providing a mechanism for abrupt block onset or breakdown on monthly timescales (Mokhov et al., 2020).

4. Stability Analysis and Oscillation Modes

Linearizing the equations about the equilibrium state, one obtains a simple oscillator

$$\ddot{x} + \omega^2 x = 0, \qquad \omega = |\Omega\sin(\theta_2-\theta_1)|,$$

with oscillation period

$T = \frac{2\pi}{\Omega |\sin(\theta_2-\theta_1)|},$

which predicts typical blocking lifetimes in the 1–3 week range. Nonlinear analysis confirms center-type, nonlinearly stable dynamics around the equilibrium for finite-amplitude perturbations. Notably, the amplitude of anticyclonic excursions exceeds that of cyclonic ones, mirroring the greater observed variability in the latitude of blocking highs versus lows (Mokhov et al., 2020).

Explicit stability conditions in the presence of additional polar vortex strength reveal sharp boundaries in parameter space for dynamically stable blocking (i.e., splitting-flow) regimes. These are expressible analytically in terms of the vortex parameters and planetary rotation, offering quantifiable criteria for block persistence or breakdown (Mokhov et al., 2016).

5. Interpretations, Application to Observed Blockings, and Model Limitations

The point vortex APV model captures the essential dynamical features of atmospheric blocking: stationary or slowly drifting dipole pairs, characteristic location and intensity ratios, and realistic temporal scales for formation, persistence, and breakdown. The model yields precise predictions for the stationary locations and mutual strengths of major blocking centers (Icelandic Low and Azores High), matches observed oscillation periods, and explains the greater morphological variability of blocking highs. Inclusion of a polar APV introduces the capacity for abrupt transitions.

Significantly, the APV framework provides an exact, closed, finite-dimensional system in contrast to the traditional infinite-dimensional PDEs inherent to geophysical fluid dynamics. It also quantifies the impact of planetary rotation and polar vortex anomalies on blocking stability (Mokhov et al., 2020, Mokhov et al., 2016).

Limitations include the idealization to point singularities (excluding finite area or filamentary vortices), omission of baroclinic and diabatic effects, and the lack of explicit stochasticity or topographic forcing. However, the APV approach is parsimonious yet physically consistent with observed atmospheric blocking phenomenology.

6. Connections to Alternative Vortex and Wave Models

The APV construction is complementary to planar point-vortex and quasigeostrophic models, such as the SQG and SQG+ point-vortex formulations with order-Rossby corrections (Lee et al., 18 Aug 2025), and hybrid vortex–wave models developed for the Charney-Hasegawa-Mima equation (Leoncini et al., 2011). These planar or quasi-planar approaches capture ageostrophic, three-dimensional, and wave–vortex coupling effects, revealing phenomena such as non-conservation of Hamiltonians, finite time vertical excursions, and chaotic, merger or splitting patterns.

In the context of atmospheric blocking, minimal “point-vortex plus wave” models explain the phase-locking, resonance, or detachment of blocking dipoles in the presence of planetary-scale Rossby waves: mode-locked states correspond to persistent blockings, while parameter excursions induce block decay and chaos (Leoncini et al., 2011). SQG+ models further bridge toward realistic three-dimensionality by introducing ageostrophic divergence and weak vertical transport, supplying a continuum from idealized vortex blockings to complex mesoscale events (Lee et al., 18 Aug 2025).

7. Synthesis and Perspective

The point-vortex model, with a focus on antipodal vortex configurations on the rotating sphere, constitutes a mathematically exact, observationally consistent reduction of atmospheric blocking dynamics to a finite-dimensional Hamiltonian system. It delineates the geometric, dynamical, and stability requirements for block initiation, stationarity, and breakdown, and can be systematically extended to include polar vortices for more complex block evolution scenarios.

By clarifying the intrinsic dynamical thresholds, conserved quantities, and topological mapping to observed atmospheric structures, the APV point-vortex framework enables both predictive parameterization and fundamental understanding of atmospheric blocking as a global dynamical phenomenon (Mokhov et al., 2020, Mokhov et al., 2016).