A Thousand-Fold Compression: Symbolic Equations That Explain Monsoon Rainfall Better Than Standard Physics Diagnostics

Symbolic distillation squeezes neural network skill into closed-form equations for South Asian monsoon precipitation, beating the standard buoyancy diagnostic by 1.5× in variance explained with ~1000× fewer parameters.

Imagine you're trying to write a recipe for a dish you've only ever watched a chef make. You have video footage of every ingredient going in, but the chef never told you the proportions. A neural network is like hiring a food scientist to taste the dish and reverse-engineer the recipe — they can tell you exactly what's in it, but their 'recipe' is a spreadsheet of 10,000 chemical measurements. Symbolic distillation is the step where you force that food scientist to hand you a recipe card that fits on an index card — and it still tastes 85% as good. The committed claim: closed-form precipitation equations derived from symbolic distillation of neural networks explain 1.5× more monsoon precipitation variance than CAPE-based buoyancy diagnostics (the current standard proxy), recovering ~85% of neural network skill with roughly a thousand times fewer parameters. The equations are not black boxes — they are readable, testable algebraic expressions. The method sits in a two-stage pipeline that the ML-for-climate community has been building toward. Stage one: train a neural network on 21 years (1999–2019) of June–August ERA5 reanalysis data to predict precipitation from vertical thermodynamic profiles (temperature, humidity, pressure across altitude). Crucially, the network learns a single scalar 'feature' for each variable — a vertical integral that isn't pinned to any particular pressure level, solving the nonlocal dependence problem. Stage two: feed those learned scalar features plus surface-state variables into PySR (a symbolic regression engine) to extract algebraic equations. The best equation uses about a thousand times fewer parameters than the neural network while retaining ~85% of its predictive skill. What makes the result interesting beyond compression is what the equations say physically. The standard diagnostic (a buoyancy measure related to CAPE) treats moisture and instability as additive — more moisture plus more instability equals more rain. The discovered equations treat them as competing controls. The heaviest precipitation occurs when the lower free troposphere is already humid, not simply when the atmosphere is unstable. Over ocean, sensible heating boosts precipitation mainly when near-surface air is both warm and moist, which also improves the predicted diurnal cycle of rainfall — a notoriously hard target. The ladder position is honest and informative. The named baseline is the widely used column-integrated buoyancy diagnostic (essentially a linearized CAPE variant). The new equation explains 1.5× as much precipitation variance. Against the full neural network, the equation recovers about 85% of skill. Classical baselines like raw CAPE or column-relative-humidity predictors are not beaten at the 'understanding why' level — they remain useful physically — but they are decisively beaten as precipitation predictors in this domain. Integrity is solid for this type of work but has the usual observational-science caveats. Training on ERA5 reanalysis means you're learning relationships baked into a data assimilation system, not direct observations. There is no independent physical experiment — the 'ground truth' is itself a model product. The 21-year dataset is split for train/validation/test, and the authors report out-of-sample metrics, but there is no pre-registration and no code availability is mentioned in the abstract. The comparison to the buoyancy diagnostic is head-to-head on the same data, which is about as fair as you can get in this setting. The successor question is revealing. The obvious next experiment is to test these equations inside a climate model — either as a convective trigger or as a diagnostic constraint on parameterized precipitation. The authors gesture toward this ('could guide data-driven precipitation models') but don't run it. Most likely read: this is being saved for the next paper, since coupling a symbolic equation into a GCM or ML weather model is a different engineering effort. The second missing experiment is out-of-distribution testing: do these equations hold for future-climate thermodynamic profiles, or do the learned features break under warming? That's the make-or-break question for climate applications and was likely beyond the scope of a methods-focused paper. The bigger field fight here is whether ML-for-weather/climate should aim for black-box skill (foundation models, GraphCast, Pangu-Weather) or interpretable hybrid approaches that preserve physical reasoning. This paper is a clear vote for interpretability — trading some skill for equations you can read, test, and argue with. If these equations hold up under climate-change thermodynamics, they become genuinely useful constraints. If they don't, they're fancy curve fits on reanalysis. That's the test to watch.