Thematic explorer

Physics AI & Causal World Models

The question I'm interested in is not whether learned models can predict physical systems; that result is increasingly established. It's what additional evidence is required before prediction can be interpreted as learned physical structure, causal understanding, or a model safe to use for counterfactual reasoning. I assembled this literature map to identify where those claims currently separate, which experiments discriminate between them, and which experiments still haven't been run.

80 papers · 8 themes

← All collections

80 papers shown

Observations, assimilation & the sensor boundary

A learned atmospheric model never sees the atmosphere. It sees instruments measuring transformations of it: radiances in satellite channels, reflectivity at radar wavelengths, station reports with inconsistent siting and quality control, radiosondes launched twice a day into a field that varies continuously. Classical data assimilation exists to solve that inverse problem and carries decades of accumulated machinery for observation operators, error covariances, and bias correction. This theme tracks what happens when learning is pushed into that layer: assimilation as generative sampling, forecast systems trained on raw observations instead of reanalysis, and neural observation operators standing in for radiative transfer. It also collects the papers that keep the sensor honest, because the truth almost every learned Earth-system model is trained and graded on is an analysis product that already contains a model. It also holds the measurements of how much of that observation stream is informative once error correlation is accounted for, which is the quantity at stake whenever raw data volume is offered as evidence of headroom.

Key threads
  • Reanalysis is not ground truth. It is a model-observation fusion, so a system trained and scored on it inherits the assimilating model's biases, and observation bias that varies with the state is not separable from model error without an anchor (Dee).
  • Learned assimilation arrived as a sampling problem: score-based and diffusion methods treat the posterior over states as something to generate rather than to optimize, and the learned forecaster supplies the adjoint that made 4D-Var expensive (Rozet; Xiao).
  • Observation-driven forecasting removed the analysis from the loop entirely. Aardvark and GraphDOP forecast from raw observations end to end, and GraphDOP's internal state recovers physical structure it was never given (Vaughan; Alexe; Lean).
  • The observation operator is now learnable, which turns a slow hand-derived map into a differentiable one and makes the surrogate's error distribution the useful artifact rather than its speed. Cross-sensor generalization is where it breaks: a model tuned to one instrument's spectral response and geometry degrades on another's, and that degradation looks exactly like model error unless the evaluation separates them (Howard; Stefanelli; Li; Jin).
  • Observation volume and assimilable information are different quantities. Past a density threshold, correlated observation error means additional data stops paying in the schemes operational centres actually run, and per-observation forecast impact is concentrated enough that some assimilated subsets measurably degrade the forecast (Liu; Cardinali).
  1. Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics

    Evensen · 1994 2922 cites

    Synthesis

    Replaces the linearized error-covariance evolution of the extended Kalman filter with a Monte Carlo ensemble, so the forecast error statistics are sampled from the nonlinear model itself rather than approximated around a trajectory.

    Why it matters The ensemble Kalman filter is the reference point every learned assimilation method is measured against, and the reason state estimation in geophysics is already a probabilistic problem.

  2. The interaction between model resolution, observation resolution and observation density in data assimilation: A one-dimensional study

    Liu · 2002 133 cites

    Synthesis

    Works out, in a one-dimensional analysis system, what an extra observation is worth once representativeness error and observation-error correlation are carried through the analysis-error formulation. Denser observations keep helping an optimal scheme, but for the suboptimal schemes operational centres run in practice, density past a threshold stops paying, and thinning observations with correlated error improves the analysis.

    Why it matters The quantitative form of the point that observation volume and assimilable information are different quantities, with the gap set by error correlation rather than by collection rate.

  3. Bias and data assimilation

    Dee · 2005 461 cites

    Synthesis

    Formalizes what happens when observations and the model are both biased: the analysis absorbs the bias into the state, variational bias correction can only work against an anchor of unbiased observations, and model bias and observation bias are not separately identifiable from the residuals alone.

    Why it matters The clearest statement of why an analysis product is not ground truth, which is the assumption almost every learned Earth-system model is trained and scored on.

  4. Monitoring the observation impact on the short-range forecast

    Cardinali · 2009 206 cites

    Synthesis

    Uses adjoint forecast sensitivity to attribute 24-hour forecast error reduction to individual observations across the full ECMWF observing system, checked against observing-system experiments over a summer and a winter month. The assimilated observations reduce forecast error in aggregate while specific subsets degrade it, traceable either to data quality or to suboptimality in the assimilation.

    Why it matters The field already measures the marginal forecast value of each observation, so an argument from raw data abundance has a per-observation measurement to answer to.

  5. Score-based Data Assimilation

    Rozet · 2023 arXiv

    Synthesis

    Trains a score model on short state trajectories, then samples the posterior over the full trajectory given arbitrary observations by conditioning at inference time, avoiding both linearization and the sequential structure of a filter.

    Why it matters Reframes assimilation as generative sampling rather than optimization, which is the template most learned-DA work has followed since.

  6. FengWu-4DVar: Coupling the Data-driven Weather Forecasting Model with 4D Variational Assimilation

    Xiao · 2023 arXiv

    Synthesis

    Couples a learned forecast model to a four-dimensional variational assimilation scheme, using the model's differentiability to get the adjoint for free and cycling the system on simulated observations.

    Why it matters Shows the learned forecaster can supply the piece 4D-Var is most expensive in, and marks the point where AI models stopped depending on someone else's analysis to run.

  7. Aardvark Weather: end-to-end data-driven weather forecasting

    Vaughan · 2024 arXiv

    Synthesis

    Replaces the entire operational pipeline (assimilation, forecast model, post-processing) with a single network trained from raw heterogeneous observations to station-level forecasts, at a fraction of the compute.

    Why it matters The first end-to-end demonstration that the analysis step is not structurally necessary, which changes what 'the model' is being evaluated on.

  8. GraphDOP: Towards skilful data-driven medium-range weather forecasts learnt and initialised directly from observations

    Alexe · 2024 arXiv

    Synthesis

    Plain-language abstract GraphDOP is an ECMWF forecast system trained and initialized exclusively from Earth-system observations, with no reanalysis inputs and no physics-based analysis anywhere in the loop. It learns relationships between observed quantities, such as brightness temperatures from polar-orbiting and geostationary satellites, and the geophysical quantities conventional instruments measure, and produces skilful predictions of weather parameters up to five days ahead.

    Motivation Every other learned forecaster is trained on reanalysis and initialized from an analysis, both produced by the physics-based assimilation system it is being compared against. That makes the usual claim structurally ambiguous: it cannot be determined how much of the skill belongs to the network and how much to the observation processing it inherited. Removing the analysis entirely is the only clean way to separate them.

    Methodology A graph model ingests observations in observation space and forms a latent representation of Earth-system state and its evolution. Training targets observed quantities rather than analyzed physical fields, so the loss never references a model-derived truth, and forecasts are issued directly in observation space. The physical state variables that other systems predict are never supplied to it.

    Results The system produces skilful forecasts of relevant weather parameters out to about five days from observations alone. The lead time is not the payoff; the representation is. A latent state built only from instrument measurements encodes coherent Earth-system dynamics, which follow-up probing work examined directly, and that turns the analysis step from a structural necessity into a design choice.

  9. Learning from nature: insights into GraphDOP's representations of the Earth System

    Lean · 2025 arXiv

    Synthesis

    Probes what GraphDOP's latent state contains and finds structure corresponding to physical fields the model was never given, recovered from observation sequences alone.

    Why it matters Turns the observation-only architecture into an argument about representation: the sensor stream carries enough information to reconstruct the physics.

  10. A Neural-Network Observation Operator for Weather Radar Reflectivity

    Stefanelli · 2025 arXiv

    Synthesis

    Learns the map from model state to radar reflectivity, replacing a forward operator whose microphysical assumptions are a known source of representativeness error in convective-scale assimilation.

    Why it matters Makes the observation operator itself a learned, differentiable component, which is where instrument-specific error is introduced and where it can now be diagnosed.

  11. Probabilistic Emulation of the Community Radiative Transfer Model

    Howard · 2025 arXiv

    Synthesis

    Emulates CRTM with calibrated uncertainty on the emulated radiances, so the surrogate reports where it disagrees with the physical forward model rather than silently substituting for it.

    Why it matters Sets the standard for surrogate observation operators: the useful artifact is the error distribution, not the speedup.

  12. WeatherReal: A Benchmark Based on In-Situ Observations for Evaluating Weather Models

    Jin · 2024 arXiv

    Synthesis

    Builds a quality-controlled global station dataset and scores AI and physical models against real in-situ measurements instead of reanalysis, documenting how much of the apparent skill ordering depends on which truth is used.

    Why it matters Makes the choice of verification truth an explicit experimental variable, which is the precondition for separating model error from analysis error.

  13. SPECTRA: Cross-Sensor Fine-Tuning for Satellite Observation Foundation Models

    Li · 2026 arXiv

    Synthesis

    Studies what breaks when a model pretrained on one instrument's spectral response and geometry is transferred to another, and what adaptation is needed to recover performance.

    Why it matters Names the failure mode that instrument scientists expect and machine-learning evaluations rarely isolate: degradation caused by the sensor, not the physics.

  14. Skillful high-resolution weather forecasting independent of physical models

    Zhao · 2026 arXiv

    Synthesis

    Builds a high-resolution forecasting system that does not depend on NWP analyses at any stage, initializing and training from observational sources.

    Why it matters Extends the observation-only line from medium-range global fields to the high-resolution regime where analysis dependence has been hardest to shed.

Learned weather & Earth-system prediction

The claim that settled between 2022 and 2024 is narrow and real: a neural network trained on reanalysis beats the operational deterministic forecast on standard headline scores at a fraction of the inference cost. FourCastNet made it plausible, Pangu-Weather and GraphCast made it undeniable under ECMWF's own verification, and GenCast moved the result into the probabilistic setting where operational value lives. What did not settle is what that skill is made of. These systems are initialized from analyses produced by physics-based assimilation, scored against those same analyses, and their advantage concentrates in the smooth, large-scale, well-observed part of the problem. This theme holds the capability line and the papers that mark its edges together, because separating fast high-skill prediction from physical understanding is the whole point.

Key threads
  • The 2022-2023 jump was architectural and then immediately commoditized: transformers on a lat-lon grid and message passing on an icosahedral mesh both beat IFS, after which the gains came from training strategy and scale rather than architecture. The scaling analysis that followed puts the return in training data rather than parameter count, and finds these models prefer width where language models prefer depth (Pathak; Bi; Lam; Yu).
  • Deterministic training produces blurry forecasts. The move to diffusion and to direct CRPS optimization is the field conceding that a mean-squared-error forecast is a conditional mean, not a forecast (Price; Lang).
  • Hybrid systems did not lose. NeuralGCM keeps a differentiable dynamical core with learned parameterizations, matches learned-model skill at medium range, and stays stable over climate-length runs; spectral nudging pushes the same convergence from the physical side (Kochkov; Polichtchouk).
  • The skill is conditional on an initialization pipeline the models do not own. Run them in an operational setting, from operational analyses, and the advantage narrows and becomes variable-specific (Ben Bouallègue; Bonavita).
  • Learned models under-represent small scales and error growth, which is how a system can win at day-5 RMSE while being wrong about how uncertainty propagates (Bonavita; Chattopadhyay).
  1. Skillful high-resolution weather forecasting independent of physical models

    Zhao · 2026 arXiv

    Synthesis

    Reports skilful high-resolution forecasts without a physics-based model anywhere in the pipeline, at resolutions where learned models previously relied on physical downscaling or analysis input.

    Why it matters Tests whether the 'AI models only work because NWP feeds them' objection survives at the scales that matter operationally.

  2. FourCastNet: A Global Data-driven High-resolution Weather Model using Adaptive Fourier Neural Operators

    Pathak · 2022 arXiv

    Synthesis

    Trains an adaptive Fourier neural operator on ERA5 to produce global forecasts at 0.25 degrees, running a full forecast in seconds and generating large ensembles at negligible marginal cost.

    Why it matters The first system to make the cost argument undeniable, which is what put learned forecasting on operational centres' agendas.

  3. Accurate medium-range global weather forecasting with 3D neural networks (Pangu-Weather)

    Bi · 2022

    Synthesis

    A 3D Earth-specific transformer with a hierarchical temporal aggregation strategy that beats the operational IFS on deterministic scores across most upper-air and surface variables.

    Why it matters The result that ended the argument about whether a purely learned model could beat the operational deterministic forecast.

  4. Learning skillful medium-range global weather forecasting (GraphCast)

    Lam · 2022 arXiv

    Synthesis

    Plain-language abstract GraphCast is a graph neural network trained directly on ERA5 reanalysis that predicts hundreds of weather variables globally at 0.25 degree resolution out to 10 days, producing a full forecast in under a minute. Against the most accurate operational deterministic system it wins on 90% of 1380 verification targets, and the same forecasts improve prediction of severe events including tropical cyclone tracks, atmospheric rivers, and temperature extremes.

    Motivation Numerical weather prediction improves by adding compute to a hand-derived model of the governing equations; it cannot absorb the historical record directly. Decades of reanalysis exist and encode how the atmosphere actually evolved, including the errors the operational system keeps repeating. The open question was whether a model trained on that record could match a system refined over forty years, and whether it would produce a usable forecast or a smoothed field that only scores well.

    Methodology An encode-process-decode architecture maps the 0.25 degree lat-lon state onto a multi-scale icosahedral mesh, runs 16 rounds of message passing over mesh edges at several resolutions so information crosses the globe within a single step, and decodes back to the grid. The model is trained autoregressively, with the number of rolled-out steps increased during training so it is optimized on the trajectory it will actually produce rather than on a single step. Input state is two consecutive time steps of five surface and six atmospheric variables across 37 pressure levels.

    Results GraphCast outperforms HRES on 90% of 1380 combinations of variable, level and lead time, with the largest margins in the upper troposphere and at longer leads. Severe-event performance is evaluated as a diagnostic on the same forecasts rather than as a trained objective, and improves on cyclone track error, atmospheric river vapor transport, and heat and cold extremes. Inference takes under a minute against roughly an hour of supercomputer time. The cost shows up in the spectrum: fields are smoother than the analysis and small-scale variance is deficient, which later operational evaluations quantified.

  5. GenCast: Diffusion-based ensemble forecasting for medium-range weather

    Price · 2023 arXiv

    Synthesis

    Plain-language abstract GenCast is a conditional diffusion model that generates ensembles of 15-day global forecasts at 0.25 degrees and 12-hour steps across more than 80 surface and atmospheric variables, producing a full ensemble in about eight minutes. It has greater skill than the ECMWF ensemble on 97.4% of 1320 evaluated targets and gives better predictions of extremes, tropical cyclone tracks, and wind power production.

    Motivation Deterministic learned models are trained on mean-squared error, whose optimum is the conditional mean, so they output a blurred best guess rather than a sample of what could happen. Nearly every weather decision that matters is made against a distribution: whether to warn, whether to commit reserve generation, whether to divert. A system that cannot generate plausible individual trajectories is not competing in the event that operational value is defined over.

    Methodology Each 12-hour step is treated as conditional generation. Given the two previous states, the model samples the next state by iterative denoising with a graph-transformer denoiser on a spherical mesh, then autoregresses to build a 15-day trajectory. Ensemble members differ only by their noise draw, so spread arises from the generative process itself rather than from a separately designed initial-condition perturbation scheme. Training is on ERA5.

    Results GenCast beats ENS by CRPS on 97.4% of 1320 targets, and the margin holds for extremes and for tropical cyclone position rather than disappearing in the tails. Generated fields carry realistic small-scale structure instead of the smoothing that deterministic training induces, which is the property the deterministic scoring regime had been hiding. The wind power evaluation is the closest this literature comes to reporting decision-relevant value, and it remains a proxy for a decision rather than a decision.

  6. Neural general circulation models for weather and climate (NeuralGCM)

    Kochkov · 2023 arXiv

    Synthesis

    Plain-language abstract NeuralGCM is a general circulation model whose dynamical core is a differentiable solver for large-scale atmospheric dynamics and whose subgrid physics is a neural network trained end to end through that solver. It matches the best learned models for 1-10 day deterministic forecasts and the ECMWF ensemble for 1-15 day probabilistic forecasts, and with prescribed sea surface temperatures it integrates stably for decades, tracking global mean temperature and producing realistic tropical cyclone frequency and trajectories at 140 km resolution.

    Motivation Purely learned forecasters had two unresolved failures: ensemble quality and long-run stability. A model that drifts or blurs after a few weeks cannot be used for climate, which is where most of the scientific value of an atmospheric model lives. Conventional GCMs have the complementary problem, since their cost and their hand-tuned parameterizations dominate both their expense and their bias structure. Each side's weakness is the other's strength, and the question was whether the coupling could be trained rather than assembled by hand.

    Methodology The dynamical core solves the primitive equations on a sphere inside a differentiable framework, so gradients propagate through many simulation steps into the learned module that supplies tendencies for unresolved processes. Training optimizes the learned physics against reanalysis along the model's own rollout rather than against offline parameterization targets, which is what keeps the learned component consistent with the solver it lives inside. Stochastic components in the learned physics generate ensemble members.

    Results Deterministic skill is competitive with GraphCast and Pangu-Weather at 1-10 days, and ensemble CRPS is competitive with ECMWF ENS through 15 days. Multi-decade integrations stay stable and track observed global mean temperature under prescribed SSTs. Behavior appears that was never trained for, including realistic tropical cyclone frequency and tracks, which is the strongest available evidence that the hybrid retains physical structure rather than pattern-matching. Compute cost is orders of magnitude below a conventional GCM at comparable resolution.

  7. Aurora: A Foundation Model of the Atmosphere

    Bodnar · 2024 arXiv

    Synthesis

    Pretrains a single atmospheric model on a heterogeneous mixture of reanalyses, forecasts and climate simulations, then fine-tunes it to tasks including air pollution and ocean waves where dedicated operational systems are the baseline.

    Why it matters The first credible demonstration that Earth-system pretraining transfers across variables and domains rather than only across lead times.

  8. AIFS - ECMWF's data-driven forecasting system

    Lang · 2024 arXiv

    Synthesis

    ECMWF's own graph-transformer forecast system, developed inside the operational verification and initialization environment rather than against a published reanalysis benchmark.

    Why it matters Marks the transition from external claims about beating IFS to an operational centre running a learned model in its own suite.

  9. Scaling Laws of Global Weather Models

    Yu · 2026 arXiv

    Synthesis

    Fits empirical scaling relations between validation loss and model size, dataset size, and compute across several learned global weather models. Aurora shows the strongest data scaling, where a tenfold larger training set cuts validation loss by up to 3.2x; under a fixed compute budget, spending on data beats spending on parameters; and these models favor width over depth, the opposite of the language-model result.

    Why it matters Puts measured exponents under the proposition that physical prediction improves with scale, and locates the return in training data rather than model size.

  10. The Rise of Data-Driven Weather Forecasting: A First Statistical Assessment of Machine Learning-Based Weather Forecasts in an Operational-Like Context

    Ben Bouallègue · 2024 131 cites

    Synthesis

    Evaluates Pangu-Weather and GraphCast inside ECMWF's operational verification framework, initialized from operational analyses and scored against the metrics and baselines the centre itself uses, finding real headline gains alongside deficiencies in near-surface fields, small scales, and extremes.

    Why it matters The independent replication that turned the published claims into a calibrated statement about where the skill is and is not.

  11. On Some Limitations of Current Machine Learning Weather Prediction Models

    Bonavita · 2024 21 cites

    Synthesis

    Documents that learned models produce spatially over-smoothed fields with deficient variance at small scales, that their apparent RMSE advantage is partly an artifact of that smoothing, and that they depend on an assimilation system they neither model nor own.

    Why it matters The most-cited statement of the gap between winning on aggregate scores and being a forecast model, written from inside the operational centre.

  12. Hybrid weather prediction using spectral nudging

    Polichtchouk · 2026 arXiv

    Synthesis

    Nudges the large scales of a physical model toward a learned forecast while leaving the small scales to the physics, so each component is used where it is strongest.

    Why it matters Represents the convergence the field arrived at after the pure-learned versus physical argument: the interesting question is the coupling, not the winner.

  13. Challenges of learning multi-scale dynamics with AI weather models: Implications for stability and one solution

    Chattopadhyay · 2023 arXiv

    Synthesis

    Analyses why data-driven emulators of multiscale chaotic systems become unstable or lose small-scale variance during long autoregressive rollouts, and what training and spectral interventions delay the failure.

    Why it matters Explains the mechanism behind the smoothing that operational evaluations later measured, rather than treating it as an implementation detail.

