Gravity inversion is the pioneer in exploring the structural characteristics of the Earth, the Moon, and other celestial bodies. Classical gravity inversion methods aim to estimate the 3D subsurface density distribution from the observed 2D surface gravity anomalies, which is an ill-posed problem. Constraints can provide vertical resolution and reduce uncertainty. However, these methods significantly increase the cost of data acquisition. This manuscript presents a novel joint inversion method to estimate subsurface density anomaly via a physics-inspired neural network. The observed signals in...