Abstract: The completion of the Bouguer gravity map for the Romanian territory, as well as the calculation of the mean gravity dataset in a 5' x 7.5' grid, enabled a 3D modelling approach for the crustal structure. The geophysical model has been built using information derived from published crustal models based on refraction seismic and borehole data, each layer getting a mean density value describing the main density contrasts. The discrepancies observed between the crustal model gravity effect and the Bouguer gravity map may lead to improvements in the geophysical model in both geometry and petrophysical features. The 3D stripped gravity map that was derived from the above mentioned gravity maps may be interpreted in terms of subcrustal density inhomogeneities, mainly related to the Trans-European Suture Zone and the Vrancea seismogenic area.