In soil magnetometry, two types of measurements are usually performed. The first type is measurements performed on the soil surface, frequently using an MS2D sensor. The second type includes measurements of magnetic susceptibility carried out in the soil profile, usually to a depth of about 30 cm. Up to now, such measurement results were analyzed separately. However, it is possible and advantageous to integrate these two types of measurements. The goal of the study was to integrate measurements of magnetic susceptibility performed on the soil surface and in the soil profile. More specifically, the goal was to obtain 3-dimensional spatial distributions of magnetic susceptibility of the topsoil horizon. Results show that it is possible to effectively integrate measurements of magnetic susceptibility performed on the soil surface and in the soil profile. Moreover, the 3-dimensional spatial distribution that is obtained shows the magnetic susceptibility of the top 20 cm of soil, which includes the soil horizons where most of the heavy metals are accumulated. The analysis of such a spatial distribution can be very helpful in delineating areas where the heightened magnetic susceptibility is a result of the influence of anthropogenic pollution from those areas where it results from lithogenic origin. It is possible to investigate where the volumes of soil with heightened magnetic susceptibility are located in the soil profile and in this way investigate which characteristic type of soil profile it is.