Spatial prediction methods have been useful to determine the variability of water quality in space and time due to difficulties in collecting spatial data across extensive spaces such as watershed. This study compares two kriging methods in predicting BOD concentration on the unmonitored sites in the Geum River Watershed and to assess its predictive performance by leave-one-out cross validation. This study has shown that cokriging method can make better predictions of BOD concentration than ordinary kriging method across the Geum River Watershed. Challenges for the application of cokriging on the spatial prediction of surface water quality involve the comparison of network-distance-based relationship and euclidean-distance-based relationship for the improvement in the predictive performance.