Foundation models & learned physical operators

Operator learning made a specific promise: learn the solution map of a family of equations, not one solution, so a trained model transfers across initial conditions, resolutions, and eventually across equations. FNO and DeepONet are the two canonical formulations, and graph-network simulators made the same argument for particle and mesh dynamics. The physics foundation model generalizes the claim again, to pretraining on many systems at once and transferring to unseen ones. Real transfer results exist and they are narrower than the framing suggests: they show up most clearly in the low-data regime, and the speedup literature they sit next to has a documented baseline problem. This theme keeps two things apart that get conflated constantly: what has been shown to transfer, and what has only been shown to fit.

Key threads
  • Resolution invariance is the reusable property. FNO gets it by parameterizing the kernel in Fourier space and DeepONet by construction; accuracy on one grid is not the claim (Li; Lu).
  • Multi-physics pretraining does transfer to held-out systems, and the honest version of the result is that it transfers in the low-data regime, where sample efficiency rather than asymptotic accuracy is the gain (McCabe; Herde).
  • 'Foundation model' names two different things: a pretrained backbone with demonstrated downstream transfer, and a large task-specific simulator with a label. Only the first is a testable claim, and Aurora is the case that argues it by beating dedicated operational systems on tasks outside forecasting (Bodnar).
  • Reported speedups over numerical solvers are systematically inflated. A review of the fluid-PDE literature found most comparisons run against baselines weaker than the domain standard, with negative results underreported (McGreivy).
  • The bottleneck moved from architecture to data, which is why the most consequential recent artifacts in this line are datasets with common interfaces rather than models (Ohana).
  1. Aurora: A Foundation Model of the Atmosphere

    Bodnar · 2024 arXiv

    Synthesis

    Uses a 1.3B-parameter 3D Swin backbone with perceiver-style encoders that accept heterogeneous inputs at different resolutions and variable sets, so downstream tasks reuse the pretrained representation rather than the architecture alone.

    Why it matters The reference point for what an Earth-system foundation model has to show: fine-tuning that beats task-specific operational systems, not just other neural networks.

  2. DeepONet: Learning nonlinear operators for identifying differential equations based on the universal approximation theorem of operators

    Lu · 2019 arXiv

    Synthesis

    Proves and instantiates an operator universal-approximation result with a branch-trunk architecture that learns mappings between function spaces from input-function and query-location pairs.

    Why it matters One of the two canonical formulations of operator learning, and the one that makes the function-space claim explicit rather than architectural.

  3. Fourier Neural Operator for Parametric Partial Differential Equations

    Li · 2020 arXiv

    Synthesis

    Parameterizes the integral kernel directly in Fourier space, giving a resolution-invariant operator that can be trained at low resolution and evaluated at high resolution, with large speedups over spectral solvers on Navier-Stokes.

    Why it matters The architecture that made 'learn the solution operator, not the solution' concrete, and the direct ancestor of FourCastNet.

  4. Learning to Simulate Complex Physics with Graph Networks

    Sanchez-Gonzalez · 2020 arXiv

    Synthesis

    Learns particle-based dynamics as message passing over a neighborhood graph, generalizing across materials, geometries, and orders of magnitude more particles than seen in training.

    Why it matters Established the inductive bias (local interactions on a graph) that GraphCast later carried onto the sphere.

  5. ClimaX: A foundation model for weather and climate

    Nguyen · 2023 arXiv

    Synthesis

    Pretrains on heterogeneous CMIP6 climate simulations with variable tokenization and aggregation, then fine-tunes to forecasting, downscaling, and climate projection tasks with different variables and resolutions.

    Why it matters The first paper to argue that Earth-system tasks share a backbone, and the one that fixed the transfer-across-tasks evaluation shape.

  6. Multiple Physics Pretraining for Physical Surrogate Models

    McCabe · 2023 arXiv

    Synthesis

    Pretrains one transformer on several distinct PDE systems simultaneously via shared normalized field embeddings, then transfers to physics not in the pretraining set, beating from-scratch training in the low-data regime.

    Why it matters The most direct test of whether pretraining across physics produces reusable structure, and the honest answer it gives is 'yes, in the small-data regime'.

  7. Poseidon: Efficient Foundation Models for PDEs

    Herde · 2024 arXiv

    Synthesis

    Scales multi-physics pretraining with an operator transformer and an all2all training scheme that exploits the semigroup structure of time-dependent PDEs, reporting large sample-efficiency gains on unseen equations.

    Why it matters The strongest quantitative transfer result in operator learning, and a demonstration that the useful pretraining signal comes from time-evolution structure.

  8. Weak baselines and reporting biases lead to overoptimism in machine learning for fluid-related partial differential equations

    McGreivy · 2024 arXiv

    Synthesis

    Plain-language abstract A systematic review of machine-learning papers claiming to beat standard numerical methods on fluid-related PDEs finds that 79% (60 of 76) compare against a weak baseline, and documents widespread outcome reporting bias and publication bias. The conclusion is that reported speedups in this literature are systematically overoptimistic.

    Motivation The value proposition of a learned PDE solver is entirely comparative: it has to be faster than a numerical method at equal accuracy. That makes the baseline the whole experiment. The baseline is chosen by authors who are usually not specialists in the numerical method they are displacing, and who have no incentive to tune it well.

    Methodology The authors define a weak baseline operationally, as a numerical method that is lower-order, unoptimized, run at an unfavorable resolution, or otherwise not the standard of practice for the problem, then screen every qualifying article that claims to outperform a numerical method on a fluid PDE. They separately look for reporting-bias signatures, including selective reporting of favorable metrics and the near-absence of published negative comparisons.

    Results 79% of qualifying claims rest on a weak baseline. Reporting biases are widespread, and both problems trace to familiar mechanisms: researcher degrees of freedom in choosing the comparison and a publication process that does not absorb negative results. The proposed remedies are cultural and structural rather than technical, which is what makes the paper generalize: the same claim structure appears wherever a learned model is said to outperform a domain-standard method, including weather.

  9. The Well: a Large-Scale Collection of Diverse Physics Simulations for Machine Learning

    Ohana · 2024 arXiv

    Synthesis

    Releases 15TB of numerical simulations across 16 physical systems, from biological flows to supernova explosions, with a common interface and baseline results.

    Why it matters Reframes the bottleneck as data rather than architecture, and supplies the diversity that makes cross-physics transfer claims testable.

Evaluation, uncertainty & physical validity

This is the theme that decides whether the others are true. Forecast verification settled its foundations long before machine learning arrived: Murphy's decomposition of what goodness means, proper scoring rules that cannot be improved by hedging, calibration and sharpness as separate axes, scale-aware spatial scores for fields whose point values are not the quantity of interest. Learned forecasting largely re-derived a headline-RMSE culture on top of that work. The papers here are the correction, and they share a structure: an aggregate score that hides the failure that matters, a benchmark that rewards the wrong behavior, a tail the training distribution barely contains, a comparison whose baseline or evaluation period is doing the work.

Key threads
  • Murphy's three kinds of goodness (consistency, quality, value) still sorts most disputes, and a headline score answers only the second, partially.
  • A deterministic score cannot certify a probabilistic system. Proper scoring rules, calibration, and sharpness are the minimum, they are not interchangeable, and models tied on RMSE can differ sharply in ensemble behavior (Gneiting; Brenowitz).
  • Aggregate skill and extreme-event skill dissociate. Models that beat IFS on mean scores fail on unseen-intensity cyclones, under-represent damaging near-surface winds in the one storm that mattered, and lose much of their margin when the evaluation is restricted to the tails (Sun; Charlton-Perez; Gabler).
  • Benchmarks reward what they measure, and a headline gain can come from the evaluation period rather than from the model: scored against a climatology the atmosphere has already drifted away from, every system trained on the past looks better than it is (Landsberg).
  • Independent standardized evaluation is now an explicit community ask, because self-reported comparisons against different baselines, initializations, and truth datasets are not comparable (Rasp; Ullrich; McGreivy).
Open gaps
  • No benchmark attributes a learned model's error between the observing system, the assimilation, and the forecast model, so 'the model is wrong' and 'the analysis is wrong here' remain indistinguishable in every reported score.
  • Extreme-event evaluation has no agreed power analysis: nobody states how many events are needed before a tail-behavior improvement is a result rather than a sample.
  1. Weak baselines and reporting biases lead to overoptimism in machine learning for fluid-related partial differential equations

    McGreivy · 2024 arXiv

    Synthesis

    Plain-language abstract A systematic review of machine-learning papers claiming to beat standard numerical methods on fluid-related PDEs finds that 79% (60 of 76) compare against a weak baseline, and documents widespread outcome reporting bias and publication bias. The conclusion is that reported speedups in this literature are systematically overoptimistic.

    Motivation The value proposition of a learned PDE solver is entirely comparative: it has to be faster than a numerical method at equal accuracy. That makes the baseline the whole experiment. The baseline is chosen by authors who are usually not specialists in the numerical method they are displacing, and who have no incentive to tune it well.

    Methodology The authors define a weak baseline operationally, as a numerical method that is lower-order, unoptimized, run at an unfavorable resolution, or otherwise not the standard of practice for the problem, then screen every qualifying article that claims to outperform a numerical method on a fluid PDE. They separately look for reporting-bias signatures, including selective reporting of favorable metrics and the near-absence of published negative comparisons.

    Results 79% of qualifying claims rest on a weak baseline. Reporting biases are widespread, and both problems trace to familiar mechanisms: researcher degrees of freedom in choosing the comparison and a publication process that does not absorb negative results. The proposed remedies are cultural and structural rather than technical, which is what makes the paper generalize: the same claim structure appears wherever a learned model is said to outperform a domain-standard method, including weather.

  2. What Is a Good Forecast? An Essay on the Nature of Goodness in Weather Forecasting

    Murphy · 1993 461 cites

    Synthesis

    Separates three kinds of forecast goodness: consistency with the forecaster's true belief, quality as correspondence with observation, and value as the benefit realized by the user, and shows they can move in opposite directions.

    Why it matters Thirty years before the AI weather literature re-derived it, this is the framework that says a headline score is one of three axes and not the one that matters most.

  3. Strictly Proper Scoring Rules, Prediction, and Estimation

    Gneiting · 2007

    Synthesis

    Develops the theory of scoring rules that are maximized only by reporting the true predictive distribution, covering CRPS, the logarithmic score, and the Brier score as instances of one framework.

    Why it matters The formal reason CRPS is the right target for ensemble systems and why a score that can be improved by hedging cannot certify a probabilistic model.

  4. Probabilistic Forecasts, Calibration and Sharpness

    Gneiting · 2007

    Synthesis

    Proposes the paradigm of maximizing sharpness subject to calibration, distinguishes probabilistic, exceedance and marginal calibration, and gives diagnostic tools including the PIT histogram.

    Why it matters Makes 'the ensemble is calibrated' and 'the ensemble is useful' two separate measurable claims, which is exactly the distinction diffusion-based forecasting is now argued over.

  5. Fuzzy verification of high-resolution gridded forecasts: a review and proposed framework

    Ebert · 2008

    Synthesis

    Reviews neighborhood-based verification methods and organizes them by the decision model each implies, showing how a high-resolution forecast can be penalized by point matching while being more useful than a smooth one.

    Why it matters The canonical answer to the double-penalty problem, and the reason a smoother field can score better while being physically worse.

  6. WeatherBench 2: A benchmark for the next generation of data-driven global weather models

    Rasp · 2024 78 cites

    Synthesis

    Defines a common evaluation framework with fixed ground truth, baselines and headline scores across deterministic and probabilistic metrics, and reports every major learned model against operational IFS on the same footing.

    Why it matters The shared scoreboard that made cross-paper comparison possible, and the artifact whose choices now shape what the field optimizes.

  7. A Practical Probabilistic Benchmark for AI Weather Models

    Brenowitz · 2024 arXiv

    Synthesis

    Shows that deterministic learned models with equal RMSE can have very different spread and skill when used to generate ensembles, and proposes a lagged-ensemble protocol that exposes the difference cheaply.

    Why it matters Demonstrates concretely that deterministic scores do not rank models the way probabilistic use does, using models the field had already declared equivalent.

  8. Do AI models produce better weather forecasts than physics-based models? A quantitative evaluation case study of Storm Ciaran

    Charlton-Perez · 2024

    Synthesis

    Compares learned and physical models on a single damaging windstorm and finds the learned models captured the synoptic evolution well while under-representing the near-surface wind extremes that caused the damage.

    Why it matters The case study that shows aggregate skill and the property a forecast user cares about can come apart on the one event that matters.

  9. Can AI weather models predict out-of-distribution gray swan tropical cyclones?

    Sun · 2024 arXiv

    Synthesis

    Plain-language abstract The authors retrain FourCastNet on 1979-2015 ERA5 with Category 3-5 tropical cyclones removed, either globally or from one ocean basin, then test the resulting models on Category 5 storms from 2018-2023. Removing the strong storms globally destroys the ability to forecast them; removing them from a single basin leaves some skill in that basin. The models generalize across basins but do not extrapolate to intensities they have never seen.

    Motivation The events that justify a forecast system are rare enough to be nearly absent from any training record. If a learned model can only interpolate within the intensities it has seen, its use for early warning on the highest-impact events rests on an untested assumption, and standard evaluation, which averages over ordinary weather, cannot detect the failure.

    Methodology Independent versions of FourCastNet are trained under three conditions: all data, Category 3-5 cyclones removed globally, and removed from the North Atlantic or Western Pacific only. Every version is evaluated both on general global weather and on 2018-2023 Category 5 storms. This is a controlled removal experiment rather than a post-hoc evaluation, so the counterfactual model actually exists.

    Results All versions score similarly on general global weather, which is the finding that matters most: the aggregate metric is blind to the ablation. The globally ablated model cannot forecast Category 5 storms. The basin-ablated models retain some skill in the withheld basin, so physical structure does transfer across regions even though regional information is implicitly encoded in the inputs. The authors expect the result to hold for other current models, which share the same learning strategy.

  10. Recommendations for Comprehensive and Independent Evaluation of Machine Learning-Based Earth System Models

    Ullrich · 2025

    Synthesis

    Sets out what an evaluation has to include before a learned Earth-system model can be trusted: physical consistency checks, out-of-sample and out-of-distribution testing, and independence between the developers and the evaluators.

    Why it matters The community's explicit statement that self-reported benchmark wins are not sufficient evidence, and the closest thing to an agreed standard.

  11. Forecasting the Future with Yesterday's Climate: Distribution Shift in AI Weather Prediction

    Landsberg · 2025 arXiv

    Synthesis

    Examines how a model trained on a historical reanalysis period behaves as the climate it was trained on drifts away from the climate it is forecasting, and how much apparent skill change is attributable to the evaluation period rather than the model.

    Why it matters Identifies a confound that grows silently over time and that no standard benchmark protocol currently controls for.

  12. Do AI weather models miss extremes?

    Gabler · 2026 arXiv

    Synthesis

    Tests learned models specifically on extreme events across variables and regions, quantifying how much of the reported skill advantage survives when the evaluation is restricted to the tails.

    Why it matters Converts a widely repeated suspicion about smoothing and extremes into a measured, variable-by-variable result.

