Geostatistical data or point-referenced data have the information on the monitoring stations of interest where the observations are measured. Practical geostatistical data are obtained from a wide variety of observational monitoring networks that are mainly operated by the Korean government. When we analyze geostatistical data and predict the expectations at unobservable locations, we can improve the reliability of the prediction by utilizing some relevant spatial data obtained from different observational monitoring networks and blend them with the measurements of our main interest. In this paper, we consider the hierarchical spatial linear model that enables us to link spatial variables from different resources but with similar patterns and guarantee the precision of the prediction. We compare the proposed model to a classical linear regression model and simple kriging in terms of some information criteria and one-leave-out cross-validation. Real application deals with Sulfur Dioxide($SO_2$) measurements from the urban air pollution monitoring network and wind speed data from the surface observation network.