An improved mathematical model of the interrelated filtration consolidation processes, heat and mass transfer as well as contact suffusion in heterogeneous soils has been built and researched. A solution of the boundary problem has been programmed with the finite elements method. The influence of the suffusion phenomena on the soil surface subsidence has been analyzed based on a model problem.