A scalable Gaussian Process (GP) framework for fusing heterogeneous spatio-temporal data sources (e.g., sensors, satellites, mobility traces) with linear-time inference via state-space methods.