System and method for surveying underground density distributions

A method for geographical surveying includes the steps of determining an observation grid comprising observation points above a solution volume in the ground; making gravitational measurements at the observation points over a survey area; compiling a matrix having elements relating a set of gravitational accelerations and gravity gradient values at the points to the mass values for volume elements in the solution volume; and calculating the three-dimensional density distribution of the volume under the survey area based on the gravitational measurements and the inverse of the matrix. By measuring the magnetic field and gradient at the observation points, the three-dimensional magnetic susceptibility distribution of the volume under the survey area can be determined.