Causal structure & mechanistic understanding

The interesting question is whether prediction pressure forces a model to learn mechanism, and the literature supports both answers depending on which claim is being made. Five things get collapsed into the word causal and they carry different burdens of proof: a predictive relationship, an internal representation that encodes a physical variable, an identified causal structure, a valid counterfactual, and a model that supports intervention. The causal-inference literature for dynamical systems supplies the standards, and it is explicit that its assumptions are strong. The interpretability work on learned weather models supplies the current evidence, and it cuts both ways: these models do develop internal structure aligned with physical fields they were never given, and they also fail dynamical tests that a physical model passes by construction. Outside weather the same question has been posed as a controlled experiment, on orbital mechanics and on simulated collisions, and both times the model that predicted well did not carry the law.

Key threads
  • Causal discovery in coupled nonlinear dynamics needs different machinery than in stochastic time series, and the methods that work depend on assumptions (causal sufficiency, stationarity, adequate sampling) that Earth-system data routinely breaks (Sugihara; Runge).
  • Invariance is the bridge that survives distribution shift: relationships that hold across environments are the ones with a claim to being causal, and the variables themselves have to be learned before that claim can be made (Peters; Schölkopf).
  • Evidence for mechanism is accumulating from the inside: probing recovers physically meaningful structure from internal activations, sensitivity analysis on a well-studied cyclone reproduces known precursors, and an observation-only model builds Earth-system representations from instrument data alone (Lean; Bano Medina; Tempest).
  • Evidence against is equally concrete: learned models do not reproduce the upscale error growth that defines atmospheric predictability, and they get some classical dynamical responses qualitatively wrong while getting others right (Selz; Hakim).
  • Being excellent at predicting a chaotic system and modeling it causally are different targets, and where the test has been run directly it returns a negative: a model trained on orbital trajectories predicts orbits accurately yet fails to apply Newtonian mechanics when adapted, and video models trained on clean 2D mechanics match the nearest training example rather than abstracting the rule (Kori; Vafa; Kang).
Open gaps
  • Which internal representations of a learned forecaster are load-bearing for its skill, and which are correlates that probing can find but ablation does not disturb?
  1. Learning from nature: insights into GraphDOP's representations of the Earth System

    Lean · 2025 arXiv

    Synthesis

    Shows a model trained only on instrument measurements develops internal representations aligned with unobserved physical variables, examined through targeted probing rather than output skill.

    Why it matters The strongest current evidence on the side that prediction pressure induces physically meaningful internal structure.

  2. Detecting Causality in Complex Ecosystems

    Sugihara · 2012 1222 cites

    Synthesis

    Introduces convergent cross mapping, which detects causal coupling in deterministic dynamical systems by testing whether one variable's attractor reconstruction can recover the other, and shows Granger causality fails in exactly this regime.

    Why it matters Establishes that causal inference in coupled nonlinear dynamics needs a different test than in stochastic time series, which is the setting every Earth-system causal claim lives in.

  3. Causal inference using invariant prediction: identification and confidence intervals

    Peters · 2015 arXiv

    Synthesis

    Exploits the fact that the conditional distribution of a target given its true causal parents is invariant across environments, turning causal discovery into a hypothesis test with confidence statements about the recovered parent set.

    Why it matters Supplies the operational definition of causal that survives distribution shift, which is the only definition that helps when a model is deployed outside its training climate.

  4. Detecting and quantifying causal associations in large nonlinear time series datasets

    Runge · 2019 436 cites

    Synthesis

    Introduces PCMCI, a two-stage conditional-independence method that controls false positives in high-dimensional autocorrelated time series where standard conditioning inflates them, with validation on nonlinear synthetic systems and climate data.

    Why it matters The method the Earth-system causal-discovery literature runs in practice, and the one whose assumptions define what its results can and cannot support.

  5. Towards Causal Representation Learning

    Schölkopf · 2021 arXiv

    Synthesis

    Argues that the variables over which causal structure is defined must themselves be learned from low-level observations, and connects that requirement to out-of-distribution generalization, transfer, and the independent-mechanisms principle.

    Why it matters The programmatic statement of why 'the model learned physics' is a representation claim, not a prediction claim, and what would be needed to establish it.

  6. Can Artificial Intelligence-Based Weather Prediction Models Simulate the Butterfly Effect?

    Selz · 2023 22 cites

    Synthesis

    Applies tiny localized perturbations to a learned model's initial state and finds the error growth does not exhibit the upscale propagation from convective scales that defines intrinsic predictability limits in the real atmosphere.

    Why it matters A direct test of whether high forecast skill implies correct dynamics, and the answer for this property is no.

  7. Dynamical Tests of a Deep Learning Weather Prediction Model

    Hakim · 2024 21 cites

    Synthesis

    Subjects a learned model to classical dynamical-meteorology experiments (steady tropical heating, extratropical cyclone development, geostrophic adjustment) and finds it reproduces some canonical responses while getting others qualitatively wrong.

    Why it matters Introduces the right evaluation instrument for mechanism claims: not held-out forecasts, but controlled experiments the physics already has known answers to.

  8. Are AI Data-Driven Weather Models Learning Atmospheric Physics? Initial Condition Sensitivity Analysis of Cyclone Xynthia

    Bano Medina · 2025

    Synthesis

    Uses gradient-based sensitivity analysis on a learned model during a well-studied explosive cyclogenesis case and compares the resulting influence patterns to the physical precursors identified in the meteorological literature.

    Why it matters Tests physical understanding through the model's own sensitivities rather than its outputs, which is the closest available proxy for asking what it is using.

  9. A Mechanistic Interpretability Tool for AI Weather Models

    Tempest · 2026 arXiv

    Synthesis

    Builds tooling to inspect internal circuits and activations of learned weather models and attribute predictions to identifiable internal computations rather than to input regions alone.

    Why it matters Moves the interpretability of physical models from saliency maps toward the circuit-level evidence that mechanism claims would require.

  10. A Unifying Perspective on Causal World Models

    Kori · 2026 arXiv

    Synthesis

    Formalizes the relationship between learned world models and structural causal models, setting out what a world model must satisfy before its rollouts can be read as interventions rather than as conditional predictions.

    Why it matters Names the exact gap between the world-model literature and the causal literature, which is the gap most 'causal AI' claims step over.

  11. What Has a Foundation Model Found? Using Inductive Bias to Probe for World Models

    Vafa · 2025 arXiv

    Synthesis

    Plain-language abstract Foundation models are built on the premise that predicting a domain well eventually amounts to understanding it, the way Kepler's accurate orbit predictions preceded Newton's laws. This paper tests that premise by adapting trained models to synthetic tasks drawn from a known world model and measuring whether their inductive bias matches it. Models trained on orbital trajectories predict orbits accurately and still fail to apply Newtonian mechanics on new physics tasks.

    Motivation Sequence prediction is the only training signal most foundation models receive, and the argument for scaling them rests on the claim that enough predictive pressure forces the underlying structure to be learned. Evaluating that claim on prediction accuracy is circular, since accuracy is what was optimized. What is needed is a test that asks what the model would do on a task it was not trained for, in a domain where the true generating structure is known in advance.

    Methodology The authors introduce an inductive bias probe: generate synthetic datasets from a postulated world model, adapt the foundation model to them, and measure whether its adaptation behaves the way a system carrying that world model would. Because the world model is chosen by the experimenter, the correct extrapolation is known exactly, which turns a question about internal understanding into a measurable comparison. The probe is run across several domains, including one where the models are trained on orbital trajectories and then adapted to force-law prediction.

    Results Models excel at their training tasks and fail to develop inductive biases aligned with the underlying world model. The orbital-mechanics case is the sharpest. Models that predict trajectories well consistently fail to apply Newtonian mechanics when adapted, and further analysis shows they behave as if they had assembled task-specific heuristics that do not generalize. Predictive accuracy on a physical system, by itself, did not yield the law generating it.

  12. How Far is Video Generation from World Model: A Physical Law Perspective

    Kang · 2024 arXiv

    Synthesis

    Trains diffusion video models on 2D mechanics with an unlimited supply of simulated data, and separates in-distribution, combinatorial, and out-of-distribution generalization. Generalization is perfect in-distribution and scales measurably for combinatorial cases but fails out of distribution, and the failure is diagnostic: the models match the nearest training example, prioritizing color over size over velocity over shape rather than abstracting a rule.

    Why it matters Scaling a generative predictor on clean, unlimited physics data did not recover the law that generated it, which bounds what observational pretraining alone can be expected to produce.

  13. Causal Confusion in Imitation Learning

    de Haan · 2019 arXiv

    Synthesis

    Behavioral cloning fits a discriminative model with no access to the causal structure of the expert-environment interaction, which produces causal misidentification: giving the policy more information can make it worse, because confounded cues that correlate with expert actions are easier to fit than the causes. Resolving it takes targeted intervention, either environment interaction or expert queries.

    Why it matters The clean demonstration that fitting behaviour is not learning what causes what, and that adding observations can move a policy in the wrong direction.

Forecast to decision

A forecast has no value until something is decided with it, and the map from score to value is not monotone. The meteorological version of this argument is old and quantitative: relative economic value curves show the ranking of two forecast systems flips depending on the user's cost-loss ratio, so better is a property of the pair, not of the forecast. The machine-learning version arrived independently as decision-focused learning, where the training loss is downstream decision regret rather than prediction error. Both literatures land on the same uncomfortable point from opposite directions: reducing prediction error in the wrong place leaves decisions unchanged, and can make them worse. Almost none of the learned-weather literature reports the downstream value of its own improvement.

Key threads
  • Relative economic value is user-specific. The same ensemble beats the deterministic forecast for some cost-loss ratios and loses for others, at the same lead time and the same headline skill (Richardson; Palmer).
  • Training on the decision objective beats training on accuracy when the optimization problem is the bottleneck, and a model with worse predictive error can produce strictly better decisions (Donti; Elmachtoub).
  • The decision-focused gain is conditional, not universal: recent work characterizes when a plain accuracy-trained predictor is already decision-optimal (Liu).
  • The problems in this literature are small and well-specified, with known cost structures. Operational weather decisions are neither, which is where the gap between the two traditions sits.
  • Where the value has been computed at industrial scale it belongs to the system, not to the forecast: the same wind-forecast improvement is worth different amounts to different grids depending on generation mix and ramp capability, and the cost-loss framework underneath most value claims assumes a risk-neutral decision maker that operational users are not (Wang; Matte).
  1. Skill and relative economic value of the ECMWF ensemble prediction system

    Richardson · 2000 323 cites

    Synthesis

    Computes the relative economic value of ensemble and deterministic forecasts across the full range of user cost-loss ratios, showing the ensemble is more valuable for most users even at lead times where its deterministic scores are unremarkable.

    Why it matters The original demonstration that forecast ranking depends on the user, and the framework that would settle most current arguments about learned ensembles if anyone applied it.

  2. The economic value of ensemble forecasts as a tool for risk assessment: From days to decades

    Palmer · 2002 173 cites

    Synthesis

    Extends the cost-loss framing across time scales from weather to climate and argues that the probabilistic forecast, not the best-guess trajectory, is the object with decision value.

    Why it matters Connects the decision argument to seasonal and climate horizons, where the deterministic forecast has no defensible meaning at all.

  3. The value of improved wind power forecasting: Grid flexibility quantification, ramp capability analysis, and impacts of electricity market operation timescales

    Wang · 2016 41 cites

    Synthesis

    Runs production-cost simulation over 480 scenarios on an IEEE 118-bus system configured to emulate the Midcontinent, California, and New England balancing authorities, and prices what a given improvement in wind power forecast accuracy is worth at each market timescale. The value tracks generation mix and available ramp capability rather than forecast error alone.

    Why it matters Converts a forecast improvement into money on a named system, and shows the same improvement is worth different amounts depending on what the system can do about it.

  4. Moving beyond the cost-loss ratio: economic assessment of streamflow forecasts for a risk-averse decision maker

    Matte · 2017 23 cites

    Synthesis

    Replaces the cost-loss ratio with a constant-absolute-risk-aversion utility when valuing ensemble streamflow forecasts, on the grounds that the standard framework assumes a risk-neutral decision maker and operational users are not.

    Why it matters The standard forecast-value metric encodes a preference assumption, so a system tuned to it can be optimizing for a decision maker who does not exist.

  5. Task-based End-to-end Model Learning in Stochastic Optimization

    Donti · 2017 arXiv

    Synthesis

    Differentiates through a stochastic optimization problem so the predictive model is trained on the resulting task cost, demonstrated on load forecasting for electricity scheduling and battery storage.

    Why it matters The first end-to-end demonstration that a model with worse predictive accuracy can produce strictly better decisions.

  6. Smart "Predict, then Optimize"

    Elmachtoub · 2017 arXiv

    Synthesis

    Defines the SPO loss, which measures decision regret rather than prediction error, and gives a convex surrogate with consistency guarantees for linear optimization problems.

    Why it matters Gives the decision-first framing a loss function with theory attached, which is why it became the reference formulation.

  7. Decision-Focused Learning: When and Why Traditional Prediction Models Fail

    Liu · 2026 arXiv

    Synthesis

    Characterizes the conditions under which a standard accuracy-trained predictor is already decision-optimal and the conditions where the decision-focused objective is necessary, rather than assuming the latter always helps.

    Why it matters Answers the question the survey left open, and prevents 'train on the decision' from becoming an unconditional recommendation.

World models, planning & physical control

Model-based control is where a learned dynamics model has to survive being acted on, and the reinforcement-learning line established the failure mode precisely: a planner will find the region where the model is wrong and optimistic. The methods that work are the ones that quantify and respect model uncertainty, not the ones that fit the dynamics best. Two demonstrations matter for large physical systems, and both are instructive for what they required rather than what they achieved: the stratospheric balloon controller and the tokamak plasma controller worked because the plant was densely instrumented, the actuation was characterized, and the loop closed fast enough to correct model error before it compounded. None of those conditions hold for a continental-scale atmosphere, which is the relevant fact whenever this literature is cited in support of intervention. The autonomous-driving stack is the other transfer argument in circulation, and its own literature is the check on it: imitation at thirty-million-example scale needed synthesized failures before the policy could drive, and policies cloned from logged behaviour are vulnerable to causal misidentification, where more observation makes them worse.

Key threads
  • Model exploitation is the central problem, and the fix is epistemic: ensembles that separate what is uncertain from what is stochastic, and rollouts short enough that the bound on policy improvement still holds (Chua).
  • Planning inside a learned model reached superhuman play without being given the rules, which is the strongest version of the argument. It also states its own preconditions: a resettable environment, cheap exact rollouts, and a scalar reward (Schrittwieser).
  • Differentiable simulation makes control of PDE-governed systems trainable end to end, with full state access and an exact solver, which is the ceiling rather than the operating point (Holl).
  • The two real deployments had tight sensing-actuation loops and bounded state spaces. Transfer to open natural systems is an assumption, and the balloon case is notable for exploiting the atmosphere rather than modifying it (Bellemare; Degrave).
  • The perceive-predict-plan stack that works in driving was not obtained by imitation alone. Thirty million expert examples needed synthesized failures and explicit penalties before the policy could drive, and cloned policies suffer causal misidentification, where adding observations degrades them (Bansal; de Haan).
  1. Causal Confusion in Imitation Learning

    de Haan · 2019 arXiv

    Synthesis

    Shows causal misidentification in benchmark control domains and in driving: a policy cloned from demonstrations degrades as its observation set grows, and only interventional data disambiguates which model is right.

    Why it matters Names the standing failure mode of any perceive-predict-plan stack learned from logged behaviour rather than from intervention.

  2. World Models

    Ha · 2018 arXiv

    Synthesis

    Learns a compressed latent representation of observations and a recurrent forward model of its dynamics, then trains a small controller entirely inside the learned model and transfers it back to the real environment.

    Why it matters The paper that made 'train the policy in the model's imagination' a concrete architecture, and the origin of the term this explorer is partly about.

  3. ChauffeurNet: Learning to Drive by Imitating the Best and Synthesizing the Worst

    Bansal · 2018 arXiv

    Synthesis

    Trains a driving policy by imitation and reports that thirty million real expert examples were not sufficient on their own, even with a perception system upstream and a controller downstream. Synthesized perturbations of expert trajectories, including collisions and road departures, together with explicit losses penalizing those events, were required before the model drove a real vehicle.

    Why it matters The perceive-predict-plan stack at industrial scale, carrying its own authors' evidence that observational expert data alone did not yield a usable policy.

  4. Mastering Atari, Go, Chess and Shogi by Planning with a Learned Model

    Schrittwieser · 2019 arXiv

    Synthesis

    MuZero learns a latent dynamics model trained only to predict reward, policy, and value, then runs tree search inside it, matching AlphaZero at Go, chess, and shogi without being told the rules and setting a new state of the art on Atari.

    Why it matters The strongest form of learn-the-dynamics-then-plan, and the clearest statement of its preconditions: a resettable environment, cheap exact rollouts, and a scalar reward, none of which a continental atmosphere supplies.

  5. Deep Reinforcement Learning in a Handful of Trials using Probabilistic Dynamics Models (PETS)

    Chua · 2018 arXiv

    Synthesis

    Separates aleatoric from epistemic uncertainty using an ensemble of probabilistic dynamics models and plans with trajectory sampling, matching model-free asymptotic performance with orders of magnitude fewer samples.

    Why it matters Establishes that what makes a learned model safe to plan against is its uncertainty representation, not its one-step accuracy.

  6. Mastering Diverse Domains through World Models (DreamerV3)

    Hafner · 2023 arXiv

    Synthesis

    Trains a single latent world model and actor-critic with fixed hyperparameters across more than 150 tasks spanning continuous control, Atari, and Minecraft, including collecting diamonds from scratch.

    Why it matters The current reference point for how far world-model planning generalizes when the environment is simulated, well-instrumented, and resettable.

  7. Learning to Control PDEs with Differentiable Physics

    Holl · 2020 arXiv

    Synthesis

    Uses a differentiable PDE solver to train a predictor-corrector control network end to end, achieving long-horizon control of fluid systems by planning through the solver rather than around it.

    Why it matters The clearest demonstration of learned control on PDE-governed dynamics, and a useful upper bound: it works with full state access and an exact simulator.

  8. Autonomous navigation of stratospheric balloons using reinforcement learning

    Bellemare · 2020 127 cites

    Synthesis

    Deploys a reinforcement-learning controller that steers stratospheric balloons by choosing altitudes to ride differing wind layers, outperforming the hand-engineered controller in flight over the Pacific under real wind uncertainty.

    Why it matters The rare case of learned control acting on a real atmospheric system, and it works by exploiting the atmosphere rather than changing it.

  9. Magnetic control of tokamak plasmas through deep reinforcement learning

    Degrave · 2022 313 cites

    Synthesis

    Trains a controller in a physics simulator of the TCV tokamak and deploys it on the real device, sustaining and shaping plasma configurations including ones that had not been demonstrated before.

    Why it matters The strongest sim-to-real result in physical control, and it depended on a trusted simulator, dense diagnostics, and millisecond actuation, none of which exist for the open atmosphere.

