We present geophysical models for the lithospheric mantle of the Siberian craton, with focus on its structureand thermo-compositional heterogeneity as constrained by various geophysical data. The latter include thermalstructure of the lithosphere based on surface heat flow data and supported by regional xenolith P-T arrays,lithosphere density heterogeneity as constrained by free-board and satellite gravity data, and the non-thermalpart of upper mantle seismic velocity heterogeneity based on joint analysis of thermal and seismic tomography data.Density structure of the cratonic lithosphere constrained independently by free-board and satellite gravityshows significant lateral variations, that are well correlated with crustal structure, surface tectonics, and regionalxenolith data. Proterozoic sutures and intracratonic basins are manifested by an increase in mantle density as comparedto light and strongly depleted lithospheric mantle of the Archean nuclei, particularly below the Anabar shield.Since we cannot identify the depth distribution of density anomalies, we complement the approach by seismicdata. An analysis of temperature-corrected seismic velocity structure indicates strong vertical and lateralheterogeneity of the cratonic lithospheric mantle, with a pronounced stratification in many Precambrian terranes.We discuss lateral and vertical heterogeneity of the cratonic lithosphere discussed in connection to regionaltectono-thermal evolution.