Spatiotemporal data consisting of measurements gathered at different times and locations is challenging to analyse due to variability and noise impact across different scales. The authors propose a statistical approach that delivers models of large-scale spatiotemporal datasets applicable to data-analysis tasks of forecasting and interpolation.