The goal of the present study is to assess the soil quality in Bulgaria using (i) an appropriate set of soil quality indicators, namely primary nutrients (C, N, P), acidity (pH), physical clay content and potentially toxic elements (PTEs: Cu, Zn, Cd, Pb, Ni, Cr, As, Hg) and (ii) respective data mining and modeling using chemometrical and geostatistical methods. It has been shown that five latent factors are responsible for the explanation of nearly 70% of the total variance of the data set available (principal components analysis) and each factor is identified in terms of its contribution to the formation of the overall soil quality-the mountain soil factor, the geogenic factor, the ore deposit factor, the low nutrition factor, and the mercury-specific factor. The obtained soil quality patterns were additionally confirmed via hierarchical cluster analysis. The spatial distribution of the patterns throughout the whole Bulgarian territory was visualized via the mapping of the factor scores for all identified latent factors. The mapping of identified soil quality patterns was used to outline regions where additional measures for the monitoring of the phytoavailability of PTEs were required. The suggested regions are located near to thermoelectric power plants and mining and metal production facilities and are characterized by intensive agricultural activity.