Intervention & the boundary of evidence

Weather modification has a seventy-year record, and the summary is that the physical chain has been observed end to end in one narrow case while operational effect sizes remain mostly unresolved. SNOWIE closed that chain for glaciogenic seeding of orographic clouds: seeding material tracked to ice initiation to snowfall reaching the ground, visible in radar, with the mass quantified. Convective seeding and hurricane modification have the opposite record, and both ended on the same two problems rather than on execution: a physical premise that did not hold, and a counterfactual that does not exist in a chaotic field. Against that background, AI weather models have produced a new genre of steering studies that find small perturbations redirecting a storm inside the model. This theme separates what has been demonstrated, what remains untestable, and what is being asserted, and treats the live question as one of experimental design: does better prediction, state estimation, and causal inference make any of these effects measurable for the first time?

Key threads
  • One case has a closed physical chain, observed rather than inferred, and the quantified effect is fractions of a millimeter of water equivalent over a limited area (French; Friedrich).
  • The failures are informative and specific: STORMFURY's premise was wrong about how much supercooled water hurricanes contain, and convective seeding assessments repeatedly could not separate effect from natural variability (Willoughby; Silverman).
  • The recurring obstacle is the counterfactual, not the intervention. The unseeded control does not exist, and the assessments that ended these programs failed on statistical power rather than on physics (NRC 2003).
  • Steering studies in AI weather models are optimization results about a surrogate. They establish that the model is sensitive to small perturbations in a targeted direction, which is a statement about the model (Huang; Liu).
  • Bounding a proposed intervention's plausible effect before testing it is a discipline that exists and is rarely applied outside the climate-intervention literature (Wood).
Open gaps
  • No AI-era steering study has been paired with a field-testable hypothesis, an actuation model that connects the perturbation to something physically deliverable, or a power analysis for detecting the claimed effect.
  • There is no validation protocol for a learned model's response to interventions absent from its training distribution, which is the exact regime every steering result occupies.
  1. Project STORMFURY: A Scientific Chronicle 1962-1983

    Willoughby · 1985 46 cites

    Synthesis

    Chronicles two decades of hurricane-seeding experiments and the observational findings that ended them: hurricanes contain too little supercooled water for the hypothesized mechanism, and the apparent successes were within the range of natural intensity fluctuation.

    Why it matters The reference case for how an intervention program fails, and it fails on the physical premise and the counterfactual, not on execution.

  2. A Critical Assessment of Glaciogenic Seeding of Convective Clouds for Rainfall Enhancement

    Silverman · 2001 36 cites

    Synthesis

    Reviews the randomized convective seeding experiments against explicit criteria for statistical and physical evidence and concludes the accumulated results do not establish an effect, with several apparently positive results failing on replication.

    Why it matters Sets the evidentiary bar the field itself adopted, and shows that decades of operational activity coexisted with no demonstrated effect.

  3. Critical Issues in Weather Modification Research

    National Research Council · 2003

    Synthesis

    Assesses the state of weather modification and concludes there is no convincing scientific proof of efficacy for intentional modification, while arguing the underlying questions remain scientifically legitimate and are limited by experimental design, statistical power, and the unobservable counterfactual.

    Why it matters The authoritative statement that separates 'not demonstrated' from 'disproven', and the framing that makes better prediction and state estimation relevant to the problem.

  4. Precipitation formation from orographic cloud seeding

    French · 2018 98 cites

    Synthesis

    Uses airborne and ground radar during the SNOWIE campaign to track the full chain from silver iodide release through ice initiation and growth to snowfall reaching the ground, with seeding lines visible in the radar field.

    Why it matters The one case where the physical chain of a weather intervention has been observed end to end rather than inferred statistically.

  5. Quantifying snowfall from orographic cloud seeding

    Friedrich · 2020 46 cites

    Synthesis

    Puts numbers on the SNOWIE cases, estimating the mass of snow attributable to seeding over specific areas and durations, on the order of fractions of a millimeter of water equivalent.

    Why it matters Converts the demonstrated mechanism into an effect size, which is what shows how far the demonstrated case is from an operationally meaningful one.

  6. Assessing the potential efficacy of marine cloud brightening for cooling Earth using a simple heuristic model

    Wood · 2021 33 cites

    Synthesis

    Estimates the achievable radiative forcing from marine cloud brightening using a transparent heuristic model, and shows the result is dominated by uncertainties in aerosol-cloud interaction rather than by deployment scale.

    Why it matters A model of how to bound an intervention's plausible effect before proposing to test it, on the intervention modality with the most active research program.

  7. Steering Tropical Cyclones Using Small Perturbations in an AI Weather Model

    Huang · 2026 arXiv

    Synthesis

    Uses gradients through a learned weather model to find small initial-condition perturbations that redirect a tropical cyclone track within the model.

    Why it matters The clearest instance of the genre this explorer exists to bound: an optimization result about a surrogate, presented in the vocabulary of intervention.

  8. Instability-Aware Steering of an Extreme Atmospheric River in an AI Weather Foundation Model

    Liu · 2026 arXiv

    Synthesis

    Targets perturbations at the dynamically unstable directions of the flow to steer an atmospheric river in a learned model, using instability structure rather than brute-force gradient search.

    Why it matters The most physically motivated of the steering studies, and still evidence about the model's sensitivity rather than about the atmosphere's.

A reading path

Start here and read in order; the path moves from foundations toward the open edge.

Companion podcast

