Research
Learning to Factorize and Adapt: A Versatile Approach Toward Universal Spatio-Temporal Foundation Models
Overview Research area: Machine learning for spatio-temporal (ST) data — specifically, spatio-temporal foundation models (STFMs) that pretrain on large multi-domain corpora and transfer to unseen fore
- arXiv
- 2601.12083
- Published
- 2026-01-17
- Authors
- Siru Zhong, Junjie Qiu, Yangyu Wu, Yiqiu Liu, Yuanpeng He, Zhongwen Rao, Bin Yang, Chenjuan Guo, Hao Xu, Yuxuan Liang
AI summary
Overview
Research area: Machine learning for spatio-temporal (ST) data — specifically, spatio-temporal foundation models (STFMs) that pretrain on large multi-domain corpora and transfer to unseen forecasting tasks.
Technical level: Advanced. The paper assumes familiarity with transformer architectures, self-supervised pretraining, graph neural networks, positional embeddings (RoPE), quantile regression, and parameter-efficient adaptation (prompting, adapters).
Scope: The paper proposes FactoST-v2, a two-stage framework that separates universal temporal pretraining from lightweight spatial adaptation, arguing this factorization is both more scalable and more transferable than joint ST pretraining. (The provided text is truncated before the experimental results section, so numerical outcomes cannot be reported here.)
What This Paper Is About
Existing spatio-temporal foundation models try to learn spatial and temporal patterns together in a single large model pretrained across many domains. The authors argue this "joint" approach is fundamentally mismatched: temporal dynamics (trends, seasonality) generalize well across domains, but spatial relationships depend on domain-specific topologies (road networks vs. power grids vs. sensor grids) that cannot be unified into one manifold. FactoST-v2 instead factorizes the problem — first learn a shared, graph-agnostic temporal backbone, then inject spatial awareness through a small adapter per domain.
Key Contributions
-
A factorized paradigm for STFMs. The framework decouples Universal Temporal Pretraining (Stage I) from Spatio-Temporal Adaptation (Stage II), circumventing the quadratic complexity of joint ST pretraining and making components architecturally plug-and-play.
-
A minimalist encoder-only backbone with arbitrary-length probabilistic modeling. The first stage replaces the earlier encoder–decoder design of FactoST-v1 with an encoder-only architecture using a learnable
[REG]token, randomized sequence masking, partial rotary position embeddings, gated attention, and a multi-quantile prediction head. This removes the fixed-horizon constraint and enables 100% weight transfer between stages. -
A streamlined adaptation module. Stage II discards FactoST-v1's heavy hierarchical alignment in favor of Spatio-Temporal Metadata Fusion, Spatio-Temporal Filtering (spatial, temporal, and time-lagged affinities), low-rank Domain-Specific Prompt Alignment, and Continual Memory Replay.
-
A broadened evaluation protocol. The journal version extends the NeurIPS 2025 conference paper with zero-shot and full-shot settings, scaling-law analysis over data and model size, interpretability and universality studies, and efficiency profiling, in addition to few-shot evaluation.
Main Findings
-
Pattern mismatch motivates factorization: The paper argues that temporal patterns share universal 1-D structure across domains (e.g., a seasonal cycle in electricity load has spectral similarities to traffic flow), whereas spatial correlations are bound to domain-specific topologies and physical distances (short-range diffusion between neighboring air-quality stations in Beijing vs. Pacific-scale teleconnections in climate indices).
-
Joint ST pretraining scales poorly: Flattening spatial and temporal dimensions into one sequence incurs complexity of roughly O(N²P²D), where N is the number of nodes and P the patch count. This produces memory and time footprints that grow quadratically and can trigger over-squashing at scale.
-
The factorized design is linear in node count: Stage I costs O(NP²D) because the backbone processes node-wise series independently; Stage II adds only O(NPD) plus O(ΔNPMD) with latent prototypes M ≪ N. Total adaptation complexity is O(N).
-
Full weight transfer is achieved: By eliminating the decoder, FactoST-v2 removes the structural gap that forced partial parameter transfer in v1, enabling the pretrained backbone to be reused intact downstream.
-
Probabilistic prediction replaces point estimation: The quantile head (optimized with pinball loss) lets the model produce prediction intervals and quantify uncertainty, which the deterministic v1 could not do. For fair comparison with existing point-forecasting baselines, the authors fall back to L1 loss on the median quantile.
-
Reported performance claims (from the abstract only): The authors state that FactoST-v2 attains state-of-the-art accuracy with "linear efficiency," outperforming existing foundation models in zero-shot and few-shot scenarios while rivaling domain-specific expert baselines. Specific benchmark numbers are not present in the supplied text.
Methodology in Plain English
The approach is a two-step recipe.
Step 1 — Learn time, ignore space. Take time series from many domains (traffic, energy, climate, etc.), normalize each window, and break it into patches. Crucially, treat each node's series independently, so no graph is involved. Train a small transformer encoder to predict future values by hiding (masking) part of the input — the amount hidden is randomized each sample, so the model sees variable-length contexts and learns a length-agnostic mapping. A special [REG] token separates the observed past from the placeholder future, letting the same head handle any horizon via truncation or rolling prediction. Instead of absolute position embeddings (which don't extrapolate), the model uses partial rotary embeddings that keep low-frequency components untouched, preserving signal amplitude. A gating mechanism suppresses noise. The output is not one number but several quantiles, trained with pinball loss so the model learns a distribution rather than a point.
Step 2 — Teach space, cheaply. For a target domain, freeze the backbone and attach a small adapter. It fuses node identifiers and calendar features into the temporal features, then reweights them through three learned affinity scores: how spatially relevant the context is, how temporally aligned it is, and how much delayed (time-lagged) influence matters — the last captured efficiently through a small set of latent prototypes rather than all-pairs attention. Global distribution shift is handled by a handful of learnable prompt tokens factored into low-rank matrices and prepended to the sequence. A memory buffer replays older samples during adaptation to prevent catastrophic forgetting.
The authors support the design with a complexity analysis (linear vs. quadratic in node count) and reference generalization-bound arguments, though the bound derivation is cut off in the provided text.
Why This Matters
Impact on research: The paper challenges the prevailing "one big joint model" assumption in spatio-temporal foundation models and offers a concrete alternative with better scaling properties. If the factorization holds up broadly, it reframes how the community thinks about ST pretraining — as two separable problems rather than one monolithic one. It also connects ST modeling to the parameter-efficient fine-tuning literature (prompts, adapters, low-rank factorization).
Real-world applications:
- Traffic forecasting on road networks, where sensor topologies differ from city to city and labels for new cities are scarce.
- Energy load forecasting on power grids, where seasonal/temporal patterns are shared but grid topology is utility-specific.
- Air quality and environmental monitoring, where sensor networks are sparse, heterogeneous, and often unlabeled in new regions.
- Climate and weather indices, where long-range teleconnections coexist with local station dynamics.
Industry relevance: A backbone with linear scaling in the number of sensors is directly deployable at city or national scale, where joint-attention models hit memory walls. Few-shot and zero-shot adaptation matters commercially because deploying a new sensor network usually means little or no historical labeled data. The lightweight adapter also means a single pretrained model can be maintained centrally and specialized cheaply per customer or region.
Future Directions
-
Validate generalization bounds empirically. The theoretical section argues factorization improves downstream accuracy via complexity arguments, but the provided text truncates before the bound is fully stated; empirical confirmation across many domains would strengthen the claim.
-
Quantify negative transfer boundaries. Factorization is presented as avoiding the negative transfer of joint pretraining, but it is unclear when spatial patterns are too complex for the lightweight adapter to capture, and whether a fallback to heavier spatial modeling is ever needed.
-
Extend beyond forecasting. The paper uses forecasting as the sole proxy task. Whether the same factorized backbone transfers to imputation, anomaly detection, and classification is an open question.
-
Test scaling limits. The scaling-law analysis covers data and model size, but behavior on truly massive graphs (millions of nodes), and how the latent-prototype count M should scale with N, remain open.
-
Assess prompt and memory sensitivity. The low-rank prompt rank r and memory replay ratio are hyperparameters; their robustness across domains and their interaction with the continual-learning setting deserve systematic study.
Target Audience
This paper is most useful to machine learning researchers and graduate students working on time-series foundation models, spatio-temporal graph neural networks, and urban computing. It is also relevant to practitioners in transportation, energy, and environmental analytics who need cross-domain forecasting with limited labeled data, and to engineers concerned with the inference and memory cost of deploying large ST models at scale. Readers without background in transformer pretraining or graph neural networks will find the methodology dense; the high-level two-stage idea is accessible, but the technical detail is not beginner-friendly.
Authors’ abstract
Spatio-Temporal (ST) Foundation Models (STFMs) promise cross-dataset generalization, yet joint ST pretraining is computationally expensive and grapples with the heterogeneity of domain-specific spatial patterns. Substantially extending our preliminary conference version, we present FactoST-v2, an enhanced factorized framework redesigned for full weight transfer and arbitrary-length generalization. FactoST-v2 decouples universal temporal learning from domain-specific spatial adaptation. The first stage pretrains a minimalist encoder-only backbone using randomized sequence masking to capture invariant temporal dynamics, enabling probabilistic quantile prediction across variable horizons. The second stage employs a streamlined adapter to rapidly inject spatial awareness via meta adaptive learning and prompting. Comprehensive evaluations across diverse domains demonstrate that FactoST-v2 achieves state-of-the-art accuracy with linear efficiency - significantly outperforming existing foundation models in zero-shot and few-shot scenarios while rivaling domain-specific expert baselines. This factorized paradigm offers a practical, scalable path toward truly universal STFMs. Code is available at https://github.com/CityMind-Lab/FactoST.