Spatial uncertainty of Busan marine clay ground, which commonly occurs during site investigation testing, data analysis and transformation modeling, has been described. In this paper geotechnical uncertainty of shear strength indicator $N_k$ has been quantified in both horizontal direction and vertical direction using geostatistical Kriging method. Most of soil data used are from 25 boring tests, 75 laboratory tests, 124 field vane tests and 25 cone penetration tests (CPT). CPT-$N_k$ data for undrained shear strength determination, which are the most important properties in geotechnical design stages, have been analysed. Comparison between cone factor from conventional CPT-based method and that of geostatistical method shows that geostatistical Kriging method is an ideal tool to quantify the spatial variability of uncertainty from self-correlation of soil property of interest, and can be recommended to identify the spatial distribution of consolidation .md shear strength of soils at any sites concerned.