The GOSAT (Greenhouse gases Observing SATellite) data provide new opportunities the most regionally complete and up-to-date assessment of $CO_2$. However, in practice, GOSAT records often suffer from missing data values mainly due to unfavorable meteorological condition in specific time periods of data acquisition. The aim of this research was to identify optimal spatial interpolation techniques to ensure the continuity of $CO_2$ from samples taken in the North East Asia. The accuracy among ordinary kriging (OK), universal kriging (UK) and simple kriging (SK) was compared based on the combined consideration of $R^2$ values, Root Mean Square Error (RMSE), Mean Error (ME) for variogram models. Cross validation for 1312 random sampling points indicate that the (UK) kriging is the best geostatistical method for spatial predictions of $CO_2$ in the East Asia region. The results from this study can be useful for selecting optimal kriging algorithm to produce $CO_2$ map of various landscapes. Also, data users may benefit from a statistical approach that would allow them to better understand the uncertainty and limitations of the GOSAT sample data.