The lateral variations of the Earth's gravitational field, deduced from orbital data of artificial satellites, indicate the existence of lateral density variations within the Earth. A density model is computed for the mantle with the following constraints: (1) the model presents perturbations to Gutenberg's earth model which are specified by spherical harmonics through the sixth order and degree; (2) the density anomalies are confined to the mantle and the crust; (3) the anomalies of the crust are determined for n = 2, .. ., 6 and m = 0, .. ., n from crustal thickness, crustal P wave velocity, and P , velocity, and those of the upper mantle for n = 2 and 3 and m = 0, ..., n are related to the lateral variations of seismic traveltime residuals; (4) the unknown density anomalies of the mantle are determined such that the total shear strain energy of the Earth is a minimum; (5) the gravitational potential of the deformed Earth (subject to the density anomalies) on its surface equals the first six orders of the spherical harmonic representation of the measured geopotential; and (6) an isotropic, elastic, and cold mantle, and a liquid core are assumed in the stress analysis.The density anomalies thus obtained exhibit a decreasing feature with depth. In the mist they are of the order OS 0 * 3 g~r n -~, in the upper mantle 0-1 g ~m -~, and in the lower mantle 0.04 g ~m -~, which are within the values deduced from independent seismic measurements. The lateral variations of the associated stress differences at shallow depths correlate with the surface feature of the Earth. This correlation disappears in the deep mantle.