Process of Earth’s density models creation leads to the solution of direct and inverse gravimetry problems. The inverse problem of gravimetry is a classic example of an ill-posed problem: in the common statement, its solution is not unique and unstably depends on input data. Therefore, it is necessary to determine solutions belonging substantial sets of correctness, choosing reasonable models of an initial approximation. In this paper the application of complex interpretation methods of seismic and gravitational data for the creation of three-dimensional models of crust and the upper mantle are presented. Original algorithms and programs were developed for implementation of these methods. They contain solution of non-linear (structural) inverse problem and the solution of the linear three-dimensional inverse problem taking note of the side sources. Coefficients of the “density-velocity” correlation formulas for a number of geo-traverses were defined. Also, we suggest a technique of tectonic maps construction, which is based on the lithostatic pressure calculation. Its idea can be applied to both two- and three-dimensional cases. In the 2D case we show the way to split the mantle to blocks with vertical boundaries. If lithostatic compensation hypothesis is adopted, the method also allows one to calculate density value for each block. Such separation of the mantle helps to diminish discrepancy between model and observed fields. In 3D case we suggest a method, which can be used to construct tectonic structure maps with information about approximate depth and height of each tectonic block.