Journal Paper Reading Notes

Literature Review, MIT EAPS, 2026

Reading notes and reproduction status for atmospheric blocking, Rossby wave breaking (RWB), and related machine-learning methods. All code, simulations, and notes live in the unified repository:

github.com/yanxingjianken/blocking_rwb_literature_review

Status legend: ✅ Reproduced   🔧 In Progress   📖 Literature Review (GIFs generated)   📋 Planned


1. Blocking Theories — 11 Canonical Mechanisms

Each subdirectory contains a Python spectral-model simulation with animated GIFs. See the blocking README for the inheritance diagram of the common barotropic solver.

#ExperimentKey ReferenceStatus
1Energy DispersionYeh (1949) — Rossby wave group velocity & energy trapping📖
2Hydraulic JumpRex (1950) — 1-D shallow-water hydraulic jump analogy📖
3Multi-Equilibria (CdV)Charney & DeVore (1979) — 3-mode ODE, S-curve bifurcation📖
4Modon / Equivalent BarotropicMcWilliams (1980), Haines & Marshall (1987) — Dipole vortex, resonant forcing in Dedalus🔧
5Wave ResonanceTung & Lindzen (1979) — Stationary wave amplitude peaks near $\bar{U}_c = 15.8$ m/s📖
6Eddy StrainingShutts (1983), Wang & Kuang (2019) — Eddy straining dipole convergence📖
7Deterministic ChaosLegras & Ghil (1985) — 25-mode strange attractor with blocking/zonal regimes📖
8Wave Accumulation & AWB/CWBSwanson (2000, 2001) — Anticyclonic vs cyclonic wave breaking diagnosis📖
9Selective AbsorptionYamazaki & Itoh (2013) — Same-sign eddies absorbed, opposite expelled📖
10NLS Soliton / EnvelopeLuo et al. (2014, 2019) — Blocking envelope as localized soliton📖
11FAWA / Traffic JamNakamura & Huang (2018) — Finite-amplitude wave activity “traffic jam”📖

Common library: blocking/common/ — pseudospectral barotropic solver (RK4) + GIF maker.


2. Rossby Wave Breaking (RWB) — Idealized Life-Cycle Experiments

Canonical baroclinic-wave experiments run with three different models (MPAS-Atmosphere, SpeedyWeather.jl, Dedalus).

ExperimentEngineReferenceStatus
Jablonowski & Williamson Dycore TestMPAS-A (x1.40962, ~240 km)JW2006, QJRMS 132, 2943–2975
Thorncroft LC1 / LC2SpeedyWeather.jl (T42, 15σ)Thorncroft et al. 1993, QJRMS 119, 17–55🔧
Polvani & Plumb Stratospheric RWBDedalusPolvani & Plumb (1992)📋
Vallis Barotropic DecaySpeedyWeather.jl (barotropic)Vallis (2004)📋
Schemm COSMO DiagnosticsCOSMO-LM archiveSchemm et al. (2013)📋

JW06 — What’s reproduced ✅

16-day hydrostatic baroclinic wave on the x1.40962 Voronoi mesh (64 ranks, no physics). Validated against the published reference: 850 hPa temperature at day 9 (cold/warm fronts), 250 hPa zonal wind (jet streak), surface pressure minimum evolution, Ertel PV on isentropic surfaces. Lat-lon regridding via convert_mpas barycentric interpolation.

Thorncroft LC1 / LC2 — In progress 🔧

LC1 (anticyclonic, equatorward wave breaking) and LC2 (cyclonic, poleward wave breaking) life cycles in SpeedyWeather.jl. LC1 uses J&W basic state with ∇⁶ hyperdiffusion; LC2 adds a barotropic cyclonic perturbation. Diagnostics in development: EKE(t), PV on θ surfaces, EP-flux cross-sections.


3. SFNO — Spectral Fourier Neural Operator

ItemDetails
PaperSpectral Fourier Neural Operator — arXiv: 2306.03838
Status🔧 Training infrastructure ready (PyTorch + torch-harmonics)
ContentsCheckpoints, literature PDF, training scripts for SFNO weather prediction

Spherical-harmonic neural operator that leverages the convolution theorem on the sphere to achieve equivariance without explicit spherical-harmonic transforms in every layer.


4. Spherical Attention — Equivariant Transformers on S²

ItemDetails
PaperAttention on the Sphere — Bonev, Rietmann, Paris, Carpentieri, Kurth (NVIDIA, 2024; arXiv: 2505.11157)
Status🔧 Training ready
Key ResultQuadrature-weighted geodesic attention reduces pole errors vs Euclidean attention on lat-lon grids
Experiments(i) Shallow-water 6-h forecast; (ii) ERA5 Z500 6-h real-data forecast, 128×256 Gaussian grid
ContentsS² vs R² Transformer implementations, checkpoints, MP4 animations, training data

5. Hilbert Transform — ERA5 Scale Decomposition & Analytic Signals

ItemDetails
PaperThe Analytic Signal and Hilbert Transform — Huang et al. (1998)
Status✅ Complete (all 15 phases of plan.md)
DataERA5 Z500, Oct–Nov 2012, 61 days, 1464 hourly snapshots
PipelineZonal wavenumber + temporal frequency filtering → 9 bands (3 spatial × 3 temporal) → Hilbert analytic signal → complex PCA
Outputs3.3 GB NetCDF, 9 PNG snapshots, 5.6 MB MP4 (9-panel animation)
ContentsDownload scripts, main pipeline, extended animation stitcher, literature PDFs