We present an innovative approach to estimate a lithosphere density distribution model based on seismic tomography and gravity data. In the studied area, the model shows that magmatic events have increased density in the middle to lower crust, which explains the observed positive gravity anomaly. We interpret the densification through crustal intrusion and magmatic underplating. The proposed method has been tested in the Alps but can be applied to other geological contexts.