8 episodes walking this map as an argument: what each study measured, where they disagree, and what that leaves open. Listen inline or read the transcript.

  • 1. The Sensor Boundary

    A learned atmospheric model has never seen the atmosphere, only instruments measuring transformations of it. What the observing system and the assimilation chain contribute to a forecast score, and why the credit is shared.

    Read transcript 15 min · 2,228 words

    A learned atmospheric model has never once seen the atmosphere. It has seen instruments measuring transformations of the atmosphere. Radiances in a satellite channel. Reflectivity at a radar wavelength. Station reports with inconsistent siting and inconsistent quality control. Radiosondes launched twice a day into a field that varies continuously in space and in time. Everything a model knows about the sky arrives through a sensor, and every sensor imposes its own geometry, its own spectral response, its own drift, and its own silence about the places it does not look.

    That sentence is the reason this series starts where it does. Welcome to an eight part walk through a literature map on physics AI and causal world models, following the eight themes of a curated explorer of about eighty papers. The question underneath all of it is not whether learned models can predict physical systems. That result is increasingly established, and we will spend the next episode on exactly how established it is. The question is what additional evidence is required before prediction can be interpreted as learned physical structure, as causal understanding, or as a model safe to use for counterfactual reasoning. Where do those claims separate? Which experiments discriminate between them? And which experiments still have not been run?

    This is episode one, the sensor boundary. Observations, data assimilation, and the layer where the physical world becomes numbers. It comes first because it is the layer that most discussions skip, and skipping it is how a claim about a model quietly becomes a claim about an entire observing system that the model did not build and does not own.

    Let us start with what the classical machinery does, because it is not a formality you can wave past on the way to the interesting part. The problem is an inverse problem. You have a state, the full three dimensional structure of the atmosphere at an instant, and you cannot measure it. You have observations, scattered, indirect, noisy, and far too few. You have a model that propagates a state forward and accumulates its own error while doing it. Data assimilation is the machinery that combines those into a best estimate.

    The reference point for the modern version is Geir Evensen in 1994, the ensemble Kalman filter. The problem he was solving is that the extended Kalman filter evolves error covariances through a linearization around a trajectory, and for a nonlinear geophysical model that linearization degrades in exactly the situations you care about. Evensen replaced it with a Monte Carlo ensemble. Run many states forward through the actual nonlinear model, and sample the forecast error statistics from the model itself rather than approximating them. That is the paper that made state estimation in geophysics an explicitly probabilistic exercise, and it is the baseline every learned assimilation method is still measured against. If you take one structural fact from it, take this one. Operational weather has been doing uncertainty quantification, at scale, in production, for thirty years. This is not a field that needs to be taught that a point estimate is insufficient.

    Now the part that matters most for machine learning, and it comes from Dick Dee in 2005, in a paper about bias and data assimilation. Suppose your observations are biased. Suppose your model is also biased. What happens? The analysis absorbs the bias into the state estimate. Variational bias correction, the standard operational tool, can only work when it is anchored by observations you have reason to believe are unbiased. And here is the sharp part. Model bias and observation bias are not separately identifiable from the residuals alone. The residual tells you they disagree. It does not tell you which one is wrong.

    Sit with what that implies for the current generation of learned weather models. Almost all of them are trained on reanalysis, and scored against reanalysis. Reanalysis is not ground truth. Reanalysis is a model observation fusion, produced by running a physics based assimilation system over the historical observing record. It contains a model. So a learned system trained on it inherits the assimilating model’s biases, and a learned system scored against it is being graded by a product that already made the choices it is being tested on. This does not make the skill fake. It makes the skill conditional, and the condition is a seventy year investment in an observing system and an assimilation pipeline that the learned model treats as a free input.

    There is a second measurement in this theme that I want to give real time to, because it answers an argument that circulates constantly right now. The argument goes like this. Language models are running into a data wall, but weather has no such wall, because the world collects an extraordinary volume of atmospheric observation every single week, more raw data than went into training a frontier language model. Therefore there is enormous headroom.

    The trouble with that argument is that it treats observation volume and assimilable information as the same quantity, and the field has measured the difference. Zhiquan Liu and Florence Rabier, in 2002, worked the problem out in a one dimensional analysis system, carrying representativeness error and observation error correlation through to the analysis error. For an optimal scheme, denser observations keep helping. For the suboptimal schemes that operational centers actually run, density past a threshold stops paying, and thinning correlated observations improves the analysis. Adding data made the estimate worse. Not through anyone’s incompetence, but because correlated error that the scheme does not model gets treated as independent information.

    Then Carla Cardinali in 2009 measured the marginal value of every observation in the ECMWF system directly, using adjoint forecast sensitivity to attribute twenty four hour forecast error reduction observation by observation, cross checked against observing system experiments over a summer month and a winter month. The assimilated observations reduce forecast error in aggregate. Specific subsets degrade it, traceable either to data quality or to suboptimality in the assimilation. So when someone argues from raw data abundance, there is already a per observation measurement waiting for them, and it has been running operationally for fifteen years. The right question is not how many bytes arrive per week. It is how much of that stream is informative after error correlation is accounted for, and nobody has published the curve showing where a learned model’s skill saturates as observing density rises.

    So that is the classical layer. Now, what happens when learning is pushed into it?

    The first move is to treat assimilation as a sampling problem instead of an optimization problem. François Rozet and Gilles Louppe, score based data assimilation, in 2023. Train a score model on short state trajectories, learn the distribution, then at inference time condition on whatever observations you happen to have and sample the posterior over the full trajectory. No linearization. No sequential filter structure. The posterior is something you generate rather than something you optimize toward. Most of the learned assimilation work since has followed that template, and you can see why it is attractive. It inherits the generative modeling toolchain, and it handles arbitrary observation configurations without rebuilding the machinery each time.

    The second move closes the loop from the other side. Xiao and colleagues, FengWu 4D-Var, also 2023. Four dimensional variational assimilation is expensive largely because it needs an adjoint, the derivative of the forecast model, and maintaining an adjoint of an operational physics model is one of the most labor intensive tasks in the field. A learned forecast model is differentiable by construction. The adjoint comes free. So they coupled a learned forecaster to a 4D-Var scheme and cycled the whole system. That is the point at which AI weather models stopped structurally depending on somebody else’s analysis in order to run.

    The third move removes the analysis entirely. Anna Vaughan and colleagues, Aardvark Weather, in 2024, replaced the full operational pipeline, assimilation and forecast model and post processing, with a single network trained from raw heterogeneous observations to station level forecasts, at a small fraction of the compute. And Mihai Alexe and colleagues, GraphDOP, trained a graph model purely on Earth system observations. No reanalysis anywhere in the pipeline. No physical state variables in the loss. It forecasts observed quantities directly in observation space, at skillful medium range, which is the cleanest existence proof I know of that useful forecast skill can be learned without ever seeing an analysis. Yumeng Zhao and colleagues in 2026 extend that line into the high resolution regime, where shedding analysis dependence has been hardest.

    Then there is a follow up on GraphDOP that belongs in this episode and also foreshadows episode five. Lean and colleagues, in 2025, probed what GraphDOP’s latent state actually contains, and found structure corresponding to physical fields the model was never given. It was trained on instrument readings. Inside, there is something that lines up with the physics. I want to be precise about what that does and does not establish, because this is exactly the kind of result that gets over read. It establishes that the observation stream carries enough information to reconstruct physical structure, and that a model under prediction pressure will build some of that structure internally. It does not establish that the structure is load bearing for the model’s skill, and it does not establish that the model can use it the way a physical model uses a conservation law. We will come back to that distinction with the tools to test it.

    The fourth move is the one closest to the instruments, and as someone who has spent time with spectroscopic data from spacecraft, it is the one I find most interesting. The observation operator is the map from model state to what the instrument would report. Historically it is hand derived physics. Radiative transfer for a satellite channel. Scattering assumptions for radar. It is also, historically, slow, and a known source of representativeness error.

    Stefanelli and colleagues, in 2025, learned the operator for weather radar reflectivity, replacing a forward operator whose microphysical assumptions are a documented weak point in convective scale assimilation. Howard and colleagues emulated the Community Radiative Transfer Model with calibrated uncertainty on the emulated radiances, and that word calibrated is doing the work. A surrogate observation operator that reports where it disagrees with the physical forward model is a scientific instrument. A surrogate that silently substitutes for it is a liability. The useful artifact is the error distribution, not the speedup.

    And then the failure mode that instrument people expect and machine learning evaluations rarely isolate. Li and colleagues, SPECTRA, in 2026, studied what breaks when a model pretrained on one satellite instrument’s spectral response and viewing geometry is transferred to another. Performance degrades, and it degrades in a way that looks exactly like model error unless your evaluation is designed to separate the two. Instrument bias, calibration drift, cross sensor generalization, provenance of the observation record. These are not administrative details. They are the mechanism by which a model can appear to misunderstand the physics when what it actually encountered was a different detector.

    The last paper here makes truth itself an experimental variable. Jin and colleagues built WeatherReal, a quality controlled global in situ station dataset, and scored AI and physical models against real measurements rather than against reanalysis. The apparent skill ordering depends on which truth you use. That is not a subtle effect at the margins. It is a reminder that when two systems are compared on reanalysis, part of what is being measured is proximity to the assimilating model.

    So what does this theme actually establish, standing back?

    It establishes that assimilation is learnable, in several distinct senses. As generative sampling of the posterior. As a differentiable coupling where the learned forecaster supplies the adjoint. And as an end to end replacement in which the analysis step disappears and the model learns from instrument readings directly. All three exist, in published, reproducible work.

    It establishes that the sensor boundary is where a large fraction of a learned model’s apparent error and apparent skill actually live. Change the truth dataset and the ranking moves. Change the instrument and the model degrades in a way that mimics physical error. Assimilate more data into a suboptimal scheme and the analysis can get worse.

    And it establishes the standard for what a claim in this area has to answer. If you assert that the physical world provides effectively unlimited training signal, you are making a claim about assimilable information, and there is a measurement for it, running in production, per observation, since 2009.

    What it does not settle is the attribution problem, and that gap runs through the whole series. No benchmark currently attributes a learned model’s error between the observing system, the assimilation, and the forecast model. So the model is wrong and the analysis is wrong right here remain indistinguishable in every score anyone reports. Everything in the next seven episodes inherits that ambiguity.

    Next episode, we go to the result that made all of this urgent. Between 2022 and 2024, neural networks trained on reanalysis beat the operational deterministic forecast on the headline scores, under the forecasting center’s own verification, at a small fraction of the inference cost. That happened. It is real, it is narrow, and the interesting work is in the words the headline leaves out. That is episode two.

  • 2. What the Forecast Models Actually Won

    The prediction argument is settled and the results are reproducible. What FourCastNet, Pangu, GraphCast, GenCast and NeuralGCM established, what operational replication added, and what all of them inherit from an analysis they cannot produce.

    Read transcript 14 min · 2,151 words

    In 2022, a neural network beat the best physics based weather forecast in the world on its own headline scores. In 2023, another one did it more convincingly, and a third did it in the probabilistic setting where operational value actually lives. By 2024 the European Centre for Medium Range Weather Forecasts was running its own learned model in its own suite. That sequence is real, it took about two years, and if you work in numerical weather prediction it was the most disorienting stretch in a generation.

    This is episode two of an eight part walk through a literature map on physics AI and causal world models. Last episode we stayed at the sensor boundary, where the physical world becomes numbers, and the takeaway was that a learned model inherits an observing system and an assimilation pipeline it did not build. Today we take the capability claim seriously and say exactly what was won, because the honest version of this result is narrower than the headline and more interesting than the headline.

    Start with cost, because cost is what changed the conversation first. Jaideep Pathak and colleagues, FourCastNet, in 2022. An adaptive Fourier neural operator trained on ERA5 reanalysis, producing global forecasts at a quarter degree resolution, running a full forecast in seconds on a single GPU. The skill was competitive rather than dominant. The economics were not competitive, they were a different category. A physics based ensemble is a supercomputer scheduling problem. A learned ensemble at that speed makes large ensembles nearly free at the margin, and once large ensembles are free you can ask questions that were previously unaffordable. That is what put learned forecasting on operational centers’ agendas.

    Then the skill argument closed, twice, within months. Kaifeng Bi and colleagues at Huawei, Pangu-Weather, a three dimensional Earth specific transformer with hierarchical temporal aggregation, beating the operational Integrated Forecast System on deterministic scores across most upper air and surface variables. And Remi Lam and colleagues at DeepMind, GraphCast, an encode process decode graph neural network on a multi mesh icosahedral grid, autoregressively rolled out during training, outperforming the operational system on the large majority of evaluated variable and lead time combinations, including several severe event proxies.

    Notice something about that pair. One is a transformer on a latitude longitude grid. The other is message passing on an icosahedral mesh. Architecturally they have almost nothing in common, and they arrived at the same place at the same time. That is the signature of a result that is about the problem rather than about the method. The atmosphere at synoptic scale, over the range where reanalysis is dense and the dynamics are large and smooth, turns out to be extremely learnable. After that pair the architecture question mostly stopped being where the gains came from, and the gains came from training strategy, from rollout schedules, and from scale.

    GraphCast also set the evaluation protocol as much as the score. It compared against the operational system on a defined variable and lead time grid, it reported severe event proxies rather than only continuous fields, and it published enough of the setup that later systems could be measured the same way. That matters more than it sounds. In a field where every group can choose its own truth dataset, its own initialization, and its own baseline, the paper that fixes a comparison protocol shapes what counts as progress for years. Much of the evaluation critique we will get to in episode four is an argument about how far that protocol goes and where it stops.

    We now have measured exponents for that scale claim, which I appreciate, because scaling arguments in this area are usually asserted rather than fit. Yu and colleagues, in 2026, fit empirical scaling relations between validation loss and model size, dataset size, and compute across several learned global weather models. Three results worth carrying. Aurora shows the strongest data scaling, where a tenfold larger training set cuts validation loss by up to a factor of three point two. Under a fixed compute budget, spending on data beats spending on parameters. And these models favor width over depth, which is the opposite of the language model result. So when someone says physical prediction will improve with scale, the measured version of that statement is that it improves with training data, and the atmosphere’s training data is bounded by the observing record, which brings us right back to episode one.

    Now the deterministic scores. There is a structural problem with them that the field walked into knowingly. If you train a model to minimize mean squared error over an ensemble of possible futures, the optimal output is the conditional mean. The conditional mean of a chaotic field is smooth. So a mean squared error trained forecast is systematically blurry, and blurriness is rewarded by the very metric being used to declare victory. A smoothed field has lower error against a noisy truth than a sharp field with the correct statistics but the wrong phase.

    Two responses. Ilan Price and colleagues, GenCast, a conditional diffusion model over fifteen day trajectories, generating sharp physically plausible ensemble members and beating the European ensemble on the large majority of continuous ranked probability score targets. That is the move that took learned forecasting out of the regime where blurring wins and into the probabilistic setting where operational users live. And a parallel line at the European Centre optimizing the probabilistic score directly rather than through a sampling procedure. Both are the field conceding, in the design of the loss function, that a mean squared error forecast is a conditional mean and not a forecast.

    Here is the part that gets flattened in popular accounts. The pure learned versus physics based framing did not survive contact with the results. Dmitrii Kochkov and colleagues, NeuralGCM, couples a differentiable dynamical core, actual fluid dynamics on the sphere, to learned parameterizations for the processes physics based models approximate anyway, and trains the whole thing end to end. It matches learned model skill at medium range. It stays numerically stable over decade long integrations. It reproduces emergent phenomena like tropical cyclone statistics that it was never explicitly trained to produce. That is not a compromise between two camps. The useful decomposition runs along a different axis, separating the parts of the system you have a reliable equation for from the parts you are approximating regardless.

    And the convergence kept going from the other direction. Polichtchouk and colleagues, in 2026, nudge the large scales of a physical model toward a learned forecast while leaving the small scales to the physics, so each component runs where it is strongest. The interesting question stopped being who wins and became how you couple them.

    Now the limits, and this is where I want to spend real time, because the limits are where the intellectual content is.

    Zied Ben Bouallègue and colleagues at the European Centre ran the independent replication in 2024. Take Pangu-Weather and GraphCast, put them inside the center’s own operational verification framework, initialize them from operational analyses rather than from the reanalysis they were trained on, and score them against the metrics and baselines the center itself uses. The headline gains are real and they survive. They also come with deficiencies in near surface fields, at small scales, and in extremes. The advantage is not uniform. It is variable specific and regime specific, and it narrows when the initialization pipeline changes.

    Massimo Bonavita, also at the European Centre, wrote the sharpest version of the critique in the same year. Learned models produce spatially over smoothed fields with deficient variance at small scales. Part of the apparent root mean square error advantage is an artifact of that smoothing rather than a reflection of better physics. And these systems depend on an assimilation system they neither model nor own. That last clause is the one to hold. In episode one we saw the pipeline being learned away, in Aardvark and GraphDOP and the high resolution follow ups, so the dependence is not permanent. But for the systems that produced the headline results, it was total.

    Ashesh Chattopadhyay and Pedram Hassanzadeh explain the mechanism underneath the smoothing rather than treating it as an implementation quirk. When you train an emulator of a multiscale chaotic system and roll it out autoregressively, the small scale variance is the first thing to go, and the failure compounds because the model then feeds its own smoothed field back to itself. They analyze why, and what training and spectral interventions delay it.

    Put those three together and you get a precise statement. These models can win at day five root mean square error while being wrong about how error grows. That combination is possible because aggregate error against a smooth truth and the propagation of uncertainty through scales are different quantities. One of them is what the leaderboard measures. The other is what a forecaster uses to decide whether to trust the forecast.

    There is one more limit that has started to dissolve, and it is worth flagging because the objection was so central. For two years the standard rebuttal to learned forecasting was that these models only work because numerical weather prediction feeds them, and that the dependence would bite hardest at high resolution, where analysis quality and physical downscaling do the most work. Zhao and colleagues, in 2026, report skillful high resolution forecasting with no physics based model anywhere in the pipeline. That does not retire the objection for the systems that produced the 2022 and 2023 headlines, which were fully dependent. It does mean the dependence is a property of particular systems rather than a structural necessity, and that the interesting version of the objection has to be restated as a claim about the observing record rather than about the analysis product.

    Two more entries round out the theme. Aurora, from Cristian Bodnar and colleagues, pretrains a single atmospheric model on a heterogeneous mixture of reanalyses, forecasts, and climate simulations, then fine tunes to tasks including air pollution and ocean waves, where dedicated operational systems are the baseline. That is the first credible demonstration that Earth system pretraining transfers across variables and domains rather than only across lead times, and it is the paper that earns the phrase foundation model in this space. We will pull on that thread properly next episode. And ECMWF’s own AIFS, developed inside the operational environment rather than against a published benchmark, marks the transition from external claims about beating the operational system to an operational center running a learned model in its own suite.

    It is also worth being concrete about what operational like means in that replication, because the phrase is doing a lot of work. Training on reanalysis and evaluating on reanalysis is a closed loop. Reanalysis is produced after the fact, with the full observing record available in both directions in time, which is a luxury no real time system has. Initializing from an operational analysis means starting from a state estimate produced under a deadline, from observations that arrived in time, with whatever gaps that implies. Models trained in the retrospective setting and run in the real time one are being asked to generalize across a distribution shift that has nothing to do with the atmosphere and everything to do with how the data was assembled. Some of the narrowing that Ben Bouallegue and colleagues measured is exactly that.

    So what did the forecast models actually win?

    They won speed by orders of magnitude, and that is not a footnote, it changes which experiments are affordable. They won the deterministic headline scores at medium range, verified independently, inside the incumbent’s own framework. They won the probabilistic setting once the loss functions were fixed. And they won a real architectural insight, which is that at synoptic scale the atmosphere is learnable to a degree very few people predicted in 2020.

    What they did not win is a claim about physical understanding, and nothing in this theme was designed to test one. The skill is concentrated in the smooth, large scale, well observed part of the problem. It is conditional on an initialization pipeline built by the systems being outperformed. It comes with under representation of small scales and of error growth. And the papers that document those limits were mostly written by the people whose model was beaten, which is worth noticing, because it means the critical literature here is unusually specific and unusually well instrumented.

    Next episode, we take the generalization claim seriously. Operator learning promised something bigger than one good weather model. Learn the solution map of a family of equations rather than one solution, and transfer across initial conditions, across resolutions, and eventually across equations entirely. That promise has real results behind it and a documented baseline problem sitting next to it. Foundation models for physical systems, and what actually transfers. That is episode three.

  • 3. What Actually Transfers

    Resolution invariance, multi-physics pretraining in the low-data regime, and the two meanings people give the phrase foundation model. Including the review that found most reported speedups over numerical solvers used weak baselines.

    Read transcript 13 min · 2,141 words

    There is a promise underneath most of the excitement about physics AI, and it is bigger than any single weather model. The promise is that you can learn the solution map of a whole family of equations rather than one solution, so a trained model transfers to new initial conditions, to new resolutions, and eventually to equations it has never seen. If that works, you do not have a fast weather model. You have something closer to a general purpose simulator for physical systems.

    This is episode three of eight, and today we test that promise against what has actually been shown. Last episode ended with Aurora, a single atmospheric model pretrained on a mixture of reanalyses and forecasts and climate simulations, then fine tuned to air pollution and ocean wave tasks where dedicated operational systems are the baseline. That is the reference point for what an Earth system foundation model has to demonstrate. Today we go back to where the idea started and follow it forward, and we keep two things apart that get conflated constantly: what has been shown to transfer, and what has only been shown to fit.

    The founding move is operator learning. A neural network normally maps vectors to vectors. An operator maps functions to functions. If you can learn the operator that takes an initial condition function to the solution function at a later time, you have learned the equation’s behavior rather than one trajectory through it.

    Two canonical formulations arrived around the same time. Lu Lu and colleagues, DeepONet, in 2019, builds on a universal approximation theorem for operators, with a branch network that reads the input function and a trunk network that reads the query location, and the two combine to evaluate the output function anywhere. What I like about DeepONet is that the function space claim is explicit and theoretical rather than architectural. It says what class of thing is being approximated and proves it can be.

    Zongyi Li and colleagues, the Fourier neural operator, in 2020, gets there differently. Parameterize the integral kernel directly in Fourier space. Because the parameters live in spectral space rather than on a grid, the learned operator is resolution invariant. You can train at low resolution and evaluate at high resolution, which is not something a convolutional network on a fixed grid does. They reported large speedups over spectral solvers on Navier Stokes. And two years later that architecture, adapted, became FourCastNet, the first global learned weather model from episode two. The lineage is direct.

    The property to hold onto from both is resolution invariance, because it is the reusable one. Accuracy on one grid is a benchmark number. An operator that is defined in function space is a different kind of object, and the difference shows up the moment you change discretization.

    It is worth being concrete about why that property is worth wanting. Numerical solvers have a hard coupling between accuracy and grid. Refine the mesh and cost rises, often steeply, and a solution computed on one mesh does not transfer to another without interpolation that introduces its own error. An operator learned in function space breaks that coupling in principle. You can train where data is cheap and evaluate where you need answers. In practice the invariance is approximate and degrades outside the spectral range the training data covered, which is a limitation worth stating plainly, because resolution invariance is often quoted as if it were unconditional. What survives is still substantial. The object being learned is defined independently of the discretization, so questions about it are questions about the physics rather than about the grid.

    The third strand is the particle and mesh line. Alvaro Sanchez-Gonzalez and colleagues, Learning to Simulate Complex Physics with Graph Networks, in 2020. Represent the system as particles, connect neighbors into a graph, learn the dynamics as message passing over that graph. It generalizes across materials, across geometries, and to orders of magnitude more particles than it saw in training. That inductive bias, local interactions on a graph, is the one GraphCast later carried onto the sphere. If you want the intellectual genealogy of the model that beat the operational forecast, it runs through fluid simulation on particle graphs, not through weather.

    That generalization result deserves a second look, because particle counts are a good stress test. A model trained on a few thousand particles and evaluated on hundreds of thousands is not interpolating within its training distribution in any ordinary sense. What makes it work is that the learned function is local. Each particle’s update depends on its neighbors, and the neighborhood statistics are similar whether the system has three thousand particles or three hundred thousand. That is a specific and checkable reason for generalization, tied to the structure of the physics rather than to model capacity, and it is the kind of explanation I want attached to every transfer claim in this field. When a model generalizes, the interesting question is what invariance it exploited.

    Now the foundation model claim proper. Tung Nguyen and colleagues, ClimaX, in 2023, pretrains on heterogeneous CMIP6 climate simulations using variable tokenization and aggregation, then fine tunes to forecasting, downscaling, and climate projection, with different variables and different resolutions at each stage. It is the first paper to argue that Earth system tasks share a backbone, and it fixed the shape of the transfer evaluation that later work uses.

    Michael McCabe and colleagues, Multiple Physics Pretraining, ran the most direct test I know of. Pretrain one transformer on several distinct partial differential equation systems at once through shared normalized field embeddings. Then transfer to physics that was not in the pretraining set. Does it beat training from scratch? Yes, in the low data regime. That is the honest form of the result, and I want to dwell on the qualifier rather than skip it, because the qualifier is the finding. Pretraining across physics buys sample efficiency. It gets you to a given accuracy with less data on the new system. It does not necessarily raise the accuracy ceiling once you have plenty of data on the new system. Those are different claims with different consequences, and only one of them is what people usually hear.

    Maximilian Herde and colleagues, Poseidon, scales that idea with an operator transformer and a training scheme that exploits the semigroup structure of time dependent equations, which is a lovely detail. If the solution operator advances time, then applying it twice should be the same as advancing twice as far, and that structure lets you generate far more training pairs from the same trajectories. They report large sample efficiency gains on unseen equations. It is the strongest quantitative transfer result in operator learning, and it suggests the useful pretraining signal is time evolution structure rather than any particular physics.

    So there is real transfer here. Now the correction, and it is a serious one.

    Nick McGreivy and Ammar Hakim, in 2024, reviewed the literature claiming machine learning speedups over numerical solvers for fluid related partial differential equations. Not one paper, the literature. What they found is that the large majority of those comparisons run against baselines weaker than the domain standard. A coarse solver, an old scheme, an implementation not tuned the way a practitioner would tune it. And separately, that reporting biases suppress negative results, so the published record is not a sample of what was tried. Their conclusion is that a large fraction of the reported speedup does not survive a fair comparison.

    I have opinions about why this happens, and none of them require anyone to be acting in bad faith. Implementing a strong numerical baseline is specialist work. It is slow, it is unglamorous, it makes your result look worse, and reviewers rarely demand it because reviewers in machine learning venues frequently cannot evaluate it. Meanwhile the incentive gradient points straight at the headline. The fix is structural rather than moral: the baseline has to be built by someone who wants it to win.

    The other thing McGreivy and Hakim’s review does is make the physics AI literature legible as a literature. When you read papers one at a time, each speedup is plausible. When you read the distribution, you can see the selection effect. That is a general lesson for reading this whole field, and it is the reason the evaluation episode is the hinge of this series.

    The last entry in this theme reframes what the bottleneck even is. Ruben Ohana and colleagues released The Well, fifteen terabytes of numerical simulations across sixteen physical systems, from biological flows to supernova explosions, with a common interface and baseline results. That is not a model. It is the substrate that makes cross physics transfer claims testable at all, and it is the most consequential artifact in this line for exactly that reason. The bottleneck moved from architecture to data, and once it does, datasets become the thing worth building.

    Which brings me to the vocabulary problem, and this one matters because it is where a lot of hand waving hides. The phrase foundation model names two different things in this area. One is a pretrained backbone with demonstrated downstream transfer, where you can point at a task the model was not trained for and show that reusing the pretrained representation beats training from scratch. The other is a large task specific simulator with a label attached. Only the first is a testable claim. The second is a description of size.

    There is a practical test you can apply while reading. Ask what the model was fine tuned to do, what the baseline on that downstream task was, and whether the downstream task differs from pretraining in kind or only in degree. Different lead times are degree. Different variables on the same grid are degree with a bit of kind. A different physical domain, evaluated against the production system that domain already has, is kind. The further a paper’s evidence sits from that last case, the more the word foundation is describing training set size rather than a demonstrated capability.

    Aurora is the case that argues the first version credibly, because the fine tuned tasks include air pollution and ocean waves, and the baselines are dedicated operational systems rather than other neural networks. Beating a specialized production system on a task outside your pretraining objective is the kind of evidence the phrase is supposed to imply. Most systems described as physics foundation models have not shown that, and many have not attempted it.

    The low data regime qualifier is also where I think the practical value of this line actually sits, and I do not mean that as a downgrade. Most physical sciences are permanently in the low data regime. A planetary spectrometer returns one observing campaign, not a benchmark corpus. A laboratory plasma experiment produces a few hundred shots. A field campaign happens once, in the weather it happened to get. If pretraining across many simulated systems buys sample efficiency on a new system with two hundred examples, that is worth more to those fields than a small accuracy gain on a system with two million examples. The framing of foundation models tends to emphasize the asymptotic claim. The measured result is about the small data end, and the small data end is where most science lives.

    So what does this theme establish?

    Resolution invariance is a real, reusable property with two clean formulations behind it, and it is the closest thing operator learning has to a general result. Multi physics pretraining transfers, measurably, in the low data regime, and the mechanism appears to be time evolution structure rather than shared physics per se. Earth system pretraining transfers across variables and domains, demonstrated once, convincingly, against operational baselines. And a large fraction of the surrounding speedup literature does not survive a fair comparison, which is a fact about the field’s evidence base rather than about any individual result.

    What it does not establish is that a model trained on many physical systems has learned physics in any sense stronger than a reusable representation. Nothing in this theme tests whether the model can apply a conservation law it was not shown, or extrapolate outside the regime it was pretrained on, or reason about a system whose governing equations differ structurally rather than parametrically. Those tests exist, and they are the subject of episode five, and they mostly come back negative.

    Before we get there, we need to talk about measurement, because every claim so far has been reported as a number, and the numbers are where the field is weakest. Next episode, evaluation, uncertainty, and physical validity. Where the headline score misleads, why aggregate skill and extreme event skill come apart, and what forecast verification worked out decades before machine learning arrived and then had to work out again. That is episode four.

  • 4. When the Score Lies

    Murphy's three kinds of goodness, proper scoring rules, the double penalty, and the experiments showing that aggregate skill and extreme-event skill come apart on the storm that mattered.

    Read transcript 13 min · 2,032 words

    In the autumn of 2023 a windstorm crossed the English Channel and did serious damage across northwest Europe. Afterwards, Andrew Charlton-Perez and colleagues ran the obvious comparison. They took the learned weather models that had by then been reported to beat the operational physics based system on aggregate scores, and they asked how those models did on that one storm. The answer was split. The learned models captured the synoptic evolution well. They knew where the low was going and roughly how deep it would get. What they under-represented was the near surface wind extremes, which is to say the specific quantity that broke things.

    That is the shape of this entire theme. An aggregate score says one thing. The property a user actually needs says another. And the gap between them is structural, with a literature going back three decades that mostly predates machine learning and that machine learning largely had to rediscover.

    This is episode four of an eight part walk through a literature map on physics AI and causal world models. Evaluation, uncertainty, and physical validity. I put this theme fourth rather than last on purpose, because it is the theme that decides whether the other seven are true. Every claim in the previous three episodes arrived as a number. This episode is about what those numbers can and cannot certify.

    Start with Allan Murphy, writing in 1993, asking what a good forecast even is. His answer separates three things that get conflated constantly. Consistency is whether the forecast matches what the forecaster actually believes. Quality is correspondence between the forecast and what was observed. Value is the benefit the forecast produces for whoever used it. Murphy’s point, and it is a sharp one, is that these three can move in opposite directions. A forecaster who knows the truth but hedges toward the safe answer sacrifices consistency and often improves quality as measured by a squared error score. A forecast can improve in quality and lose value if the improvement lands in a regime nobody makes decisions in. Thirty years later the learned forecasting literature reports headline scores, which is one of Murphy’s three axes, measured partially. Almost every argument I have watched about whether an AI weather model is better than the operational system turns out on inspection to be an argument about which of Murphy’s three goodnesses the two sides are talking about.

    Then there is the question of what a score is even allowed to do. Tilmann Gneiting and Adrian Raftery worked out the formal theory in 2007. A scoring rule is strictly proper when it is optimized only by reporting your true predictive distribution, which means you cannot improve your score by shading your forecast toward something you do not believe. The continuous ranked probability score, the logarithmic score, and the Brier score all fall out of that one framework as instances. This matters because the alternative is a score that rewards hedging, and a system optimized against a hedgeable score will hedge. Root mean square error is the canonical case. It is minimized by the conditional mean, so a deterministic model trained to minimize it learns to produce something like the average of the possible futures, which is smoother than any real atmospheric state and which is exactly the behavior that shows up later as missing extremes.

    The same year, Gneiting, Fadoua Balabdaoui, and Raftery published the companion idea. Maximize sharpness subject to calibration. Calibration is whether the stated probabilities are honest, so that events you call thirty percent happen about thirty percent of the time. Sharpness is how concentrated the forecast is, how much it actually commits. A climatological forecast is perfectly calibrated and useless. A confident wrong forecast is sharp and dangerous. You want the most concentrated distribution that remains honest, and the two properties have to be measured separately, with tools like the probability integral transform histogram, because a single number cannot tell you which one you lost. That distinction is not a technicality. It is the axis along which the current arguments about diffusion based ensemble forecasting are actually being fought, and it is invisible to anyone reading a leaderboard.

    Now the spatial problem. Beth Ebert’s 2008 review of neighborhood verification is the canonical treatment of what people call the double penalty. Imagine a forecast that predicts a squall line perfectly in structure and intensity but displaced by thirty kilometers. Point by point matching penalizes it twice, once for predicting rain where none fell and once for missing rain where it did. A blurred forecast that puts weak rain everywhere avoids both penalties and scores better while being less useful to anyone deciding whether to close a road. Ebert’s contribution is to organize the neighborhood based methods by the decision model each one implies, which reframes verification as a question about the user rather than about the field. It is also the reason you should be suspicious when a learned model’s advantage over a physical model is largest in exactly the variables where smoothing pays.

    Which brings us to the papers that measure the gap directly. Noah Brenowitz and colleagues took learned models that the field had already declared roughly equivalent, because their deterministic scores matched, and asked what happens when you use them the way forecasts are actually used, as the seed of an ensemble. Models with equal root mean square error turned out to have substantially different spread and different probabilistic skill. They propose a lagged ensemble protocol that exposes this cheaply, without retraining anything. The result to hold onto is that deterministic scores do not rank models the way probabilistic use ranks them, demonstrated on the specific models whose equivalence the field had accepted.

    Then the tails. Sun and colleagues ran what I think is the cleanest experiment in this theme. Take the most intense tropical cyclones out of the training set entirely, then test whether learned models can forecast storms of that intensity. They cannot. The models fail to extrapolate to unseen intensity even though the governing physics has not changed at all between the storms they saw and the storms they did not. That converts a suspicion into a mechanism. Extreme event failure is not vaguely about smoothing or about loss functions. It is out of distribution generalization failure, isolated by experiment. And it has a direct implication for the counterfactual question this whole series is about, because a model that cannot extrapolate to an unseen intensity of a phenomenon it has seen thousands of times is not a model you can interrogate about a world that differs from the one it was trained on.

    Gabler and colleagues in 2026 do the systematic version, testing learned models on extremes across variables and regions and quantifying how much of the reported skill advantage survives when you restrict the evaluation to the tails. Much of it does not. The value of that paper is that it converts a widely repeated intuition into a variable by variable number, which is the form a claim has to be in before you can act on it.

    Now the subtlest one, and my favorite in this theme. Landsberg and colleagues ask what happens when the climate a model was trained on drifts away from the climate it is forecasting. Learned models are trained on historical reanalysis. The atmosphere is not stationary. So the evaluation period differs from the training period in distribution, and part of any measured skill change is attributable to that drift rather than to the model. The consequence is uncomfortable. It means a benchmark number is not a fixed property of a model, it is a property of a model and an evaluation window jointly, and the confound grows silently over time. No standard protocol currently controls for it. When you see two systems compared across different evaluation years, some unknown fraction of the difference belongs to the calendar.

    Against that background, WeatherBench 2, from Stephan Rasp and colleagues, is the piece of infrastructure that made the field legible at all. Fixed ground truth, fixed baselines, headline deterministic and probabilistic scores, every major learned model reported against the operational Integrated Forecast System on the same footing. Before it, cross paper comparison was mostly impossible, because everyone chose their own baseline and their own truth dataset. After it, comparison is routine. It is also, and this is the part worth sitting with, an artifact whose choices now shape what the field optimizes. Benchmarks reward what they measure. A scoreboard that leads with deterministic scores at fixed lead times produces models that are excellent at deterministic scores at fixed lead times. That is not a criticism of the benchmark, which is careful and explicit about its own limits. It is an observation about what happens when a community gets a shared objective.

    Paul Ullrich and colleagues wrote the 2025 statement of what a real evaluation would have to include before a learned Earth system model can be trusted. Physical consistency checks, out of sample and out of distribution testing, and independence between the people who built the model and the people who evaluate it. That last requirement is the one I would underline. The field’s evidence base is largely self reported, against baselines and initializations and truth datasets chosen by the developers. Ullrich is the community saying out loud that this is not sufficient, which is a normal stage for a field to pass through and a healthy one.

    Nick McGreivy and Ammar Hakim make the same argument with numbers, in the fluid dynamics setting rather than the weather setting. They reviewed the literature claiming machine learning speedups over numerical solvers for partial differential equations and found that the large majority compare against baselines weaker than the state of the art in the target field, and that reporting biases suppress the negative results that would balance the record. Their finding generalizes past fluids to any claim of the form a learned model beats the domain standard method, which is the claim structure this entire explorer is built around.

    So what does the theme leave open? Two gaps, and both of them are the kind that would take real work rather than a clever idea.

    The first is attribution. There is no benchmark that decomposes a learned model’s error into the part contributed by the observing system, the part contributed by the assimilation, and the part contributed by the forecast model. Which means the model is wrong and the analysis is wrong in this region remain indistinguishable in every score anyone reports. That is the loop closing back to episode one. The sensor boundary is not just a philosophical point about what a model has seen. It is a live measurement problem that makes every reported error number partially uninterpretable.

    The second is statistical. Extreme event evaluation has no agreed power analysis. Nobody states how many events are required before an improvement in tail behavior counts as a result rather than as a sample. We have a field that has correctly identified extremes as the place where learned models are weakest, evaluating on event counts that are small by construction, without a convention for how small is too small. That one is fixable, and someone should fix it.

    Here is where that leaves the argument. We have three episodes of evidence that learned models predict physical systems well, and one episode establishing that the measurements behind that evidence are weaker than they look, particularly in the tails, particularly under distribution shift, and particularly when the evaluator and the developer are the same person. None of this says the models are bad. Several of them are genuinely excellent. It says the evidence supports a narrower claim than the one usually made.

    And it sets up the question the rest of the series is about. Everything so far has been about correspondence between a prediction and an observation. Next episode we ask whether any of it implies that the model has learned the physics. Probes of learned representations, models that predict orbits perfectly and cannot apply Newtonian mechanics, symbolic recovery of governing equations, and the difference between a model that fits a system and a model that understands it. That is episode five.

  • 5. The Kepler to Newton Problem

    Five different claims get collapsed into the word causal. Probing finds physical structure inside learned forecasters, dynamical tests find the wrong error growth, and the direct test of whether prediction delivers the law comes back negative.

    Read transcript 14 min · 2,175 words

    There is a story people tell about machine learning and physics, and it goes like this. Johannes Kepler stared at Tycho Brahe’s planetary observations and extracted three empirical laws describing how planets move. Kepler could predict. He could not explain. Then Isaac Newton came along and produced the law underneath, gravitation, from which Kepler’s rules fall out as consequences, and which also explains falling apples and tides and the orbits Kepler never measured. The hopeful version of the story says that today’s predictive models are at the Kepler stage, and that scale or the right probe will take them to Newton.

    Keyon Vafa and colleagues ran that story as an experiment. They trained models on orbital trajectories. The models predict orbits accurately. Then they adapted those models to new physics tasks that require actually applying Newtonian mechanics, and measured what inductive bias transferred. The answer was that the models behaved as though they had learned task specific heuristics rather than the law. Excellent at the trajectories they were trained on. Unable to carry the mechanics forward.

    This is episode five of an eight part walk through a literature map on physics AI and causal world models. Causal structure and mechanistic understanding. The theme where the central question of the series gets asked directly. Does prediction pressure force a model to learn mechanism?

    The literature supports both answers, and the reason it supports both is that the word causal is being asked to carry five different claims with five different burdens of proof. It is worth separating them before going further, because almost every confused argument in this space is two people holding different items from this list.

    The first claim is that there is a predictive relationship. Variable A helps forecast variable B. That is a statement about correlation structure in a distribution and it is the easiest to establish.

    The second is that the model contains an internal representation encoding a physical variable. Somewhere in the activations there is something that tracks temperature, or vorticity, or a quantity nobody put in the training data. That is a representation claim, established by probing, and it is not the same as the first.

    The third is that a causal structure has been identified. This variable causes that one, in a form you could write down as a graph and test.

    The fourth is that the model supports a valid counterfactual. If this had been different, that would have followed.

    And the fifth is that the model supports intervention. You can reach in, set something, and read out what actually happens rather than what typically co-occurs.

    Those five are ordered roughly by strength, and evidence for one does not transfer to the next. Most claims that a model has learned physics are claims of type two being reported as though they established types four and five.

    That ladder is not academic bookkeeping. Consider what each rung licenses in practice. Rung one lets you forecast. Rung two lets you say something about how the model works internally, which is useful for debugging and for trust but not for planning. Rung three lets you attribute, which is what an attribution study of a heatwave needs. Rung four lets you answer what would have happened without the aerosol forcing, or without the reservoir, which is what most climate policy questions actually are. Rung five lets you plan, because planning means evaluating actions you have not taken. When a company or a paper says its model supports counterfactual reasoning about physical systems, the useful response is to ask which rung the evidence sits on, and the answer is usually two.

    Start with the inference machinery, because the causal literature for dynamical systems already exists and is explicit about its own assumptions. George Sugihara and colleagues introduced convergent cross mapping in 2012, which detects causal coupling in deterministic dynamical systems by testing whether the reconstructed attractor of one variable can recover the other. The important part of that paper for our purposes is the demonstration that Granger causality, the workhorse of time series causal inference, fails in exactly this regime. Coupled nonlinear dynamics needs different machinery than stochastic time series, and Earth system data is coupled nonlinear dynamics.

    Jakob Runge and colleagues built the method the Earth system community actually runs, called P C M C I, a two stage conditional independence approach that controls false positives in high dimensional autocorrelated data where naive conditioning inflates them. It works, and it comes with assumptions: causal sufficiency, meaning no unobserved common driver, stationarity, and adequate sampling. Earth system data routinely violates all three. That is not a criticism of the method, which states its conditions honestly. It is a warning about what a discovered graph from that machinery licenses you to say.

    Jonas Peters, Peter Bühlmann, and Nicolai Meinshausen supply what I think is the most useful operational definition available. Causal inference using invariant prediction. The idea is that the conditional distribution of a target given its true causal parents stays invariant across environments, while spurious relationships shift. So you turn causal discovery into a hypothesis test across environments, and you get confidence statements about the recovered parent set. Invariance is the property that survives distribution shift, which makes it the only definition that helps when a model is deployed outside the climate it was trained on. Bernhard Schölkopf and colleagues then make the programmatic argument in Towards Causal Representation Learning: the variables over which causal structure is defined have to themselves be learned from low level observations before any of this applies. That paper is the clearest statement of why the model learned physics is a representation claim rather than a prediction claim, and what it would take to establish it.

    Now the evidence, which cuts both ways, and I want to give both sides properly because the field tends to cite whichever half suits the argument.

    On the side of mechanism, the strongest current result comes from Lean and colleagues probing GraphDOP. GraphDOP, recall from episode one, is trained only on instrument measurements, with no reanalysis and no physical model in the loop. Probing its internal activations recovers structure aligned with physical variables it was never given and never directly observed. A model that has only ever seen radiances carries something inside it that tracks the fields underneath. That is a genuine result and it is the best evidence available that prediction pressure induces physically meaningful internal structure.

    Baño Medina and colleagues take a different route, running gradient based sensitivity analysis on a learned model during Cyclone Xynthia, a well studied explosive cyclogenesis case, and comparing the influence patterns the model reveals against the physical precursors the meteorological literature already identified. Asking the model what it is using, rather than checking what it output. And Tempest and colleagues in 2026 build the tooling to push this further, inspecting internal circuits and attributing predictions to identifiable internal computations rather than to input regions. That is the move from saliency maps toward circuit level evidence, which is the kind of evidence a mechanism claim would actually require.

    On the other side, the negative results are equally concrete and in some ways sharper, because they test properties that physical models satisfy by construction.

    Tobias Selz and George Craig asked whether a learned weather model exhibits the butterfly effect. Apply a tiny localized perturbation to the initial state and watch how error grows. In the real atmosphere, and in physical models that resolve the relevant scales, error injected at convective scales propagates upscale, and that upscale propagation is what sets the intrinsic predictability limit. The learned model does not do it. Forecast skill is high. The error growth dynamics are wrong.

    Gregory Hakim and Sanjit Masanam do the systematic version, subjecting a learned model to classical dynamical meteorology experiments where the answer is already known: steady tropical heating, extratropical cyclone development, geostrophic adjustment. The model reproduces some canonical responses and gets others qualitatively wrong. What I value most about that paper is methodological. It establishes the right instrument for mechanism claims, which is not held out forecast skill but controlled experiments the physics already has answers to. If you want to know whether a model knows something, do not test it on more of what it was trained on. Perturb it in a way the training distribution never contained and check the response against theory.

    Then the direct tests, run outside weather, where the experiment can be clean.

    We opened with Vafa. The orbital mechanics case is the sharpest because the ground truth is unambiguous and the model’s prediction quality is not in dispute. Bingyi Kang and colleagues run the complementary experiment on video generation. Train diffusion video models on two dimensional mechanics with an unlimited supply of simulated data, so data quantity is explicitly not the bottleneck, and separate three kinds of generalization: within the training distribution, combinatorial recombination of seen elements, and genuinely out of distribution. In distribution generalization is perfect. Combinatorial generalization improves measurably with scale. Out of distribution generalization fails. And the failure mode is diagnostic rather than mysterious. The models retrieve the nearest training example, prioritizing color first, then size, then velocity, then shape, which is a ranking of surface attributes and not a physical hierarchy at all. Unlimited clean data on simple mechanics did not produce the rule.

    Pim de Haan, Dinesh Jayaraman, and Sergey Levine add the result from imitation learning, and it is the one I find most unsettling. Behavioral cloning fits a discriminative model with no access to the causal structure of how the expert and the environment interacted. The consequence is causal misidentification: giving the policy more information can make it worse, because confounded cues that correlate with expert actions are often easier to fit than the actual causes. More observations, worse policy. Resolving it requires targeted intervention, either interacting with the environment or querying the expert. Passive data will not do it, no matter how much of it there is.

    Kori and colleagues, in 2026, write down the formal version of the gap. A unifying perspective on causal world models, setting out what a world model has to satisfy before its rollouts can be read as interventions rather than as conditional predictions. That distinction is the whole ballgame. A rollout answers the question, given what I have seen, what typically follows. An intervention answers, if I set this variable to this value, what follows. Those coincide only under conditions that a learned model has no reason to satisfy and that nobody currently checks.

    It is also worth noticing what the Kepler to Newton story quietly assumes. Kepler did not become Newton. A different person, working eighty years later, with different tools and a different question, produced the law. The empirical regularities were necessary input and they were not sufficient, and the step from one to the other was not more data or more staring. It was a change in what was being asked. The optimistic reading of the analogy treats the transition as something a predictive system undergoes on its own given enough scale. Nothing in the history supports that, and nothing in the experiments does either. The more defensible reading is that learned models are useful Kepler stage instruments, producing regularities at a rate no human could, and that the Newton step still belongs to whoever designs the experiment that discriminates between explanations.

    So where does the theme land? It leaves one gap open that I think is the most tractable important experiment in this entire explorer. We now have probing results showing that learned forecasters contain internal representations of physical variables. We do not know which of those representations are load bearing for the model’s skill and which are correlates that probing can find but ablation would not disturb. Probing tells you information is present and decodable. It does not tell you the model uses it. The experiment is ablation and causal mediation on internal representations of a learned weather model, and it is a matter of engineering effort rather than of new theory.

    Put the two halves together and the honest summary is narrower than either camp’s version. Learned models develop internal structure aligned with physics nobody handed them, which is real and interesting. They fail dynamical tests that physical models pass by construction. And where the Kepler to Newton premise has been tested directly, on orbits and on simulated collisions, the model that predicted well did not carry the law. Prediction pressure produces representation. There is no current evidence that it produces mechanism in the sense required for counterfactuals.

    Which raises the practical question. If a model cannot be trusted for counterfactual reasoning, what can it be trusted for? Quite a lot, as it turns out, and the answer runs through the people who actually use forecasts. Next episode, forecast to decision. Value versus skill, the users who need probabilities rather than point predictions, and the surprising finding that the link between a better forecast and a better outcome is weaker and stranger than anyone assumes. That is episode six.

  • 6. Better for Whom

    Forecast value is a property of a forecast and a user together, and the ranking of two systems flips along the range of users. The relative-economic-value curve has existed since 2000 and nobody has run it on a learned ensemble.

    Read transcript 13 min · 2,105 words

    Suppose you grow fruit, and a frost overnight will destroy the crop unless you run the wind machines and the heaters. Running them costs you something. Call that cost C. Losing the crop costs you rather more. Call that loss L. The only number that matters to your decision is the ratio between them. If protection costs one tenth of what the crop is worth, you should protect whenever the forecast probability of frost exceeds one tenth, because over a season that rule loses you less than any other rule. If protection costs nearly as much as the crop, your threshold is high and you protect only when the forecast is nearly certain.

    Now here is the part that matters for this series. Two forecast systems can be compared, and which one is better for you depends on where your ratio sits. Not on which one has the lower error. On your ratio.

    David Richardson published the quantitative version of that argument in 2000, computing the relative economic value of the European ensemble prediction system against the deterministic forecast across the full range of cost to loss ratios. The ensemble is more valuable for most users, including at lead times where its deterministic scores are unremarkable. Which means the ranking of two forecast systems is not a property of the forecasts. It is a property of the pair consisting of a forecast and a user, and it flips as you walk along the range of users.

    This is episode six of an eight part walk through a literature map on physics AI and causal world models. Forecast to decision. A forecast has no value until something is decided with it, and the map from score to value is not monotone. You can reduce prediction error and change no decision at all. You can reduce prediction error and make decisions worse.

    It is worth being precise about what that curve is, because it is simpler than its reputation. For a given user, you take three numbers. The expense of always acting on climatology, which is the naive strategy. The expense of acting on a perfect forecast, which is the floor. And the expense of acting on the forecast system under test, following the optimal threshold rule for that user’s ratio. Relative economic value is where the system sits between the naive strategy and the floor, scaled so that climatology is zero and perfection is one. Sweep the ratio from near zero to near one and you get a curve rather than a number. Two systems produce two curves, and the curves cross. Everything contentious about comparing forecast systems is visible in where they cross and how wide each one’s useful range is, and none of it is visible in a single score.

    Tim Palmer extended the argument in 2002 across time scales, from weather out to seasonal and climate horizons, and made the stronger claim that the probabilistic forecast is the object carrying decision value. At a two week horizon the deterministic trajectory is a convenient fiction. At a seasonal horizon it has no defensible meaning at all, and the only honest product is a distribution. That is worth holding next to episode four, where we saw that the field’s headline scores are largely deterministic. The evaluation culture and the decision literature have been pointing in different directions for twenty years.

    The value framing has been run at industrial scale, and the results are instructive about where value actually lives. Wang and colleagues in 2016 ran production cost simulation over four hundred and eighty scenarios on a standard one hundred and eighteen bus test system, configured to emulate the Midcontinent, California, and New England balancing authorities, and priced what a given improvement in wind power forecast accuracy is worth at each market timescale. The finding is that the value tracks generation mix and available ramp capability rather than forecast error alone. Improve the forecast for a grid that cannot ramp fast enough to exploit the information and you have created very little. The same improvement delivered to a system with flexible generation is worth substantially more. So the value of a forecast improvement belongs to the system it lands in, not to the forecast.

    Matte and colleagues in 2017 attack a different assumption inside the same framework. The cost to loss ratio treats the decision maker as risk neutral, meaning they are indifferent between a certain loss and a gamble with the same expected loss. Operational users are not risk neutral. A reservoir operator does not weigh a small chance of catastrophic downstream flooding the way expected value says they should, and they are right not to. Matte and colleagues value ensemble streamflow forecasts under a constant absolute risk aversion utility instead, and the valuations change. The point generalizes past hydrology. The standard forecast value metric has a preference assumption baked into it, so a system tuned to maximize it is optimizing for a decision maker who may not exist.

    That is the meteorological tradition. Machine learning arrived at the same uncomfortable place from the other direction, independently, and gave it a different name.

    Priya Donti, Brandon Amos, and Zico Kolter published task based end to end model learning in stochastic optimization in 2017. The setup is that you have a prediction feeding an optimization problem, which is the structure of most operational decisions. The conventional approach trains the predictor on prediction error and then hands its output to the optimizer. Donti and colleagues differentiate through the optimization problem itself, so the predictive model is trained on the resulting task cost. Demonstrated on load forecasting for electricity scheduling and for battery storage. And the headline result is the one that should unsettle anyone who reports benchmark scores: a model with worse predictive accuracy can produce strictly better decisions. Not equal. Better. Because the errors it makes are in places the optimizer does not care about, and the errors the accurate model makes are in places it does.

    Adam Elmachtoub and Paul Grigas gave the framing a loss function with theory behind it in the same year. Their smart predict then optimize loss measures decision regret rather than prediction error, and they supply a convex surrogate with consistency guarantees for linear optimization problems. That is why it became the reference formulation. It is a tractable objective rather than an argument.

    Now the correction, because a decade of enthusiasm for decision focused learning needed one. Liu and colleagues in 2026 characterize when a plain accuracy trained predictor is already decision optimal, and when the decision focused objective is genuinely necessary. That work matters because train on the decision had started to circulate as an unconditional recommendation, and it is not one. The gain is conditional. Knowing the conditions is what turns it from a slogan into a method.

    Here is the tension that this theme leaves standing, and I do not think anyone has resolved it.

    The decision focused learning literature works on problems that are small and well specified. A scheduling problem with a known cost structure. A shortest path problem with a known objective. Under those conditions you can differentiate through the decision because you can write the decision down. Operational weather decisions are neither small nor well specified. There is no single objective function for a national meteorological service. Its outputs feed aviation, agriculture, energy trading, flood response, emergency management, and the public, and those users have different cost to loss ratios, different risk attitudes, different lead time requirements, and objectives that conflict with each other. Which decision do you train on?

    That question has no clean answer, and I want to say plainly that it may not have one. But it does have a partial answer that the field is not using, and this is the practical thing I would take from this theme. The relative economic value framework does not require you to pick a single user. It sweeps the entire range of cost to loss ratios and reports value as a curve. You can compute that curve for a learned model exactly as easily as for a physical one. Richardson did it in 2000, with the computing resources of 2000. It would settle a large fraction of the current arguments about whether learned ensembles are better than operational ones, because it would show for which users they are and for which users they are not.

    Almost none of the learned weather literature reports the downstream value of its own improvement. That absence is the gap in this theme, and it is not a research gap in the sense of requiring an invention. The instrument exists, it is twenty five years old, it is well understood, and nobody is picking it up. Reported skill improvements sit at some unknown place on a curve that nobody plots.

    There is a second observation worth making, which connects this theme back to episode four. Recall the double penalty, where a blurred forecast scores better than a sharply placed one that is displaced. From a decision perspective the blurred forecast is often worse, because a user acting on a spatially specific threshold gets nothing from a field that has been smoothed toward its mean. The verification literature knew this. Beth Ebert’s neighborhood methods are explicitly organized by the decision model each one implies. So the tools for connecting scores to decisions were built inside the verification community, and they are the tools that got left out when the benchmark culture formed around aggregate deterministic error.

    And there is a third, which is about what improvement even means at this point. A forecast that is more accurate by half a percent on a global average, for a user whose threshold sits in the tail of the distribution, may be worse for them than the previous system if the improvement came from smoothing. That is not a hypothetical. It is the mechanism connecting the missing extremes result from episode four to the value question here. Aggregate improvement and user relevant improvement can have opposite signs, and only one of the two gets reported.

    So what does this theme establish? That forecast value is user specific and computable, with a framework in place since 2000. That the value of an improvement belongs to the system that acts on it, measured in money on named electricity systems. That the standard value metric encodes a risk neutrality assumption operational users do not satisfy. That training on the decision objective can beat training on accuracy, demonstrated end to end, and that the gain is conditional in ways now characterized. And that the learned weather literature has adopted none of it.

    The experiment I would like to see is not complicated. Take two systems already evaluated head to head on a public benchmark, pick a handful of variables with genuine operational users behind them, ten meter wind, two meter temperature, precipitation, and compute relative economic value curves for both across the full ratio range at several lead times. Report where the curves cross. Report the ratio ranges where the learned system wins, where it loses, and where the difference is inside the noise. That is a paper someone could write this quarter from artifacts that already exist, and it would convert a long-running argument conducted in aggregate scores into a statement about which users are affected and by how much. The reason it has not been written, as far as I can tell, is that the two communities that would have to write it together do not share a venue.

    The through line for this series is that skill is not evidence of value, in the same way that prediction is not evidence of mechanism. Each of the last three episodes has been a version of the same structure. A number gets reported. The number supports a narrower claim than the one made. And the instrument that would support the wider claim already exists in a neighboring literature that the field has not read.

    Next episode we go a step further, to the case where the model is not just informing a decision but is the thing being planned inside. World models, planning, and physical control. Model based reinforcement learning worked out the failure mode precisely, which is that a planner will find the exact region where your model is wrong and optimistic and it will go there. We will look at what makes a learned model safe to act inside, at the two real deployments on physical systems, a stratospheric balloon and a tokamak plasma, and at what conditions those deployments needed that a continental atmosphere does not supply. That is episode seven.

  • 7. Acting Inside a Learned Model

    A planner will find the region where the model is wrong and optimistic. What the stratospheric balloon and tokamak deployments required, why rollout length is a safety parameter, and which conditions the open atmosphere fails.

    Read transcript 13 min · 2,123 words

    Let an optimizer plan inside a learned dynamics model and watch what it does. The optimizer is not searching for good actions in the world. It is searching for good actions according to the model, and those are different objectives. Anywhere the model is wrong in the optimistic direction, predicting an outcome better than reality would deliver, the planner will find it and go there, because that region is exactly where the objective looks best. The failure is not a bug in any particular architecture. It is what optimization does to an imperfect model, and the model based reinforcement learning community named it early and precisely. They call it model exploitation.

    That single dynamic organizes this entire theme, and it is the reason this episode is not about which model fits the dynamics best.

    This is episode seven of an eight part walk through a literature map on physics AI and causal world models. World models, planning, and physical control. The place where a learned model stops informing a decision and starts being the thing decisions are made inside.

    Start with the fix, because it is not what you would guess. Kurtland Chua, Roberto Calandra, Rowan McAllister, and Sergey Levine published the method usually called PETS in 2018, and the central move is epistemic. They use an ensemble of probabilistic dynamics models to separate two kinds of uncertainty that get conflated. Aleatoric uncertainty is noise inherent in the system, which more data will not remove. Epistemic uncertainty is ignorance, which more data will remove, and which tells you where the model does not know what it is talking about. Plan with trajectory sampling that respects that distinction and you match the asymptotic performance of model free methods with orders of magnitude fewer samples. The lesson worth carrying out of that paper is that what makes a learned model safe to plan against is its uncertainty representation rather than its one step accuracy. A model that is slightly worse at prediction and honest about where it is ignorant is more useful for control than a model that is slightly better and uniformly confident.

    There is a corollary that took the field a while to accept. Rollout length is a safety parameter. Every step you roll forward inside a learned model compounds its error, and the theoretical bounds on how much a policy improved inside the model improves in the real world degrade with horizon. So the practical recipe that emerged is short rollouts from real states, repeatedly, rather than long imagined trajectories from a single starting point. That is a strange conclusion if you came in thinking of a world model as a simulator. A simulator is valuable in proportion to how far it can run. A learned dynamics model used for control is trustworthy in inverse proportion to how far you run it, and the engineering discipline is knowing where to cut.

    The architectural origin point for this line is David Ha and Jürgen Schmidhuber’s World Models in 2018, which is where the phrase this explorer is partly named after entered circulation. Compress observations into a latent representation, learn a recurrent forward model of the latent dynamics, then train a small controller entirely inside that learned model and transfer it back to the real environment. Train the policy in the model’s imagination. It is a clean idea and it made the subsequent literature possible.

    Then the strongest version of the argument. Julian Schrittwieser and colleagues built MuZero, which learns a latent dynamics model trained only to predict reward, policy, and value, and then runs tree search inside it. It matches AlphaZero at Go, chess, and shogi without ever being told the rules of those games, and sets a new state of the art on Atari at the same time. If you want a demonstration that planning inside a learned model can reach superhuman performance in a domain whose mechanics were never supplied, that is the paper.

    It is also, read carefully, a statement of its own preconditions, and this is the part that gets skipped when it is cited in physical settings. MuZero requires a resettable environment, so you can return to a state and try something else. It requires cheap exact rollouts, so you can afford millions of them. And it requires a scalar reward, so the objective is unambiguous. Go supplies all three. Danijar Hafner and colleagues push the same architecture much further with the third version of Dreamer, a single world model and actor critic with fixed hyperparameters across more than one hundred and fifty tasks, spanning continuous control, Atari, and Minecraft, including collecting diamonds from scratch, which had resisted everything before it. That is the current reference point for how far world model planning generalizes when the environment is simulated, well instrumented, and resettable.

    A continental scale atmosphere supplies none of those three conditions. It is not resettable. Its rollouts are neither cheap nor exact. And there is no scalar reward, as episode six established at length, because the objective is a family of conflicting user specific objectives. So when the reinforcement learning results are cited as evidence that we are close to planning inside learned physical world models, the honest reading is that they establish what is possible under conditions that the target domain does not meet.

    Philipp Holl, Vladlen Koltun, and Nils Thuerey provide the cleanest result for systems governed by partial differential equations. They use a differentiable solver to train a predictor corrector control network end to end, planning through the solver rather than around it, and achieve long horizon control of fluid systems. That is a genuine demonstration of learned control on the kind of dynamics we care about. It is also a useful upper bound rather than an operating point, because it works with full state access and an exact simulator. Both of those are things the open atmosphere does not give you, and episode one was largely about why.

    Now the two real deployments, and I want to spend time on what each of them required rather than on what each achieved, because the requirements are the transferable content.

    Marc Bellemare and colleagues published autonomous navigation of stratospheric balloons using reinforcement learning in Nature in 2020. The controller steers a balloon by choosing altitude, riding wind layers that move in different directions, and it outperformed the hand engineered controller in flight over the Pacific under real wind uncertainty. That is learned control acting on a real atmospheric system, which makes it the closest thing in this literature to the thing people imagine when they talk about AI controlling weather. Notice how it works. The controller does not modify the atmosphere in any way. It exploits an atmospheric structure that already exists, selecting among wind layers rather than changing them. The action space is a single scalar, altitude. The state is a balloon. Success is measured by staying near a location. Every one of those properties is what made the problem tractable.

    Jonas Degrave and colleagues published magnetic control of tokamak plasmas through deep reinforcement learning in 2022. Train a controller in a physics simulator of the TCV tokamak, deploy it on the real device, and sustain and shape plasma configurations including ones that had not been demonstrated before. As a sim to real result on a physical system this is the strongest in the literature and I do not want to undersell it. It depended on three things. A trusted simulator, meaning decades of plasma physics encoded in a model good enough to train against. Dense diagnostics, meaning the plant is instrumented to a degree that few physical systems are. And millisecond actuation, meaning the control loop closes fast enough to correct model error before it compounds.

    That last one is the general principle and it deserves stating on its own. The two successful deployments both closed their loop faster than model error could accumulate. Model exploitation is bounded when you can observe the consequence of an action and correct quickly. It is unbounded when you act and then wait. For a continental atmosphere, the observation is partial, the actuation is undefined, and the feedback horizon is days. That combination inverts every condition that made the balloon and the tokamak work.

    The third transfer argument in circulation comes from autonomous driving, where the perceive predict plan stack demonstrably works at scale, and its own literature is the best check on it. Mayank Bansal, Alex Krizhevsky, and Abhijit Ogale reported in ChauffeurNet that thirty million real expert driving examples were not sufficient on their own to produce a policy that could drive a real vehicle, even with a perception system upstream and a controller downstream. What was required was synthesizing failures the expert never made, including collisions and road departures, and adding explicit losses penalizing them. Expert data alone taught the model what good driving looks like and did not teach it what to do when things go wrong, because the expert never went wrong. The fix was to manufacture counterfactual experience.

    Which is the same finding as the imitation learning result we met in episode five. Pim de Haan, Dinesh Jayaraman, and Sergey Levine showed that a policy cloned from demonstrations suffers causal misidentification, and that it can get worse as its observation set grows, because confounded cues correlating with expert actions are easier to fit than the actual causes. Only interventional data disambiguates. Two independent lines, one from industrial engineering practice and one from causal analysis, arriving at the identical conclusion: passive observation of expert behavior does not yield a policy, and the missing ingredient is intervention.

    Put the theme together and it makes an argument that is more coherent than any individual paper in it.

    Planning inside a learned model is real, it is powerful, and it has explicit preconditions. The preconditions are a resettable or fast correcting environment, a bounded and characterized action space, dense instrumentation, and a clear objective. Where those hold, results range from superhuman to operationally deployed. Where they do not hold, the central failure mode is a planner exploiting model error, and the mitigation is representing uncertainty honestly rather than fitting the dynamics better.

    The open atmosphere fails the preconditions on every axis, and it is worth saying which axis is worst. It is not model quality. Learned atmospheric models are good, as episode two established. It is that there is no characterized actuation. Nobody has a transfer function from an intervention to an atmospheric response with error bars on it. The balloon case is instructive precisely because it sidesteps this: it acts on a balloon, not on the sky.

    So the gap I would flag here is narrower and more useful than a general statement that we are far away. Every result in this theme that touches a physical system had a characterized actuator and a fast loop. The research question is what a control problem looks like when the actuation is uncharacterized and the feedback is slow, and whether uncertainty aware planning of the kind PETS established can be extended to that regime with any guarantee at all. That question is not being worked on inside the weather machine learning community, and the reinforcement learning community is not aiming at it either.

    There is one more thing worth extracting from the two deployments, which is what they imply about evaluation. Neither the balloon controller nor the tokamak controller was validated by showing that its internal dynamics model predicted well. They were validated by running the closed loop on the real plant and measuring what happened. That is a different and much stronger form of evidence than anything in the first six episodes, because the environment gets a vote. A forecast benchmark cannot push back on a wrong model. A plasma can. If there is a single methodological transfer available from the control literature to the physics machine learning literature, it is that closing the loop on a real system is the evaluation that cannot be gamed, and that where you cannot close the loop, you should be correspondingly less confident in what your model knows.

    Which brings us to the last episode, and to the claim that sits at the top of this whole literature map. If a learned model can be interrogated about interventions, then it is an instrument for asking what if, which is a different kind of object from a forecast system. Companies are being built on that premise right now. Next episode, intervention and the boundary of evidence. What the strongest claims in this space actually assert, what evidence would settle them, which experiments the literature has already run, and which ones nobody has run yet. That is episode eight, and it is where the whole series has been going.

  • 8. The Boundary of Evidence

    Seventy years of weather modification, one physical chain observed end to end and measured in fractions of a millimeter, and what today's AI steering studies establish about the model rather than the atmosphere.

    Read transcript 14 min · 2,109 words

    In 2026 a group published a paper showing that a tropical cyclone can be steered by a small perturbation to the initial conditions. Change a modest amount of the atmospheric state in the right places, run the model forward, and the track moves. A second group did a more sophisticated version, targeting the perturbation at the dynamically unstable directions of the flow rather than searching by brute force gradient descent, and steered an extreme atmospheric river.

    Both results are real. Both are entirely about the model. What they establish is that a learned weather model is differentiably sensitive to small changes in a targeted direction, which is a property of a function. Whether the atmosphere is sensitive in the same way, and whether any physically deliverable action could produce the perturbation in question, are separate questions that neither paper asks and that no experiment in the literature currently answers.

    This is episode eight, the last of a walk through a literature map on physics AI and causal world models. Intervention and the boundary of evidence. The theme where the claims are largest and the evidence is oldest.

    Because there is a seventy year record here, and the AI weather community has mostly not read it.

    Project STORMFURY ran hurricane seeding experiments from 1962 to 1983. Hugh Willoughby and colleagues wrote the chronicle, and what killed the program is worth stating precisely, because it is not what people assume. It did not fail on execution or on funding or on nerve. It failed on two things. The physical premise turned out to be wrong: hurricanes contain far too little supercooled water for the hypothesized seeding mechanism to operate. And the apparent successes fell within the range of natural intensity fluctuation, so there was no way to attribute them.

    Bernard Silverman’s 2001 assessment of glaciogenic seeding of convective clouds reaches a parallel conclusion by a different route. He reviews the randomized experiments against explicit criteria for statistical and physical evidence and concludes the accumulated results do not establish an effect, with several apparently positive results failing on replication. Decades of operational seeding activity coexisted with no demonstrated effect. Not with a demonstrated absence of effect. With no demonstration either way.

    The National Research Council’s 2003 assessment is the authoritative version and it draws the distinction that matters more than any other in this theme. There is no convincing scientific proof of efficacy for intentional weather modification. The underlying questions remain scientifically legitimate. What limits them is experimental design, statistical power, and the counterfactual that does not exist, because you cannot run the unseeded version of the same storm.

    Hold onto that. Not demonstrated is a different claim from disproven, and the reason for the difference is a measurement problem rather than a physics problem. Which is exactly the kind of problem that better prediction, better state estimation, and better causal inference could in principle change. That is the whole reason this theme belongs in a literature map about learned physical models.

    There is one case where the chain has been closed. Jeffrey French and colleagues, working in the SNOWIE campaign, used airborne and ground based radar to track a glaciogenic seeding event end to end. Silver iodide released into an orographic cloud, ice initiation, particle growth, snowfall reaching the ground, with the seeding lines visible in the radar field as they propagated. Observed, not inferred from a statistical comparison across storms. It is a genuinely beautiful piece of measurement.

    Katja Friedrich and colleagues then quantified it, estimating the mass of snow attributable to seeding over specific areas and durations. The number is on the order of fractions of a millimeter of water equivalent. That is the effect size for the one intervention whose mechanism has been observed working. Anyone reasoning about weather intervention should have both halves of that in mind at once: the mechanism is real and the demonstrated magnitude is very small, in the narrowest and most favorable configuration anyone has managed to instrument.

    Robert Wood’s 2021 analysis of marine cloud brightening is the methodological model I would hold up here. Before proposing to test an intervention, he estimates its achievable radiative forcing with a transparent heuristic model, and shows the answer is dominated by uncertainties in aerosol cloud interaction rather than by deployment scale. Bound the plausible effect first. Find out what dominates the uncertainty. Then decide whether the experiment is worth designing. That discipline exists, it is standard in the climate intervention literature, and it is absent from the AI steering genre entirely.

    So put the steering studies against that background. Huang and colleagues, and Liu and colleagues, produce optimization results inside a surrogate that was trained exclusively on unperturbed reanalysis. Which means the intervention regime is precisely the regime with no training support. Recall the finding from episode five, that models trained on trajectories fail to apply the mechanics underneath, and from episode four, that learned models fail to extrapolate to unseen intensities of phenomena they have seen thousands of times. A model with those documented properties is being asked what happens under a state it has never seen, and its answer is being reported as an intervention.

    There is a piece of history that makes the steering genre look different, and it belongs here. Finding the perturbation that most efficiently changes a forecast is not a new capability. Operational centers have been computing singular vectors and running adjoint sensitivity since the 1990s, precisely to identify the fastest growing directions in the flow, and they use them for targeted observation, meaning they decide where to send an aircraft to sample because that is where uncertainty matters most. Episode one covered Carla Cardinali’s use of the same adjoint machinery to price individual observations. So the atmosphere’s sensitivity structure has been studied with physics based tools for thirty years, and the field that studied it drew the conclusion that this structure tells you where to measure. The learned model version rediscovers the sensitivity through a differentiable surrogate and draws the conclusion that it tells you where to act. The mathematics is similar. The inference is not, and the difference is that a targeted observation is an action whose effect you can verify while a targeted perturbation is not.

    The two gaps in this theme follow directly and both are stated as experiments rather than as complaints.

    No AI era steering study has been paired with a field testable hypothesis, an actuation model connecting the perturbation to something physically deliverable, or a power analysis for detecting the claimed effect. Those three are what would turn a gradient into a proposal. The perturbation the optimizer finds is a vector in state space. To be an intervention it has to be reachable by some physical action, and nobody has written down the map.

    And there is no validation protocol for a learned model’s response to interventions absent from its training distribution. That one is embarrassingly tractable. Take the perturbation a learned model recommends. Apply it inside a physics based model. Compare the trajectories. Characterize how the disagreement grows with perturbation size, with lead time, and with distance from the training distribution. That experiment requires no new theory and no new data, only the willingness to check. Until somebody runs it, every steering result is a claim about a function.

    It is worth asking what learned models could genuinely contribute to the counterfactual problem, because there is an answer and it is not steering. The reason STORMFURY and the convective seeding programs could not resolve their effect sizes is that natural variability swamped the signal and there was no control storm. Everything that reduces the variance of the expected outcome makes that detection problem easier. A better analysis reduces initial condition uncertainty. A larger and cheaper ensemble characterizes the natural spread more finely, and learned models made large ensembles cheap, which is not a small thing. A calibrated probabilistic forecast supplies the distribution of what would have happened, which is the closest available substitute for a control. None of that is intervention. All of it is the statistical infrastructure an intervention experiment would need, and it is the part of this literature that actually bears on the question the NRC identified. The contribution runs through improving the counterfactual estimate rather than through finding perturbations inside a surrogate.

    Which is the right point to close the series, because the same structure has appeared in all eight episodes.

    We started at the sensor boundary, with the fact that a learned atmospheric model has never seen the atmosphere, only instruments, and that credit for its skill is shared with an observing system and an assimilation chain it did not build and cannot reproduce. Then the prediction results, which are genuine: the learned models won that argument, reproducibly, and are now running operationally inside the centers they were once benchmarked against. Operator learning gave us one clean general property, resolution invariance, plus measurable transfer in the low data regime, alongside the finding that much of the surrounding speedup literature does not survive a fair comparison.

    Episode four was the turn. The measurements behind all of it are weaker than they look, in the tails, under distribution shift, and when the developer and the evaluator are the same person. Episode five asked whether any of it implies mechanism, and found that prediction pressure does produce internal representations aligned with physics nobody supplied, while the direct tests of the Kepler to Newton premise come back negative. Episode six pointed out that forecast value is user specific, computable, and uncomputed. Episode seven found that planning inside a learned model works under conditions including a fast correcting loop and a characterized actuator, neither of which the open atmosphere supplies. And intervention, this episode, has a long empirical record whose limiting factor was always the counterfactual, which the current AI results have not engaged.

    The pattern is consistent. At every step, a real and impressive capability is demonstrated, and then a claim is made one rung higher on the ladder than the evidence reaches. That ladder, from episode five, is worth repeating because it is the most portable thing in this map. Predictive relationship. Internal representation. Identified causal structure. Valid counterfactual. Support for intervention. Most claims about learned physical understanding are rung two evidence carrying rung four or rung five language.

    I want to be clear that this is not a case for pessimism. Every one of those gaps is an experiment somebody can run, and several of them are cheap. Hold a model fixed and vary the observing system, dropping a satellite channel or degrading station density to a historical level, and you get the first attribution of forecast error between the sensors, the analysis, and the model. Build surrogate fields with matched autocorrelation and seasonality but no causal coupling, run them through the standard causal discovery pipeline, and you get the empirical false positive rate that the Earth system causal literature currently reports without. Probe a learned forecaster’s internal representations and then ablate them, and you find out which ones are load bearing rather than merely decodable. Compute a relative economic value curve for a modern ensemble system and a named user, and you convert a score improvement into a statement about outcomes. Publish a calibrated error bound over a reachable state set for one variable, one region, one lead time, and you produce the first honest interface between the learned modeling literature and the control theory that any intervention proposal would need.

    None of those require a new architecture. Several of them are a few months of work on artifacts that already exist. What they have in common is that each one converts a claim into a measurement, and the field’s bottleneck right now is measurement rather than capability.

    The last question is the one the NRC assessment left open in 2003 and that nobody has revisited with modern tools. Given current observing density and a learned model’s demonstrated sensitivity, what randomization design, what sample size, and what detection threshold would make a seeding effect of the SNOWIE magnitude measurable in the field? That is an answerable question. It might have a discouraging answer. Either way it would be the first time in twenty years that anyone checked whether the terms of the problem have changed.

    That is the map. Eight themes, about eighty papers, and a set of experiments that would move claims from asserted to established. The explorer this series follows is on the site, with every paper linked, each with a takeaway and a note on why it belongs. If you want to argue with any of this, the papers are the place to start.

