The optimal choice of sites to make spatial prediction is critical for a better understanding of really spatio-temporal data. It is important to obtain the essential spatio-temporal variability of the process in determining optimal design, because these data tend to exhibit both spatial and temporal variability. Two new methods of prediction for spatially correlated functional data are considered. The first method models spatial dependency by fitting variogram to empirical variogram, similar to ordinary kriging (univariate approach). The second method models spatial dependency by linear model co-regionalization (multivariate approach). The variance of prediction method was chosen as the optimization design criterion. An application to CO concentration forecasting was conducted to examine possible differences between the design and the optimal design without considering temporal structure.