With the development of remote sensing and geostatistical technology, complex environmental variables are increasingly easily quantified and applied in modelling soil organic carbon (SOC). However, this emphasizes data redundancy and multicollinearity problems adding to the difficulty in selecting dominant influential auxiliary variables and uncertainty in estimating SOC stocks. The current paper considers the spatial characteristics of SOC density (SOCD) to construct prediction models of SOCD on the basis of reducing the data dimensionality and complexity using the principal component analysis (PCA) method. A total of 260 topsoil samples were collected from Chahe town, China. Eight environmental variables (elevation, aspect, slope, normalized difference vegetation index, normalized difference moisture index, nearest distance to construction area and road, and land use degree comprehensive index) were pre-analysed by PCA and then extracted as the main principal component variables to construct prediction models. Two geostatistical approaches (ordinary kriging and ordinary co-kriging) and two regression approaches (ordinary least squares and geographically weighted regression (GWR)) were used to estimate SOCD. Results showed that PCA played an important role in reducing the redundancy and multicollinearity of the auxiliary variables and GWR achieved the highest prediction accuracy in these four models. GWR considered not only the spatial characteristics of SOCD but also the related valuable information of the auxiliary attributes. In summary, PCA-GWR is a promising spatial method used here to predict SOC stocks.