Open problems

Where the literature is thin and the next contribution could land.

  1. Attribute forecast error between the observing system, the analysis, and the model

    Every reported score folds three error sources into one number. WeatherReal and the Ullrich recommendations both push evaluation toward in-situ observations, and neither supplies a decomposition. The experiment is available: hold the model fixed and vary the observing system (drop a satellite channel, inject a realistic calibration drift, degrade station coverage to a historical density), then measure how much of the score change is attributable to each. Instrument-level provenance exists in the assimilation systems; nothing currently connects it to learned-model verification. The same design answers the marginal-information question the field keeps assuming away: operational assimilation already measures forecast impact per observation, and no published curve shows where a learned model's skill saturates as observing density rises.

  2. Test whether physics foundation models encode reusable structure or only fit the target family

    Transfer is currently argued from held-out-PDE accuracy in the low-data regime. That is consistent with a model learning a good initialization rather than physical structure. A representation-level protocol would test the properties that would have to be there: conserved quantities recoverable by linear probe, correct dimensional scaling under nondimensionalization, and composition of learned operators matching composition of the underlying ones. The Well supplies the system diversity to run it; the protocol does not exist.

  3. Establish negative controls for spatiotemporal causal discovery

    PCMCI-family methods are applied at scale to Earth-system data whose properties violate their assumptions, and the results are reported without an empirical false-positive rate. What is missing is the equivalent of a scrambled-label control: surrogate fields with matched autocorrelation, spatial structure, and seasonality but no causal coupling, run through the same pipeline, so a discovered network can be compared against what the method finds in the absence of causation.

  4. Validate counterfactuals from learned models against a physical simulator

    AI steering studies optimize perturbations inside a model trained only on unperturbed reanalysis, so the intervention is exactly the regime with no training support. The validation is straightforward to specify and has not been run: apply the perturbation the learned model recommends inside a physics-based model, compare the trajectories, and characterize how the disagreement grows with perturbation size, lead time, and how far the perturbed state sits from the training distribution.

  5. Give extreme-event evaluation a power analysis

    Extreme-event benchmarks are proliferating without a stated answer to how many events are needed to detect a given improvement in tail behavior. Large ensembles help with sampling and do not resolve the question, because the reanalysis truth for the rarest events is itself model-dependent. A worked power analysis, with an explicit treatment of truth uncertainty in the tail, would tell the field which of its extreme-event comparisons are currently informative.

  6. Report the value curve, not just the score

    The relative-economic-value framework is twenty-five years old, takes an ensemble and a cost-loss ratio as input, and has not been applied publicly to GenCast-class systems for any named user. Running it for a specific decision (a grid operator's reserve commitment, an airline's fuel loading) would convert a CRPS improvement into a statement about outcomes, and would show which parts of the learned-forecasting advantage a real user can access.

  7. Produce an error bound a controller can consume

    Safe learning-based control requires a bound on model error over the reachable state set. No learned Earth-system model publishes anything of that form, which means every proposal to act on one is missing the input its own control theory requires. Even a restricted result (a calibrated bound for one variable, one region, one lead time, one class of perturbations) would be the first honest interface between this literature and the intervention literature.

  8. Design a weather intervention experiment that could actually detect its effect

    The programs that ended did so on statistical power, not physics: the unseeded counterfactual does not exist, and natural variability swamped the signal. Better state estimation and observation-driven modeling change the terms of that problem for the first time since the NRC assessment. The open question is concrete: given current observing density and a learned model's demonstrated sensitivity, what randomization design, sample size, and detection threshold would make a seeding effect of the SNOWIE magnitude measurable, and is that number reachable?