Soil bulk density is one of the main direct indicators of soil health, and is an important aspect of models for determining agroecosystem services potential. By way of applying multi-regression methods, we have created a distributed prediction of soil bulk density used subsequently for topsoil carbon stock estimation. The soil data used for this study were from the Slovakian partial monitoring system-soil database. In our work, two models of soil bulk density in an equilibrium state, with different combinations of input parameters (soil particle size distribution and soil organic carbon content in %), have been created, and subsequently validated using a data set from 15 principal sampling sites of Slovakian partial monitoring system-soil, that were different from those used to generate the bulk density equations. We have made a comparison of measured bulk density data and data calculated by the pedotransfer equations against soil bulk density calculated according to equations recommended by Joint Research Centre Sustainable Resources for Europe. The differences between measured soil bulk density and the model values vary from -0.144 to 0.135 g cm^-3 in the verification data set. Furthermore, all models based on pedotransfer functions give moderately lower values. The soil bulk density model was then applied to generate a first approximation of soil bulk density map for Slovakia using texture information from 17 523 sampling sites, and was subsequently utilised for topsoil organic carbon estimation.