This study aimed to produce a soil organic carbon (SOC) content map with high accuracy and spatial resolution using the most effective factors in the model. The spatial SOC estimation success of Inverse Distance Weighting (IDW), Ordinary Kriging (OK), Empirical Bayesian Kriging (EBK), Multi-Layered Perception Network (MLP) and MLP-OK Hybrid models were compared to obtain the most reliable model in estimating the SOC content. The study area was located in Besni district in the Southeastern Anatolia Region of Turkey. Total of 132 surface (0-30 cm) soil samples were collected from the covers 1330 km2 land and analyzed for SOC, lime, clay and sand content and soil reaction included in the estimation models. Mean annual precipitation and temperature, elevation, compound topographic index, enhanced vegetation and normalized difference vegetation index, were also used as the inputs in the modelling. The spatial distribution of SOC was determined using a MLP and a two-stage ensemble model (MLP-OK) combining the estimation of OK residuals. Soil surveys and covariates were used to train and validate the MLP-OK hybrid model. The MLP-OK model provided a more accurate estimation of SOC content with minimal estimation errors (ME: -0.028, 45 MAE: 0.042, RMSE: 0.066) for validation points compared to the other models. The MLP-OK model outperformed other models by 75.09 to 77.92%. The MLP-OK model estimated the lower and upper limits of the estimated and the measured values in a consistent manner compared to the other models. The spatial distribution map of SOC content obtained by ANN-kriging approach was significantly affected by ancillary variables, and revealed more detail than other interpolation methods in the northern, central, southwestern and southeastern parts of the study area. The results revealed that the assembling of MLP with OK model can contribute to obtain more reliable regional, national and global spatial soil information.