The tool, called METEORv1.6, is the latest evolution of a climate emulation framework first introduced by the same group. At its heart lies a technique known as pattern scaling, which decomposes a climate model’s response to forcing into spatial patterns that evolve over characteristic timescales. Where traditional pattern scaling assumes that the same spatial fingerprint of warming applies regardless of what caused it, METEOR takes a more sophisticated view. It represents the climate response as a sum of impulse response patterns, each with its own emergence timescale, driven by the forcing trajectory of each agent separately. Greenhouse gas patterns are derived from an idealized experiment in which carbon dioxide concentrations are abruptly quadrupled, while aerosol responses are estimated from residual signals in concatenated historical and future simulations. The result is a model that can capture hysteresis effects and time-evolving responses, which are essential for scenarios involving temperature overshoots or shifting mixes of short-lived and long-lived forcers.

Version 1.6 adds two capabilities that the original framework lacked. The first is a monthly climate variability model, dubbed METEOR-NOISE, which generates realistic sub-annual climate sequences complete with seasonal cycles and year-to-year fluctuations. The second is a modular impact assessment framework that translates those climate projections directly into sector-relevant metrics. The initial implementation focuses on heating and cooling degree days, quantities that drive estimates of energy demand for buildings and that cannot be meaningfully computed from annual mean temperatures alone. Together, these additions turn METEOR from an annual-resolution pattern scaler into an end-to-end pipeline that runs from emission pathways to probabilistic impact estimates at the monthly resolution that sectoral applications actually require.

The technical machinery behind METEOR-NOISE is a two-stage statistical construction. First, the model characterizes how the seasonal cycle itself changes as the planet warms. At every grid point, a harmonic regression captures the annual and semi-annual cycles using sine and cosine functions, but crucially the amplitudes and phases of those harmonics are allowed to modulate linearly with smoothed global mean temperature. This temperature-dependent parameterization lets the model represent genuinely non-stationary seasonal behavior, such as the well-documented pattern in which Arctic winters warm far faster than summers as sea ice retreats, shifting both the amplitude and the timing of the seasonal cycle. The regression uses just nine features per location, a remarkably compact description that nevertheless captures spatially varying changes in seasonality across the globe.

What remains after the seasonal signal is removed are the monthly anomalies: the chaotic, weather-driven fluctuations that no deterministic model can predict individually but whose statistical structure can be learned. METEOR-NOISE applies Principal Component Analysis to these anomalies, identifying by default forty spatial patterns that together explain the bulk of the variance. For surface temperature, the seasonal model accounts for roughly 91 percent of total variance, and the forty principal components capture 75 to 80 percent of what remains, for a combined total near 98 percent. Precipitation, being inherently noisier and strongly skewed, is harder: the seasonal model explains about 32 percent of its variance, and the components add another 29 to 45 percent. The leading patterns often project onto recognizable modes of variability such as the El Niño-Southern Oscillation or the North Atlantic Oscillation, though the authors are careful to note that the components are treated as a statistical basis for reproducing covariance rather than as physically labelled processes.

The temporal evolution of those forty pattern amplitudes is then modeled with a Vector Autoregressive process with exogenous variables, or VARX. Each month’s value of every component depends on the previous two months’ values of all components, plus global temperature and a seasonal signal as external predictors, plus a random shock drawn from a multivariate Gaussian distribution whose covariance matrix is learned directly from the training data. This design captures three things at once: persistence, since anomalies tend to linger for months; teleconnections, since activity in one spatial pattern can influence others with a lag; and modulation by the warming trend itself. It is a meaningful departure from approaches that apply independent autoregressive processes at each grid point, because the cross-mode dependencies preserve the spatial coherence that makes regional climate statistics realistic.

Training is strikingly frugal. The model learns from a single ensemble member of a CMIP6 Earth system model, concatenating the historical simulation with a mid-range future scenario to build a record of roughly 250 years, supplemented by the carbon dioxide quadrupling experiment for the core pattern scaling. Once trained, generation is fast: producing one hundred ensemble members of regional or global mean temperature spanning 350 years at monthly resolution takes about ninety seconds on a standard laptop, and point-based extraction for a single city takes roughly ten seconds, meaning thousand-member ensembles can be generated in under two minutes. Training itself takes between thirty and sixty seconds. The trick is that regional means can be computed directly from the forty principal component time series without ever reconstructing full spatial fields, reducing memory demands by two to three orders of magnitude.

Validation against CMIP6 output suggests the framework earns its keep. When trained on the SSP2-4.5 scenario of the Norwegian Earth System Model, the emulator reproduces the distribution, spread, and seasonal structure of the source model’s output across global, regional, and city scales, and captures long-term trends through 2100. In city-level tests, the emulator reproduces the shrinking seasonal temperature cycle in Delhi under warming and the near-total absence of cooling demand in Oslo until the middle of the century, a threshold behavior it was never explicitly trained to represent. More demanding still, when trained only on SSP2-4.5 and then asked to project scenarios it had never seen, including a strong mitigation pathway, a high-emissions pathway, and a temperature overshoot scenario, the framework tracked the non-monotonic temperature trajectories of five different Earth system models across all spatial scales. Temperature skill ratios comparing emulated members against the actual model output cluster near unity, indicating the realizations are statistically indistinguishable from genuine ensemble members.

The framework is candid about its limits. Precipitation changes under the most extreme emissions scenario are somewhat underestimated, reflecting the stronger nonlinearity of hydrological responses, and some regional precipitation ensembles prove slightly underdispersive. The Gaussian noise assumption likely understates the probability of extremes, with diagnostics showing heavier-than-Gaussian tails, and the authors suggest future extensions using t-distributed innovations or empirical resampling. Nonlinear catastrophes such as a collapse of the Atlantic Meridional Overturning Circulation, and novel interventions like solar radiation management, remain beyond the reach of any pattern scaling approach. Daily and sub-daily generation, needed for extreme event statistics, is on the roadmap.

What makes METEORv1.6 potentially transformative is less any single statistical trick than the combination: a physically motivated, forcing-aware core paired with a statistically calibrated variability model and an extensible impact layer, all open-source and fast enough for interactive use. The authors envision it embedded within integrated assessment models, supplying spatially resolved, probabilistic climate feedbacks in real time as users explore alternative emission pathways. In a field where the distance between an emission choice and its regional consequences has traditionally been measured in months of supercomputing time, a tool that makes that relationship explorable on demand represents a genuine shift in what climate risk assessment can look like.

Subject of Research: Development and validation of the METEORv1.6 climate emulator for monthly climate variability and integrated impact assessment

Article Title: METEORv1.6: Spatial climate variability and integrated impact emulation

Article References: Sanderson, B. M., Sandstad, M., Dewey, M., Steinert, N. J., & Fjeldså, J. L. (2026). METEORv1.6: Spatial climate variability and integrated impact emulation. Geoscientific Model Development, 19(19), 9263-9288. https://doi.org/10.5194/gmd-19-9263-2026

Image Credits: AI Generated

DOI: 10.5194/gmd-19-9263-2026

Keywords: climate emulation, pattern scaling, METEOR, CMIP6, seasonal cycle, internal variability, PCA, VARX, degree days, energy demand, climate risk assessment, Geoscientific Model Development