Probabilistic forecasting

We have a time series dataset (spatiotemporal, but not an image/video). The dataset is in 3D, where each (x,y,t) coordinate has a numeric value (such as the sea temperature at that location and at that specific point in time). So we can think of it as a matrix with a temporal component. The dataset is similar to this but with just one channel:


We need to predict/forecast the future (next few time steps) values for the whole region (i.e., all x,y coordinates in the dataset) along with the uncertainty.

Can you all suggest any architecture/approach that would suit my purpose well? Thanks!

Not an expert in the area, but the example in the Forecasting III: hierarchical models tutorial may be a good start.