# Toward GPU-Resident Climate Models: A Feasibility Study on Lossy Compression for the Spherical Harmonic Transform's Communication Bottleneck

Preprint · Lorenzo Breschi, Flavio Vella

- Archive id: arxiv:2609.24294
- Published: 2026-09-21
- arXiv: https://arxiv.org/abs/2609.24294
- Source URL: https://arxiv.org/abs/2609.24294
- PDF: https://arxiv.org/pdf/2609.24294
- Provenance: arxiv

## Abstract

Operational pseudospectral atmospheric models such as the ECMWF Integrated Forecasting System (IFS) run today almost exclusively on CPUs; GPU ports are under active development but not yet used in production. These models rely on the Spherical Harmonic Transform (SHT). Each time-step requires forward and inverse SHTs, and both passes depend on global pencil transposition that redistribute multi-dimensional arrays across compute nodes. At large node counts these global collectives dominate wall-clock time. We investigate GPU-resident lossy compression, using representative fields from the DYAMOND high-resolution operational dataset as input, and combining measured GPU compression throughput with SimGrid network simulation, we show that ZFP at 16 bits per value (rate-16), about the same storage budget as float16, matches the communication-time reduction of float16 truncation while delivering approximately $\\mathbf{1600\\times}$ lower mean relative error. ZFP at 8 bits per value (rate-8) achieves approximately $\\mathbf{1.93\\times}$ the speedup of float16 while retaining $\\mathbf{4\\times}$ lower mean relative error.

_Abstract reproduced under CC0-1.0 (arXiv descriptive metadata)._

---
Space Kitty Archive · https://space-kitty.com · API: https://space-kitty.com/api/v1 · MCP: https://space-kitty.com/mcp · Premium: $0.20/request via x402 (Base Sepolia (testnet)). Metadata: arXiv (CC0), OpenAlex (CC0), Europe PMC (CC BY/CC0 works), ClinicalTrials.gov.